Program to solve 2nd order ODE to simulate simple pendulum

Objective: Simple Pendulum Program to solve 2nd order ODE: import math import matplotlib.pyplot as plt import numpy as np from scipy.integrate import odeint #Function that returns dz/dt def model(theta,t,g,b,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,…