Simple Pendulum Oscillations

Objective: clear allclose allclc%inputsb=0.05;l=1;g=9.81;m=1;theta_0 = [0;3];t_span=linspace(0,20,600);%function[t,Results]=ode45(@(time,theta) ode_func(time,theta,b,g,l,m), t_span, theta_0);%for loopfor i=1:length(Results(:,1))%displacement & velocity plotsubplot(5,4,[13,14,17,18]);plot(t,Results(:,1),'color','b');hold on;plot(t,Results(:,2),'color','r');hold…