Solving 2nd order ODE using Python

Objective: Objective: To write a program to simulate the transient behavior of a simple pendulum. Program: #ODE import math import matplotlib.pyplot as plt import numpy as np from scipy.integrate import odeint b = 0.05 #damping co-eff g = 9.81 #gravity l = 1 #length of pendulum m = 1 #mass of the pendulum ct = 1 theta_0…