/**/ use R ::= QQ[x,y,z];
/**/ f := x^3*z+x*y+x*z+y+2*z;
/**/ Cx := CoefficientsWRT(f, x); -- same as...
/**/ Cx := CoefficientsWRT(f, [x]);
/**/ indent(Cx);
[
record[PP := x^3, coeff := z],
record[PP := x, coeff := y +z],
record[PP := 1, coeff := y +2*z]
]
/**/ f = sum([M.coeff * M.PP | M in Cx]);
true
/**/ Foreach M in Cx Do Print " +(", M.coeff, ")*", M.PP; EndForeach;
+(y +2*z)*1 +(y +z)*x +(z)*x^3
/**/ Cxz := CoefficientsWRT(f, [x,z]);
/**/ indent(Cxz);
[
record[PP := x^3*z, coeff := 1],
record[PP := x*z, coeff := 1],
record[PP := x, coeff := y],
record[PP := z, coeff := 2],
record[PP := 1, coeff := y]
]
|