How to make bode plot of transfer function?

Hello, i am trying to make a bode plot of the transfer function of a twin-t notch filter, that i am analyzing. I was able to produce the transfer function, and the bode plot by hand, but i am struggling to do it in Matlab, here is what i have so far: r=320; %Resistance c=100*10^-9; %Capactiance p=(1/(r.*c)); %Beta a=0.95; %Signma, adjusts the wiper of a potentiometer f=(5000); %Target Frequency s= (1i*f); h=((s^2 + p^2) ./ (s^2+4*p*s*(1-a)+p^2)); %Transfer function mh=abs(h); %Magitude of h ah=angle(h); %Phase angle of h, in rads. I am trying to plot the body function, and here is one i made for the transfer function H=2/(s+1):   % Example Transfer Function: g(s) = 2/(s+1) % Numerator num = [2]; % Denominator den = [1 1]; % Transfer Function G = tf(num, den) % Plot Frequency Response bode(G), grid   My question is, how would i represent the numerator and the denominator, since they are both powers of of s, which from my understanding the computer understands as s=jw, and thus i solved the transfer function as such. I have tried to understand how exactly the bode function takes it inputs, but i am not understanding that...