
(* Realize lambda_6 *)

Print["\n\n             SALEM 6"]

<<cox.m

(*-----------------Find Salem action --------------*)
salem = salem6
qs = 3 QEuler[salem]
fs = Mat[salem]
polys = -(1+x)
qs = QTwist[qs,fs,polys]

(*-----------------Find twist of QA[12] with suitable signature--------------*)
qc = QA[12]
fc = Cox[qc]
polyc = 1+x-x^3
qc = QTwist[qc,fc,polyc]

{qd,fd} = GlueAuto[qs,fs,qc,fc,3]

(*-----------------Find twist of QA[2] with suitable signature--------------*)
q3 = QA[2]
f3 = Cox[q3]
q3 = 2 q3

(*-----------------Identity factor --------------*)
q1 = {{2,1},{1,20}}
f1 = Id[q1]
{qx,fx} = GlueAuto[q3,f3,q1,f1,3]

{q,f} = GlueAuto[qd,fd,qx,fx,2]

{q,f} = GlueSelf[q,f,13]

DiK3[q,f]
