solving second order ODE -simp le pendulum

Objective: simple pendulum program clear all close all clc % Solving second order ODE % Let theta = y and angular velocity = v t=linspace(0,20,1050); theta=[0 3]; y=lsode(\"pendulum_fun\",theta,t); figure(1) subplot(3,1,1) plot(t,y(:,1)) xlabel(\'Time\',\'FontWeight\',\'bold\',\'FontSize\',10) ylabel(\'angular displacement\',\'FontWeight\',\'bold\',\'FontSize\',10)…