krac2 - continuous SISO system limit feedback gain
krac2(sys) computes the gains g such that the system sys fedback by g (sys/.g) has 2 real equal poles.
h=syslin('c',352*poly(-5,'s')/poly([0,0,2000,200,25,1],'s','c'));
xbasc();evans(h,100)
g=krac2(h)
hf1=h/.g(1);roots(denom(hf1))
hf2=h/.g(2);roots(denom(hf2))