
(* Realize lambda_10(2) *)

Print["\n\n             SALEM, DEGREE 10 (2)\n"]
<<cox.m

salem = SmallSalems10[[2]] 
prime=43

(*-----------------Salem part --------------*)
qs = QEuler[salem]
fs = Mat[salem]
ps = PrimeTwists[qs,fs,prime,1] 
ps = Select[ps,AbGp[QTwist[qs,fs,#],prime]=={prime,prime}&] 
polys = ps[[1]]
qs = QTwist[qs,fs,polys]
charfs = GlueChar[qs,fs,prime]

(*-----------------Cyclotomic = Transcendental part --------------*)
cyclo = Cyclotomic[42,x]
qc = QEuler[cyclo]
fc = Mat[cyclo]
ps = PrimeTwists[qc,fc,prime,2] 
ps = Select[ps,GlueChar[QTwist[qc,fc,#],fc,prime]==charfs&]
polyc = ps[[1]]
qc = QTwist[qc,fc,polyc]

(*-----------------Glue --------------*)
{q,f} = GlueAuto[qs,fs,qc,fc,prime]

DiK3[q,f,polys]
