Modified Bessel function of the first kind, order 0: $I_0(x)$.

Composite implementation via the Abramowitz & Stegun polynomial approximation (Table 9.8.1):

• |x| <= 3.75 — polynomial series in $(x / 3.75)^2$.
• |x| > 3.75 — polynomial series in $3.75 / |x|$, multiplied by $e^{|x|} / \sqrt{|x|}$ to recover the exponential growth.

Differentiable through both branches (engine primitives only).