

(* Cannot realize lambda_20 *)
<<cox.m
Print["\n\n\n             SALEM 20\n"]

DiStrat[salem20]

qpic = QEuler[salem20]
fpic = Mat[salem20]
p = FixSig[qpic,fpic,1][[1]]
qpic = QTwist[qpic,fpic,p]

r = GeomOb[qpic,fpic]
Print["Obstruction to untwisted f:  ",r]
Print["Verified... ",VerifyOb[qpic,fpic,r]]

{a1,a2} = PrimeTwists[qpic,fpic,11]

DiGlue[QTwist[qpic,fpic,a1],fpic]
DiGlue[QTwist[qpic,fpic,a2],fpic]
