R = 4.5000000000000000000000000000000000000000000000000000000000000000000000000000000000000000,  M_0 <= [348.41120457649094528661630450571098660276359846610104712098000000000000000000000000000000 +/- 2.18e-88]
  p=2 term            <= [3.1191623125197538923775445465617945563397506046211486435430600427202712988636262202263850 +/- 2.58e-89]
  p=3 term (cancelled)<= [1.6479184330021645370928678553837885569712358367341241776020415004562414398279134503104236 +/- 3.39e-89]
  p>=5 sum (to 1000000) <= [26.443102927049308558809611161863886646746386979609063293693664904357686900320112185290135 +/- 8.62e-89]
  p>1000000 tail        <= [0.00073336777261923156815334346204125616965752714212274557388189838658736906507417368663009117 +/- 5.40e-94]
  bracket bound B     <= [31.210917040345469199162491410810570954548452184535397487625707396727493323544188942669025 +/- 3.30e-88]

M_1 <= [10874.233201964654502103163152649953469825451610683212985790292437627046646516885106790819 +/- 7.31e-86]

HEAD (i <= 16, s1_head.log) <= [453.08379944484570437816698792067755090294197229672161723333000000000000000000000000000000 +/- 2.23e-88]
TAIL (i >  16, Cauchy)      <= 33.110710144042968750000000000000000000000000000000000000000000000000000000000000000000000
S_1 = HEAD + TAIL              <= [486.19450958888912787551787438479650842247322229672161723333000000000000000000000000000000 +/- 1.12e-88]