function [ngc,dgc]=lag2(ng0,dg0,w,KK,Pm)
[mu,pu]=bode(KK*ng0,dg0,w);
wgc=spline(pu,w,Pm+5-180);
ngv=polyval(KK*ng0,j*wgc);
dgv=polyval(dg0,j*wgc);
g=ngv/dgv;
alph=abs(1/g);
T=10/alph*wgc;
ngc=[alph*T,1];
dgc=[T,1];