import numpy as np
import matplotlib.pyplot as plt
from scipy import signal

# numH = []
# denH =
# sysH = signal.TransferFunction(numH, denH)
# impulse = signal.impulse(sysH)


num = 1
den = [1, 2, 2]
sys = signal.TransferFunction(num, den)
t = [0.0001 * t for t in range(60000)]
y = signal.step(sys, t)


plt.plot(t, y)
plt.xlabel('t (sec)')
plt.ylabel('y(t)')
plt.title('MAE 171A - Homework 3.3\n Step Response')
plt.show()