The
Macaulay2 package
SumsOfSquares decomposes polynomials as sums of squares.
It is based on methods to rationalize sums-of-squares decompositions due to Parrilo
and Peyrl. The package features a data type for sums-of-squares polynomials,
support for external semidefinite programming solvers, and optimization over
varieties.