clear num=[0,0,10] den=[1,2,10] p=roots(den) [u,w]=solve( w^2=10 , 2*w*u=2 , u,w ) [y,x,t]=step(num,den) plot(t,y) [yss,n]=max(y) finalvalue=dcgain(num,den) percentovershoot=100*(yss-finalvalue)/finalvalue timetopeak=t(n