import matplotlib.pyplot as plt
import numpy as np
import skued

times = np.concatenate(
    (np.arange(-10, -2, step=1), np.arange(-2, 2, step=0.05), np.arange(2, 10, step=1))
)
data = skued.exponential(times, 0, 1, 1)
conv = skued.with_irf(2)(skued.exponential)(times, 0, 1, 1)

plt.figure()
plt.plot(times, data, ".k", label='No IRF')
plt.plot(times, conv, "-r", label='With IRF')
plt.legend()
plt.show()