rtdpy ncstr
import matplotlib.pyplot as plt
import rtdpy
for n in [1, 2, 5, 10, 100]:
a = rtdpy.Ncstr(tau=1, n=n, dt=.001, time_end=5)
plt.plot(a.time, a.exitage, label="n={}".format(n))
plt.legend()
plt.xlabel('Time')
plt.ylabel('Exit Age Function')
plt.title('Impulse Responses')
Proud Parrot