main(d,prc)=\ prc=46000;\ e=initell([0,0,0,-1/d^2,1/4/d^3]);\ level=globalred(e);\ e=chell(e,compo(level,2));\ e=initell(e);\ level=compo(level,1);\ print(e,factor(level));\ f = anell(e,prc);\ q=exp(-2*pi/sqrt(level));\ s=2*sum(0,n=1,prc,f[n]/n*q^n);\ s/e[15]; \\ p=phitate(koly)