Journal article
Algorithm to Calculate a Large Number of Roots of the Cross-Product of Bessel Functions
Lab. d'Electromagn. et d'Acoust., (LEMA), Ecole Polytech. Fed. de Lausanne (EPFL), Lausanne, Switzerland1
This paper describes an algorithm to calculate a large number of roots of the cross-product of Bessel functions and of their first derivatives. The algorithm initially finds the roots of the zeroth order using an auxiliary function that exhibits the same roots as the original cross-products but with better behavior for numerical root search with the Newton-Raphson algorithm.
In order to find the roots for higher orders, the algorithm follows a pyramidal scheme using the interlacing property of the cross-product of Bessel functions. The algorithm shows globally convergent behavior for a large range of values of the argument and of the order of the Bessel functions. The roots can be computed to any precision, limited only by the computer implementation, and the convergence is attained in six iterations per root in average, showing a much better performance than previous works for the calculation of these roots.
Language: | English |
---|---|
Publisher: | IEEE |
Year: | 2013 |
Pages: | 2180-2187 |
ISSN: | 15582221 and 0018926x |
Types: | Journal article |
DOI: | 10.1109/TAP.2012.2231929 |
Algorithm Bessel functions Convergence Electromagnetic scattering Electromagnetic waveguides Electromagnetics Equations Indexes Mathematical model McMahon's expansion Newton-Raphson algorithm Newton-Raphson method Newton–Raphson auxiliary function convergence cross-product cross-product interlacing property cross-product roots globally convergent behavior interlacing properties numerical root search behavior pyramidal scheme