257/CH9/EX9.30/example_9_30.sce


1
2
3
4
5
6
7
8
9
10
11

s=%s
sys=syslin('c',(k)/((s)*(s^2+8*s+17)))
clf
evans(sys)

//damping ratio 0.5

f=(-1.15+%i*2)
disp("k for damping ratio 0.5 is")
disp(abs(((f)*(f^2+8*f+17))))