Solving 2nd order ODE and creating animation in MATLAB

Objective: Program for solving 2nd order ODE and simulating the simple pendulum motion: clear all close all clc %inputs b = 0.05; %damping coefficient m = 1; %mass g = 9.81; %acceleration due to gravity l = 1; %length of simple pendulum %initial condition theta_0 = [0;3]; %interval of integration t_span = linspace(0,20,500);…