Bessel Function Calculator — J₀, J₁, J₂(x)
Calculate Bessel functions of the first kind J_n(x) for orders n = 0, 1 and 2 at any argument x. Widely used in wave propagation, heat conduction, quantum mechanics and signal processing.
Order n
First zero of J0: x ≈ 2.4048
How does this calculator work?
J_n(x) = Σ (−1)^m / [m!(m+n)!] · (x/2)^{2m+n}. Enter order n (0, 1 or 2) and argument x to compute the Bessel function of the first kind. First zeros: J₀ at 2.4048, J₁ at 3.8317, J₂ at 5.1356. The calculator plots all three curves from x = 0 to 15.
Formula
How this is calculated
Bessel functions of the first kind J_n(x) arise as solutions to Bessel's differential equation: x² y″ + x y′ + (x² − n²) y = 0. They appear wherever a problem has cylindrical symmetry — vibrating circular membranes, heat flow in cylinders, electromagnetic waves in waveguides, quantum mechanical wave functions and the diffraction of light.
The series expansion converges for all finite x and any non-negative integer order n. Each term alternates sign, so the function oscillates around zero with slowly decreasing amplitude as x grows — qualitatively similar to a damped sinusoid. The functions are not periodic: successive zeros are not equally spaced, though the gap approaches π as x → ∞. The first zeros are: J₀ at x ≈ 2.4048, J₁ at x ≈ 3.8317 and J₂ at x ≈ 5.1356.
This calculator uses the direct power series with 40 terms, which provides full double-precision accuracy (≈ 15 significant digits) for |x| ≤ 20 and n ≤ 2. For larger x or higher orders, asymptotic expansions or the Miller backward recurrence algorithm are preferred. The recurrence relation J_{n+1}(x) = (2n/x) J_n(x) − J_{n−1}(x) is shown as a cross-check: when applied forward from J₀ and J₁ it is numerically unstable for large x, but the values agree closely in the valid range.
Frequently asked questions
They are the natural eigenfunctions of problems with cylindrical or spherical geometry. Examples: the normal modes of a circular drum head use J₀ and J₁; the TE and TM modes in a circular waveguide use J_n and its derivative; heat conduction in a cylinder uses J₀; Fraunhofer diffraction through a circular aperture gives an Airy pattern described by J₁.
The series converges for all finite x, but for large x (say x > 20) the terms first grow then cancel — catastrophic cancellation in floating-point arithmetic reduces precision. For |x| ≤ 20 and n ≤ 2, 40 terms gives 15-digit accuracy. For larger x, the asymptotic expansion J_n(x) ≈ √(2/πx) · cos(x − nπ/2 − π/4) is more numerically stable.
The three-term recurrence J_{n+1}(x) = (2n/x) J_n(x) − J_{n−1}(x) is exact but numerically unstable when applied forward (increasing n) for large x, because small errors in J_n grow. The backward recurrence (Miller's method, starting from large n) is stable. For n ≤ 2 and |x| ≤ 20 the forward recurrence is accurate enough for a cross-check.
Also known as
TG we-Calculate Editorial Team. (2026). Bessel Function Calculator — J₀, J₁, J₂(x) [Online calculator]. TG we-Calculate. https://we-calculate.com/calculator/bessel-function-calculator
TG we-Calculate Editorial Team. "Bessel Function Calculator — J₀, J₁, J₂(x)." TG we-Calculate. 2026. https://we-calculate.com/calculator/bessel-function-calculator.
TG we-Calculate Editorial Team, "Bessel Function Calculator — J₀, J₁, J₂(x)," TG we-Calculate, 2026. [Online]. Available: https://we-calculate.com/calculator/bessel-function-calculator
@misc{wecalculate_bessel_function_calculator, title = {Bessel Function Calculator — J₀, J₁, J₂(x)}, author = {{TG we-Calculate Editorial Team}}, howpublished = {\url{https://we-calculate.com/calculator/bessel-function-calculator}}, year = {2026}, note = {TG we-Calculate} }
Did this calculator help you?