m := 4;
n := 4;
g := 13;
F<i>:=CyclotomicField(4);
R<u,v>:=PolynomialRing(F,2);
X := (-u^8 - 4*u^7 - 6*u^6 - 4*u^5 - u^4)*v^8 + (16*u^7 + 48*u^6 + 48*u^5 + 16*u^4)*v^7 + (i*u^11 + 2*i*u^10 + i*u^9 + 8*u^8 + 16*u^7 - 96*u^6 - 208*u^5 - 104*u^4 - i*u^3 - 2*i*u^2 - i*u)*v^6 + (-12*i*u^10 - 12*i*u^9 - 96*u^7 - 96*u^6 + 352*u^5 + 352*u^4 + 12*i*u^2 + 12*i*u)*v^5 + (-6*i*u^11 + 54*i*u^9 - 24*u^8 + 432*u^6 - 664*u^4 + 6*i*u^3 - 54*i*u)*v^4 + (48*i*u^10 - 48*i*u^9 - 64*i*u^8 + (64*i + 192)*u^7 + (-64*i - 192)*u^6 + (64*i - 704)*u^5 + (-64*i + 704)*u^4 + 64*i*u^3 - 112*i*u^2 + 112*i*u)*v^3 + (10*i*u^11 - 20*i*u^10 - 90*i*u^9 + (200*i + 32)*u^8 + (-196*i - 64)*u^7 + (192*i - 384)*u^6 + (-188*i + 832)*u^5 + (184*i - 416)*u^4 - 198*i*u^3 + 212*i*u^2 - 106*i*u)*v^2 + (-40*i*u^10 + 120*i*u^9 - 144*i*u^8 + (112*i - 128)*u^7 + (-96*i + 384)*u^6 + (96*i - 384)*u^5 + (-112*i + 128)*u^4 + 144*i*u^3 - 120*i*u^2 + 40*i*u)*v + u^12 + (-4*i - 4)*u^11 + (16*i + 10)*u^10 + (-20*i - 20)*u^9 + 15*u^8 + (24*i + 24)*u^7 + (-32*i - 52)*u^6 + (24*i + 24)*u^5 + 15*u^4 + (-20*i - 20)*u^3 + (16*i + 10)*u^2 + (-4*i - 4)*u + 1;
q := (2-u*v-v)/2;
t := -q/u;
E := [R!1,(q^2-1)*(t^2-1)/16,(q^2-1)*(t^2-1)/16,0,0];
P := [R!0,0];
Q := [(q+1)*(t^2-1)/8,(q+1)^2*(t-1)^2*(t+1)/32];