1 2 3 4 5 6
function y=locus6(z) s=%s H=syslin('c',1/(s*(s+3)*(s^2+(2*s)+2))) evans(H,z) y=0; endfunction
