(1) max |J0_hankel - besselj| at test points: 8.0e-32
(2) m=10..3000: max (2 pi m)|z_m - F_m| = 0.202365; args strictly decreasing; min m^2 (arg z_m - arg z_(m+1)) = 0.429442
(3) t=10: finite difference -0.00387154384014  formula -0.00387154384014
(3) t=37.5: finite difference -0.00035429964331  formula -0.00035429964331
(3) t=1000: finite difference -7.62450451192e-7  formula -7.62450451192e-7
