Hi Shashwat,

Your code took a bit a time of run, so I wonder if you're trying to do much at once. It seems you are looking to model dichloromethane (D), benzene (B), toluene (T), ethyl benzene (E), and m-xylene simultaneously, which is fine perhaps if you have compiled code, but running a script as it may not be optimal. Have you tried to build a simpler model, for instance emulating the results as shown in Figure 7A?

Hi Shashwat,

I ran you code by decreasing dt to shorten the simulation time and to me the time course profile looks in similar trend with the published result. The error in your model I find is the scaling of blood flow and alveolar ventilation rate, as you are muliplying with body weight, where it should be BW^0.75. if you want to scale with body weight then you should use this value for QCC 4.8 L/h/kg instead of 18, which is basically L/h/BW^0.75.

Thanks,

Raju

Hi Raju ,

Thanks for your input. It has started woking now.

Sahay