m := 2; n := 2; g := 0; F := Rationals(); R:=PolynomialRing(F,2); X := u; b := (16*v^2-1)/16; E := [1,-b,-b,0,0]; P := [-(4*v+1)/8,(16*v^2+8*v+1)/32]; Q := [R!0,0];