Pendulum Motion

Objective: pendulum_main.m clear all close all clc t = linspace(0,20,500); theta0 = [0 3]; g = 9.81; b = 0.05; m = 1; L = 1; ANON = @(theta0, t) ode_second_order_pendulum(theta0,t, g,b,m,L); [theta] = lsode(ANON, theta0, t); %Plotting variation of Angular displacement and velocity w.r.t. time figure(1) subplot(2,1,1) graphs = plot(t,…