These functions expand a polynomial into its homogeneous terms
and return a model formula.
sympoly(x, n) creates a formula whose right-hand side represents the
polynomial of degree n in the variable x. Each
homogeneous term x^k is a separate term in the formula.
sympoly(x, y, n) creates a formula representing the
polynomial of degree n in the two variables x and
y.
If f is a formula containing a term of the form
polynom(…) then expand.polynom(f) replaces this term
by its expansion as a sum of homogeneous terms.