simple Pendulum using ODE

Objective: https://www.youtube.com/watch?v=5NVw9bxcNls close all clear all clc %inputs b = 0.05; g = 9.81; l = 1; m = 0.1; %initial condition% theta_0 = [0;5]; %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)) hold on plot(t,results(:,2))…