lit
    JINRLIB

COUL2 -

C326

: J.H.Gunn,K.S.Koelbig :

popaa c oo oa c326_1,c326_2
p poo o ro c326_3,c326_4,
a a ao sigma 0.1<=eta<=30, 5<=ro<=30 L=0,1,2,...,Lmax; Lmax<=100.
c326_1 c326_2 opc a a ac p paoo pa (c.[1],.14 [2])

c326_5
c aco:
c326_6
c326_7
c326_8

:

: - SUBROUTINE
: - COUL2

:

CALL COUL2(F,FD,G,GD,SIG,RHO,ETA,LMAX), :

F - acc papoc 101, copa c a FL, FL→F(L+1);
FD - acc papoc 101, copa c a F'L, F'L→FD(L+1);
G - acc papoc 101, copa c a GL, GL→G(L+1);
GD - acc papoc 101, copa c a G'L, G'L→GD(L+1);
SIG - acc papoc 101, copa c a GL, GL→SIG(L+1);
RHO - aao a ro;
ETA - aao a eta;
LMAX - aao a Lmax (oo pac a oo ).

:

Kpp ooc oca [3] [4]. c oa o 6-7 p c ao. a oaa pp, a aa ooc oop ca ocac.

:

B ca, oa ro eta a aao oac, a a ac coo:
COUL2 ... RHO OR ETA NOT IN RANGE pa oaac pa .

:

Ec ooo c oo FL F'L , o oc cooa popa C0UL1 (316).

:

  1. Cpao o ca o pa M.paoa .a.
    oca, "aa", 1979.
  2. Gunn J.H. Algorithm 300, Coulomb wave functions, Comm. ACM,10, 1967, 244-245.
  3. Kolbig K.S. Crtification of Algorithm 300, Coulomb wav functions, Comm. ACM,12, 1969, No.5.
  4. Kolbig K.S. Rmark on Algorithm 300, Coulomb wav functions, CERN/DD/CO/69/3.
:
       . . .
       IMPLICIT REAL*8 (A-H,O-Z)
       DIMENSION F(101),FD(101),G(101),GD(101),SIG(101)
       RH0=5.0D0
       ETA=10.0D0
       LMAX=20
       CALL COUL2(F,FD,G,GD,SIG,RH0,ETA,LMAX)
       . . .
:
       F(1)=     .1720746435D-05     F(21)=     .4042841079D-17
       FD(1)=    .3097601327D-05     FD(21)=    .1837788672D-16
       G(1)=     .1676374624D+06     G(21)=     .2778705490D+17
       GD(1)=   -.2793706069D+06     GD(21)=   -.1210368262D+18
       SIG(1)=   .1380291297D+02     SIG(21)=   .3057615005D+02


home up e-mail