Solving second order ODE

Objective: clear all close all clc % Input Variables b=0.05; %damping coefficient g=9.81; %acceleration due to gravity l=1; %length of the string m=1; %mass of the pendulum bob %Initial_position theta_0 = [0; 3]; %Time t_span = linspace(0 ,10, 500); %Solve ODE [t,results] = ode45(@(t,theta) ode_func(t, theta, b, g, l ,m), t_span,…