Understanding the Bessel function calculator — concept and principle
This calculator evaluates the Bessel functions of the first kind Jν(x) and the second kind Yν(x) for an order ν from 0 to 100, whole or fractional, and an argument x. Bessel functions appear in vibrating drumheads, heat flow in cylinders, electromagnetic waves in round waveguides and signal processing.
How it is calculated
Bessel's equation: x² w″ + x w′ + (x² − ν²) w = 0
Series: Jν(x) = (x/2)^ν Σ (−x²/4)^k ÷ (k! Γ(ν + k + 1))
Integer order: Jₙ(x) = (1/π) ∫₀^π cos(x sin θ − nθ) dθ
Yν(x) = (1/π) ∫₀^π sin(x sin θ − νθ) dθ − (1/π) ∫₀^∞ (e^(νt) + e^(−νt) cos νπ) e^(−x sinh t) dt, for x > 0
- The page uses the series for small x and the integral forms otherwise, with the three-term recurrence for high orders, and shows 6 significant digits.
- Yν(x) → −∞ as x → 0, so it is shown as — for x ≤ 0. For whole n, Jₙ(−x) = (−1)ⁿ Jₙ(x).
The principle
Bessel's equation appears whenever a wave or diffusion problem is written in cylindrical coordinates. Jν is the solution that stays finite at the centre, and Yν is a second, independent solution that blows up there. For large x both behave like slowly decaying cosine waves with amplitude about √(2/(πx)), which is why the ripples of a drum or a pond fade outward. Neighbouring orders are linked by the recurrence J_(ν−1)(x) + J_(ν+1)(x) = (2ν/x) Jν(x).
Worked example
Default inputs: ν = 2, x = 5.
- Known values: J₀(5) = −0.177597 and J₁(5) = −0.327579.
- Recurrence: J₂(5) = (2 × 1 ÷ 5) J₁(5) − J₀(5) = 0.4 × (−0.327579) + 0.177597 = 0.0465651.
- Likewise Y₀(5) = −0.308518 and Y₁(5) = 0.147863, so Y₂(5) = 0.4 × 0.147863 + 0.308518 = 0.367663.
- Order 2, argument 5: an integer order (written Jₙ, Yₙ).
| Order ν (0 to 100, may be fractional) | Jν(x) — Bessel function of the first kind |
|---|---|
| 0.5 | -0.342168 |
| 1 | -0.327579 |
| 1.5 | -0.169651 |
| 2 | 0.0465651 |
| 2.5 | 0.240377 |
| 3 | 0.364831 |
| 4 | 0.391232 |
How to read the result
Jν(x) oscillates and changes sign; its zeros give, for example, the vibration frequencies of a circular membrane. Small values like J₂(5) ≈ 0.047 simply mean x = 5 is close to a zero of J₂, the first of which is near x = 5.14. For large orders and small x, Jν is extremely small and Yν extremely large and negative.
Common mistakes
- Confusing Yν (also written Nν, Neumann function) with the modified functions Iν and Kν.
- Expecting Yν to exist at x = 0.
- Mixing up order and argument.
- Using too few digits when values are near a zero.
Frequently asked questions
What is a Bessel function used for?
Solving wave, heat and diffusion problems in round geometries, such as drums, pipes and antennas.
What is the difference between J and Y?
J is finite at x = 0; Y is a second independent solution that tends to −∞ there.
Can the order be fractional?
Yes. Half-integer orders give spherical Bessel functions, which can be written with sine and cosine.
Sources: NIST Digital Library of Mathematical Functions (DLMF), §10.2 Definitions and §10.9 Integral Representations