Let $P \in \Bbb R[X]$. Show that if $\forall x \in [0,1], P(x) \ge 0$ therefore $\exists A,B,C,D \in \Bbb R[X], P=A^2+XB^2+(1-X)C^2+X(1-X)D^2$.
I showed that this is true for $\deg P=2$ but I don't see how to generalize my proof. I tried to introduce coefficients for the polynomials, also tried Lagrange interpolation without success. Does someone have any hint?
If we consider $P \ge 0$ on $\Bbb R$ a classical problem shows that $P=A^2+XB^2$. Adapting the proof of this result seems difficult.
Edit: I found this reference (a very similar problem with solution): George Pólya and Gábor Szegő, Problems and Theorems in Analysis II (problem $47$, p.$78$).