
(* Realize lambda_18 *)
                Print["\n\n             SALEM, DEGREE 18 \n"]

<<cox.m
prime = 13

(*--------------------------Salem part---------------*)
salem = salem18
qs0 = QEuler[salem]
fs = Mat[salem]
polys = (3*x - x^3)*(2 - 4*x^2 + x^4)*(2 - x - 4*x^2 + x^4)
qs = QTwist[qs0,fs,polys]

(*--------------------------Cyclotomic part---------------*)
cyclo = Cyclotomic[12,x]
real  = CxToReal[cyclo]
qc0 = QEuler[cyclo]
fc  = Mat[cyclo]
(* Factors of 13 in Z[x]/real:  4 +/- x *)
polyc = 4+x
qc = QTwist[qc0,fc,polyc]

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

DiK3[q,f,polys]

DiK3Salem[salem]
