Program showing the pendulum motion

Objective: clear allclose allclc b=0.5;g=9.81;l=1;m=1; %initial conditiontheta_0=[0;3]; %time windowt_span=linspace(0,20,500); % solving ode[t,results]=ode45(@(t,theta) ode_func(t,theta,b,g,l,m), t_span, theta_0); %storing t & results in different matricesstore_0 = results(:,1); ct = 1; for i= 1:length(results(:,1)) result=results(i);…