Solving a second order ODE

Objective: %{ program to solve a second order ODE. Method applied is that we convert it to n number of first order ODEs and then solve. created: 12/06/2018 eqn: (d^2(theta)/dt^2) + (b/m)(d(theta)/dt) + (g/L)sin(theta) = 0 %} clc clear all close all %inputs b = 0.5; g = 9.81; L = 1; m = 1; %initial condition theta_0 = [0;5]; %initial…