1 2 3 4 5 6 7
function y1=locus2(z) s=%s D=s*(s+1)*(s+2) H=syslin('c',1,D) evans(H,z) y1=0; endfunction
