Simple Harmonic Motion

Objective: YouTubLink : https://youtu.be/AnfUQa9Oork/ clc clear all close all %inputs l=1.00; m =1.0; g = 9.81; k = (g/(l*m))^2; % ODE funtion f = @(t,y) [y(2); -k*sin(y(1))]; % initial values for Disp y_ini= [pi/2;0]; %time points t_span = linspace(0,10,500); %Solve ODE [t,y] = ode45(f,t_span,y_ini); %Plot of Results figure(201)…