Journal of East China Normal University(Natural Sc ›› 2019, Vol. 2019 ›› Issue (1): 76-82,92.doi: 10.3969/j.issn.1000-5641.2019.01.009

• Physics and Electronics • Previous Articles     Next Articles

Program for calculating the integer order of Bessel functions with complex arguments

REN Hong-hong1, GUO Ying-chun1, WANG Bing-bing2   

  1. 1. School of Physics and Materials Science, East China Normal University, Shanghai 200241, China;
    2. Laboratory of Optical Physics, Beijing National Laboratory of Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
  • Received:2017-12-20 Online:2019-01-25 Published:2019-01-24

RichHTML

18

PDF

413

Abstract: Fortran source code for calculating the integer order of Bessel functions of the first kind with complex arguments is presented. The method is based on the backward recurrence relation of Bessel functions. Values of the Bessel function generated by our program are in-agreement with the values generated by Matlab to at least 12 significant digits. We use the program to calculate the integer order of Bessel functions of the first kind with pure imaginary arguments, provided by Xu Shiliang's Fortran algorithm assembly. The results show that the first 6 significant digits are accurate. We also analyze the algorithm for calculating Bessel functions with complex arguments, which use the infinite sum of the product of the real arguments of the Bessel function and the pure imaginary arguments of the Bessel function, provided by Xu Shiliang's algorithm assembly. The results show that this algorithm does not always get accurate values for Bessel functions with complex arguments. The reason lies with the fact that the term in the sum larger than the function value causes the loss of significant digits.

Key words: Bessel function with complex arguments, Fortran source code, backward recurrence relation of Bessel function

CLC Number: