MOTION OF PENDULUM

Objective: PENDULUM MAIN CODE..... clear all close all clc b=0.5; g=9.81; %in m/s^2 l=1; % in meter m=1; % in kg %initial condition theta_0 = [0;3]; %time points t_span = linspace(0,20,300); %solve ODE [t, result ] = ode45(@(t,theta) ode_func(t,theta,b,g,l,m) , t_span ,theta_0); plot(t , result(:,1),\'linewidth\',2) hold on plot(t…