I am not exactly certain what you are doing or what your pharmacokinetic model is. Not knowing more, I would do something like this (assuming ‘t’ is the time vector and the integrated values are stored in matrix ‘y’):
Use first set of initial conditions ‘ic’
tspan = linspace(0, 10, 10);
Run ode15s with your ODE function
ic = y(end,:);
t0 = t(end);
ic(appropriate_compartment) = dose(1);
Use new set of initial conditions ‘ic’
tspan = linspace(t0, t0+60, 60);
Run ode15s with your ODE function
ic = y(end,:);
t0 = t(end);
ic(appropriate_compartment) = dose(2);
... continuing for as many steps as necessary. Adjust the number of elements in ‘tspan’ with the third argument of linspace. You will need to vertically concatenate the ‘t’ vectors and ‘y’ matrices from each step (loop iteration) to create the complete simulation.
I’m simply sketching this, corresponding to similar problems I’ve worked on.
Apparently the SimBiology has a set of functions that will make alll this straightforward, although I have no experience with it. (I don’t do such simulations frequently enought to justify buying it.)
