Solving second order ODEs for oscillating pendulum using MATLAB

Objective: clear allclose allclc %inputsg=9.81;l=1;m=1;b=0.05; %initialconditiontheta_0=[0;3]; %time pointst_span=linspace(0,20,500); %solving ode[t,results] = ode45(@(t,theta) sec_order_ode_func(t,theta,b,g,l,m),t_span,theta_0); figure(1)plot(t,results(:,1))hold onplot(t,results(:,2))grid onxlabel(\'time\')ylabel(\'plot\') figure(2)plot(results(:,1),results(:,2));grid…