Week 3 - Solving second order ODEs

Objective: clear all close all clc b=0.05; g=9.81; l=1; m=0.1; t_span=linspace(0,10,500); theta_0=[0;5]; [t, results]=ode45(@(t,theta) ode_func(t,theta,b,g,l,m),t_span,theta_0); subplot(2,2,1) plot(t,results(:,1)) xlabel(\'time\') ylabel(\'Angular Displacement\') title(\'Angluar Displacement vs time\') subplot(2,2,3) plot(t,results(:,2))…