A General Analytic Approximation Technique for the Modified Bessel Functions with Fractional Order Iv (x), Applied to I1/6(x) and I1/7(x)

Abstract

A general technique to obtain simple analytic approximations for the first kind of modified Bessel functions. The general procedure is shown, and the parameter determination is explained through the applications to this particular case I 1/6 ( x ) and I 1/7 ( x ) . In this way, it shows how to apply the technique to any particular order ν , in order to obtain an approximation valid for any positive value of the variable x. In the present method power series and asymptotic expansion are simultaneously used. The technique is an extension of the multipoint quasirational approximation method, MPQA. The main idea is to look for a bridge function between the power and asymptotic expansion of the I 1/6 ( x ) , and similar procedure for I 1/7 ( x ) . To perform this, rational functions are combined with hyperbolic ones and fractional powers. The number of parameters to be determined for each case is four. The maximum relative errors are 0.0049 for ν=1/6 , and 0.0047 for ν=7 . However, these relative errors decrease outside of the small region of the variables, wherein the maximum relative errors are reached. There is a clear advantage of this procedure compared with any other ones.

Share and Cite:

Olivares, J. , Martín, P. and Veliz, M. (2024) A General Analytic Approximation Technique for the Modified Bessel Functions with Fractional Order Iv (x), Applied to I1/6(x) and I1/7(x). Advances in Pure Mathematics, 14, 768-777. doi: 10.4236/apm.2024.149041.

1. Introduction

Bessel functions have become very important in Physics, Electrodynamics and other areas of science [1]-[7]. In order to simplify its calculations several procedures can be used. The power series and asymptotic expansions are well known. However, they are mainly useful for small or large values of the variable. In recent works, analytic approximations valid for all values of the variable have been determined with high accuracy for the modified Bessel functions I 0 ( x ) and I 1 ( x ) , using the multipoint quasirational approximations technique, MPQA [8]-[11]. Approximations of polynomial type can also be found as [12] well as other interesting [13] [14]. Now a general expression to determine precise approximations for the modified Bessel functions of fractional order will be shown, and the procedure to determine the parameters explained by its applications to the particular case of the main I 1/6 ( x ) and I 1/7 ( x ) will be considered. The fact of working with not entire order involves some changes, which will be explained in Section 2, where the theoretical treatment will be exposed. The main idea of the MPQA technique is to determine a bridge function connecting the power and asymptotic expansion in such a way that part of both expansions can be reproduced. In order to do that, it is necessary to use rational functions combined with hyperbolic ones, as well as fractional powers. The new simple approximations functions here presented will be shown to be good for all real values of the variable x with good accuracy and easy to calculate. The parameter determination as well as the relative errors of the approximations will be considered in Section 3 for I 1/6 ( x ) , and in Section 4 for I 1/7 ( x ) . Finally, in the final section or Conclusion, the main results will be shown.

2. Theoretical Treatment

In this section, it will be first discussed in detail, how to determine the structure and parameters of the modified Bessel function of order ν, as a general case, and later the application to the particular case will be considered I 1/6 ( x ) . The power series I ν ( x ) is well known [9], and it is given by

I ν ( x )= ( x 2 ) ν k=0 1 k!Γ( k+1+ν ) ( x 2 ) 2k (1)

The second series important in the kind of treatment here considered is the asymptotic expansion, given by

I ν ( x )~ e x 2πx ( 1 4 ν 2 1 8x + ) (2)

In the present technique a bridge function connecting both series power and asymptotic, MPQA technique, must be defined. The way to determine this particular function is performed using several unknown parameters to be determined later. There are several ways to do this, depending on how many parameters are used, which also leads to approximations with high accuracy. Here we will consider a simple approximation producing an accuracy adequate to most of the applications. With this idea in mind, the bridge function will be defined as

I ν ( x )= | x | ν cosh( x ) 2 ν Γ( ν+1 ) ( 1+ λ 2 x 2 ) ( 2ν+1 )/4 p 0ν + p 1ν x 2 1+ q ν x 2 (3)

where the parameters p 0ν , p 1ν , q ν and λ will be determined using the power series and the asymptotic expansions already known. It is interesting to point out, that all the functions in Equation (3) are even functions since in this way the symmetry of the function I ν ( x ) present and the approximation will be good for positive values of x, and also for negative values.

The procedure to be followed is first to show the way to determine the equations for the parameters determination. This will be performed using the power series and the asymptotic expansions. Later we will proceed to the calculation of the parameters and the way to do this avoiding the usual problems of the defects in Padé technique, that is, a zero in the denominator with a nearby zero at the numerator [11]. These themes will be better explained in relation to the examples here treated in detail.

