p(t) = p_i + \int_{-\infty}^t \, p_u(t-\tau) \, \delta q(\tau) \, d \tau

p(t) = p_i + \sum_k \, p_u(t-\tau_k) \, q(\tau_k)