A numerical method to solve linear integro-differential equations is
presented. This method has been used to solve the QCD Altarelli-Parisi
evolution equations within the H1 Collaboration at DESY-Hamburg. Mathematical
aspects and numerical approximations are described. The precision of the method
is discussed.Comment: 18 pages, 4 figure