Calculate the HRF from a time vector and parameters. Optionally compute the
first or second derivative of the HRF instead. Form of HRF is similar to SPM
but here the response and undershoot are scaled so the difference of the HRFs
peaks at 1 and -c
Usage
HRF_main(t, a1 = 6, b1 = 1, a2 = NULL, b2 = NULL, c = 1/6, o = 0)