These routines provides access to variable history at lagged times.
The lagged time t must not be less than t0, nor
should it be greater than the current time of gradient calculation.
The routine cannot be directly called by a user, and will only work
during the integration process as triggered by the
Access history at time
Vector of variable history at time