Add an option that makes it possible to ensure that the last point is exactly for Tmax; make sure that we never exceed Tmax

But isn't the new version more correct?  In our case, this parameter is used to determine the intrinsic acceleration.  Imagine that ``t0 = 0 s``, ``h = 1 s``, and that we have a burst starting at ``0.5 s``.  We will come here with ``t1 + c_[0] * h = 1 s`` and will consider the intrinsic acceleration in the computation of ``f``.  Seems like the right thing to do.  Unfortunately, the Wolfram page referenced above is botched at the moment (formulas showing up as red rectangles) so it's hard to know for sure.