I have a set of 3 differential equations which I solve with ODE. However I would like to add a constaint to c.MOAQ, which value needs to be below 0.001. How can I implement this into the ODE block?
Apart from that, it should be simple to find c.MOX, by integration of the rhs of the first equation. I assume that r.MO, r.EX and c.x are known functions.
Then you fill c.MOX into the second equation and solve that differential equation. Again, D must be a known function of t. This gives you c.MOAQ.
Finally c.MODD is found after filling c.MOAQ in the rhs of the third equation and integrating that.