syms s 
f =((s+1)^3)*(s^2+1);
% f1 = expand(f); expand，collect皆可
f1=collect(f);
% f1=maple('sort',f1,s); 可有可无
N = maple('degree',f1,s);
deg=[];
for i=0:N 
    deg=[deg,double(maple('coeff',f1,s,N-i))]; % 最高阶在前
end

或者
syms s 
f =((s+1)^3)*(s^2+1);
% f1 = expand(f); 
f1=collect(f);
deg=sym2poly(f1);