This algorithm uses a high-order, variable step Runge-Kutta like method in the region where the potential term dominates, and an exponential or Bessel fitted method in the asymptotic region. This approach can be used to compute scattering phase shifts in an efficient and reliable manner. A Fortran p
FORTRAN program for a numerical solution of the nonsinglet Altarelli-Parisi equation
✍ Scribed by R. Kobayashi; M. Konuma; S. Kumano
- Publisher
- Elsevier Science
- Year
- 1995
- Tongue
- English
- Weight
- 760 KB
- Volume
- 86
- Category
- Article
- ISSN
- 0010-4655
No coin nor oath required. For personal study only.
✦ Synopsis
We investigate a numerical solution of the flavor-nonsinglet Altarelli-Parisi equation with next-to-leading-order as 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 (Fb 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-4000/500. Accurate evolution results are obtained by taking approximately twenty Laguerre polynomials.
📜 SIMILAR VOLUMES
DIFEQ with boundary conditions ~p(a)= /~(b) = 0 are solved. The main purpose of this paper is to present a perturbative proce-Catalogue number: ACCC dure for the calculations of approximate eigenvalues of the Schrodinger equation. Program obtainable from: CPC Program Library, Queen's University of
dq N dp N ϭ const, (3a) A numerical method for the time evolution of systems described by Liouville-type equations is derived. The algorithm uses a lattice of numerical markers, which follow exactly Hamiltonian trajectories, to represent the operator d/dt in moving (i.e., Lagrangian) coordinates. H
The Chebyshevian Multi-step Theory of Lyche has been developed and applied to the numerical solution of the radial form of the Schrodinger equation. Significant improvements over previously reported approaches are found.