Pêndulo Forçado

g=978. (* in cm/s^2 *)

Solução da equação diferencial do exemplo 7.5 do Marion

solucao=ParametricNDSolve[{θ''[t]-a ω^2 Cos[θ[t]-ω t]/b+g Sin[θ[t]]==0,θ[0]==0,θ'[0]==0},θ,{t,0,10},{a,ω,b}]

Fixando raio do disco giratório em 3 cm. 
ω é a frequência de rotação (rd) e b comprimento do pêndulo, em cm

Manipulate[ Plot[Evaluate[ θ[3,ω,b][t]/.solucao ],{t,0,10},AxesLabel->{"t (s)", "θ (rd)"}],{{ω,10},1,100},{{b,30},10,100}]