The particular case of
n = 2 was already solved by Hilbert in 1893. The general problem was solved in the affirmative, in 1927, by
Emil Artin, for positive semidefinite functions over the reals or more generally
real-closed fields. An algorithmic solution was found by
Charles Delzell in 1984. A result of
Albrecht Pfister shows that a positive semidefinite form in
n variables can be expressed as a sum of 2
n squares. Dubois showed in 1967 that the answer is negative in general for
ordered fields. In this case one can say that a positive polynomial is a sum of weighted squares of rational functions with positive coefficients. McKenna showed in 1975 that all positive semidefinite polynomials with coefficients in an ordered field are sums of weighted squares of rational functions with positive coefficients only if the field is dense in its real closure in the sense that any interval with endpoints in the real closure contains elements from the original field. A generalization to the matrix case (matrices with polynomial function entries that are always positive semidefinite can be expressed as sum of squares of symmetric matrices with rational function entries) was given by Gondard,
Ribenboim and Procesi, Schacher, with an
elementary proof given by Hillar and Nie. In
complex analysis and
complex geometry, the Hermitian analogue, requiring the squares to be squared norms of holomorphic polynomials, was proven for strictly positive polynomials by
Quillen using techniques based on
elliptic partial differential equations. The sum-of-squares representation is unique if it exists, which was first observed by Putinar in the context of optimization. == Minimum number of square rational terms ==