ODE- Pendulum Animation

Objective: clear all;close all;clc;b=0.5; %damping coeffg=9.81; %gav(m/s^2)l=1; %metersm=1;%Kg %initial conditiontheta_0=[0;3]; %time spant_span=linspace(0,20,500); %solve ODE[t,results]=ode45(@(t,theta) ode_fun(t,theta,b,g,l,m),t_span,theta_0); angular_disp = results(:,1);ct=1; for i=1:length(t_span) %Coordinates for pendulum x0=0;…