d = conv([1 1],[1 1]);
d = conv(d,[1 1]);

Gp = tf(1e4*[1 0 0],d);

margin(Gp*Gd)
xlim([1e-3 1e8]);