I am a 5-year student studying Doctor of Pharmacy program in Faculty of Pharmaceutical Sciences. I am trying to write MCSim code for IV bolus PBPK model to describe a level of drug in plasma. In the fact, the level of drug should decrease over time but when I run my model, the trend of the level of drug is increase. Could you correct my code, please?
[Model file]
States = {M_liv, # mass liver
M_kid, # mass kidney
M_rapf, # mass rapidly-perfused
M_slpf, # mass slowly-perfused
M_plsm, # mass venous blood
M_Ke};
Outputs = {C_liv, C_kid, C_rapf, C_slpf, C_plsm}; # concentration
Inputs = {M_iv}; # Input siR mass
Dynamics {dt (M_plsm) = M_iv + Cout_kid * Q_kid + Cout_liv * Q_liv + Cout_rapf * Q_rapf + Cout_slpf * Q_slpf - (Q_kid + Q_liv + Q_rapf + Q_slpf) * C_plsm;}
[Input file]
Integrate (Lsodes, 1e-5, 1e-7, 1);
Simulation { M_iv = 0.224; # Initial Value of M_iv
PrintStep (C_plsm, 0,72, 12); # Concentration in plasma}