Hello Dave: do this
n=[1,0,4];
d=[ 4.0000 14.5245 27.3703 33.2917 27.3703 14.5245 4.0000];
f=logspace(-3,1,300);
w=2*pi*f;
s=sqrt(-1)*w;
close all
semilogx(f,20*log10(abs(polyval(n,s) ./ polyval(d,s)))); grid on; xlabel('Hz'); ylabel('db'); title('Bode plot');