I need help with the Matlab code

Consider the frictionless rod, i.e. =0=0. The equation of motion becomes

(d^2 r)/(dt^2 )-m^2 r=-mg sin(t)

with g=9.81 m/s^2 and a constant angular speed .

The rod is initially horizontal, and the initial conditions for the bead are r(0)=r_0 and r (0)=v_0.

Suppose =2, graph the solutions, (), for the initial conditions given here: _0=0 and initial velocities of _0=2.40, 2.45, 2.50, and the initial velocity (0) you found in part B). Use 05