Solving second order ODEs with animation

Objective: animation : https://youtu.be/1kDMFahAp8U clear all close all clc %constants b = 0.05; g = 9.81; l = 1; m = 0.1; %initial condition, distance and velocity theta_0 = [0 ; 5]; %time span tspan = linspace(0,20,500); %SOLVING THE ODE [t, results] = ode45(@(t,theta) odefun(t,theta,b,g,l,m), tspan, theta_0) ; figure(1) hold…