I want to reproduce an (old) biochemistry paper of Gilman and Ross, i.e. " Genetic algorithm selecetion of a regulatory structure that directs flux in a simple metabolic model." ( The following link leeds to the paper: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1236362/) I know this is an old paper, but maybe someone here remembers it. I really need the answers of the following questions because the enzymatic model explained in the paper is very important for my master thesis.
First: The author defines a system of two ODEs but the initial conditions of these ODEs are not printed in the paper. The author just writes: "Initial conditions for the integration are determined by letting the network relax to steady state. " (see p. 1324 on the lower right) . But how is this meant? In this paper the author wants to optimize the mentioned ODEs , so how can I find initial conditions, when there are not all parameters of the ODEs given?
Second: On the top of page 1324 are two needstate functions defined, called $\xi_T$ and $\xi_F$. When I plot $\xi(F)$, the plot looks completely different from the one depicted in figure 3b in the paper. So the function defined in figure 3b) must be wrong. Did anyone realize that error and has a workaround for that? Of course I could define my own functions but I want to stay as close to the model in the paper as possible.
Thank you very much for your help.