from mpmath import mp, besseljzero, gamma, pi, mpf
mp.dps = 120
n = mpf(21)
lam = besseljzero(n / 2, 1) ** 2
kappa = pi ** (n / 2) / gamma(1 + n / 2)
2 * (1 + n / 2) ** (1 + 2 / n) / (n * lam * kappa ** (2 / n))The generator isolates $j_{n/2,1}$ by sign changes using Sage's arb-backed RealBallField. For integer orders it evaluates the Bessel function by a rigorously tailed alternating power series and recurrence; for half-integer orders it uses the spherical-Bessel recurrence from sine and cosine.
The one-dimensional value is checked against $27/(16\pi^2)$.