/* CFT.gp */

pol = y^3+4*y-6;
bnf = bnfinit(pol,1);
p = 737840520470346748252167770839;
dec = idealprimedec(bnf,p); pr1 = dec[1]; pr2 = dec[2];
modf = [pr1,1;pr2,1]; modoo = [0];
bnr = bnrinit(bnf,[modf,modoo],,3);
bnrinit(bnf,[modf,modoo]);
bnrclassfield(bnr)
pol2 = bnrclassfield(bnr,,2)
bnrdisc(bnr)
nfdisc(pol2)