257/CH9/EX9.28/example_9_28.sce


1
2
3
4
5
6
7
8
9

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

//stability

[Ki,s]=kpure(sys)
disp("k marginal is")
disp(Ki)