ai = s * (s^2+4*s+4);
b = s+1;

%unstable p
[p,q] = axbyc(ai,b, (s+3)^4)
[p,q] = axbyc(ai,b, (s+3)^5)

%stable p
[p,q] = axbyc(ai,b, (s+1.5)^4)
[p,q] = axbyc(ai,b, (s+1.5)^5)

