多项式展开的逆过程的MATLAB实现

在maple软件中有convert()函数，可以让例如a^2-2ab+b^2形式的多项式写成(a-b)^2形式，在MATLAB里这个功能并没有明确给出，有时候会出现很长的式子看不明白，例如式子：

H=27*s^3*x1^3 - 27*s^3*x1^2 + 9*s^3*x1 - s^3 + 27*s^2*x1^3 - 18*s^2*x1^2 + 3*s^2*x1 + 9*s*x1^3 - 3*s*x1^2 + x1^3

这个式子看起来非常复杂，项数多，而且还是双变量的形式，所以把它写成简单形式显得尤为重要(单变量也适合)；上述式子其实是

(x1 - s + 3*s*x1)^3

这个式子看起来就简单多了。本人是用MATLAB自带的factor()函数实现，代码如下：

function convert=convert(f)
n=length(factor(expand(f)));
R=factor(f);
g=1;
for i=n:-1:1
    g=g*R(i);
end
convert=g;
end

例如上述例子：

再如：G=s^3+3*s^2*t+3*s*t^2+t^3