lcfa.sph: Array of logarithmic derivative and ratio of Spherical Bessel...

Description Usage Arguments Details See Also Examples

Description

Array of logarithmic derivative and ratio of Spherical Bessel functions.

Usage

1
lcfa.sph(nmax, x, code = "C", NMAX = 200)

Arguments

nmax

Maximun value of n.

x

The argument of the functions. Can be complex.

code

If the C code or R built in.

NMAX

Maximum number of iterations on the Lentz algorithm.

Details

Calculate by downward recurrence the derivative of Spherical Bessel functions c_n=j_n'(x)/j_n(x) and simultaneously the ratio ρ_n=j_n(x)/j_{n+1}(x), from 0 to nmax. The starting values are calculated by Lentz continued fraction method.

See Also

lcfe.sbl, lcfe.sbd, lcfe.sbi.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
nmax<-10
x<-5
u.c<-lcfa.sph(nmax,x,code="C")
u.r<-lcfa.sph(nmax,x,code="R")
u<-data.frame(
   # Logarithmic Derivatives
   C.LogDev=u.c$cn,R.LogDev=u.r$cn,
   # Ratio between Spherical Bessel functions
   C.SphRat=u.c$rn,R.SphRat=u.r$rn)
print(u)

wendellopes/rvswf documentation built on May 4, 2019, 4:19 a.m.