"Source: Text/Geometry/CrvG3.text";
"Line: 97";
"Date: Thu Sep 25 14:47:45 2025";
"Main: Fri Sep 26 14:29:19 2025";
// original file: Text/Geometry/CrvG3.text, line: 97
// Example: H136E1 ()
print "Example: H136E1";
ei := GetEchoInput();
SetEchoInput(true);
P<x,y,z> := PolynomialRing(Rationals(), 3);
PP := ProjectiveSpace(P);
f1 := x^3*y + y^3*z + z^3*x;
f2 := x^4 + 7*x^3*z + 3*x^2*y^2 - 3*x^2*z^2 - 6*x*y*z^2 - 5*x*z^3 +
2*y^3*z + 3*y^2*z^2 + 2*y*z^3 - 4*z^4;
C1 := Curve(PP, f1); DO1 := DixmierOhnoInvariants(C1 : normalize := true);
C2 := Curve(PP, f2); DO2 := DixmierOhnoInvariants(C2 : normalize := true);
DO1 eq DO2;
IsIsomorphicPlaneQuartics(C1, C2);
SetEchoInput(ei);
