g=tf([1 2],[1 4 5])
[mag,phase, w]=bode(g)
margin(mag,phase,w)