The procedure to determine the parameter will be better explained by application to particular cases, which in this work will be by particular case of ν=1/6 and ν=1/7 .

In the case of ν=1/6 , the power series will be written as

I 1/6 ( x )= x 1/6 2 1/6 Γ( 7 6 ) ( 1+ 3 14 x 2 + 9 728 x 4 + ) (4)

and in the case of the asymptotic expansion, this will be

I 1/6 ( x )~ e x 2πx ( 1+ 1 9x + 5 81 x 2 + ) (5)

The idea is how to find a bridge function linking both expansions. It is clear that the rational function must be multiplied by fractional powers of x, in order to reproduce the structure of the power series. The power series shows also that after this factor, the series has only even powers. This kind of symmetry has to be also preserved in the rational function, and also new factors will be considered, in order to obtain the correct asymptotic structure. In this case, there is also the exponential factor e x , but this factor is not even as required in the preceding lines, therefore this must be replaced by cosh( x ) , which has the exponential behavior at the infinite, and the symmetry of I 1/6 ( x ) is preserved. There is also the fractional power 1/ x , in the asymptotic expansion, which has to be considered, as well as the fractional power x 1/6 coming from the power series considerations. In order to accomplish all the restrictions coming from the power and asymptotic expansions, the structure of the analytic approximations I 1/6 ( x ) will be written as

I ˜ 1/6 ( x )= x 1/6 cosh( x ) 2 1/6 Γ( 7 6 ) ( 1+ λ 2 x 2 ) 1/3 ( p 0 + p 1 x 2 1+q x 2 ) (6)

where the important point is that

x 1/6 x 2/3 = 1 x 1/2 (7)

as it is required. Now the parameters p 0 , p 1 ,q and λ have to be determined, as it will be shown in the next section.

In the case of I 1/7 ( x ) , the power series will be

I 1/7 ( x )= x 1/7 2 1/7 Γ( 8 7 ) ( 1+ 7 32 x 2 + 49 3840 x 4 + ) (8)

and the asymptotic expansion

I 1/7 ( x )~ e x 2πx ( 1+ 45 392x + ) (9)

Now, using an analysis in some ways similar to the previous one, the analytic approximation will have the structure

I ˜ 1/7 ( x )= x 1/7 cosh( x ) 2 1/7 Γ( 8 7 ) ( 1+ λ ˜ 2 x 2 ) 9/ 28 ( p ˜ 0 + p ˜ 1 x 2 1+ q ˜ x 2 ) (10)

where now

2 9 28 1 7 = 92 14 = 1 2 (11)

3. Parameters Determination and Results

In the case of the I 1/6 ( x ) , two equations from the power series and one equation for the asymptotic expansion will be used. The fourth equation for the parameters will come from numerical calculations, as will be explained below. Before equalizing power terms from the series for I 1/6 ( x ) and I ˜ 1/6 ( x ) , it is necessary to rationalize the equation, and in this way, nonlinear equations will be avoided. Thus the product will be,

( 1+ λ 2 x 2 ) 1/3 ( 1+q x 2 ) I 1/6 ( x ) =( 1+ 1 3 λ 2 x 2 + )( 1+q x 2 )( 1+ 3 14 x 2 + ) =1+( 1 3 λ 2 +q+ 3 14 ) x 2 + ( p 0 + p 1 x 2 )( 1+ 1 2 x 2 + ) (12)

Therefore

p 0 =1 (13)

1 2 p 0 + p 1 = 1 3 λ 2 +q+ 3 14 (14)

The third equation for this case, it is obtained from the asymptotic expansion as

p 1 2× 2 1/6 Γ( 7 6 ) λ 2/3 q = 1 2π (15)

From these three equations, all the parameters are obtained as a function of λ . Thus

q= 1 3 λ 2 + 2 7 1 2 1/6 Γ( 7 6 ) 2 π λ 2/3 (16)

p 1 = 2 1/6 Γ( 7 6 ) 2 π λ 2/3 ( 1 3 λ 2 + 2 7 1 2 1/6 Γ( 7 6 ) 2 π λ 2/3 ) (17)

Now we have to avoid the so-called de “deffects” in the Padé method, that is a zero at the denominator with another zero at the numerator near the first one, but not exactly equal. Our procedure to avoid this problem is to make a plot of q as a function of λ , and to avoid the negative values of q. It is clear that if q is negative, there real zeros at the denominator one positive and anothers negative. This means that there is the so-called deffect in the approximation which are want to avoid. Thus, Figure 1 shows the values of λ , where q is positive.

