Solving square polynomial systems : a practical method using Bezout matrices

Let $f$ be a polynomial system consisting of $n$ polynomials $f_1,\cdots, f_n$in$n$variables$x_1,\cdots, x_n$,withcoefficientsin$$\mathbb{Q}$$ and let $⟨f⟩$ be the ideal generated by $f$. Such a polynomial system, which has as many equations as variables is called a square system. It may be zero-dimensional, i.e the system of equations $f=0$ has finitely many complex solutions, or equivalently the dimension of the quotient algebra $A = $\mathbb{Q}$[x]/\langle f\rangle$ is finite. In this case, the companion matrices of $f$ are defined as the matrices of the endomorphisms of $A$, called multiplication maps, $x_j : \left\vert \begin{array}{c} h \mapsto x_jh \end{array} \right.$,writteninsomebasisof$A$. We present a practical and efficient method to compute the companion matrices of $f$ in the case when the system is zero-dimensional. When it is not zero-dimensional, then the method works as well and still produces matrices having properties similar to the zero-dimensional case. The whole method consists in matrix calculations. An experiment illustrates the method’s effectiveness.