Precise Special Functions

Precise Special Functions computes high-precision values of specialized mathematical functions: Bessel (J, Y, I, K), Riemann zeta, Gauss hypergeometric ₂F₁, complete elliptic integrals (K and E), and the Gamma function. Backed by arbitrary-precision arithmetic via mpmath — up to 1000 decimal digits. Suited for physics simulations, number theory, engineering calculations, and any context where a standard float or an inline code approximation is not precise enough. Every result is deterministic and cross-checkable against reference implementations. Ask naturally: "what is J₀(7.2) to 30 digits" and you get an exact value.