\p38
setrand(1429412696);bnfinit(x^2+29051222508*x-12).clgp
setrand(1); bnfinit(x^8 + 12*x^6 + 30*x^4 + 24*x^2 + 4).reg
setrand(1); bnfinit(x^4 - 3*x^2 + 49).reg

nfinit(factor(polzagier(9,5))[2,1],2).disc
nfinit(Pol([1,0,42,112,728,3248,14224,3392,289478,-804944,2966908,-11015200,17342836,-108601584,381107816,-1679988352,6252186465,-14812800240,28868620970,-27997506768,-33428758132,98285772160,-51592356424,-39975211584,55983352320,-24670808064,5337884160,-733917184,87744512])).disc

setrand(3);bnfinit(x^2-(1130481^2+4)).clgp
setrand(2);bnfinit(x^4 - x^3 + 63*x^2 - 22*x + 1004).clgp
setrand(1);bnfinit(x^8 - 8*x^6 + 38*x^4 - 143*x^2 + 121).clgp
bnfcertify(bnfinit(x^2-40!));
zetakinit(2*x-1);

nf=nfinit(y^5-4*y^3+2*y+11);
v = [4/3, -1, y^2+y+1, [1,2,3,4,5]~];
for (i=1, #v, print( nfelttrace(nf,v[i]) ))
for (i=1, #v, print( nfeltnorm(nf,v[i]) ))

funs = [nfeltadd, nfeltdiv, nfeltdiveuc, nfeltdivrem, nfeltmod, nfeltmul];

try(f) = for (i=1, #v, for (j=1,#v, print( f(nf, v[i],v[j])) ))
for (i = 1, #funs, try(funs[i]))

nfisincl(nfinit(x-1),y)
