Pendulum

Objective: close all clear all clc %inputs b=0.05; m=0.1; g=9.81; l=1; %initial condition theta_0=[0;3]; %time points t_span=linspace(0,20,400); %solve ode [t,results]=ode45(@(t,theta_1) ode_func(t,theta_1,b,g,l,m),t_span,theta_0); figure(1) plot(t,results(:,1)) hold on plot(t,results(:,2)) ylabel(\'plot\') xlabel(\'time\') %alternative…