besj(3f) - [M_bessel] compute the J Bessel function for a given argument and order
Synopsis
Description
Options
Example
subroutine besj(x,n,bj,d,ier)
N must be greater than or equal to zero, but it must be less than
20+10*X-X** 2/3 For X less than or equal to 15 90+X/2 For X greater than 15
none
Recurrence relation technique described by H. Goldstein and R.M. Thaler,"Recurrence Techniques for the Calculation of Bessel Functions",M.T.A.C.,V.13,PP.102-108 and I.A. Stegun and M. Abramowitz,"Generation of Bessel Functions on High Speed Computers",M.T.A.C.,V.11,1957,PP.255-257
Description of Parameters
X The argument of the J Bessel function desired N The order of the J Bessel function desired BJ The resultant J Bessel function D Required accuracy IER Resultant error code whereIER=0 No error IER=1 N is negative IER=2 X is negative or zero IER=3 Required accuracy not obtained IER=4 Range of N compared to X not correct (see remarks)
Nemo Release 3.1 | besj (3) | February 23, 2025 |