Poisson mean absolute deviation, closed form vs direct sum on 899 points in (0, 24.3): max error 1.07e-50
series identity for log(peak_(m+1)/peak_m), m = 0..59: max error 6.99e-50; all negative: True
peak_1 = 2 (3/2)^(3/2) e^(-3/2) in [0.819832557884, 0.819832557884] < 0.82: True; certified log(peak_(m+1)/peak_m) < 0 for m = 1..2000 (interval arithmetic): True
2 sqrt(l) e^-l at l = 0.32: [0.8215438532, ...], at l = 0.74: [0.8208578181, ...]; both > 0.82: True (log of it is concave, so > 0.82 on all of [0.32, 0.74])
s <= 0.42 iff q <= 1 - e^-0.42 = [0.3429531802, ...] (text: q <= 0.34): True
RESULT: ALL CHECKS PASSED
