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
)
,
{\displaystyle J_{\nu }\left(xe^{\frac {3\pi i}{4}}\right),\,}
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
)
,
{\displaystyle e^{-\nu \pi i/2}K_{\nu }\left(xe^{\frac {\pi i}{4}}\right),\,}
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 .
kei(x )
kei(x ) for x between 0 and 14.
k
e
i
(
x
)
e
x
/
2
{\displaystyle \mathrm {kei} (x)e^{x/{\sqrt {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
.
{\displaystyle {\begin{aligned}&\mathrm {kei} _{n}(x)=-\ln \left({\frac {x}{2}}\right)\mathrm {bei} _{n}(x)-{\frac {\pi }{4}}\mathrm {ber} _{n}(x)\\&-{\frac {1}{2}}\left({\frac {x}{2}}\right)^{-n}\sum _{k=0}^{n-1}\sin \left[\left({\frac {3n}{4}}+{\frac {k}{2}}\right)\pi \right]{\frac {(n-k-1)!}{k!}}\left({\frac {x^{2}}{4}}\right)^{k}\\&+{\frac {1}{2}}\left({\frac {x}{2}}\right)^{n}\sum _{k\geq 0}\sin \left[\left({\frac {3n}{4}}+{\frac {k}{2}}\right)\pi \right]{\frac {\psi (k+1)+\psi (n+k+1)}{k!(n+k)!}}\left({\frac {x^{2}}{4}}\right)^{k}.\end{aligned}}}
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
{\displaystyle \mathrm {kei} (x)=-\ln \left({\frac {x}{2}}\right)\mathrm {bei} (x)-{\frac {\pi }{4}}\mathrm {ber} (x)+\sum _{k\geq 0}(-1)^{k}{\frac {\psi (2k+2)}{[(2k+1)!]^{2}}}\left({\frac {x^{2}}{4}}\right)^{2k+1}}
and the asymptotic series
k
e
i
(
x
)
∼
−
π
2
x
e
−
x
2
[
f
2
(
x
)
sin
β
+
g
2
(
x
)
cos
β
]
,
{\displaystyle \mathrm {kei} (x)\sim -{\sqrt {\frac {\pi }{2x}}}e^{-{\frac {x}{\sqrt {2}}}}[f_{2}(x)\sin \beta +g_{2}(x)\cos \beta ],}
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: