Omega =1; a = 0.02; A = NDSolve[{Theta'[t] == Omega + a Sin[Theta[t]], Theta[0] == 0}, Theta, {t, 0, 200}]; Plot[Evaluate[Mod[Theta[t],2Pi] /. A], {t, 0, 100}]
