Verification of the polynomials:

/*
Invariants for D=2
J1=G2^2/H4
J3=G2*H6/H4^2
 (J2=J1/J3=G2*H4/H6)

Invariants for D=5
J2=G2^5/H10
J4=H6*G2^2/H10
 (J1=J4/J2=H6/G2^3)
*/


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


l=3
l=5
l=11

*************************************************
ch="(...)/Gundlach-Q2/Inert/";


DD=read(Str(ch,"PolQ2Gundlach-",l,"/DenQ2Gund",l));
D(xx,yy)=substvec(DD,[x,y],[xx,yy]);

J1='J1;J2='J2;j1='j1;j2='j2;j3='j3;
M1=read("logJ1");
M2=read("logJ2");

[J1,J3]=[1+3*I,5+7*I]*1.;
Z=taufromgundNiv2(J1,J3);

L1=vector(l^2+1,i,read(Str(ch,"PolQ2Gundlach-",l,"/NumPhi/numphi",i)));
L2=vector(l^2+1,i,read(Str(ch,"PolQ2Gundlach-",l,"/NumPsi/numpsi",i)));

M=read(Str(ch,"PolQ2Gundlach-",l,"/classeGD2G0",l));

Phi=X^(l^2+1)+sum(i=1,l^2+1,substvec(L1[i],[x,y],[J1,J3])/D(J1,J3)^2*X^(i-1));
Phip=Phi';
Psi=sum(i=1,l^2+1,substvec(L2[i],[x,y],[J1,J3])/D(J1,J3)^4*X^(i-1))/Phip;

for(i=1,l^2+1,z=1/l*mtau2(M[i],Z);[J1p,J3p]=JGundlachNiv2([z[1],z[2]],M1,M2);print(i"  "round(10^1000*subst(Phi,X,J1p))"  "round(10^1000*(J3p-subst(Psi,X,J1p)))));