L=[8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 40, 41, 44];for(i=1,13,k=zetakinit(x^2-L[i]);print("k=zetakinit(x^2-",L[i],")");for(j=16,18,a=round(zetak(k,1-2*j)*32!)/32!;print("z(1-",2*j,")=",a, " ",factor(denominator(a)))))