Hello all,

I have noticed that the time in `integrate_ode_bdf(dz_dt, z0, 0, ts, theta, x_r, x_i);`

may go beoynd the limits. For example, if `ts=1:52`

, then time within `integrate_ode_bdf`

could be like 54-something. I assume it could be expected, but it makes also a bit difficult to integrate then the time-varied covariates. For example, if x_r gives a covariate over the time period (0,52), then it may become a bit surprising that integrate_ode_bdf goes to times beyond that.

Am I correct that the main recommendation would be to use `integrate_ode_rk45`

in this case?

Thank you in advance