clc; clear; close all;

s = tf('s');                        % define the variable s as a transfer function of value s
Q = 2;                              % define Q factor
w0 = 2*pi*1000;                     % define corner frequency
G = 10/(1+s/(w0*Q)+(s/w0)^2);       % Define G as a transfer function as shown
h = bodeplot(G);                    % Plot the Bode plot of G(s)
setoptions(h, 'FreqUnits', 'Hz');   % change frequency scale from rad/sec to Hz
set(findall(gcf,'type','line'),'linewidth',2)
grid on;