Second order differential equation

Objective: calling functions for the main prgram function [dtheta_dt] = func_ode(t,theta,b,g,l,m) 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 main program clear all close all clc b = 0.05; g = 9.81; l = 1; m = 1; %intial conddition…