w0 = 2*pi*1e6;
Q  = 5;

s = tf('s');
H = (s*w0)^2 / (1 + s/w0/Q + (s/w0)^2);

bode(H);