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)

k=length(t)

while (y(k)>0.98*finalvalue)&(y(k)<1.02*finalvalue)

k=k-1

end

settlingtime=t(k)

p

w=vpa(w,4)

u=vpa(u,4)

yss

timetopeak

finalvalue

settlingtime

-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)　k = length (t)　while (y (k)>　0.98* finalvalue) &　(y (k) < 1.02* finalvalue) k = k-1　end settlingtime = t (k)　p w = vpa (w, 4) u = vpa (u, 4) yss timetopeak finalvalue settlingtime