J0(0) = 1
At the center of a circular membrane or waveguide, the zeroth-order Bessel function of the first kind is unity. This is the only non-singular Bessel function value at the origin for the J family.
J0(0) = 1.000000
Mission control · cylindrical Bessel J
Bessel Function Calculator
Evaluate cylindrical Bessel functions Jv(x) and related Hankel forms in your browser, then read the full guide: what they are (#what-is-bessel-function), the governing equation (#bessel-function-formula), how to compute them (#how-to-calculate-bessel-functions), worked examples (#bessel-function-examples), and where they appear in physics and engineering (#applications). Open the on-page tool at #calculator when you are ready to try your own v and x.
Bessel function formulas
Start from the differential equation, then name the standard solutions. The calculator on this page evaluates Jv, Yv, and Hankel forms numerically; modified Iv and Kv are summarized below for reference.
Bessel's differential equation (cylindrical, order v):
x2 y'' + x y′ + (x2 - v2) y = 0
First kind: Jv(x) (finite at x = 0 for integer v)
Second kind: Yv(x) (singular at x = 0)
Modified equation (replace x2 with -x2):
Iv(x), Kv(x) (exponential growth/decay instead of oscillation)
For integer n >= 0, a series definition is Jn(x) = sumk=0^infty (-1)^k / (k!(n+k)!) * (x/2)(n+2k). Non-integer v uses Gamma functions in place of factorials. In practice, libraries switch between series, continued fractions, and asymptotic expansions so values stay accurate for large |x|.
Hankel functions used in scattering are Hv(1)(x) = Jv(x) + i Yv(x) and Hv(2)(x) = Jv(x) - i Yv(x). The on-page calculator reports both, along with derivatives.
How to calculate Bessel functions
Hand calculation is practical only for small |x| or the lowest orders. For design work you will use recurrence relations, published tables, or software. The steps below mirror what a robust library does conceptually.
- State the order and argument: Write v and x in the same scaling as your PDE (for example k r after separation). Check whether you need Jv, Yv, Hankel, or modified Iv/Kv.
- Use series for small |x|: For |x| below roughly 1, truncate the power series for Jn or Jv with Gamma coefficients. Watch rounding when many terms are needed.
- Apply recurrence for nearby orders: Relations such as Jv-1(x) + Jv+1(x) = (2v/x) Jv(x) step from a known value to neighboring orders once one Jv is available.
- Switch to asymptotic forms for large |x|: For large x, use trigonometric approximations (Airy-type envelopes with cos(x - pi/4 - v pi/2)) so you do not sum hundreds of series terms.
- Numerical libraries and CAS: Production code calls AMOS, Cephes, SciPy (jv, yv), MATLAB besselj/bessely, or Mathematica BesselJ. Match branch cuts when x is complex.
- Verify with this calculator: Spot-check Jv(x), Yv(x), derivatives, and Hankel values at your v and x before embedding numbers in a report or spreadsheet.
Bessel function examples
These worked values are typical checks against tables or software. Plug the same v and x into the calculator above to confirm.
J1(3.8317...) ≈ 0
The first positive zero of J0 is about 3.8317. At that x, J1(x) is also a commonly tabulated root because J0'(x) = -J1(x), so the derivative of J0 vanishes there.
J1(3.83170597) ≈ 0 (first zero of J0)
Y0(x) near x = 0
The second kind has a logarithmic singularity: Y0(x) ~ (2/pi) ln(x) as x -> 0+. Any model that must stay finite at the origin must use Jv, not Yv, in that region.
Y0(0.1) ≈ -1.45 (illustrative; use calculator for your x)
Waveguide cutoff (engineering)
A circular waveguide TE mode involves zeros of Jm. If radius a = 1 cm and the first J1 zero is 3.8317, a rough scale for the lowest cutoff frequency involves c times that zero divided by 2 pi a (with mu and epsilon of the fill).
f_c ~ (c / 2 pi a) * 3.8317 for the dominant TE family (order-of-magnitude check)
Bessel function of the first kind, Jv(x)
Jv(x) is the solution chosen to be finite at x = 0 for integer orders. Physically it often represents a standing radial mode inside a disk, fiber, or cavity.
For integer n, Jn is an even function when n is even and odd when n is odd. Zeros of Jn set resonant frequencies: each zero is a radial node in a circular drumhead or a cutoff condition in a round waveguide.
Derivatives appear in boundary conditions involving slope. Remember J0'(x) = -J1(x) and the general recurrence Jv-1(x) - Jv+1(x) = 2 Jv'(x) / x.
Bessel function of the second kind, Yv(x)
Yv(x) is the second linearly independent solution on the positive real axis. It blows up logarithmically at x = 0, so it is excluded from problems that include the origin.
In radiation and scattering, Yv appears inside Hankel functions Hv(1) and Hv(2), which represent outgoing and incoming cylindrical waves when combined with Jv.
Near zeros of Jv, Yv has poles in the quotient definition Yv = (Jv cos(v pi) - J-v) / sin(v pi) for non-integer v. Numerical libraries handle these limits stably.
Modified Bessel functions
When the Laplacian in cylindrical coordinates leads to imaginary separation constants, the ordinary Bessel equation becomes the modified Bessel equation with a minus sign on x2. Solutions grow or decay exponentially instead of oscillating.
Iv(x) is the analog of Jv and is finite at the origin for integer v. Kv(x) is the analog of Yv and decays for large x on the real axis.
Modified Bessel equation:
x2 y'' + x y′ - (x2 + v2) y = 0
Iv(x) : modified first kind
Kv(x) : modified second kind
Bessel function graphs
On the real axis, Jv(x) looks like a damped sinusoid for large x: amplitude falls roughly as 1/sqrt(x). Yv(x) has the same envelope but is phase-shifted and singular at zero.
The calculator plot at the top shows Jv(x) and Yv(x) for your chosen order on 0 <= x <= 10. Compare shapes as you change v: higher order pushes the first peak rightward.
Modified Iv and Kv are monotone (for real positive x and many orders): Iv grows while Kv decays, which matches heat diffusion in an infinite cylinder.
J0 and J1
J0 starts at 1 and oscillates with decreasing amplitude. J1 starts at 0, rises, then oscillates. Their zero crossings set mode shapes in round geometry.
Yv singularity
Plots must start at x > 0 for Yv. The logarithmic spike near the origin is why only Jv is used for fields that must stay bounded at r = 0.
Complex argument
For complex z, magnitude and phase plots replace a single xy curve. Use complex mode in the calculator table when z = x + iy is not on the real line.
Hankel as J + iY
Real and imaginary parts of Hv(1) combine Jv and Yv traces. In the far field they approximate traveling cylindrical waves.
How to use this Bessel function calculator
The tool at the top of the page (#calculator) evaluates cylindrical Bessel functions for order v between -99 and 99. Choose real x for a full table plus a Jv and Yv plot, or complex z = Re(x) + i Im(x) for tabulated complex values.
Results include Jv, J′v, Yv, Y′v, Hankel Hv(1), Hv(2), and their derivatives at your point. Nothing is sent to a server: arithmetic runs locally after the page loads.
Use it to sanity-check homework, confirm a MATLAB or Python call, or capture numbers for a lab report. For step-by-step theory, see the sections on formulas, examples, and applications below; this guide does not repeat the input widgets.
Common Bessel function mistakes
Most errors come from mixing families (J vs Y vs I vs K), using the wrong scaling for x, or ignoring singularities at the origin.
- Using Yv(x) at or very near x = 0 when the physics requires a finite field.
- Forgetting that x in the equation is dimensionless: compare k r, not r alone, to tabulated zeros.
- Confusing modified Iv with ordinary Jv: they solve different sign equations.
- Picking the wrong Hankel function for propagation direction (H(1) vs H(2)).
- Assuming Jv is always even: only integer orders have clean parity; half-integer v needs care.
- Evaluating non-integer v at negative x without specifying the branch cut (especially in complex mode).
- Reporting Yv where Jv zeros were intended (cutoff formulas use zeros of Jm, not Ym).
Bessel functions vs Fourier series
Fourier series expand a periodic function on a line segment: sines and cosines are eigenfunctions of d2/dx2 on an interval with periodic boundaries.
Bessel functions are the radial eigenfunctions on a disk or cylinder: Jm(k r) replaces sin(k x) when the geometry is round and the PDE is separated in polar coordinates.
You often use both in one problem: Fourier or Laplace modes in z or time, and Bessel modes in the transverse plane. If the domain is rectangular, use Fourier in x and y; if it is circular, use Bessel in r and cos(m theta) in angle.
Bessel functions in physics and engineering
Whenever a field lives on a circle, cylinder, or sphere, radial modes are almost always Bessel functions. The list below focuses on physics and engineering settings where Jv and Yv are central.
Electromagnetic waveguides
TE and TM modes in circular metal pipes satisfy boundary conditions on Jm and its derivatives at the wall radius.
Antennas and scattering
Cylindrical and wire antennas use Hankel functions (J + iY) to match radiation conditions at infinity.
Structural vibrations
Circular plates and membranes expand displacement in Jm(alpha r) with alpha fixed by the edge (clamped, simply supported, or free).
Acoustics and ultrasonics
Piston sources on a circular baffle and room modes in cylindrical tanks involve Bessel radiation integrals.
Heat and mass diffusion
Transient conduction in a disk or cylinder uses J0 and J1 in series solutions; steady semi-infinite problems often bring in modified Kv.
Quantum and plasma physics
Radial wavefunctions in cylindrical traps and certain plasma dispersion relations use Jv and Yv with complex arguments.
Frequently asked questions
What is the difference between Jv and Yv?
Both solve the same Bessel equation. Jv is finite at x = 0 for integer v and is used when the origin is part of the domain. Yv is singular at x = 0 and appears in radiation problems, often inside Hankel functions.
What is the difference between real and complex mode in the calculator?
Real mode evaluates at a point on the number line and draws Jv(t) and Yv(t) for t from 0 to 10. Complex mode sets z = Re(x) + i Im(x) and returns complex values in the table only.
What input ranges does the calculator allow?
Order v must be between -99 and 99. Real x must lie in [-20, 20]. In complex mode, both Re(x) and Im(x) must each be in [-20, 20]. Values outside these ranges show a clear error message.
How are Hankel functions related to J and Y?
Hv(1)(x) = Jv(x) + i Yv(x) and Hv(2)(x) = Jv(x) - i Yv(x). They package outgoing and incoming cylindrical waves in scattering and waveguide problems.
When should I use modified Bessel functions Iv and Kv?
Use them when the governing equation has a minus sign on x2 (evanescent fields, steady heat in semi-infinite media, some potentials). They grow or decay exponentially rather than oscillate like Jv.
Why does the plot start slightly above x = 0?
Yv has a singularity at x = 0. The graph samples from a small positive x so Yv stays finite on the canvas while still showing the Jv curve from the origin.
Do you store my inputs?
No. All calculations run locally in your browser after the page loads. Your v and x never leave your device unless you copy them elsewhere.
How do Bessel zeros relate to resonant frequencies?
Boundary conditions on a circular domain pick discrete zeros of Jm (or derivatives). Each zero is a radial eigenvalue; combined with angular order m and boundary conditions in z or time, it sets a resonant frequency or cutoff.
Can v be non-integer?
Yes. The calculator accepts non-integer v in the allowed range. Series use Gamma functions, and libraries handle half-integer and fractional orders with stable recurrences.
How accurate are the numbers compared to MATLAB or SciPy?
The page uses standard numerical libraries in the browser. For publication or safety-critical work, cross-check one point in your preferred CAS and note the precision shown in the table.