second Order ODEs

Objective: %Second ODE clear all close all clc %inputs b = 0.05; %Damping Coefficient g = 9.81; %Gravity m/s^2 m = 1; %Mass of the Pendulum in kg l = 1; %Length of the string attached to pendulum in meters %initial conditions theta_0 = [0;3]; %[angular displacement=0;angular velocity=3] %time points t_span = linspace(0,20,500); %pendulum…