Solving ODE of motion of Simple Pendulum with damping

Objective: Given below is the ODE describing the transient behavior of the Pendulum system. `(d^2theta)/dt^2+b/m(dθ)/dt+g/Lsintheta=0` Second order ODE can be solved using 'ode45' predefined MATLAB function by resolving the second order ODE into two first-order ODE's. let us rewrite the equation as `(d^2theta)/dt^2 =…