Undamped Pendulum Motion solved by ODE

Objective: %solving ODE for the undamped pendulum motion clear all close all clc %inputs g= 9.81; l= 1; m=1; %initial conditions init= [0;8]; %time points t_span = linspace(0,5,200); %solving ODE [t,x]= ode45(@(t,theta) ode_function(t,theta,g,l,m),t_span,init); % plot the solution figure(1) plot (t,x(:,1)); grid on; hold on; plot(t,x(:,2));…