Analytical solution: \[Cp={Dose*k_a \over V_d*(k_a-k_e)}*(e^{-k_e*t} - e^{-k_a*t})\] \[k_e={V_d \over t_{1/2}}\]