Simple Pendulum simulation using second order ODE

Objective: Simple Pendulum Problem Youtube Link: https://youtu.be/ADz-aGrM5Gk clear all close all clc b = 0.5; g = 9.81; l = 1; m = 1; % Initial Condition theta_0 = [0;3]; % Time points t_span = linspace(0,20,500); % Solve ODE [t,results] = ode45(@(t,theta) ode_func(t,theta,b,g,l,m),t_span,theta_0); figure(1) plot(t,results(:,1))…