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

ch="(...)/Theta-Q5/Split/";

D=5;

MM=read("actionG24");
M=read(Str(ch,"../classeSLDG24-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;

*********

l=5;
beta=3-(1+sqrt(5))/2;
betac=3-(1-sqrt(5))/2;

Den(x,y)=16*x^8*y^8 - 8*x^4*y^4 + x^4 + y^4;
Tab=vector(l+1,i,2);Tab[l+1]=1;

*********

l=11;
beta=4-(1+sqrt(5))/2;
betac=4-(1-sqrt(5))/2;

*********

l=19;
beta=5-(1+sqrt(5))/2;
betac=5-(1-sqrt(5))/2;

************************************

l=29;
beta=5+(1+sqrt(5))/2;
betac=5+(1-sqrt(5))/2;

************************************

l=31;
beta=5+2*(1+sqrt(5))/2;
betac=5+2*(1-sqrt(5))/2;

************************************

l=41;
beta=6+(1+sqrt(5))/2;
betac=6+(1-sqrt(5))/2;

************************************

l=59;
beta=7+2*(1+sqrt(5))/2;
betac=7+2*(1-sqrt(5))/2;

************************************

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

MCl=read(Str(ch,"PolQ5Theta-",l,"/classeSLDG24-D",D,"G",l));
Tab=vector(l+1,i,1);

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

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

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

for(i=1,length(MCl),z=mtauD(MCl[i],Z,D);[b1p,b2p,b3p]=EvalBjp(TauToOmega([z[1]/beta,z[2]/betac],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)))));