BesselK0

BesselK0 or \(K_0\) is holomorphic function, solution \(f\) of equation
 * \( \!\!\!\!\!\!\!\!\!\!\ (1) ~ ~ ~ \displaystyle

f''(z)+f(z)/z=f(z) \) that has asymptotic
 * \( \!\!\!\!\!\!\!\!\!\!\ (2) ~ ~ ~ \displaystyle

K_0(z) = \exp(-z)\sqrt{\frac{\pi}{2z}} ~ \Big( 1+ O(1/z)\Big) \) at large values of \(\Re(z)\). The same equation (1) holds also for \(f(z)=\mathrm{BesselY0}(\mathrm i z)\), where BesselY0 is Neumann function.

Complex map of \(f=\mathrm{BesselK0}(x+\mathrm i y)\) is shown in upper figure with lines \(u\!=\!\Re(f)\!=\!\rm const\) and lines \(v\!=\!\Im(f)\!=\!\rm const\) in the \(x\),\(y\) plane.

Integral representation

 * \( \!\!\!\!\!\!\!\!\!\!\ (o) ~ ~ ~ \displaystyle

K_0(z) = \int_0^\infty \exp\!\Big( -z \cosh(t) \Big) ~ \mathrm d t \)

Relation to other cylindric functions

 * \( \!\!\!\!\!\!\!\!\!\!\ (o) ~ ~ ~ \displaystyle

K_0(z)=\frac{\mathrm i \pi}{2} H_0(\mathrm i z) \) where \(H_0\) is the Hankel function of zero order, id est, BeesselH0.

Behavior along the real axis
Esplicit plot of BesselK0 is shown in figure at right, \(y=\mathrm{BesselK0}(x)=K_0(x)\). For comparison, the two asymptotics are shown,
 * \( \!\!\!\!\!\!\!\!\!\!\ (\mathrm A1) ~ ~ ~ \displaystyle y= \exp(-z)\sqrt{\frac{\pi}{2z}} ~ ~ ~ \) and
 * \( \!\!\!\!\!\!\!\!\!\!\ (\mathrm A2) ~ ~ ~ \displaystyle y= \ln \frac{\pi}{2 x} - \mathrm{EulerGamma} ~ ~ ~ \)

Here EulerGamma is the Euler constant,
 * \( \!\!\!\!\!\!\!\!\!\!\ (4) ~ ~ ~ \displaystyle

\mathrm{EulerGamma}=\lim_{n \rightarrow \infty} \left( \sum_{k=1}^{n} \frac{1}{k} - \ln(n) \right) \approx 0.577215664901532861 \) The first of them approximates the function at large values of the argument, and the second does for the small values of the argument.

Expansions
At small values of the argument, BesselK0 expands as follows:
 * \( \!\!\!\!\!\!\!\!\!\!\ (5) ~ ~ ~ \displaystyle

K_0(z) = \ln(2/z)  - \mathrm{EulerGamma} + \frac{1}{4}\left( \ln(2/z) + 1 - \mathrm{EulerGamma} \right)~ z^2 + ... \frac{1}{128}\left( \ln(2/z) + 3 - 2\mathrm{EulerGamma} \right)~ z^4 + ... \)

At large values of the argument BesselK0 expands as follows:
 * \( \!\!\!\!\!\!\!\!\!\!\ (6) ~ ~ ~ \displaystyle

K_0(z) = \exp(-z) \sqrt{\frac{\pi}{2z}} \left( 1-\frac{1}{8z} +\frac{9}{128 z^2} -\frac{75}{1024 z^3}+ .. \right)\) The Mathematica easy calculates some tens of terms in the expansions above. The leading terms of these expansions are shown in the explicit plot at right. These expansions are used in the numerical implementation.