Bessel boundary equations solved directly

← Verification · Cylinders · Plane wave
Reference
Bessel boundary equations solved directly
Kind
analytic: a closed form, limit or conservation law: independent of any code
Calculated
coefficients, cross widths
Agreement
1e-12
A cylinder of n = 1.5 in air at 600 nm, normal incidence: the coefficients from the Bessel boundary equations Each order's coefficient solved directly from the two boundary conditions in SciPy's Bessel functions (Bohren & Huffman §8.4), summed into Q_sca = (2/x) Σ|t_m|², against the cylinder transfer matrices.
100200300400500600012345radius [nm]Qsca
E along the axis, PyStratify H along the axis, PyStratify E along the axis, Bessel equations H along the axis, Bessel equations
Hover to read values.
  • largest relative difference: 2.21×10⁻¹⁵
  • PyStratify the lines, the reference the circles

Run it yourself: unified_checks.py · test_unified.py · figures.py · UNIFIED_CHECKS.md on GitHub, PyStratify 97e52a9