Second order ODE using MATLAB

Objective: Main Program : clear all close all clc % Inputs l = 1 % lenth of the string in metre m = 1 % mass of the ball in kg b = 0.05 % damping coefficient g = 9.81 % acceleration due to gravity in m/s^2 % Initial condition t_span = linspace(0,20,500); % time theta_value = [0;3]; % angular displacement & angular velocity %…