Figure 1. Plot of the curve q( λ ) , in the regions ( 0,0.92 )( 1.38,+ ) the values of q( λ ) are positive.

Now, it is obtained that for λ in the intervals ( 0,0.92 )( 1.38,+ ) , the values of q verify q>0 . In order to determine the optimum value of λ , the procedure will be to give values to λ inside the right intervals, to determine the parameters for each value, and to calculate the maximum relative error ϵ max ( λ ) for that value. The optimum λ will be that with the lowest ϵ max ( λ ) . In this case, this values is λ 0 =0.3675 , and the approximation will be

I ˜ 1/6 = cosh( x ) x 1/6 2 1/6 Γ( 7 6 ) ( 1+ λ 0 2 x 2 ) 1/3 ( 1+ 2 1/6 Γ( 7 6 ) 2 π λ 0 2/3 x 2 ( 1 3 λ 0 2 + 2 7 1 2 1/6 Γ( 7 6 ) 2 π λ 0 2/3 ) 1+( 1 3 λ 0 2 + 2 7 1 2 1/6 Γ( 7 6 ) 2 π λ 0 2/3 ) x 2 ) (18)

The relative error of this approximation is given in Figure 2. The maximum values of the errors are ϵ max ( 2.4 )=0.0049 for x=2.4 , and ϵ max ( 11.1 )=0.004 . However relative errors of the approximation out of the small regions around these values, are smaller than 0.002.

Figure 2. Relative errors of the approximation of I ˜ 1/6 ( x ) as a function of the variable x.

The same procedure applied to I 1/7 ( x ) , leads to the equations

( 1+ λ ˜ 2 x 2 ) 9/ 28 ( 1+ q ˜ x 2 ) I 1/7 ( x ) =( 1+ 9 28 λ ˜ 2 x 2 + )( 1+ q ˜ x 2 )( 1+ 7 32 x 2 + ) =1+( 9 28 λ ˜ 2 + q ˜ + 7 32 ) x 2 + ( p ˜ 0 + p ˜ 1 x 2 )( 1+ 1 2 x 2 + ) (19)

The equation to determine the parameters is now

p ˜ 0 =1 (20)

1 2 p ˜ 0 + p ˜ 1 = 9 28 λ ˜ 2 + q ˜ + 7 32 (21)

p ˜ 1 2× 2 1/7 Γ( 8 7 ) 2 π λ ˜ 9/ 14 q ˜ = 1 2π (22)

From these equations, it is obtained

q ˜ = 9 28 λ ˜ 2 + 9 32 1 2 1/7 Γ( 8 7 ) 2 π λ ˜ 9/ 14 (23)

p ˜ 1 = 2 1/7 Γ( 8 7 ) 2 π λ ˜ 9/ 14 ( 9 28 λ ˜ 2 + 9 32 1 2 1/7 Γ( 8 7 ) 2 π λ ˜ 9/ 14 ) (24)

The plot of q ˜ as a function of λ ˜ is shown in Figure 3. The q ˜ is positive for λ ˜ ( 0,0.93 )( 1.35,+ ) . Now looking for λ ˜ optimum, it is found λ ˜ =0.37 . With this value, the parameters are determined and the result will be,

I ˜ 1/7 ( x )= cosh( x ) x 1/7 2 1/7 Γ( 8 7 ) ( 1+ λ ˜ 0 2 x 2 ) 9/ 28 ( 1+ 2 1/7 Γ( 8 7 ) 2 π λ ˜ 0 9/ 14 x 2 ( 9 28 λ ˜ 0 2 + 9 32 1 2 1/7 Γ( 8 7 ) 2 π λ ˜ 0 9/ 14 ) 1+( 9 28 λ ˜ 0 2 + 9 32 1 2 1/7 Γ( 8 7 ) 2 π λ ˜ 0 9/ 14 ) x 2 ) (25)

The relative errors for this approximation are shown in Figure 4. The maximum values of the relative errors will be in this case ϵ max ( 2.3 )=0.0047 and ϵ max ( 10.8 )=0.005 . As in the case of I 1/6 ( x ) , the relative errors are much smaller outside of the nearby regions of x=2.3 and x=10.8 .

4. Conclusion

In this work a technique to obtain quasi-rational analytic approximations for the

Figure 3. Plot of the curve q ˜ ( λ ˜ ) , in the regions of λ ˜ ( 0,0.93 )( 1.35,+ ) the values of q ˜ ( λ ˜ )>0 .

Figure 4. Relative errors of the approximation I ˜ 1/6 ( x ) as a function of the variable λ .

modified Bessel functions of general order ν is shown, and the structure of the function in the general case is shown. However, the procedure to determine the parameters of the approximation is shown by application to the particular cases of ν=1/6 and ν=1/7 . The analytic approximations of I 1/6 ( x ) and I 1/7 ( x ) have been presented using rational functions with second order polynomials combined with the hyperbolic functions cosh( x ) and the fractional powers x 1/6 and x 1/7 , corresponding to each case. The numbers of parameters to be determined are four for each Bessel function. However, the accuracy is about as low as 0.005 for each case. The analytic approximations are very simple compared with their accuracy. The method used is an extension of the MPQA technique, and this is a clear advantage compared with other ones, since a unique approximation is valid for any value of the variable x, and the symmetry is preserved. The procedure here presented for I 1/6 ( x ) and I 1/7 ( x ) is also very general, and a general structure for any modified Bessel function of order ν is also described, therefore new approximations functions for any ν can be also determined, in agreement with the needs of other authors.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

References

[1] Watson, G.N. (1966) A Treatise on the Theory of Bessel Functions. 2nd Edition, Cambridge University Press, 85-224.
[2] Jackson, J.D. (1998) Classical Electrodynamics. 3rd Edition, John Willey & Sons, Inc., 112-119.
[3] Rothwell, E.J. (2009) Exponential Approximation of the Bessel Functions I0(x), J0(x), I1(x), J1(x), Y0(x) and , with Applications to Electromagnetic Scattering, Radiation and Difraction. [EM Programer’s Notebooks]. IEEE Antennas and Propagation Magazine, 51, 138-147.
https://doi.org/10.1109/MAP.2009.5251217
[4] Khosravian-Arab, H., Dehghan, M. and Eslahchi, M.R. (2017) Generelized Bessel Functions: Theory and their Applications. Math. Mathematical Methods in the Applied Sciences, 40, 6389-6410.
https://doi.org/10.1002/mma.4463
[5] El Halba, E.M., Boustimi, M., Ez-zariy, L. and Belafhal, A. (2017) Focusing Properties of Radially Polarized Modified Bessel-Modulated Gaussian Beam by a High Numerical Aperture Objective. Optics & Laser Technology, 88, 40-53.
https://doi.org/10.1016/j.optlastec.2016.08.008
[6] Lyons, M.E.G., Joy Solomi, R. and Rajendran, L. (2023) Formal Analysis of Isothermal Reaction/Diffusion in a Defined General Geometry. International Journal of Electrochemical Science, 18, Article 100296.
https://doi.org/10.1016/j.ijoes.2023.100296
[7] Segura, J. (2023) Simple Bounds with Best Possible Accuracy for Ratios of Modified Bessel Functions. Journal of Mathematical Analysis and Applications, 526, Article 127211.
https://doi.org/10.1016/j.jmaa.2023.127211
[8] Martin, P., Castro, E., Paz, J.L. and De Freitas, A. (2009) Multi-Point Quasi-Rational Approximants in Quantum Chemistry, In: Paz, J.L. and Hernandez, A., Eds., New Developments in Quantum Chemistry, Transworld Research Network, 55-78.
[9] Martin P., Olivares J., Maass F. and Valero E. (2018) Analytic Approximations for Special Functions, Applied to the Modified Bessel Functions and. Results in Physics, 11, 1028-1033.
https://doi.org/10.1016/j.rinp.2018.10.029
[10] Martin, P., Olivares, J. and Maass, F. (2018) High Accuracy Analytic Approximations for Modified Bessel Functions. International Research Conference Proceedings, New York, April 2018, 1307-6892.
[11] Martin P., Rojas E., Olivares J. and Sotomayor A. (2021) Quasi-Rational Analytic Approximation for the Modified Bessel Function with High Accuracy, Symmetry, 13, Article 741.
https://doi.org/10.3390/sym13050741
[12] Abramoritz, M. and Stegun, I.A. (1970) Handbook of Mathematical Functions, Ninth Printing. Dover Books, 355-379.
[13] Alyazidi A. (2021) On the Uniform and Simultaneous Approximations of Functions. Advances in Pure Mathematics, 11, 785-790.
https://doi.org/10.4236/apm.2021.1110052
[14] Namin, F.A. (2022) Analytical Solutions to Definite Integrals for Combinations of Legendre, Bessel and Trigonometric Functions Encountered in Propagation and Scattering Problems in Spherical Coordinates. Journal of Applied Mathematics and Physics, 10, 2690-2697.
https://doi.org/10.4236/jamp.2022.109179

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.