2nd order ode

Objective: clear all close all clc %Inputs l = 1; m = 1; b = 0.05; g = 9.81; %Initial Condition theta_0 = [0;3]; %Time points t_span = linspace(0,20,500); %solve ODE [t, results] = ode45(@(t,theta) ode_func(t, theta, b,g,l,m),t_span, theta_0); plot(t,results(:,1)) hold on plot(t,results(:,2)) ylabel(\'Plot\') xlabel(\'time\') …