M2 = I_0(5) upper bound = 27.239871823604446
[right side x=60, y:0->5] certified steps: 31
[top side y=5, x:60->0] certified steps: 256
lifted quadrant index at last interior sample Q=39; endpoint argument = 2*pi*10
continuous change of arg g on right+top sides = 19 * pi
=> number of zeros of J0-1 in B minus the unit quarter disc: N = 9
z_1 ~ 6.947506242533062 + 1.880440653100628 i   |g(zh)| <= 2.74e-79, |g'(zh)| >= 0.942216, Rouche disc radius 2^-150: CERTIFIED
z_2 ~ 13.276697958613063 + 2.205702959317610 i   |g(zh)| <= 2.09e-79, |g'(zh)| >= 0.971538, Rouche disc radius 2^-150: CERTIFIED
z_3 ~ 19.579253925689088 + 2.401190563939702 i   |g(zh)| <= 1.14e-79, |g'(zh)| >= 0.981361, Rouche disc radius 2^-150: CERTIFIED
z_4 ~ 25.873373935357950 + 2.541420474435272 i   |g(zh)| <= 5.43e-80, |g'(zh)| >= 0.986207, Rouche disc radius 2^-150: CERTIFIED
z_5 ~ 32.163666085292618 + 2.650840397868881 i   |g(zh)| <= 2.22e-79, |g'(zh)| >= 0.989077, Rouche disc radius 2^-150: CERTIFIED
z_6 ~ 38.451875827821837 + 2.740576973117263 i   |g(zh)| <= 1.89e-79, |g'(zh)| >= 0.990968, Rouche disc radius 2^-150: CERTIFIED
z_7 ~ 44.738818789431924 + 2.816641954352064 i   |g(zh)| <= 1.06e-78, |g'(zh)| >= 0.992307, Rouche disc radius 2^-150: CERTIFIED
z_8 ~ 51.024930018195782 + 2.882654594138808 i   |g(zh)| <= 4.67e-76, |g'(zh)| >= 0.993303, Rouche disc radius 2^-150: CERTIFIED
z_9 ~ 57.310463905466626 + 2.940963696531522 i   |g(zh)| <= 2.19e-73, |g'(zh)| >= 0.994072, Rouche disc radius 2^-150: CERTIFIED
9 pairwise disjoint certified discs, each with exactly one (simple) zero; count matches N = 9: these are all zeros of J0-1 in B \ {0}.
arg z_1 > arg z_2: Im(z_1 conj z_2) = 9.641908 > 1.7e-44
arg z_2 > arg z_3: Im(z_2 conj z_3) = 11.306136 > 2.6e-44
arg z_3 > arg z_4: Im(z_3 conj z_4) = 12.367785 > 3.5e-44
arg z_4 > arg z_5: Im(z_4 conj z_5) = 13.155215 > 4.4e-44
arg z_5 > arg z_6: Im(z_5 conj z_6) = 13.782783 > 5.3e-44
arg z_6 > arg z_7: Im(z_6 conj z_7) = 14.305010 > 6.2e-44
arg z_7 > arg z_8: Im(z_7 conj z_8) = 14.752397 > 7.1e-44
arg z_8 > arg z_9: Im(z_8 conj z_9) = 15.143805 > 8.0e-44
arg z_9 > sup arg D_10: Im(z_9 conj b*) = 10.784692 > 0  (b* = 62.8159 + 3.0353 i)
elapsed 1.3s
