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