L = [7, x^3 - 5*x^2 + 6*x - 1, 13, x^3 - 9*x^2 + 14*x - 1, 19, x^3 - 354315*x^2 + 39645147398*x - 1424861898171643];for(i=1,3,k=zetakinit(L[2*i]);print("k=zetakinit(",L[2*i],")");for(j=1,18,a=round(zetak(k,1-2*j)*32!)/32!;print("z(1-",2*j,")=",a, " ",factor(denominator(a)))))