P1 = x^4-5; polgalois(P1) P2 = x^4-x^3-7*x^2+2*x+9; polgalois(P2) P3 = x^4-x^3-3*x^2+x-1; polgalois(P3) Q1 = nfsplitting(P1) Q2 = nfsplitting(P2) Q3 = nfsplitting(P3) Q3 = polredbest(Q3) gal = galoisinit(Q3); gal.gen ord = gal.orders prod(i=1,#ord,ord[i]) galoisidentify(gal) L = galoissubgroups(gal); #L R1 = galoisfixedfield(gal,L[25])[1]; polgalois(R1) R2 = galoisfixedfield(gal,L[28])[1]; polgalois(R2) nf = nfinit(Q3); factor(nf.disc) dec3 = idealprimedec(nf,3); pr3 = dec3[1]; [#dec3, pr3.f, pr3.e] ram3 = idealramgroups(nf,gal,pr3); #ram3 galoisidentify(ram3[1]) galoisisabelian(ram3[1]) galoisidentify(ram3[2]) galoisidentify(ram3[3]) dec11 = idealprimedec(nf,11); pr11 = dec11[1]; [#dec11, pr11.f, pr11.e] ram11 = idealramgroups(nf,gal,pr11); #ram11 galoisidentify(ram11[1]) galoisidentify(ram11[2]) dec2 = idealprimedec(nf,2); pr2 = dec2[1]; [#dec2, pr2.f, pr2.e] frob2 = idealfrobenius(nf,gal,pr2); permorder(frob2) N = 7*13*19; L1 = polsubcyclo(N,3); L2 = [P | P <- L1, #factor(nfinit(P).disc)[,1] == 3] G = znstar(N) H = mathnfmodid([1,0;-1,1;0,-1],3); pol = galoissubcyclo(G,H) factor(nfinit(pol).disc) bnf = bnfinit(a^2-a+50); bnf.cyc R = bnrclassfield(bnf)[1] [cond,bnr,subg] = rnfconductor(bnf,R); cond subg R2 = bnrclassfield(bnf,,2) bnr = bnrinit(bnf,12); bnr.cyc [deg,r1,D] = bnrdisc(bnr); [deg,r1] D bnrclassfield(bnr) bnrclassfield(bnr,,1) bnr = bnrinit(bnf,7); bnr.cyc bnrclassfield(bnr,3) pr41 = idealprimedec(bnf,41)[1]; bnrisprincipal(bnr,pr41,0) bnr = bnrinit(bnf,[102709,43512;0,1],1); bnr.cyc bnrclassfield(bnr,[9,3;0,1]) bnf=bnfinit(a^2-217); bnf.cyc bnrinit(bnf,1).cyc bnrinit(bnf,[1,[1,1]]).cyc quadhilbert(-31) lift(quadray(13,7)) bnf = bnfinit(x^2+2*3*5*7*11); bnf.cyc bnr = bnrinit(bnf,1,1); gal = galoisinit(bnf); m = bnrgaloismatrix(bnr,gal)[1]