2nd Order ODE

Objective: clear all close all clc b=0.5; g=9.81; l=1; m=1; %initial conditions theta_0=[0;3]; %time points t_span=linspace(0,20,100); [t,theta]= ode45('ode_func1', t_span, theta_0); figure(1) plot(t,theta(:,1)) hold on plot(t,theta(:,2),'red') legend('displacement vs time','velocity vs time') figure(2) frame=1; for i=1:length(t_span)…