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

`// NOTE(phl): This may lead to convergence as bad as (1.5 Δt)^5 rather than Δt^5`