Second order ODEs & pendulam motion in MATLAB

Objective: %PENDULAM MOTION https://youtu.be/-hjV3fAfiqg %PLOT 1 %MAIN PROGRAM clear all close all clc b=0.5; g=9.81; %in m/s^2 l=1; % in meter m=1; % in kg %initial condition theta_0 = [0;3]; %time points t_span = linspace(0,20,300); %solve ODE [t, result ] = ode45(@(t,theta) ode_func(t,theta,b,g,l,m) , t_span…