Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aczc_v1_0.gz(5 Kbytes)
Manuscript Title: On numerical Bessel transformation.
Authors: B. Sommer, J.G. Zabolitzky
Program title: FTRANS
Catalogue identifier: ACZC_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 16(1979)383
Programming language: Fortran.
Computer: CDC 7600.
Operating system: SCOPE 3.4.4.
RAM: 7K words
Word size: 60
Keywords: General purpose, Fourier transform, Bessel transform, Spherical bessel, Filon integration, Numerical integration.
Classification: 4.6, 4.7.

Nature of problem:
We present a computer program to calculate the three dimensional Fourier or Bessel transforms and definite integrals with Bessel functions. Numerical integration of systems containing Bessel functions occurs in many physical problems, e.g. electromagnetic form factor of nuclei, all transitions involving multipole expansions at high momenta.

Solution method:
Filon's integration rule is extended to spherical Bessel functions. The numerical error is of the order of the Simpson error term of the function which has to be transformed. Thus one gets a stable integral even at large arguments of the transformed function.

Running time:
The test run compiled in 1.37 s and executed in 0.45 s on the CDC 7600.