Solving second order ODEs

Objective: AIM= soving second order odes Main Code import numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt import math def model(theta,t,b,g,l,m): theta1=theta[0] theta2=theta[1] dtheta1_dt=theta2 dtheta2_dt=-(b/m)*theta2-(g/l)*math.sin(theta1) dtheta_dt=[dtheta1_dt,dtheta2_dt] return dtheta_dt b=0.05…