syms s
z=(s^4+10*s^2+9)/(s^3+4*s)
f1=maple('convert',z,'parfrac',s);
pretty(f1)
f2=maple('convert',1/z,'parfrac',s);
pretty(f2)