Solving Second order non-Stiff ODE of Underdamped Oscillation using MATLAB

Objective: The following is the code of the function used for solving the ODE: function [dtheta_dt] = ODESOLVER(t,theta,c,g,L,m) theta1 = theta(1); theta2 = theta(2); dtheta1_dt = theta2; dtheta2_dt = (((-c)/m)*theta2) - ((g/L)*sin(theta1)); dtheta_dt = [dtheta1_dt;dtheta2_dt]; end The follwing is the code used to plot and animate:…