Pendulum animation

Objective: function [dtheta_dt] = ode_funct(t ,theta,b,g,l,m) theta1 = theta(1); theta2 = theta(2); dtheta1_dt = theta2; dtheta2_dt = -(b/m)*theta2-(g/(l*m))^2*sin(theta1); dtheta_dt = [dtheta1_dt;dtheta2_dt]; end clear all close all clc b=0.5; l = 1; m = 1; g = 9.81; %initial conditions theta_0 = [0;3]; %time points t_span = linspace(0,10,100);…