Procedure from library finvar.lib (see finvar_lib).
Usage:
power_products(dv,d);
dv: an <intvec> giving the degrees of homogeneous polynomials, d: the
degree of the desired power products
Return:
a size(dv)*m <intmat> where each column ought to be interpreted as
containing the exponents of the corresponding polynomials. The product
of the powers is then homogeneous of degree d.