Title ofprogram.' CLMINT
wavelength much smaller than the integration interval. When we choose R 2 = i~ii,the convergence of the integral is very slow.
Catalogue number.' ACCM
Method of solution
Program obtainable from: CPC Program Library, Queen's Uni-If we split up the integrand into a rapidly and slowly varying versity of Belfast, N. Ireland (see application form in this issue) function of the integration variable r, it is possible to construct an asymptotic series for the integral with the fast oscillating Computer: Burroughs 7700 at Eindhoven University of Tech-integrand. The remaining integral is easily obtained by Gausnology, Eindhoven, The Netherlands; CDC 175/100 at Utrecht sian quadrature. Because we are integrating over large r values University, Utrecht, The Netherlands (compared to the classical turning points), it is possible to use the familiar WKB-approximation of Coulomb wave functions Programming language used: FORTRAN 77 in the step-by-step integration. The subroutine has been set up to acl~ievea relative accuracy of i0~.
Operating system: MCP (Burroughs); NOS/BE (CDC) Restriction on the complexity of the problem Peripherals used: line printer, card reader The subroutine has only been tested for significant 1, i~and k values i.e. 0 ~1< 2000, 0 <~< 1000, 0 <k < 50. The interval No. of lines in the combined program and test deck: 1544 [R1, R2] can be any part of [0, 00). When R1 does not exceed the lowest turning point (especially R1 0), only the integrals No. of bits in a word: 52 (Burroughs); 60 (CDC)
with the regular Coulomb wave functions F1(~, p) are reliable. Furthermore the case A = 0 is not included.