Second Order Differential Equation

Objective: MAIN PROGRAMME clear all close all clc % inputs b = 0.02; g = 9.81; l = 1; m =0.1; % initial condition theta_0 = [0;5]; % time points t_span = linspace(0,20,800) % solving ode ct=1 for i=1:length(t_span) [t,results] = ode45(@(t,theta) ode_func(t,theta,b,g,l,m),t_span,theta_0); x1=0; y1=0; x2=l*sin(results(i,1)); y2=-l*cos(results(i,1));…