Please solve it with matlab. Thank you.

Numerically integrate and plot the response of an underdamped system determined by m = 100 kg, k = 20,000 N/m, and c = 200 kg/s, subject to the initial conditions of x_0 = 0.01 m and v_0 = 0.1 m/s, and the applied force F (t) = 150 cos 5t. Then plot the exact response as computed by equation (2.33). Compare the plot of the exact solution to the numerical simulation, equation (2.33) (- omega^2 A_s + 2 zeta omega_n omega B_s + omega^2 _n A_s - f_0) cos omega t + (- omega^2 B_s - 2 zeta omega_n omega A_s + omega^2 _n B_s) sin omega t = 0