Simple Pendulum Oscillation

Objective: 1. Main Program: clear all close all clc b = 0.5; g = 9.81; l = 1; m = 1; ct = 1; x0 = 0; y0 = 0; % Initial Conditions theta_0 = [0;3]; % Time Points t_span = linspace(0, 20, 500); % Solve ODE [t, results] = ode45(@(t,theta) pend_func(t,theta,b,g,l,m), t_span, theta_0); figure (1) plot(t, results(:,1), \'color\', \'b\')…