fg = 1000;
s = 1i * 2 * pi * fg;
omg = 2 * pi * fg;

for n = 1:10

    Hd(n) = ( 1 / ( (s/omg)^2 + sqrt(2) * (s/omg) + 1 ) );
    Hg(n) = ( 1 / ( (omg/s)^2 + sqrt(2) * (omg/s) + 1 ) );
    X(n) = fg * n * 0.1;

end

plot( X, abs(Hd) );
figure;
plot( X, abs(Hg) );