d = conv([0.1 1],[0.1 -1]);
d = conv(d,[1e-5 1]);

Gp = tf(sqrt(10)*1e4,d);
margin(Gp*Ge)