Simulation of the transient behaviour of a simple pendulum

Objective: import math from scipy.integrate import odeint import matplotlib.pyplot as plt import numpy as np def mod(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 g=9.81 l=1 m=1 #initial condition theta_0=[0,3]…