Addition-subtraction chain

From formulasearchengine
Revision as of 17:57, 19 March 2012 by en>CRGreathouse (narrow cat)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

In applied mathematics, the Kelvin functions berν(x) and beiν(x) are the real and imaginary parts, respectively, of

Jν(xe3π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 Kν(xeπ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 x ei φ, φ ∈ [0, 2π). With the exception of Bern(x) and Bein(x) for integral n, the Kelvin functions have a branch point at x = 0.

ber(x)

ber(x) for x between 0 and 10.
ber(x)/ex/2 for x between 0 and 100.

For integers n, bern(x) has the series expansion

bern(x)=(x2)nk0cos[(3n4+k2)π]k!Γ(n+k+1)(x24)k

where Γ(z) is the Gamma function. The special case ber0(x), commonly denoted as just ber(x), has the series expansion

ber(x)=1+k1(1)k(x/2)4k[(2k)!]2

and asymptotic series

ber(x)ex22πx[f1(x)cosα+g1(x)sinα]kei(x)π,

where α=x/2π/8, and

f1(x)=1+k1cos(kπ/4)k!(8x)kl=1k(2l1)2
g1(x)=k1sin(kπ/4)k!(8x)kl=1k(2l1)2

bei(x)

bei(x) for x between 0 and 10.
bei(x)/ex/2 for x between 0 and 100.

For integers n, bein(x) has the series expansion

bein(x)=(x2)nk0sin[(3n4+k2)π]k!Γ(n+k+1)(x24)k

where Γ(z) is the Gamma function. The special case bei0(x), commonly denoted as just bei(x), has the series expansion

bei(x)=k0(1)k(x/2)4k+2[(2k+1)!]2

and asymptotic series

bei(x)ex22πx[f1(x)sinαg1(x)cosα]ker(x)π,

where α, f1(x), and g1(x) are defined as for ber(x).


ker(x)

For integers n, kern(x) has the (complicated) series expansion

kern(x)=12(x2)nk=0n1cos[(3n4+k2)π](nk1)!k!(x24)kln(x2)bern(x)+π4bein(x)+12(x2)nk0cos[(3n4+k2)π]ψ(k+1)+ψ(n+k+1)k!(n+k)!(x24)k
ker(x) for x between 0 and 10.
ker(x)ex/2 for x between 0 and 100.

where ψ(z) is the Digamma function. The special case ker 0(x), commonly denoted as just ker(x), has the series expansion

ker(x)=ln(x2)ber(x)+π4bei(x)+k0(1)kψ(2k+1)[(2k)!]2(x24)2k

and the asymptotic series

ker(x)π2xex2[f2(x)cosβ+g2(x)sinβ],

where β=x/2+π/8, and

f2(x)=1+k1(1)kcos(kπ/4)k!(8x)kl=1k(2l1)2
g2(x)=k1(1)ksin(kπ/4)k!(8x)kl=1k(2l1)2.


kei(x)

For integers n, kein(x) has the (complicated) series expansion

kein(x)=12(x2)nk=0n1sin[(3n4+k2)π](nk1)!k!(x24)kln(x2)bein(x)π4bern(x)+12(x2)nk0sin[(3n4+k2)π]ψ(k+1)+ψ(n+k+1)k!(n+k)!(x24)k
kei(x) for x between 0 and 10.
kei(x)ex/2 for x between 0 and 100.

where ψ(z) is the Digamma function. The special case kei 0(x), commonly denoted as just kei(x), has the series expansion

kei(x)=ln(x2)bei(x)π4ber(x)+k0(1)kψ(2k+2)[(2k+1)!]2(x24)2k+1

and the asymptotic series

kei(x)π2xex2[f2(x)sinβ+g2(x)cosβ],

where β, f2(x), and g2(x) are defined as for ker(x).


See also

References

External links

  • Weisstein, Eric W. "Kelvin Functions." From MathWorld—A Wolfram Web Resource. [1]
  • GPL-licensed C/C++ source code for calculating Kelvin functions at codecogs.com: [2]