h(p, t-τJ) is an array. When you subtract a scalar from it you get an error. Did you need intend for it to be an array?
h(p, t-τJ)