m := 2;
n := 21;
g := 61;
F := Rationals();
R<u,v>:=PolynomialRing(F,2);
X := (u^16 + 2*u^15 + u^14)*v^30 + (14*u^17 + 34*u^16 - 2*u^15 - 71*u^14 - 84*u^13 - 56*u^12 - 28*u^11 - 8*u^10 - u^9)*v^29 + (91*u^18 + 262*u^17 - 128*u^16 - 1273*u^15 - 1856*u^14 - 1295*u^13 - 518*u^12 - 41*u^11 + 112*u^10 + 56*u^9 + 9*u^8)*v^28 + (364*u^19 + 1222*u^18 - 1030*u^17 - 10159*u^16 - 18276*u^15 - 14793*u^14 - 5033*u^13 + 1180*u^12 + 2675*u^11 + 1431*u^10 + 51*u^9 - 153*u^8 - 36*u^7)*v^27 + (1001*u^20 + 3874*u^19 - 4468*u^18 - 48529*u^17 - 105968*u^16 - 104375*u^15 - 38454*u^14 + 17089*u^13 + 28689*u^12 + 15846*u^11 + 359*u^10 - 4004*u^9 - 1096*u^8 + 168*u^7 + 84*u^6)*v^26 + (2002*u^21 + 8866*u^20 - 12770*u^19 - 156403*u^18 - 407096*u^17 - 490345*u^16 - 231519*u^15 + 96315*u^14 + 192905*u^13 + 110463*u^12 + 6494*u^11 - 37109*u^10 - 17841*u^9 + 3160*u^8 + 2696*u^7 + 84*u^6 - 126*u^5)*v^25 + (3003*u^22 + 15158*u^21 - 26048*u^20 - 361933*u^19 - 1105584*u^18 - 1608395*u^17 - 1003620*u^16 + 264787*u^15 + 877403*u^14 + 576453*u^13 + 72591*u^12 - 181570*u^11 - 136391*u^10 + 5220*u^9 + 35777*u^8 + 4290*u^7 - 3035*u^6 - 503*u^5 + 126*u^4)*v^24 + (3432*u^23 + 19734*u^22 - 39446*u^21 - 623051*u^20 - 2206556*u^19 - 3813525*u^18 - 3062229*u^17 + 254169*u^16 + 2778560*u^15 + 2321345*u^14 + 541613*u^13 - 532919*u^12 - 582618*u^11 - 99010*u^10 + 183102*u^9 + 75758*u^8 - 26772*u^7 - 11638*u^6 + 1283*u^5 + 709*u^4 - 84*u^3)*v^23 + (3003*u^24 + 19734*u^23 - 45210*u^22 - 813813*u^21 - 3312672*u^20 - 6696107*u^19 - 6666794*u^18 - 598727*u^17 + 6329897*u^16 + 7042995*u^15 + 2701656*u^14 - 832931*u^13 - 1572965*u^12 - 661170*u^11 + 375960*u^10 + 416004*u^9 - 33822*u^8 - 110493*u^7 - 2782*u^6 + 10410*u^5 + 769*u^4 - 542*u^3 + 36*u^2)*v^22 + (2002*u^25 + 15158*u^24 - 39446*u^23 - 813813*u^22 - 3788136*u^21 - 8840645*u^20 - 10545671*u^19 - 2521901*u^18 + 10808856*u^17 + 15945543*u^16 + 9025223*u^15 + 292559*u^14 - 2861448*u^13 - 2050200*u^12 - 87150*u^11 + 929392*u^10 + 269300*u^9 - 342496*u^8 - 130790*u^7 + 61463*u^6 + 18698*u^5 - 3450*u^4 - 1294*u^3 + 242*u^2 - 9*u)*v^21 + (1001*u^26 + 8866*u^25 - 26048*u^24 - 623051*u^23 - 3312672*u^22 - 8840645*u^21 - 12270188*u^20 - 4320415*u^19 + 14489083*u^18 + 27170877*u^17 + 20749128*u^16 + 5158635*u^15 - 3460695*u^14 - 4057456*u^13 - 2055500*u^12 + 359404*u^11 + 930875*u^10 - 216635*u^9 - 484292*u^8 + 46742*u^7 + 135490*u^6 - 1586*u^5 - 11844*u^4 - 939*u^3 + 699*u^2 - 59*u + 1)*v^20 + (364*u^27 + 3874*u^26 - 12770*u^25 - 361933*u^24 - 2206556*u^23 - 6696107*u^22 - 10545671*u^21 - 4320415*u^20 + 15884378*u^19 + 35313167*u^18 + 33770491*u^17 + 13919143*u^16 - 2320034*u^15 - 5936254*u^14 - 4706682*u^13 - 2299916*u^12 + 638292*u^11 + 842095*u^10 - 391663*u^9 - 240353*u^8 + 251422*u^7 + 104460*u^6 - 48217*u^5 - 13282*u^4 + 1627*u^3 + 1117*u^2 - 172*u + 6)*v^19 + (91*u^28 + 1222*u^27 - 4468*u^26 - 156403*u^25 - 1105584*u^24 - 3813525*u^23 - 6666794*u^22 - 2521901*u^21 + 14489083*u^20 + 35313167*u^19 + 39640258*u^18 + 21515373*u^17 + 59339*u^16 - 7164534*u^15 - 5678188*u^14 - 4392848*u^13 - 1378159*u^12 + 1682477*u^11 + 959700*u^10 - 338305*u^9 - 11561*u^8 + 233970*u^7 - 8246*u^6 - 67505*u^5 - 976*u^4 + 4094*u^3 + 1044*u^2 - 300*u + 15)*v^18 + (14*u^29 + 262*u^28 - 1030*u^27 - 48529*u^26 - 407096*u^25 - 1608395*u^24 - 3062229*u^23 - 598727*u^22 + 10808856*u^21 + 27170877*u^20 + 33770491*u^19 + 21515373*u^18 + 1340338*u^17 - 7720692*u^16 - 4419550*u^15 - 1992248*u^14 - 1995606*u^13 + 764263*u^12 + 2171667*u^11 + 304949*u^10 - 530548*u^9 + 16436*u^8 + 114847*u^7 - 61201*u^6 - 39583*u^5 + 11270*u^4 + 3215*u^3 + 566*u^2 - 355*u + 20)*v^17 + (u^30 + 34*u^29 - 128*u^28 - 10159*u^27 - 105968*u^26 - 490345*u^25 - 1003620*u^24 + 254169*u^23 + 6329897*u^22 + 15945543*u^21 + 20749128*u^20 + 13919143*u^19 + 59339*u^18 - 7720692*u^17 - 3411434*u^16 + 2533744*u^15 + 1233709*u^14 - 496789*u^13 + 1083430*u^12 + 750477*u^11 - 783131*u^10 - 511324*u^9 + 122668*u^8 + 78023*u^7 - 40853*u^6 - 9377*u^5 + 10459*u^4 + 334*u^3 + 270*u^2 - 299*u + 15)*v^16 + (2*u^30 - 2*u^29 - 1273*u^28 - 18276*u^27 - 104375*u^26 - 231519*u^25 + 264787*u^24 + 2778560*u^23 + 7042995*u^22 + 9025223*u^21 + 5158635*u^20 - 2320034*u^19 - 7164534*u^18 - 4419550*u^17 + 2533744*u^16 + 3728960*u^15 - 512841*u^14 - 1359681*u^13 + 144121*u^12 - 431018*u^11 - 817218*u^10 - 42299*u^9 + 239109*u^8 + 53952*u^7 - 23633*u^6 - 1223*u^5 + 3724*u^4 - 1046*u^3 + 341*u^2 - 174*u + 6)*v^15 + (u^30 - 71*u^29 - 1856*u^28 - 14793*u^27 - 38454*u^26 + 96315*u^25 + 877403*u^24 + 2321345*u^23 + 2701656*u^22 + 292559*u^21 - 3460695*u^20 - 5936254*u^19 - 5678188*u^18 - 1992248*u^17 + 1233709*u^16 - 512841*u^15 - 2530732*u^14 - 633717*u^13 + 662529*u^12 - 154974*u^11 - 227854*u^10 + 202117*u^9 + 147793*u^8 + 285*u^7 - 17348*u^6 - 38*u^5 + 137*u^4 - 857*u^3 + 380*u^2 - 62*u + 1)*v^14 - (84*u^29 + 1295*u^28 + 5033*u^27 - 17089*u^26 - 192905*u^25 - 576453*u^24 - 541613*u^23 + 832931*u^22 + 2861448*u^21 + 4057456*u^20 + 4706682*u^19 + 4392848*u^18 + 1995606*u^17 + 496789*u^16 + 1359681*u^15 + 633717*u^14 - 1391252*u^13 - 1134228*u^12 + 23037*u^11 + 90805*u^10 - 61896*u^9 - 3149*u^8 + 20301*u^7 + 6122*u^6 - 2016*u^5 + 677*u^4 + 602*u^3 - 206*u^2 + 10*u)*v^13 - (56*u^29 + 518*u^28 - 1180*u^27 - 28689*u^26 - 110463*u^25 - 72591*u^24 + 532919*u^23 + 1572965*u^22 + 2050200*u^21 + 2055500*u^20 + 2299916*u^19 + 1378159*u^18 - 764263*u^17 - 1083430*u^16 - 144121*u^15 - 662529*u^14 - 1134228*u^13 - 250480*u^12 + 285759*u^11 + 106663*u^10 + 23031*u^9 + 27644*u^8 + 7908*u^7 - 1915*u^6 - 3081*u^5 + 340*u^4 + 374*u^3 - 43*u^2)*v^12 - (28*u^29 + 41*u^28 - 2675*u^27 - 15846*u^26 - 6494*u^25 + 181570*u^24 + 582618*u^23 + 661170*u^22 + 87150*u^21 - 359404*u^20 - 638292*u^19 - 1682477*u^18 - 2171667*u^17 - 750477*u^16 + 431018*u^15 + 154974*u^14 + 23037*u^13 + 285759*u^12 + 155766*u^11 - 36185*u^10 - 14385*u^9 + 7356*u^8 + 1162*u^7 - 1581*u^6 - 2338*u^5 - 158*u^4 + 100*u^3)*v^11 - (8*u^29 - 112*u^28 - 1431*u^27 - 359*u^26 + 37109*u^25 + 136391*u^24 + 99010*u^23 - 375960*u^22 - 929392*u^21 - 930875*u^20 - 842095*u^19 - 959700*u^18 - 304949*u^17 + 783131*u^16 + 817218*u^15 + 227854*u^14 + 90805*u^13 + 106663*u^12 - 36185*u^11 - 75168*u^10 - 17094*u^9 + 2337*u^8 + 3113*u^7 + 1276*u^6 - 682*u^5 - 121*u^4)*v^10 - (u^29 - 56*u^28 - 51*u^27 + 4004*u^26 + 17841*u^25 - 5220*u^24 - 183102*u^23 - 416004*u^22 - 269300*u^21 + 216635*u^20 + 391663*u^19 + 338305*u^18 + 530548*u^17 + 511324*u^16 + 42299*u^15 - 202117*u^14 - 61896*u^13 + 23031*u^12 - 14385*u^11 - 17094*u^10 - 1094*u^9 + 77*u^8 + 1562*u^7 + 1210*u^6 + 22*u^5)*v^9 + (9*u^28 - 153*u^27 - 1096*u^26 + 3160*u^25 + 35777*u^24 + 75758*u^23 - 33822*u^22 - 342496*u^21 - 484292*u^20 - 240353*u^19 - 11561*u^18 + 16436*u^17 + 122668*u^16 + 239109*u^15 + 147793*u^14 + 3149*u^13 - 27644*u^12 - 7356*u^11 - 2337*u^10 - 77*u^9 + 2442*u^8 + 594*u^7 - 165*u^6)*v^8 - (36*u^27 - 168*u^26 - 2696*u^25 - 4290*u^24 + 26772*u^23 + 110493*u^22 + 130790*u^21 - 46742*u^20 - 251422*u^19 - 233970*u^18 - 114847*u^17 - 78023*u^16 - 53952*u^15 - 285*u^14 + 20301*u^13 + 7908*u^12 + 1162*u^11 + 3113*u^10 + 1562*u^9 - 594*u^8 - 264*u^7)*v^7 + (84*u^26 + 84*u^25 - 3035*u^24 - 11638*u^23 - 2782*u^22 + 61463*u^21 + 135490*u^20 + 104460*u^19 - 8246*u^18 - 61201*u^17 - 40853*u^16 - 23633*u^15 - 17348*u^14 - 6122*u^13 + 1915*u^12 + 1581*u^11 - 1276*u^10 - 1210*u^9 - 165*u^8)*v^6 - (126*u^25 + 503*u^24 - 1283*u^23 - 10410*u^22 - 18698*u^21 + 1586*u^20 + 48217*u^19 + 67505*u^18 + 39583*u^17 + 9377*u^16 + 1223*u^15 + 38*u^14 - 2016*u^13 - 3081*u^12 - 2338*u^11 - 682*u^10 + 22*u^9)*v^5 + (126*u^24 + 709*u^23 + 769*u^22 - 3450*u^21 - 11844*u^20 - 13282*u^19 - 976*u^18 + 11270*u^17 + 10459*u^16 + 3724*u^15 + 137*u^14 - 677*u^13 - 340*u^12 + 158*u^11 + 121*u^10)*v^4 - (84*u^23 + 542*u^22 + 1294*u^21 + 939*u^20 - 1627*u^19 - 4094*u^18 - 3215*u^17 - 334*u^16 + 1046*u^15 + 857*u^14 + 602*u^13 + 374*u^12 + 100*u^11)*v^3 + (36*u^22 + 242*u^21 + 699*u^20 + 1117*u^19 + 1044*u^18 + 566*u^17 + 270*u^16 + 341*u^15 + 380*u^14 + 206*u^13 + 43*u^12)*v^2 - (9*u^21 + 59*u^20 + 172*u^19 + 300*u^18 + 355*u^17 + 299*u^16 + 174*u^15 + 62*u^14 + 10*u^13)*v + u^20 + 6*u^19 + 15*u^18 + 20*u^17 + 15*u^16 + 6*u^15 + u^14;
q := (u+v+2)/(v-u);
t := (u-v)*(u+v)*(u+v+2)/(u^3+u^2*v+2*u^2+u*v^2+2*u*v+v^3+2*v^2-2);
E := [0,t^2-2*q*t-2,0,-(t^2-1)*(q*t+1)^2,0];
P := [R!0,0];
Q := [(t+1)*(q*t+1),t*(q*t+1)*(t+1)];