pendulum ode

Objective: %Calling a function pendulum function [dtheta_dt] = pendulum_ode(t, theta,g,l,m) w = g/(l*m); theta1 = theta(1); theta2 = theta(2); dtheta1_dt = theta2; dtheta2_dt = -w*w*sin(theta1); dtheta_dt = [dtheta1_dt;dtheta2_dt]; end % main program clear all close all clc % inputs g = 9.81; l = 1.5; m = 0.5;…