d = conv([1 0 0],[1e-5 1]);
Gp = tf(sqrt(10)*1e6,d);
margin(Gp*Ge)