SOLVING 2nd Order ODE in OCTAVE

Objective: function[dtheta_dt]=ode_funco(theta,t) b=0.05; g=9.81; l=1; m=0.1; theta1=theta(1) theta2=theta(2) dtheta1_dt=theta2; dtheta2_dt=-(b/m)*theta2-(g/l)*sin(theta1); dtheta_dt=[dtheta1_dt;dtheta2_dt]; end clear all close all clc b=0.05; g=9.81; l=1; m=0.1; %Initial Conditions theta_0=[0;3]; %Time Points t_span=linspace(0,20,500);…