lit
    JINRLIB

DBESI0, DBESK0, DEBSI0, DEBSK0 - I0, I1, K0, K1

C313

: K.S.Koelbig :

popaa ce ecce I0(X),I1(X),K0(X),K1(X) coo apa X (X - ooo K0(X),K1(X))

:

: - FUNCTION
: - DBESI0,DBESI1,DBESK0,DBESK1,
DEBSI0,DEBSI1,DEBSK0,DEBSK1

:

o apco pa:

DBESI0(X) a a I0(X)
DBESI1(X) a a I1(X)
DBESK0(X) a a K0(X)
DBESK1(X) a a K1(X)
DEBSI0(X) a a EXP(-ABS(X))*I0(X)
DEBSI1(X) a a EXP(-ABS(X))*I1(X)
DEBSK0(X) a a EXP(ABS(X))*K0(X)
DEBSK1(X) a a EXP(ABS(X))*K1(X).

a ap ocac opaopo DOUBLE PRECISION.

:

coec apoca cee pa eea.

:

- 14 ao, - .

Opae:

1. X<=0 DBESK0, DBESK1, DEBSK0, DEBSK1 oaac pa , coo:'DBESIK...NON-POSITIVE ARGUMENT X= ', pa a popa.
2. DEBSI0,DEBSI1,DEBSK0,DEBSK1 p p po ao , a DEXP.

:

  1. C.W.Clenshaw. Chebyshev series for mathematical functions.
    ( National Physical Lab., Math. Tables, vol.5, London, 1962).
  2. Y.L.Luke, Mathematical functions and their approximations
    ( Academic Press, New York,1975).
:
       . . .
       Y=DBESK1(0.1D0)
       . . .
:
       Y=9.853844780870537


home up e-mail