m := 3; n := 3; g := 1; R:=PolynomialRing(Rationals(),2); X := 3*v^3 + (3*u^3 + 9*u^2)*v^2 + (u^6 + 6*u^5 + 9*u^4)*v + 4*u^6 + 33*u^5 + 108*u^4 + 186*u^3 + 189*u^2 + 108*u + 27; // Elliptic curve 27a3 t := (u+1)/(u+2); x := (u+1)^2*(u-v)/((v+1)*(u+2)^3); E := [-(t^3-t^2-1),-t^2*(t-1)*(t^2-t+1),-t^2*(t-1)*(t^2-t+1),0,0]; P := [x]; Q := [R!0,0];