function [Gc,KK]=rlag(G,zeta,ew,Tz);
G=tf(G);[r,k]=rlocus(G);
za=zeta/sqrt(1-zeta^2);ri=r(1,find(imag(r(1,:))>0));
ra=imag(ri)./real(ri);KK=spline(ra,k(find(imag(r(1,:))>0)),1/za);
syms x Gn Gd;Gn=poly2sym(G.num{1});Gd=poly2sym(G.den{1});ess=limit(Gn*KK/Gd*x);
beta=round(1/sym2poly(ess)/ew);Tp=Tz/beta;Gc=tf([1 Tz],[1 Tp]);