2nd order ODE

Objective: b=0.5; g=9.81; l=1; m=1; %Initial condition theta_0=[0;3]; %Time points t_span=linspace(0,20,500); %solve ODE [t,results]=ode(@(t,theta)ode_func(t,theta,b,g,l,m),t_span,theta_0); function[dtheta_dt]=ode_func(t.theta,b,g,l,m) theta1=theta(1) theta2=theta(2) dtheta1_dt=theta2; dtheta2_dt=-(b/m)*theta2-(g/l)*sin(theta1);…