FORTRAN program for a numerical solution of the nonsinglet Altarelli-Parisi equation Original Research Article

We investigate a numerical solution of the flavor-nonsinglet Altarelli-Parisi equation with next-to-leading-order αs corrections by using Laguerre polynomials. Expanding a structure function (or a quark distribution) and a splitting function by the Laguerre polynomials, we reduce an integrodifferential equation to a summation of finite number of Laguerre coefficients. We provide a FORTRAN program for Q2 evolution of nonsinglet structure functions (F1, F2, and F3) and nonsinglet quark distributions. This is a very effective program with typical running time of a few seconds on SUN-IPX or on VAX-4000500. Accurate evolution results are obtained by taking approximately twenty Laguerre polynomials.