I don’t know is this is the right place but while trying to implement multi dose IV I happened to look at torsten codes and was confused by IV rate implementation. It applies both to 1cmp and 2cmp models.

In the Pred1_one function there is a statement which updates concentration in central compartment:

if ((init[1] != 0) || (rate[1] != 0)) {

a[0] = 1;

pred(0, 1) += PolyExp(x=dt, dose=init[1], rate=0, xinf=0, tau=0, ss=false, a, alpha, n=1) +

PolyExp(dt, 0, rate[1], dt, 0, false, a, alpha, 1);

}

which translates into:

bolusResult = init[1] * exp (-k10 * dt);

result = (1 - exp (-k10 * dt)) / k10;

pred(0, 1) = bolusResult + rate[1] * result;

However units of rate in mg/h and concentration in mg/L become inconsistent. However, if we replace k10 with CL (L/h) then units become consistent.

Please pardon for my ignorance.