I came across an error in the bode plot function. When you specify the frequency limits of the bode plot using the 'omega'-statement, the plot seems to scale linear instead of logarithmic. Hence, the plot is wrong.

To reproduce this try the following:
sys = control.tf([1], [1,1])
control.bode(sys)
control.bode(sys, omega=[0.1, 10])