| Geometric resolution |
1.Use from the interpreter
Mmx]
use "geomsolvex"
Mmx]
X == coordinate ('x); Y == coordinate ('y);
Mmx]
x == polynomial_dag (1 :> Rational, X);
y == polynomial_dag (1 :> Rational, Y);
Mmx]
f1 == x^2 + y^2 - 1;
f2 == x^2 + x * y - 2;
Mmx]
probable_geometric_solve_reduced_regular ([f1, f2])
Mmx]
p: Integer == 101;
Mmx]
x == polynomial_dag (1 :> Integer, X);
y == polynomial_dag (1 :> Integer, Y);
Mmx]
f1 == x^2 + y^2 - 1;
f2 == x^2 + x * y - 2;
Mmx]
probable_geometric_solve_reduced_regular ([f1 mod p, f2
mod p])
For long computations, verbosity can be enabled by setting geomsolvex_verbose.