I think its because of your initial conditions for theta and omega. I did the same as you did but without the range in the initial conditions and I got a plot. Try changing the initial condition for theta to Pi/2 and omega for 1 to see if the error disappears.