The system solver behaves slightly differently than the scalar solver. For linear systems of equations, see solve[linear].

Multiple solutions are reported only once. Solutions of irreducible polynomials of degree 2, 3, or 4 are represented implicitly by RootOfs.

Inequalities ($\mathrm{expr1}\ne \mathrm{expr2}$) are either part of the input or are generated from denominators, singularities, and so on. That the solutions will satisfy the equations and inequalities is guaranteed.

If there are more unknowns than equations, the solver will select the variables for which the system is easiest to solve.

If the global variable _MaxSols is set to an integer value, the solver will stop after _MaxSols solutions are found. By default this value is set to 100.

$\mathrm{solve}\left(\left\{x+2y=3\,y+\frac{1}{x}=1\right\}\,\left[x\,y\right]\right)$

$\left[\left[x=\mathrm{-1}\,y=2\right]\,\left[x=2\,y=\frac{1}{2}\right]\right]$

$\mathrm{solve}\left(a{x}^2+bx+c\,\left[a\,b\,c\,x\right]\right)$

$\left[\left[a=a\,b=b\,c=-a{x}^2-bx\,x=x\right]\right]$

$\mathrm{solve}\left(\left\{x+y^2\,x^2-y^2-y\,x\ne 0\right\}\right)$

$\left\{x=-{\mathrm{RootOf}\left({\mathrm{\_Z}}^3-\mathrm{\_Z}-1\right)}^2\,y=\mathrm{RootOf}\left({\mathrm{\_Z}}^3-\mathrm{\_Z}-1\right)\right\}$