Kelvin functions
Plot of the Kelvin function ber(z) in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1 function ComplexPlot3D
In applied mathematics, the Kelvin functions berν (x ) and beiν (x ) are the real and imaginary parts , respectively, of
J ν ( x e 3 π i 4 ) ,
where x is real, and Jν (z ) , is the ν th order Bessel function of the first kind. Similarly, the functions kerν (x ) and keiν (x ) are the real and imaginary parts, respectively, of
e − ν π i / 2 K ν ( x e π i 4 ) ,
where Kν (z ) is the ν th order modified Bessel function of the second kind.
These functions are named after William Thomson, 1st Baron Kelvin .
While the Kelvin functions are defined as the real and imaginary parts of Bessel functions with x taken to be real, the functions can be analytically continued for complex arguments xe iφ , 0 ≤ φ < 2π . With the exception of bern (x ) and bein (x ) for integral n , the Kelvin functions have a branch point at x = 0.
Below, Γ(z ) is the gamma function and ψ (z ) is the digamma function .
ber(x )
ber(x ) for x between 0 and 20.
b e r ( x ) / e x / 2 for x between 0 and 50.
For integers n , bern (x ) has the series expansion
b e r n ( x ) = ( x 2 ) n ∑ k ≥ 0 cos [ ( 3 n 4 + k 2 ) π ] k ! Γ ( n + k + 1 ) ( x 2 4 ) k ,
where Γ(z ) is the gamma function . The special case ber0 (x ), commonly denoted as just ber(x ), has the series expansion
b e r ( x ) = 1 + ∑ k ≥ 1 ( − 1 ) k [ ( 2 k ) ! ] 2 ( x 2 ) 4 k
and asymptotic series
b e r ( x ) ∼ e x 2 2 π x ( f 1 ( x ) cos α + g 1 ( x ) sin α ) − k e i ( x ) π ,
where
α = x 2 − π 8 ,
f 1 ( x ) = 1 + ∑ k ≥ 1 cos ( k π / 4 ) k ! ( 8 x ) k ∏ l = 1 k ( 2 l − 1 ) 2
g 1 ( x ) = ∑ k ≥ 1 sin ( k π / 4 ) k ! ( 8 x ) k ∏ l = 1 k ( 2 l − 1 ) 2 .
Plot of the Kelvin function bei(z) in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1 function ComplexPlot3D
bei(x )
bei(x ) for x between 0 and 20.
b e i ( x ) / e x / 2 for x between 0 and 50.
For integers n , bein (x ) has the series expansion
b e i n ( x ) = ( x 2 ) n ∑ k ≥ 0 sin [ ( 3 n 4 + k 2 ) π ] k ! Γ ( n + k + 1 ) ( x 2 4 ) k .
The special case bei0 (x ), commonly denoted as just bei(x ), has the series expansion
Plot of the Kelvin function ker(z) in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1 function ComplexPlot3D b e i ( x ) = ∑ k ≥ 0 ( − 1 ) k [ ( 2 k + 1 ) ! ] 2 ( x 2 ) 4 k + 2
and asymptotic series
b e i ( x ) ∼ e x 2 2 π x [ f 1 ( x ) sin α − g 1 ( x ) cos α ] − k e r ( x ) π ,
where α, f 1 ( x ) , and g 1 ( x ) are defined as for ber(x ).
ker(x )
ker(x ) for x between 0 and 14.
k e r ( x ) e x / 2 for x between 0 and 50.
For integers n , kern (x ) has the (complicated) series expansion
k e r n ( x ) = − ln ( x 2 ) b e r n ( x ) + π 4 b e i n ( x ) + 1 2 ( x 2 ) − n ∑ k = 0 n − 1 cos [ ( 3 n 4 + k 2 ) π ] ( n − k − 1 ) ! k ! ( x 2 4 ) k + 1 2 ( x 2 ) n ∑ k ≥ 0 cos [ ( 3 n 4 + k 2 ) π ] ψ ( k + 1 ) + ψ ( n + k + 1 ) k ! ( n + k ) ! ( x 2 4 ) k . Plot of the Kelvin function kei(z) in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1 function ComplexPlot3D
The special case ker0 (x ), commonly denoted as just ker(x ), has the series expansion
k e r ( x ) = − ln ( x 2 ) b e r ( x ) + π 4 b e i ( x ) + ∑ k ≥ 0 ( − 1 ) k ψ ( 2 k + 1 ) [ ( 2 k ) ! ] 2 ( x 2 4 ) 2 k
and the asymptotic series
k e r ( x ) ∼ π 2 x e − x 2 [ f 2 ( x ) cos β + g 2 ( x ) sin β ] ,
where
β = x 2 + π 8 ,
f 2 ( x ) = 1 + ∑ k ≥ 1 ( − 1 ) k cos ( k π / 4 ) k ! ( 8 x ) k ∏ l = 1 k ( 2 l − 1 ) 2
g 2 ( x ) = ∑ k ≥ 1 ( − 1 ) k sin ( k π / 4 ) k ! ( 8 x ) k ∏ l = 1 k ( 2 l − 1 ) 2 .
kei(x )
kei(x ) for x between 0 and 14.
k e i ( x ) e x / 2 for x between 0 and 50.
For integer n , kein (x ) has the series expansion
k e i n ( x ) = − ln ( x 2 ) b e i n ( x ) − π 4 b e r n ( x ) − 1 2 ( x 2 ) − n ∑ k = 0 n − 1 sin [ ( 3 n 4 + k 2 ) π ] ( n − k − 1 ) ! k ! ( x 2 4 ) k + 1 2 ( x 2 ) n ∑ k ≥ 0 sin [ ( 3 n 4 + k 2 ) π ] ψ ( k + 1 ) + ψ ( n + k + 1 ) k ! ( n + k ) ! ( x 2 4 ) k .
The special case kei0 (x ), commonly denoted as just kei(x ), has the series expansion
k e i ( x ) = − ln ( x 2 ) b e i ( x ) − π 4 b e r ( x ) + ∑ k ≥ 0 ( − 1 ) k ψ ( 2 k + 2 ) [ ( 2 k + 1 ) ! ] 2 ( x 2 4 ) 2 k + 1
and the asymptotic series
k e i ( x ) ∼ − π 2 x e − x 2 [ f 2 ( x ) sin β + g 2 ( x ) cos β ] ,
where β , f 2 (x ), and g 2 (x ) are defined as for ker(x ).
References
Abramowitz, Milton ; Stegun, Irene Ann , eds. (1983) [June 1964]. "Chapter 9" . Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables . Applied Mathematics Series. Vol. 55 (Ninth reprint with additional corrections of tenth original printing with corrections (December 1972); first ed.). Washington D.C.; New York: United States Department of Commerce, National Bureau of Standards; Dover Publications. p. 379. ISBN 978-0-486-61272-0 . LCCN 64-60036 . MR 0167642 . LCCN 65-12253 .
Olver, F. W. J.; Maximon, L. C. (2010), "Bessel functions" , in Olver, Frank W. J. ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions , Cambridge University Press, ISBN 978-0-521-19225-5 , MR 2723248 .
External links
Weisstein, Eric W. "Kelvin Functions." From MathWorld—A Wolfram Web Resource.
GPL-licensed C/C++ source code for calculating Kelvin functions at codecogs.com: