Simple Pendulum Motion

Objective: clear all close all clc b = 0.05 g = 9.81 l = 1 m = 1 theta_0 = [0;5] t_span = linspace(0,20,500) [t,results] = ode45( @(t,theta) ode_func1(t,theta,b,g,l,m),t_span, theta_0); subplot(5,4,[13,14,17,18]); plot(t,results(:,1),\'color\',\'b\'); hold on plot(t,results(:,2),\'color\',\'r\'); hold off subplot(5,4,[15,16,19,20]);…