import matplotlib.pyplot as pltimport numpy as npimport control as ctrlR = 22e3C = 470e-6circuit = ctrl.ss(-1/(R*C), 1/(R*C), 1, 0)t, resp = ctrl.step_response(circuit)plt.plot(t, resp)plt.title('Step Response of the RC Circuit')plt.grid()print("The eigenvalue is: ",-1/(R*C))print("numpy gives the eigenvalue, and eigen vector respectively as: ")np.linalg.eig(circuit.A)t = np.arange(0, 120, 0.01)u = 3.3*148/256-3.3*40/256*np.heaviside(t-60,1)plt.plot(t,u)x0 = 3.3*128/256t, vc = ctrl.input_output_response(circuit, t, u, x0)plt.plot(t,vc)plt.grid()plt.title("Response of RC Circuit")plt.legend(["u","$v_c$"])plt.xlabel("time (seconds)")plt.ylabel("Voltage (Volts)")plt.show()