s = tf('s');
Gs = (s^2 + 4*s + 2) / (s^3*(s^2+4)*(s^2+4*s));