Sorry, here is a "better" version of "coefficient_list": coefficient_list(p,s) := block( [n:floor(hipow(expand(p),s))], map(lambda([u],coeff(expand(p),s,u)),makelist(i,i,0,n)) )$ Wilhelm