Macaulay2 » Documentation
Packages » NumericalAlgebraicGeometry :: solveSystem(List)
next | previous | forward | backward | up | index | toc

solveSystem(List) -- solve a system of polynomial equations



Solve a system of polynomial equations using homotopy continuation methods.

i1 : R = CC[x,y];
i2 : F = {x^2+y^2-1, x*y};
i3 : solveSystem F 

o3 = {{3.02251e-20+2.05294e-20*ii, -1}, {-3.02251e-20-2.05294e-20*ii, 1}, {1,
     8.30314e-21+4.44918e-20*ii}, {-1, -8.30314e-21-4.44918e-20*ii}}

o3 : List
i4 : R = CC[x,y];
i5 : F = {x^2+y^2-1, x*y, x*(y+1)};
i6 : solveSystem F 

o6 = {{-6.73915e-17-4.6016e-18*ii, -1}, {4.27328e-17+5.26917e-17*ii, 1}}

o6 : List

The system is assumed to have finitely many solutions. If it is not square (number of equations = number of variables), squareUp is applied and solutions to the original system are then picked out from the resulting (larger) set of solutions.

The output (produced by track with default options) contains all points obtained at the end of homotopy paths when tracking starting at the totalDegreeStartSystem. In particular, this means that solving a system that has fewer than Bezout bound many solutions will produce points that are not marked as regular. See track for detailed examples.

Ways to use this method: