One way to compute the gradients of such a pulse program is via the stochastic parameter shift rule.
However, this involves a Monte Carlo Integration, which, on the other hand, necessitates executing quantum circuits for many intermediate times τ.
4/11