hyperellratpoints(x*(x+1)*(2*x+1)*6,[10^7,1]) sum(k=1,24,k^2) V=matkerint(Mat(concat([2..24],70))); M=matrix(23,23,i,j,V[,i]~*V[,j]-2*V[24,i]*V[24,j]); Q=qfminim(M);Q[1..2] G=qfauto(M); G[1] #G[2] #qforbits(G,Q) [mf, F, v] = mffromqf(2*M); mfparams(mf) mfcoefs(F,10) L1=lfunqf(matdiagonal([1,1,1])); L2=lfunqf(matdiagonal([4,1,1])); L3=lfunqf(matdiagonal([4,4,1])); 3*lfun(L1,1/2)+6*lfun(L2,1/2)-12*lfun(L3,1/2)