$\begin{cases}\mathrm{IPR\_3\_flux}-\frac{\mathrm{V\_p}P\mathrm{k\_p}}{\mathrm{k\_p}+P} & \text{if $\mathrm{time\_}\le 15$}\\ \frac{-\mathrm{V\_p}P\mathrm{k\_p}}{\mathrm{k\_p}+P} & \text{otherwise}\end{cases}$ $\mathrm{J\_flux}-\mathrm{J\_pump}+\mathrm{J\_leak}$ $\mathrm{k\_flux}\mathrm{mu}h(b+\frac{(1-b)c}{\mathrm{k\_1}+c})$ $\frac{\mathrm{gamma}c^2}{\mathrm{k\_gamma}^2+c^2}$ $\mathrm{beta}$ $\frac{P^3}{\mathrm{k\_mu}^3+P^3}$ $\frac{\frac{\mathrm{k\_2}^2}{\mathrm{k\_2}^2+c^2}-h}{\mathrm{tau\_h}}$