Solving second order ODE - Simple pendulum

Objective: https://www.youtube.com/watch?v=laCkS4meztk clear all close all clc b = 0.5; g = 9.81; l = 1; m = 1; % initial condition theta_0 = [0;3]; % time t_span = linspace(0,20,500); %solution [t,results] = ode45(@(t,theta) ode_func(t, theta ,b,g,l,m),t_span,theta_0); theta_1 = results(:,1); figure(1) plot(t,results(:,1)) hold…