Simulating the transient behavior a Simple pendulum using Python and creating an animation of its motion

Objective: import numpy as np from scipy.integrate import odeint import math import matplotlib.pyplot as plt def model(theta,t,b,l,g,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 #Inputs l=1 m=1 b=0.05 g=9.81 #Initial…