lit
    JINRLIB

COUL1 -

C316

: K.S.Kolbig :

opopaa paaa ao a po ooo oa
c361 c aao ooc cpoa eta>=<0, ro>p.=0 (.ooc),
L=0,1,2,...,Lmax<=100. c316 c paoo pa (c.[1], .14):

c316

:

: - SUBROUTINE
: - COUL1

:

CALL COUL1(ETA,RO,LMAX,N,F), :
  ETA,RO,LMAX cooc op a eta, ro, Lmax.

N - (INTEGER) co a aa p, N<13;
F - (REAL*8) acc papoc >=LMAX+1;
a F0,F1,...,FL acac opopao
F(1),F(2),...,F(LMAX+1).

:

ro>=0 (), 0<=LMAX<=100. po oa ooc o o, oa ro>=100 () |eta|>=100 ().

:

pop ooc aopa c.[3]. o cc, o o N aa p, oa apoaa o coa / ro, |eta| o.
a ro L>>1 (app, ro<10-4 L=50, ro<10-1 L=100, eta=1 oo ca) o o ca, oa ooo op oo N aa p paa. a oaac pa c eta, c ro=0, .. c316=0 c eta.

:

  1. B ca, oa pocc c coc, aac aoca:
    COUL1...CONVERGENCE DIFFICULTY IN THE GENERATION OF THE COEFICIENTS LAMBDA SUB L. COUL1...CONVERGENCE DIFFICULTY.
  2. Ec ro<0, aac aoca:
    COUL1...RO IS LESS THEN ZERO.

:

opopaa pca coo oppa po oo pop, acao W.Gautschi [2].
aac o pao oca [2].

:

  1. Cpao o ca o a .paoa .a.
    oca, "Haa", 1979.
  2. Gautschi W., Algorithm 292, Rgular Coulomb Wav Functions, Comm. ACM 9, 1966, p.793-795.
  3. Kolbig K.S., Crtification of Algorithm 292, CERN-DD/CO, 67/9.
:
       . . .
       IMPLICIT REAL*8
       DIMENSION F(3)
       ETA=15.0D0
       RO=9.0D0
       LMAX=2
       N=11
       CALL COUL1(ETA,RO,LMAX,N,F)
       . . .

:

       F(1)= .455113609922D-07
       F(2)= .410013566987D-07
       F(3)= .333068257494D-07


home up e-mail