allocatemem(2^33)
\p 3000
\r Hilbert.gp


l=3
l=7


ch="(...)/Theta-Q5/Inert/";
D=5;
PB3(x,y,z)=-1/2*(x^4+y^4+z^4+x^4*y^4+x^4*z^4+y^4*z^4)+x*y*z+x^5*y*z+x*y^5*z+x*y*z^5-x^2*y^2*z^2;


MM=read("actionG24");
M=read(Str(ch,"../classeSLDG24-",5));

[b1,b2]=[3+7*I,-5+2*I];
b3=polroots(PB3(b1,b2,x))[1];
Z=taufromthetNivD(b1,b2,M,MM,D)[1];

MCl=read(Str(ch,"PolQ5Theta-",l,"/classeSLDG24-D5G",l));

DD=read(Str(ch,"PolQ5Theta-",l,"/DenQ5Theta",l));
Den(xx,yy)=substvec(DD,[x,y],[xx,yy]);

L1=vector(l^2+1,i,read(Str(ch,"PolQ5Theta-",l,"/NumPhi/numphi-",i)));
L2=vector(l^2+1,i,read(Str(ch,"PolQ5Theta-",l,"/NumPsi2/numpsi2-",i)));
L3=vector(l^2+1,i,read(Str(ch,"PolQ5Theta-",l,"/NumPsi3/numpsi3-",i)));

Phi=X^(l^2+1)+sum(i=1,l^2+1,sum(k=1,5,substvec(L1[i][k],[x,y],[b1,b2])*b3^(k-1))*X^(i-1)/Den(b1,b2));
Phip=Phi';
Psi2=sum(i=1,l^2+1,sum(k=1,5,substvec(L2[i][k],[x,y],[b1,b2])*b3^(k-1))/Den(b1,b2)*X^(i-1))/Phip;
Psi3=sum(i=1,l^2+1,sum(k=1,5,substvec(L3[i][k],[x,y],[b1,b2])*b3^(k-1))/Den(b1,b2)*X^(i-1))/Phip;


for(i=1,length(MCl),z=mtauD(MCl[i],Z,D);[b1p,b2p,b3p]=EvalBjp(TauToOmega([z[1]/l,z[2]/l],D));print(i"  "round(10^1000*subst(Phi,X,b1p))"  "round(10^1000*(b2p-subst(Psi2,X,b1p)))"  "round(10^1000*(b3p-subst(Psi3,X,b1p)))));