w = logspace(-2, 4, 200)

%

db = 20 * log10(sqrt((10 + 5 * w.^2).^2 + (49 * w).^2) ./ (100 + w.^2))

%

deg = atand(49 * w ./ (10 + 5 * w.^2))

%

semilogx(w, db)

%

semilogx(w, deg)

%

H = tf([5 1], [1 10])

%

bode(H)