Hello, I am trying to reproduce a model for glycolytic oscillations in mathematica. The model has just three differential equations and many parameters. However there are some relations that are assumed to hold that I do not know how to implement. I am using NDSolve.
The state variables of my system are ATP, F6P and FBP and they change over time (t). The assumptions include, among others, these:
Being Ka=1 and Atot=500; ADP and AMP are parameters.
I do not know how to incorporate these relations in the system. The values are needed to solve numerically but as they depend on each other I cannot define them with an expression.
I will appreciate any help.