Question: i need help with the matlab code consider the frictionless...

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 0≤𝑡≤5