# Tail t >= 300: analytic part on s in [0,1/300]; non-analytic terms checked at t=300..4800 (each is c t^m e^{-bt}, m/b <= 36).
import sys; sys.path.insert(0,'.')
import engine as E
from engine import A, iv
S = iv(A(0), A(1)/300)
P1 = E.P1_eval(S); M1 = E.Mcomb_eval(S); rz = 2*E.BETA*(1+6*S)/E.GAMMA1
an = P1 + M1 + iv(-rz.upper(), rz.upper())
print("analytic part on s in [0,1/300]:", an.lower(), an.upper())
for t in [300, 600, 1200, 2400, 4800]:
    o = E.cell(t, t + 1e-6)
    rest = sum((v for k, v in o.items() if k not in ('U1_P1','M1_main','U6_zeros','TOTAL')), A(0))
    print(t, "sum of non-analytic terms:", float(rest.lower()), float(rest.upper()))