clear all
s=tf('s');
gs1=100/(s^3+11*s^2+10*s+1)
gs0=feedback(gs1,1)
ltiview('step',gs1,'-',gs0,'-r')
