G = 200000 * tf([1, 1], [1, 110, 11000, 100000])

%

bode(G)

%

w = [1, 10, 1000]

%

Hjw = zeros(size(w));
for i = 1:length(w)
    Hjw(i) = evalfr(G, w(i) * 1i);
end
Hjw

%

Mag = abs(Hjw)

%

Phase = angle(Hjw)

%

20 * log10(Mag)
rad2deg(Phase)