
(*-----------------Show lambda_18(2) cannot be realized------------------*)
<<cox.m

s=SmallSalems18[[2]]

Print["\n	Salem 18 (2)\n"]

prime=5
qs = QEuler[s]
fs = Mat[s]
ps = PrimeTwists[qs,fs,prime]
ps = Select[ps,AbGp[QTwist[qs,fs,#]]=={prime,prime}&]
polys = ps[[1]]
qs = QTwist[qs,fs,polys]

DiGlue[qs,fs]

fixes = FixSig[qs,fs,1]
If[Length[fixes]==0,Print["Signature cannot be fixed"]]

qx = QTwist[qs,fs,polys]
DiGlue[qx,fs]
Print["Further twists have period too large:  ",GluePeriod[qx,fs]];
