—   order ($\nu\geq0$)

—   argument ($x\geq0$)

$\operatorname{bei}_\nu(x)=\Im\!\left(e^{\nu\pi i/2}I_\nu(xe^{\pi i/4})\right)$, where $I_\nu$ is the modified Bessel function of the first kind [1].

For an integer $n\geq0$, $\operatorname{bei}_n(x)=\left(\dfrac{x}{2}\right)^n\sum_{k=0}^{\infty}\dfrac{\sin(3n\pi/4+k\pi/2)}{k!(n+k)!}\left(\dfrac{x^2}{4}\right)^k$ [3].

$\operatorname{bei}_{\nu-1}(x)+\operatorname{bei}_{\nu+1}(x)=-\dfrac{\nu\sqrt{2}}{x}\bigl(\operatorname{bei}_\nu(x)+\operatorname{ber}_\nu(x)\bigr)$ for $x>0$ and $\nu\geq1$, where $\operatorname{ber}_\nu(x)$ is the companion Kelvin function [2].

For an integer $n\geq0$, $\operatorname{bei}_{-n}(x)=(-1)^n\operatorname{bei}_n(x)$ and $\operatorname{bei}_n(-x)=(-1)^n\operatorname{bei}_n(x)$ [1].

$J_\nu$ denotes the principal branch of the Bessel function of the first kind. The rotated argument has modulus $x$ and angle $3\pi/4$ [1].

The Kelvin function $\operatorname{ber}_\nu(x)$ is the real part of the same rotated Bessel value. The notation $\operatorname{bei}(x)$ without an order means $\operatorname{bei}_0(x)$ [1].

The continuous value at the origin is $\operatorname{bei}_\nu(0)=0$ for $\nu\geq0$, as follows from the Bessel power series [3].

At integer order the defining Bessel function is entire, so no branch choice remains. (4) extends the listed integer-order values to negative orders and negative real arguments.

It preserves the existing order-zero extremal heights and finite initial value where present. Order zero gains the sparse references $x=1/2,1,2,3,5,8$ and $x=\sqrt2,2\sqrt2,4\sqrt2$. The arguments $x=\sqrt2,2\sqrt2,4\sqrt2$ give equal integer-sized real and imaginary components in the rotated Bessel argument; this is a normalization-based reference, not a universal physical scale. Orders $1,2,3,1/2$ are sampled only at $x=1$, with order $1$ also at $x=2$, to cover adjacent integer-order relations and the elementary half-order reduction. These finite exploratory slices avoid a Cartesian grid and make no demand claim for each radius. [7] [8] [9]

The reference inputs and finite cutoffs are explicit choices rather than a completeness or application-demand claim. The attached balanced-selection.json specifies every argument exactly. Existing entries retain their identifiers, decimal strings and annotations.

For $\nu=0$, $c_1,\ldots,c_6$ are the unique zeros of the derivative in $(3863/1024,483/128)$, $(8479/1024,265/32)$, $(13047/1024,1631/128)$, $(8803/512,17607/1024)$, $(1385/64,22161/1024)$, $(26713/1024,13357/512)$, respectively. These labels identify arguments, not the values at them.

The retained pre-extension subfamily has the following argument definitions: These are intrinsic extremal heights, with the finite origin limit where it exists. The exact argument is the unique zero of the derivative in the explicitly listed rational interval. The intervals identify roots; their endpoints are not evaluation arguments. Six isolated critical arguments are retained per selected order. Numbering orders this finite selection, without claiming a global count of all preceding zeros. The order-zero Kelvin functions are the standard ber, bei, ker, kei functions. [5]. [6]

The attached generate.py forwards to balanced_generate.py, the standalone generator for the entire current selection. Its default offline check recomputes both new and retained values. The previous entry point is preserved byte-for-byte as pre-balanced-generator.py and in revision history. Earlier supporting sources and checks are preserved as provenance for their original subfamilies; balanced-provenance.json identifies their source hashes and the retained website revision.

Python

import mpmath as mp
mp.mp.dps = 150
left = mp.mpf('3863') / 1024
right = mp.mpf('483') / 128
argument = mp.findroot(lambda point: mp.diff(
    lambda value: mp.bei(0, value), point), (left, right))
print(mp.nstr(mp.bei(0, argument), 100))

DLMF 10.61: Definitions and Basic Properties

DLMF 10.63: Recurrence Relations and Derivatives

DLMF 10.65: Power Series

DLMF 10.9: Integral Representations

Argument selection: 10.61

Argument selection: 10.63

DLMF 10.61

DLMF 10.63

DLMF 10.73

T490 —   $\operatorname{ber}_\nu$ and $\operatorname{bei}_\nu$ are the real and imaginary parts of the same rotated Bessel value

Values_of_Bessel_functions_of_the_first_kind —   $\operatorname{bei}_\nu$ is the imaginary part of $J_\nu$ after rotation of its argument by $3\pi/4$

Values_of_modified_Bessel_functions_of_the_first_kind —   the rotated $I_\nu$ representation in (1) gives the same Kelvin function

Values_of_Bessel_functions_of_the_second_kind —   $Y_\nu$ is the second standard solution of the Bessel equation whose rotated solutions give the Kelvin functions

Values_of_modified_Bessel_functions_of_the_second_kind —   rotating $K_\nu$ gives the Kelvin pair $\operatorname{ker}_\nu,\operatorname{kei}_\nu$, while rotating $I_\nu$ gives $\operatorname{ber}_\nu,\operatorname{bei}_\nu$

Entries are of type: real number

digits: 100 (Unknown key)

How well the digits are known: proven

Table is complete: no (it preserves the existing order-zero extremal heights and finite initial value where present. Order zero gains the sparse references $x=1/2,1,2,3,5,8$ and $x=\sqrt2,2\sqrt2,4\sqrt2$. The arguments $x=\sqrt2,2\sqrt2,4\sqrt2$ give equal integer-sized real and imaginary components in the rotated Bessel argument; this is a normalization-based reference, not a universal physical scale. Orders $1,2,3,1/2$ are sampled only at $x=1$, with order $1$ also at $x=2$, to cover adjacent integer-order relations and the elementary half-order reduction. These finite exploratory slices avoid a Cartesian grid and make no demand claim for each radius)

How they were obtained:

Retained-entry provenance (unchanged numerical data): Arb ball arithmetic encloses the function and its argument uncertainty at 773 bits for 100 requested significant digits. Critical arguments are recomputed from their defining rational intervals, with opposite endpoint signs, a Gronwall bound proving the second derivative has a fixed sign throughout the interval, and a separately certified narrow bracket. Turning points use exact integers and enclosed square roots. Exact origin values remain exact integers. Every stored interval contains an independently recomputed 900-bit enclosure and agrees with mpmath at 160 and 210 decimal digits. Kelvin and Struve values additionally agree with numerical quadrature of elementary integral representations; these agreement checks do not establish the proven bounds. Coulomb checks include the unit Wronskian, the vanishing turning-point coefficient, and the free-case Bessel identities. Nonzero controls and a deliberately displaced value test the comparisons. See checks-evidence.json for per-entry measurements. Earlier grid-specific checks and sources remain historical evidence and do not certify this selection.

more

Balanced extension: balanced_generate.py recomputes the entire attached selection, including retained values. New values use Arb ball arithmetic with exact rational inputs or enclosed algebraic and hyperbolic inputs at 773 working bits for 100 significant digits. Exact origin anchors are stored as integers. The attached critical-argument source is preserved and hash-identified; its certified argument enclosures are propagated through the evaluation. All old numerical strings, higher precision and manual entry annotations are preserved. balanced-checks.json records full independent higher-precision Arb recomputation and independent mpmath evaluations at 160 and 210 digits for every new value, with finite-width controls and family-specific identities or integral checks. The generator does not claim authorship of retained data; their prior supporting attachments remain unchanged, and the previous entry point is preserved as pre-balanced-generator.py.