Evolution of the system variables
P_ic
DP_ic = k_E*P_ic*(2.0*G_1*(P_a-P_1) - 2.0*G_2*(P_2-P_c) + G_f*(P_c-P_ic)*heavi(P_c-P_ic) + G_pv*(P_c-P_v) - G_vs*(P_v-P_vs) - G_o*(P_ic-P_vs)*heavi(P_ic-P_vs) + I_i)
******
P_vs
DP_vs = 1.0/C_ve*(G_vs*(P_v-P_vs) + G_o*(P_ic-P_vs)*heavi(P_ic-P_vs) - G_ve*(P_vs-P_cv))
******
P_v
DP_v -1.0DP_ic = 1.0/C_vi*(G_pv*(P_c-P_v) - G_vs*(P_v-P_vs))
******
r_1
Dr_1 = (G_1*(P_a-P_1) - G_1*G_2/(G_1+G_2)*(P_1-P_2))/K_V1/r_1
******
r_2
Dr_2 = (G_1*G_2/(G_1+G_2)*(P_1-P_2) - G_2*(P_2-P_c))/K_V2/r_2
******
_gluc
D_gluc = - Vmax_glucN*_gluc*pow(_NAD, 2.0)*pow(_ADP, 2.0)*pow(_Phos, 2.0)/((Km_glucNgl + _gluc)*(pow(Km_glucNN, 2.0) + pow(_NAD, 2.0))*(pow(Km_glucNA, 2.0) + pow(_ADP, 2.0))*(pow(Km_glucNP, 2.0) + pow(_Phos, 2.0)))+ k_ecgluc*(_egluc - _gluc)/Vol_exm
******
_NAD
D_NAD = - 2.0*Vmax_glucN*_gluc*pow(_NAD, 2.0)*pow(_ADP, 2.0)*pow(_Phos, 2.0)/((Km_glucNgl + _gluc)*(pow(Km_glucNN, 2.0) + pow(_NAD, 2.0))*(pow(Km_glucNA, 2.0) + pow(_ADP, 2.0))*(pow(Km_glucNP, 2.0) + pow(_Phos, 2.0)))+ Vmax_PyL*(1.0/(1.0 + Km_PyLPy/_Py0 + Km_PyLNH/_NADH + Km_PyLPN/_Py0/_NADH + Km_PyLH/_Hy) - VR_PyL/(1.0 + Km_PyLL/_L0 + Km_PyLN/_NAD + Km_PyLLN/_L0/_NAD))+(k_MAshut*_NADH*_Hy*_mNAD - k_nMAshut*_mNADH*_mHy*_NAD)/Vol_exm
******
_ADP
D_ADP = - 2.0*Vmax_glucN*_gluc*pow(_NAD, 2.0)*pow(_ADP, 2.0)*pow(_Phos, 2.0)/((Km_glucNgl + _gluc)*(pow(Km_glucNN, 2.0) + pow(_NAD, 2.0))*(pow(Km_glucNA, 2.0) + pow(_ADP, 2.0))*(pow(Km_glucNP, 2.0) + pow(_Phos, 2.0)))-(k_PCrATP*_PCr*_ADP*_Hy - k_nPCrATP*_Cr*_ATP)+ Vmax_ATPuse*_ATP/(Km_ATPuse + _ATP)- 2.0*(k_ADPATP*pow(_ADP, 2.0) - k_nADPATP*_ATP*_AMP)- Vmax_ANT*((1.0-k1_ANT*_ATP*_mADP/(_ADP*_mATP))/(1.0 + k2_ANT*exp(-0.5*PD_mit/26.73)*_ATP/_ADP)/(1.0 + k3_ANT*_mADP/_mATP))/Vol_exm+ Vmax_KATP*(_ATP/(Km_KATPA + _ATP)*(_eK/_K)/(Km_KATPK + (_eK/_K))*(_Na/_eNa)/(Km_KATPNa + (_Na/_eNa)))/Vol_exm
******
_Phos
D_Phos = - 2.0*Vmax_glucN*_gluc*pow(_NAD, 2.0)*pow(_ADP, 2.0)*pow(_Phos, 2.0)/((Km_glucNgl + _gluc)*(pow(Km_glucNN, 2.0) + pow(_NAD, 2.0))*(pow(Km_glucNA, 2.0) + pow(_ADP, 2.0))*(pow(Km_glucNP, 2.0) + pow(_Phos, 2.0)))+ Vmax_ATPuse*_ATP/(Km_ATPuse + _ATP)- Vmax_Phosshut*((_Phos*_Hy - _mPhos*_mHy)/((_Phos*_Hy + K_Phosshut*(1.0 + p_Phosshut))*(_mPhos*_mHy + K_Phosshut*(1.0 + p_Phosshut)) - pow(K_Phosshut, 2)*pow(p_Phosshut, 2)))/Vol_exm+ Vmax_KATP*(_ATP/(Km_KATPA + _ATP)*(_eK/_K)/(Km_KATPK + (_eK/_K))*(_Na/_eNa)/(Km_KATPNa + (_Na/_eNa)))/Vol_exm
******
_NADH
D_NADH = + 2.0*Vmax_glucN*_gluc*pow(_NAD, 2.0)*pow(_ADP, 2.0)*pow(_Phos, 2.0)/((Km_glucNgl + _gluc)*(pow(Km_glucNN, 2.0) + pow(_NAD, 2.0))*(pow(Km_glucNA, 2.0) + pow(_ADP, 2.0))*(pow(Km_glucNP, 2.0) + pow(_Phos, 2.0)))- Vmax_PyL*(1.0/(1.0 + Km_PyLPy/_Py0 + Km_PyLNH/_NADH + Km_PyLPN/_Py0/_NADH + Km_PyLH/_Hy) - VR_PyL/(1.0 + Km_PyLL/_L0 + Km_PyLN/_NAD + Km_PyLLN/_L0/_NAD))-(k_MAshut*_NADH*_Hy*_mNAD - k_nMAshut*_mNADH*_mHy*_NAD)/Vol_exm
******
_Py0
D_Py0 = + 2.0*Vmax_glucN*_gluc*pow(_NAD, 2.0)*pow(_ADP, 2.0)*pow(_Phos, 2.0)/((Km_glucNgl + _gluc)*(pow(Km_glucNN, 2.0) + pow(_NAD, 2.0))*(pow(Km_glucNA, 2.0) + pow(_ADP, 2.0))*(pow(Km_glucNP, 2.0) + pow(_Phos, 2.0)))- Vmax_PyL*(1.0/(1.0 + Km_PyLPy/_Py0 + Km_PyLNH/_NADH + Km_PyLPN/_Py0/_NADH + Km_PyLH/_Hy) - VR_PyL/(1.0 + Km_PyLL/_L0 + Km_PyLN/_NAD + Km_PyLLN/_L0/_NAD))- Vmax_Pyshut*((_Py0*_Hy - _mPy0*_mHy)/((_Py0*_Hy + K_Pyshut*(1.0 + p_Pyshut))*(_mPy0*_mHy + K_Pyshut*(1.0 + p_Pyshut)) - pow(K_Pyshut, 2)*pow(p_Pyshut, 2)))/Vol_exm+ k_ecPy0*(_ePy0 - _Py0)/Vol_exm
******
_Hy
D_Hy = + 4.0*Vmax_glucN*_gluc*pow(_NAD, 2.0)*pow(_ADP, 2.0)*pow(_Phos, 2.0)/((Km_glucNgl + _gluc)*(pow(Km_glucNN, 2.0) + pow(_NAD, 2.0))*(pow(Km_glucNA, 2.0) + pow(_ADP, 2.0))*(pow(Km_glucNP, 2.0) + pow(_Phos, 2.0)))+ N1_H*Vmax_NQ*((_mNADH/_mNAD)/((_mNADH/_mNAD) + (Km_NQN))*_mO2/(_mO2 + Km_NQO)*(_mHy/_Hy)/((_mHy/_Hy) + Km_NQH))/Vol_exm+ N2_H*Vmax_FQ*((_mFADH/_mFAD)/((_mFADH/_mFAD) + (Km_FQF))*_mO2/(_mO2 + Km_FQO)*(_mHy/_Hy)/((_mHy/_Hy) + Km_FQH))/Vol_exm- (k1_lk*(exp(k2_lk*PMF) - 1.0))/Vol_exm- 3.0*VmaxC_oxp*(((100.0*p_a_C + p_c1_C*exp(3.0*0.0374*50.0))*(K_F1*_mATP/_mADP/_mPhos) - exp(3.0*0.0374*PMF)*(p_a_C + p_c2_C*(K_F1*_mATP/_mADP/_mPhos)))/((1.0 + p_1_C*(K_F1*_mATP/_mADP/_mPhos))*exp(3.0*0.0374*50.0) + (p_2_C + p_3_C*(K_F1*_mATP/_mADP/_mPhos))*exp(3.0*0.0374*PMF)))/Vol_exm- Vmax_PyL*(1.0/(1.0 + Km_PyLPy/_Py0 + Km_PyLNH/_NADH + Km_PyLPN/_Py0/_NADH + Km_PyLH/_Hy) - VR_PyL/(1.0 + Km_PyLL/_L0 + Km_PyLN/_NAD + Km_PyLLN/_L0/_NAD))-(k_PCrATP*_PCr*_ADP*_Hy - k_nPCrATP*_Cr*_ATP)+(k_LH*_lac - k_nLH*_L0*_Hy)+(k_CHi*_CO2 - k_nCHi*_Hy*_BiC)+(k_Pbuf*PbufH - k_nPbuf*Pbuf*_Hy)-(k_MAshut*_NADH*_Hy*_mNAD - k_nMAshut*_mNADH*_mHy*_NAD)/Vol_exm- Vmax_Pyshut*((_Py0*_Hy - _mPy0*_mHy)/((_Py0*_Hy + K_Pyshut*(1.0 + p_Pyshut))*(_mPy0*_mHy + K_Pyshut*(1.0 + p_Pyshut)) - pow(K_Pyshut, 2)*pow(p_Pyshut, 2)))/Vol_exm- Vmax_Phosshut*((_Phos*_Hy - _mPhos*_mHy)/((_Phos*_Hy + K_Phosshut*(1.0 + p_Phosshut))*(_mPhos*_mHy + K_Phosshut*(1.0 + p_Phosshut)) - pow(K_Phosshut, 2)*pow(p_Phosshut, 2)))/Vol_exm+ Vmax_ecL0H*(_eL0*_eHy - _L0*_Hy)/((_eL0*_eHy + K_ecL0H*(1.0 + p_ecL0H))*(_L0*_Hy + K_ecL0H*(1.0 + p_ecL0H)) - pow(K_ecL0H, 2)*pow(p_ecL0H, 2))/Vol_exm
******
_ATP
D_ATP = + 2.0*Vmax_glucN*_gluc*pow(_NAD, 2.0)*pow(_ADP, 2.0)*pow(_Phos, 2.0)/((Km_glucNgl + _gluc)*(pow(Km_glucNN, 2.0) + pow(_NAD, 2.0))*(pow(Km_glucNA, 2.0) + pow(_ADP, 2.0))*(pow(Km_glucNP, 2.0) + pow(_Phos, 2.0)))+(k_PCrATP*_PCr*_ADP*_Hy - k_nPCrATP*_Cr*_ATP)- Vmax_ATPuse*_ATP/(Km_ATPuse + _ATP)+(k_ADPATP*pow(_ADP, 2.0) - k_nADPATP*_ATP*_AMP)+ Vmax_ANT*((1.0-k1_ANT*_ATP*_mADP/(_ADP*_mATP))/(1.0 + k2_ANT*exp(-0.5*PD_mit/26.73)*_ATP/_ADP)/(1.0 + k3_ANT*_mADP/_mATP))/Vol_exm- Vmax_KATP*(_ATP/(Km_KATPA + _ATP)*(_eK/_K)/(Km_KATPK + (_eK/_K))*(_Na/_eNa)/(Km_KATPNa + (_Na/_eNa)))/Vol_exm
******
_mPy0
D_mPy0 = - Vmax_PyN*(_mPy0/(_mPy0 + Km_PyNP)*_Ox/(_Ox + Km_PyNO)/(1.0 + Ka_PyNA/_mADP)*pow(_mNAD, 2.0)/(pow(_mNAD, 2.0) + pow(Km_PyNN, 2.0)))+ Vmax_Pyshut*((_Py0*_Hy - _mPy0*_mHy)/((_Py0*_Hy + K_Pyshut*(1.0 + p_Pyshut))*(_mPy0*_mHy + K_Pyshut*(1.0 + p_Pyshut)) - pow(K_Pyshut, 2)*pow(p_Pyshut, 2)))/Vol_inm
******
_Ox
D_Ox = - Vmax_PyN*(_mPy0/(_mPy0 + Km_PyNP)*_Ox/(_Ox + Km_PyNO)/(1.0 + Ka_PyNA/_mADP)*pow(_mNAD, 2.0)/(pow(_mNAD, 2.0) + pow(Km_PyNN, 2.0)))+ Vmax_sucN*_suc*_mNAD*_mFAD/((Km_sucNs + _suc)*(Km_sucNN + _mNAD)*(Km_sucNF + _mFAD))
******
_mNAD
D_mNAD = - 2.0*Vmax_PyN*(_mPy0/(_mPy0 + Km_PyNP)*_Ox/(_Ox + Km_PyNO)/(1.0 + Ka_PyNA/_mADP)*pow(_mNAD, 2.0)/(pow(_mNAD, 2.0) + pow(Km_PyNN, 2.0)))- Vmax_AKGN*_AKG*_mNAD*_mADP*_mPhos/((Km_AKGNA + _AKG)*(Km_AKGNN + _mNAD)*(Km_AKGNAD + _mADP)*(Km_AKGNP + _mPhos))- Vmax_sucN*_suc*_mNAD*_mFAD/((Km_sucNs + _suc)*(Km_sucNN + _mNAD)*(Km_sucNF + _mFAD))+ 2.0*Vmax_NQ*((_mNADH/_mNAD)/((_mNADH/_mNAD) + (Km_NQN))*_mO2/(_mO2 + Km_NQO)*(_mHy/_Hy)/((_mHy/_Hy) + Km_NQH))/Vol_inm-(k_MAshut*_NADH*_Hy*_mNAD - k_nMAshut*_mNADH*_mHy*_NAD)/Vol_inm
******
_AKG
D_AKG = + Vmax_PyN*(_mPy0/(_mPy0 + Km_PyNP)*_Ox/(_Ox + Km_PyNO)/(1.0 + Ka_PyNA/_mADP)*pow(_mNAD, 2.0)/(pow(_mNAD, 2.0) + pow(Km_PyNN, 2.0)))- Vmax_AKGN*_AKG*_mNAD*_mADP*_mPhos/((Km_AKGNA + _AKG)*(Km_AKGNN + _mNAD)*(Km_AKGNAD + _mADP)*(Km_AKGNP + _mPhos))
******
_mCO2
D_mCO2 = + 2.0*Vmax_PyN*(_mPy0/(_mPy0 + Km_PyNP)*_Ox/(_Ox + Km_PyNO)/(1.0 + Ka_PyNA/_mADP)*pow(_mNAD, 2.0)/(pow(_mNAD, 2.0) + pow(Km_PyNN, 2.0)))+ Vmax_AKGN*_AKG*_mNAD*_mADP*_mPhos/((Km_AKGNA + _AKG)*(Km_AKGNN + _mNAD)*(Km_AKGNAD + _mADP)*(Km_AKGNP + _mPhos))-(k_CHm*_mCO2 - k_nCHm*_mHy*_mBiC)+ k_CO2shut*(_CO2 - _mCO2)/Vol_inm
******
_mNADH
D_mNADH = + 2.0*Vmax_PyN*(_mPy0/(_mPy0 + Km_PyNP)*_Ox/(_Ox + Km_PyNO)/(1.0 + Ka_PyNA/_mADP)*pow(_mNAD, 2.0)/(pow(_mNAD, 2.0) + pow(Km_PyNN, 2.0)))+ Vmax_AKGN*_AKG*_mNAD*_mADP*_mPhos/((Km_AKGNA + _AKG)*(Km_AKGNN + _mNAD)*(Km_AKGNAD + _mADP)*(Km_AKGNP + _mPhos))+ Vmax_sucN*_suc*_mNAD*_mFAD/((Km_sucNs + _suc)*(Km_sucNN + _mNAD)*(Km_sucNF + _mFAD))- 2.0*Vmax_NQ*((_mNADH/_mNAD)/((_mNADH/_mNAD) + (Km_NQN))*_mO2/(_mO2 + Km_NQO)*(_mHy/_Hy)/((_mHy/_Hy) + Km_NQH))/Vol_inm+(k_MAshut*_NADH*_Hy*_mNAD - k_nMAshut*_mNADH*_mHy*_NAD)/Vol_inm
******
_mHy
D_mHy = + Vmax_PyN*(_mPy0/(_mPy0 + Km_PyNP)*_Ox/(_Ox + Km_PyNO)/(1.0 + Ka_PyNA/_mADP)*pow(_mNAD, 2.0)/(pow(_mNAD, 2.0) + pow(Km_PyNN, 2.0)))+ Vmax_AKGN*_AKG*_mNAD*_mADP*_mPhos/((Km_AKGNA + _AKG)*(Km_AKGNN + _mNAD)*(Km_AKGNAD + _mADP)*(Km_AKGNP + _mPhos))+ Vmax_sucN*_suc*_mNAD*_mFAD/((Km_sucNs + _suc)*(Km_sucNN + _mNAD)*(Km_sucNF + _mFAD))- (N1_H+2.0)*Vmax_NQ*((_mNADH/_mNAD)/((_mNADH/_mNAD) + (Km_NQN))*_mO2/(_mO2 + Km_NQO)*(_mHy/_Hy)/((_mHy/_Hy) + Km_NQH))/Vol_inm- N2_H*Vmax_FQ*((_mFADH/_mFAD)/((_mFADH/_mFAD) + (Km_FQF))*_mO2/(_mO2 + Km_FQO)*(_mHy/_Hy)/((_mHy/_Hy) + Km_FQH))/Vol_inm+ (k1_lk*(exp(k2_lk*PMF) - 1.0))/Vol_inm+ 3.0*VmaxC_oxp*(((100.0*p_a_C + p_c1_C*exp(3.0*0.0374*50.0))*(K_F1*_mATP/_mADP/_mPhos) - exp(3.0*0.0374*PMF)*(p_a_C + p_c2_C*(K_F1*_mATP/_mADP/_mPhos)))/((1.0 + p_1_C*(K_F1*_mATP/_mADP/_mPhos))*exp(3.0*0.0374*50.0) + (p_2_C + p_3_C*(K_F1*_mATP/_mADP/_mPhos))*exp(3.0*0.0374*PMF)))/Vol_inm+(k_CHm*_mCO2 - k_nCHm*_mHy*_mBiC)+(k_MAshut*_NADH*_Hy*_mNAD - k_nMAshut*_mNADH*_mHy*_NAD)/Vol_inm+ Vmax_Pyshut*((_Py0*_Hy - _mPy0*_mHy)/((_Py0*_Hy + K_Pyshut*(1.0 + p_Pyshut))*(_mPy0*_mHy + K_Pyshut*(1.0 + p_Pyshut)) - pow(K_Pyshut, 2)*pow(p_Pyshut, 2)))/Vol_inm+ Vmax_Phosshut*((_Phos*_Hy - _mPhos*_mHy)/((_Phos*_Hy + K_Phosshut*(1.0 + p_Phosshut))*(_mPhos*_mHy + K_Phosshut*(1.0 + p_Phosshut)) - pow(K_Phosshut, 2)*pow(p_Phosshut, 2)))/Vol_inm
******
_mADP
D_mADP = - Vmax_AKGN*_AKG*_mNAD*_mADP*_mPhos/((Km_AKGNA + _AKG)*(Km_AKGNN + _mNAD)*(Km_AKGNAD + _mADP)*(Km_AKGNP + _mPhos))- VmaxC_oxp*(((100.0*p_a_C + p_c1_C*exp(3.0*0.0374*50.0))*(K_F1*_mATP/_mADP/_mPhos) - exp(3.0*0.0374*PMF)*(p_a_C + p_c2_C*(K_F1*_mATP/_mADP/_mPhos)))/((1.0 + p_1_C*(K_F1*_mATP/_mADP/_mPhos))*exp(3.0*0.0374*50.0) + (p_2_C + p_3_C*(K_F1*_mATP/_mADP/_mPhos))*exp(3.0*0.0374*PMF)))/Vol_inm+ Vmax_ANT*((1.0-k1_ANT*_ATP*_mADP/(_ADP*_mATP))/(1.0 + k2_ANT*exp(-0.5*PD_mit/26.73)*_ATP/_ADP)/(1.0 + k3_ANT*_mADP/_mATP))/Vol_inm
******
_mPhos
D_mPhos = - Vmax_AKGN*_AKG*_mNAD*_mADP*_mPhos/((Km_AKGNA + _AKG)*(Km_AKGNN + _mNAD)*(Km_AKGNAD + _mADP)*(Km_AKGNP + _mPhos))- VmaxC_oxp*(((100.0*p_a_C + p_c1_C*exp(3.0*0.0374*50.0))*(K_F1*_mATP/_mADP/_mPhos) - exp(3.0*0.0374*PMF)*(p_a_C + p_c2_C*(K_F1*_mATP/_mADP/_mPhos)))/((1.0 + p_1_C*(K_F1*_mATP/_mADP/_mPhos))*exp(3.0*0.0374*50.0) + (p_2_C + p_3_C*(K_F1*_mATP/_mADP/_mPhos))*exp(3.0*0.0374*PMF)))/Vol_inm+ Vmax_Phosshut*((_Phos*_Hy - _mPhos*_mHy)/((_Phos*_Hy + K_Phosshut*(1.0 + p_Phosshut))*(_mPhos*_mHy + K_Phosshut*(1.0 + p_Phosshut)) - pow(K_Phosshut, 2)*pow(p_Phosshut, 2)))/Vol_inm
******
_suc
D_suc = + Vmax_AKGN*_AKG*_mNAD*_mADP*_mPhos/((Km_AKGNA + _AKG)*(Km_AKGNN + _mNAD)*(Km_AKGNAD + _mADP)*(Km_AKGNP + _mPhos))- Vmax_sucN*_suc*_mNAD*_mFAD/((Km_sucNs + _suc)*(Km_sucNN + _mNAD)*(Km_sucNF + _mFAD))
******
_mATP
D_mATP = + Vmax_AKGN*_AKG*_mNAD*_mADP*_mPhos/((Km_AKGNA + _AKG)*(Km_AKGNN + _mNAD)*(Km_AKGNAD + _mADP)*(Km_AKGNP + _mPhos))+ VmaxC_oxp*(((100.0*p_a_C + p_c1_C*exp(3.0*0.0374*50.0))*(K_F1*_mATP/_mADP/_mPhos) - exp(3.0*0.0374*PMF)*(p_a_C + p_c2_C*(K_F1*_mATP/_mADP/_mPhos)))/((1.0 + p_1_C*(K_F1*_mATP/_mADP/_mPhos))*exp(3.0*0.0374*50.0) + (p_2_C + p_3_C*(K_F1*_mATP/_mADP/_mPhos))*exp(3.0*0.0374*PMF)))/Vol_inm- Vmax_ANT*((1.0-k1_ANT*_ATP*_mADP/(_ADP*_mATP))/(1.0 + k2_ANT*exp(-0.5*PD_mit/26.73)*_ATP/_ADP)/(1.0 + k3_ANT*_mADP/_mATP))/Vol_inm
******
_mFAD
D_mFAD = - Vmax_sucN*_suc*_mNAD*_mFAD/((Km_sucNs + _suc)*(Km_sucNN + _mNAD)*(Km_sucNF + _mFAD))+ 2.0*Vmax_FQ*((_mFADH/_mFAD)/((_mFADH/_mFAD) + (Km_FQF))*_mO2/(_mO2 + Km_FQO)*(_mHy/_Hy)/((_mHy/_Hy) + Km_FQH))/Vol_inm
******
_mFADH
D_mFADH = + Vmax_sucN*_suc*_mNAD*_mFAD/((Km_sucNs + _suc)*(Km_sucNN + _mNAD)*(Km_sucNF + _mFAD))- 2.0*Vmax_FQ*((_mFADH/_mFAD)/((_mFADH/_mFAD) + (Km_FQF))*_mO2/(_mO2 + Km_FQO)*(_mHy/_Hy)/((_mHy/_Hy) + Km_FQH))/Vol_inm
******
_mO2
D_mO2 = - Vmax_NQ*((_mNADH/_mNAD)/((_mNADH/_mNAD) + (Km_NQN))*_mO2/(_mO2 + Km_NQO)*(_mHy/_Hy)/((_mHy/_Hy) + Km_NQH))/Vol_inm- Vmax_FQ*((_mFADH/_mFAD)/((_mFADH/_mFAD) + (Km_FQF))*_mO2/(_mO2 + Km_FQO)*(_mHy/_Hy)/((_mHy/_Hy) + Km_FQH))/Vol_inm+ k_O2shut*(_O2 - _mO2)/Vol_inm
******
_mBiC
D_mBiC = +(k_CHm*_mCO2 - k_nCHm*_mHy*_mBiC)+(k_BiCshut*_BiC - k_nBiCshut*_mBiC)/Vol_inm
******
_L0
D_L0 = + Vmax_PyL*(1.0/(1.0 + Km_PyLPy/_Py0 + Km_PyLNH/_NADH + Km_PyLPN/_Py0/_NADH + Km_PyLH/_Hy) - VR_PyL/(1.0 + Km_PyLL/_L0 + Km_PyLN/_NAD + Km_PyLLN/_L0/_NAD))+(k_LH*_lac - k_nLH*_L0*_Hy)+ Vmax_ecL0H*(_eL0*_eHy - _L0*_Hy)/((_eL0*_eHy + K_ecL0H*(1.0 + p_ecL0H))*(_L0*_Hy + K_ecL0H*(1.0 + p_ecL0H)) - pow(K_ecL0H, 2)*pow(p_ecL0H, 2))/Vol_exm
******
_PCr
D_PCr = -(k_PCrATP*_PCr*_ADP*_Hy - k_nPCrATP*_Cr*_ATP)
******
_Cr
D_Cr = +(k_PCrATP*_PCr*_ADP*_Hy - k_nPCrATP*_Cr*_ATP)
******
_AMP
D_AMP = +(k_ADPATP*pow(_ADP, 2.0) - k_nADPATP*_ATP*_AMP)-(k_AMPAd*_AMP - k_nAMPAd*_Ad)
******
_Ad
D_Ad = +(k_AMPAd*_AMP - k_nAMPAd*_Ad)+ k_ecAd*(_eAd - _Ad)/Vol_exm
******
_lac
D_lac = -(k_LH*_lac - k_nLH*_L0*_Hy)+ k_ecL*(_eL - _lac)/Vol_exm
******
_CO2
D_CO2 = -(k_CHi*_CO2 - k_nCHi*_Hy*_BiC)- k_CO2shut*(_CO2 - _mCO2)/Vol_exm+ k_ecCC*(_eCO2 - _CO2)/Vol_exm
******
_BiC
D_BiC = +(k_CHi*_CO2 - k_nCHi*_Hy*_BiC)-(k_BiCshut*_BiC - k_nBiCshut*_mBiC)/Vol_exm
******
PbufH
DPbufH = -(k_Pbuf*PbufH - k_nPbuf*Pbuf*_Hy)
******
Pbuf
DPbuf = +(k_Pbuf*PbufH - k_nPbuf*Pbuf*_Hy)
******
_eL
D_eL = -(k_LH*_eL - k_nLH*_eL0*_eHy)- k_ecL*(_eL - _lac)/Vol_ecs+ k_LL*(L - _eL)/Vol_ecs- _eL/Vol_ecs*(G_o*(P_ic-P_vs)*heavi(P_ic-P_vs))
******
_eL0
D_eL0 = +(k_LH*_eL - k_nLH*_eL0*_eHy)- Vmax_ecL0H*(_eL0*_eHy - _L0*_Hy)/((_eL0*_eHy + K_ecL0H*(1.0 + p_ecL0H))*(_L0*_Hy + K_ecL0H*(1.0 + p_ecL0H)) - pow(K_ecL0H, 2)*pow(p_ecL0H, 2))/Vol_ecs- Vmax_L0H*(_eL0*_eHy - L0*Hy)/((_eL0*_eHy + K_L0H*(1.0 + p_L0H))*(L0*Hy + K_L0H*(1.0 + p_L0H)) - pow(K_L0H, 2)*pow(p_L0H, 2))/Vol_ecs- _eL0/Vol_ecs*(G_o*(P_ic-P_vs)*heavi(P_ic-P_vs))
******
_eHy
D_eHy = +(k_LH*_eL - k_nLH*_eL0*_eHy)+(k_CHe*_eCO2 - k_nCHe*_eHy*_eBiC)- Vmax_ecL0H*(_eL0*_eHy - _L0*_Hy)/((_eL0*_eHy + K_ecL0H*(1.0 + p_ecL0H))*(_L0*_Hy + K_ecL0H*(1.0 + p_ecL0H)) - pow(K_ecL0H, 2)*pow(p_ecL0H, 2))/Vol_ecs- Vmax_L0H*(_eL0*_eHy - L0*Hy)/((_eL0*_eHy + K_L0H*(1.0 + p_L0H))*(L0*Hy + K_L0H*(1.0 + p_L0H)) - pow(K_L0H, 2)*pow(p_L0H, 2))/Vol_ecs
******
_eCO2
D_eCO2 = -(k_CHe*_eCO2 - k_nCHe*_eHy*_eBiC)- k_ecCC*(_eCO2 - _CO2)/Vol_ecs+ k_CC*(CO2 - _eCO2)/Vol_ecs- _eCO2/Vol_ecs*(G_o*(P_ic-P_vs)*heavi(P_ic-P_vs))
******
_eBiC
D_eBiC = +(k_CHe*_eCO2 - k_nCHe*_eHy*_eBiC)+ Vmax1_BiC*BiC/(Km_BiC + BiC)/Vol_ecs- Vmax2_BiC*_eBiC/(Km_BiC + _eBiC)/Vol_ecs- _eBiC/Vol_ecs*(G_o*(P_ic-P_vs)*heavi(P_ic-P_vs))
******
_eK
D_eK = + k_Kinj- 2.0*Vmax_KATP*(_ATP/(Km_KATPA + _ATP)*(_eK/_K)/(Km_KATPK + (_eK/_K))*(_Na/_eNa)/(Km_KATPNa + (_Na/_eNa)))/Vol_ecs- k_ecKK*(_eK - _K)/Vol_ecs+ Vmax1_KK*K/(Km_KK + K)/Vol_ecs- Vmax2_KK*_eK/(Km_KK + _eK)/Vol_ecs- _eK/Vol_ecs*(G_o*(P_ic-P_vs)*heavi(P_ic-P_vs))
******
X
DX = -(k_XH*X*Hy - k_nXH*X_H)-(k_XO4*X*pow(O2, 2.5) - k_nXO4*X_O4) + min(2.0*q*(X_a - X), q*X_a)/Vol_c
******
Hy
DHy = - n_H*(k_XH*X*Hy - k_nXH*X_H)- n_H*(k_XO4H*X_O4*Hy - k_nXO4H*X_O4_H)+(k_CH*CO2 - k_nCH*Hy*BiC)+(k_LH*L - k_nLH*L0*Hy)+ Vmax_L0H*(_eL0*_eHy - L0*Hy)/((_eL0*_eHy + K_L0H*(1.0 + p_L0H))*(L0*Hy + K_L0H*(1.0 + p_L0H)) - pow(K_L0H, 2)*pow(p_L0H, 2))/Vol_c + min(2.0*q*(Hy_a - Hy), q*Hy_a)/Vol_c
******
X_H
DX_H = +(k_XH*X*Hy - k_nXH*X_H)-(k_XHO*X_H*pow(O2, 2.5) - k_nXHO*X_O4_H) + min(2.0*q*(X_H_a - X_H), q*X_H_a)/Vol_c
******
X_O4
DX_O4 = -(k_XO4H*X_O4*Hy - k_nXO4H*X_O4_H)+(k_XO4*X*pow(O2, 2.5) - k_nXO4*X_O4) + min(2.0*q*(X_O4_a - X_O4), q*X_O4_a)/Vol_c
******
X_O4_H
DX_O4_H = +(k_XO4H*X_O4*Hy - k_nXO4H*X_O4_H)+(k_XHO*X_H*pow(O2, 2.5) - k_nXHO*X_O4_H) + min(2.0*q*(X_O4_H_a - X_O4_H), q*X_O4_H_a)/Vol_c
******
O2
DO2 = - 4.0*(k_XO4*X*pow(O2, 2.5) - k_nXO4*X_O4)- 4.0*(k_XHO*X_H*pow(O2, 2.5) - k_nXHO*X_O4_H)- k_OO*(O2 - _eO2)/Vol_c + min(2.0*q*(O2_a - O2), q*O2_a)/Vol_c
******
CO2
DCO2 = -(k_CH*CO2 - k_nCH*Hy*BiC)-(k_CY*CO2*Y - k_nCY*CY)- k_CC*(CO2 - _eCO2)/Vol_c + min(2.0*q*(CO2_a - CO2), q*CO2_a)/Vol_c
******
BiC
DBiC = +(k_CH*CO2 - k_nCH*Hy*BiC)- Vmax1_BiC*BiC/(Km_BiC + BiC)/Vol_c+ Vmax2_BiC*_eBiC/(Km_BiC + _eBiC)/Vol_c + min(2.0*q*(BiC_a - BiC), q*BiC_a)/Vol_c
******
L
DL = -(k_LH*L - k_nLH*L0*Hy)- k_LL*(L - _eL)/Vol_c + min(2.0*q*(L_a - L), q*L_a)/Vol_c
******
L0
DL0 = +(k_LH*L - k_nLH*L0*Hy)+ Vmax_L0H*(_eL0*_eHy - L0*Hy)/((_eL0*_eHy + K_L0H*(1.0 + p_L0H))*(L0*Hy + K_L0H*(1.0 + p_L0H)) - pow(K_L0H, 2)*pow(p_L0H, 2))/Vol_c + min(2.0*q*(L0_a - L0), q*L0_a)/Vol_c
******
Y
DY = -(k_CY*CO2*Y - k_nCY*CY) + min(2.0*q*(Y_a - Y), q*Y_a)/Vol_c
******
CY
DCY = +(k_CY*CO2*Y - k_nCY*CY) + min(2.0*q*(CY_a - CY), q*CY_a)/Vol_c
******
X_v
DX_v = -(k_XH*X_v*Hy_v - k_nXH*X_H_v)-(k_XO4*X_v*pow(O2_v, 2.5) - k_nXO4*X_O4_v) + q*(max(2.0*X - X_a, 0.0) - X_v)/Vol_ven
******
Hy_v
DHy_v = - n_H*(k_XH*X_v*Hy_v - k_nXH*X_H_v)- n_H*(k_XO4H*X_O4_v*Hy_v - k_nXO4H*X_O4_H_v)+(k_CH*CO2_v - k_nCH*Hy_v*BiC_v)+(k_LH*L_v - k_nLH*L0_v*Hy_v) + q*(max(2.0*Hy - Hy_a, 0.0) - Hy_v)/Vol_ven
******
X_H_v
DX_H_v = +(k_XH*X_v*Hy_v - k_nXH*X_H_v)-(k_XHO*X_H_v*pow(O2_v, 2.5) - k_nXHO*X_O4_H_v) + q*(max(2.0*X_H - X_H_a, 0.0) - X_H_v)/Vol_ven
******
X_O4_v
DX_O4_v = -(k_XO4H*X_O4_v*Hy_v - k_nXO4H*X_O4_H_v)+(k_XO4*X_v*pow(O2_v, 2.5) - k_nXO4*X_O4_v) + q*(max(2.0*X_O4 - X_O4_a, 0.0) - X_O4_v)/Vol_ven
******
X_O4_H_v
DX_O4_H_v = +(k_XO4H*X_O4_v*Hy_v - k_nXO4H*X_O4_H_v)+(k_XHO*X_H_v*pow(O2_v, 2.5) - k_nXHO*X_O4_H_v) + q*(max(2.0*X_O4_H - X_O4_H_a, 0.0) - X_O4_H_v)/Vol_ven
******
O2_v
DO2_v = - 4.0*(k_XO4*X_v*pow(O2_v, 2.5) - k_nXO4*X_O4_v)- 4.0*(k_XHO*X_H_v*pow(O2_v, 2.5) - k_nXHO*X_O4_H_v) + q*(max(2.0*O2 - O2_a, 0.0) - O2_v)/Vol_ven
******
CO2_v
DCO2_v = -(k_CH*CO2_v - k_nCH*Hy_v*BiC_v)-(k_CY*CO2_v*Y_v - k_nCY*CY_v) + q*(max(2.0*CO2 - CO2_a, 0.0) - CO2_v)/Vol_ven
******
BiC_v
DBiC_v = +(k_CH*CO2_v - k_nCH*Hy_v*BiC_v) + q*(max(2.0*BiC - BiC_a, 0.0) - BiC_v)/Vol_ven
******
L_v
DL_v = -(k_LH*L_v - k_nLH*L0_v*Hy_v) + q*(max(2.0*L - L_a, 0.0) - L_v)/Vol_ven
******
L0_v
DL0_v = +(k_LH*L_v - k_nLH*L0_v*Hy_v) + q*(max(2.0*L0 - L0_a, 0.0) - L0_v)/Vol_ven
******
Y_v
DY_v = -(k_CY*CO2_v*Y_v - k_nCY*CY_v) + q*(max(2.0*Y - Y_a, 0.0) - Y_v)/Vol_ven
******
CY_v
DCY_v = +(k_CY*CO2_v*Y_v - k_nCY*CY_v) + q*(max(2.0*CY - CY_a, 0.0) - CY_v)/Vol_ven
******
_O2
D_O2 = - k_O2shut*(_O2 - _mO2)/Vol_exm+ k_ecOO*(_eO2 - _O2)/Vol_exm
******
_Na
D_Na = - 3.0*Vmax_KATP*(_ATP/(Km_KATPA + _ATP)*(_eK/_K)/(Km_KATPK + (_eK/_K))*(_Na/_eNa)/(Km_KATPNa + (_Na/_eNa)))/Vol_exm+ k_ecNa*(_eNa - _Na)/Vol_exm
******
_K
D_K = + 2.0*Vmax_KATP*(_ATP/(Km_KATPA + _ATP)*(_eK/_K)/(Km_KATPK + (_eK/_K))*(_Na/_eNa)/(Km_KATPNa + (_Na/_eNa)))/Vol_exm+ k_ecKK*(_eK - _K)/Vol_exm
******
_eNa
D_eNa = + 3.0*Vmax_KATP*(_ATP/(Km_KATPA + _ATP)*(_eK/_K)/(Km_KATPK + (_eK/_K))*(_Na/_eNa)/(Km_KATPNa + (_Na/_eNa)))/Vol_ecs- k_ecNa*(_eNa - _Na)/Vol_ecs+ k_Na*(Na - _eNa)/Vol_ecs- _eNa/Vol_ecs*(G_o*(P_ic-P_vs)*heavi(P_ic-P_vs))
******
_egluc
D_egluc = - k_ecgluc*(_egluc - _gluc)/Vol_ecs+ Vmax_gluc*(gluc - _egluc)/((gluc + K_gluc*(1.0 + p_gluc))*(_egluc + K_gluc*(1.0 + p_gluc)) - pow(K_gluc, 2)*pow(p_gluc, 2))/Vol_ecs
******
_eO2
D_eO2 = - k_ecOO*(_eO2 - _O2)/Vol_ecs+ k_OO*(O2 - _eO2)/Vol_ecs- _eO2/Vol_ecs*(G_o*(P_ic-P_vs)*heavi(P_ic-P_vs))
******
_ePy0
D_ePy0 = - k_ecPy0*(_ePy0 - _Py0)/Vol_ecs
******
_eAd
D_eAd = - k_ecAd*(_eAd - _Ad)/Vol_ecs
******
gluc
Dgluc = - Vmax_gluc*(gluc - _egluc)/((gluc + K_gluc*(1.0 + p_gluc))*(_egluc + K_gluc*(1.0 + p_gluc)) - pow(K_gluc, 2)*pow(p_gluc, 2))/Vol_c + min(2.0*q*(gluc_a - gluc), q*gluc_a)/Vol_c
******
K
DK = - Vmax1_KK*K/(Km_KK + K)/Vol_c+ Vmax2_KK*_eK/(Km_KK + _eK)/Vol_c + min(2.0*q*(K_a - K), q*K_a)/Vol_c
******
Na
DNa = - k_Na*(Na - _eNa)/Vol_c + min(2.0*q*(Na_a - Na), q*Na_a)/Vol_c
******
Ca_i1
DCa_i1 = + k_Cain1- Vmax_Caout1*Ca_i1/(Km_Caout + Ca_i1)
******
Ca_i2
DCa_i2 = + k_Cain2- Vmax_Caout2*Ca_i2/(Km_Caout + Ca_i2)
******
NO1
DNO1 = +(k_NOprod1 - k_NOout1*NO1)
******
NO2
DNO2 = +(k_NOprod2 - k_NOout2*NO2)
******
MLC1
DMLC1 = -(k_MLC1ac*MLC1 - k_MLC1inac*MLCp1)
******
MLCp1
DMLCp1 = +(k_MLC1ac*MLC1 - k_MLC1inac*MLCp1)
******
MLC2
DMLC2 = -(k_MLC2ac*MLC2 - k_MLC2inac*MLCp2)
******
MLCp2
DMLCp2 = +(k_MLC2ac*MLC2 - k_MLC2inac*MLCp2)
******
K_v
DK_v = + q*(max(2.0*K - K_a, 0.0) - K_v)/Vol_ven
******
gluc_v
Dgluc_v = + q*(max(2.0*gluc - gluc_a, 0.0) - gluc_v)/Vol_ven
******
Na_v
DNa_v = + q*(max(2.0*Na - Na_a, 0.0) - Na_v)/Vol_ven
******
q
2.0*G_2*(P_2-P_c) - q = 0
******
P_c
G_f*(P_c-P_ic)*heavi(P_c-P_ic) - 2.0*G_2*(P_2-P_c) + G_pv*(P_c-P_v) = 0
******
sigma_v1
2.0*(G_1*(P_a-P_1) - G_1*G_2/(G_1+G_2)*(P_1-P_2)) - 2.0*K_V1*r_1*r_01*sigma_v1/eta_1 = 0
******
sigma_v2
2.0*(G_1*G_2/(G_1+G_2)*(P_1-P_2) - G_2*(P_2-P_c)) - 2.0*K_V2*r_2*r_02*sigma_v2/eta_2 = 0
******
P_1
T_e1 + T_m1 + T_v1 - P_1*r_1 + P_ic*(r_1+h_1) = 0
******
P_2
T_e2 + T_m2 + T_v2 - P_2*r_2 + P_ic*(r_2+h_2) = 0
******
sigma_e1
sigma_e01*(exp(K_sigma1*(r_1-r_01)/r_01) - 1.0) - sigma_coll1 - sigma_e1 = 0
******
sigma_e2
sigma_e02*(exp(K_sigma2*(r_2-r_02)/r_02) - 1.0) - sigma_coll2 - sigma_e2 = 0
******
T_e1
sigma_e1*h_1 - T_e1 = 0
******
T_e2
sigma_e2*h_2 - T_e2 = 0
******
T_m1
T_max1*exp(-pow(fabs((r_1-r_m1)/(r_t1-r_m1)), n_m1)) - T_m1 = 0
******
T_m2
T_max2*exp(-pow(fabs((r_2-r_m2)/(r_t2-r_m2)), n_m2)) - T_m2 = 0
******
T_v1
sigma_v1*h_1 - T_v1 = 0
******
T_v2
sigma_v2*h_2 - T_v2 = 0
******
h_1
-r_1 + sqrt(r_1*r_1 + 2.0*r_01*h_01 + h_01*h_01) - h_1 = 0
******
h_2
-r_2 + sqrt(r_2*r_2 + 2.0*r_02*h_02 + h_02*h_02) - h_2 = 0
******
G_1
K_G1*pow(r_1, 4) - G_1 = 0
******
G_2
K_G2*pow(r_2, 4) - G_2 = 0
******
V_1
K_V1*r_1*r_1 - V_1 = 0
******
V_2
K_V2*r_2*r_2 - V_2 = 0
******
C_vi
1.0/(k_ven*(P_v-P_ic-P_vl)) - C_vi = 0
******
G_vs
G1_vs*(P_v-P_ic)/(max(P_v-P_vs, 0.001)) - G_vs = 0
******
r_mca
r_mcan*(1.0/K_mca*log((P_a-P_ic)/(P_an-P_icn))+1.0) - r_mca = 0
******
V_mca
******
T_max1
T_max01*(k_aut1*M_1 + (1.0 - k_aut1)) - T_max1 = 0
******
T_max2
T_max02*(k_aut2*M_2 + (1.0 - k_aut2)) - T_max2 = 0
******
M_1
K1_MLC*pow((MLCp1/MLCp1n), n_MLC)/(pow(K2_MLC, n_MLC) + pow((MLCp1/MLCp1n), n_MLC)) - M_1 = 0
******
M_2
K1_MLC*pow((MLCp2/MLCp2n), n_MLC)/(pow(K2_MLC, n_MLC) + pow((MLCp2/MLCp2n), n_MLC)) - M_2 = 0
******
Vmax_glucN
Vmax_glucNn*(Inh_glucN + 1.0)/(1.0 + Inh_glucN*(_ATP/_ATPn)*(_AMPn/_AMP)) - Vmax_glucN = 0
******
PMF
PD_mit + C_elec*log(_Hy/_mHy) - PMF = 0
******
PD_mit
PD_mitn/pH_mc_diff*(log10(_Hy/_mHy)) - PD_mit = 0
******
gKpot1
gKpot1tot*(m_K*pow((_eK/_eKn), n_K) + c_K)/(1.0 + m_K*pow((_eK/_eKn), n_K) + c_K) - gKpot1 = 0
******
gKpot2
gKpot2tot*(m_K*pow((_eK/_eKn), n_K) + c_K)/(1.0 + m_K*pow((_eK/_eKn), n_K) + c_K) - gKpot2 = 0
******
gKV1
gKV1tot*(1.0/(1.0 + exp(-(V_inact - PD_mem1)/kV_inact)))*(1.0/(1.0 + exp((V_act - PD_mem1)/kV_act))) - gKV1 = 0
******
gKV2
gKV2tot*(1.0/(1.0 + exp(-(V_inact - PD_mem2)/kV_inact)))*(1.0/(1.0 + exp((V_act - PD_mem2)/kV_act))) - gKV2 = 0
******
gKATP1
gKATP1tot*(m_pH*_eHy/_eHyn + m_Ad*_eAd/_eAdn + c_pHAd)/(1.0 + m_pH*_eHy/_eHyn + m_Ad*_eAd/_eAdn + c_pHAd) - gKATP1 = 0
******
gKATP2
gKATP2tot*(m_pH*_eHy/_eHyn + m_Ad*_eAd/_eAdn + c_pHAd)/(1.0 + m_pH*_eHy/_eHyn + m_Ad*_eAd/_eAdn + c_pHAd) - gKATP2 = 0
******
gKpres1
gKpres1tot*(m_NO*NO1/NO1n + m_Ca*Ca_i1/Ca_i1n + c_V*(exp((PD_mem1 - PD_mem1n)/kV_gKCa)))/(c_pres*exp((P_1 - P_ic - (P_1n - P_icn))/m_pres) + m_NO*NO1/NO1n + m_Ca*Ca_i1/Ca_i1n + c_V*(exp((PD_mem1 - PD_mem1n)/kV_gKCa))) - gKpres1 = 0
******
gKpres2
gKpres2tot*(m_NO*NO2/NO2n + m_Ca*Ca_i2/Ca_i2n + c_V*(exp((PD_mem2 - PD_mem2n)/kV_gKCa)))/(c_pres*exp((P_2 - P_ic - (P_2n - P_icn))/m_pres) + m_NO*NO2/NO2n + m_Ca*Ca_i2/Ca_i2n + c_V*(exp((PD_mem2 - PD_mem2n)/kV_gKCa))) - gKpres2 = 0
******
gK1
gKpot1 + gKpres1 + gKATP1 + gKV1 - gK1 = 0
******
gK2
gKpot2 + gKpres2 + gKATP2 + gKV2 - gK2 = 0
******
PD_mem1
(PD_K + gNa/gK1*PD_Na + gCl/gK1*PD_Cl)/(1.0 + gNa/gK1 + gCl/gK1) - PD_mem1 = 0
******
PD_mem2
(PD_K + gNa/gK2*PD_Na + gCl/gK2*PD_Cl)/(1.0 + gNa/gK2 + gCl/gK2) - PD_mem2 = 0
******
gCa1
gCa1_tot/(1.0 + exp((V_Ca - PD_mem1)/kV_Ca)) - gCa1 = 0
******
gCa2
gCa2_tot/(1.0 + exp((V_Ca - PD_mem2)/kV_Ca)) - gCa2 = 0
******
k_Cain1
gCa1*(PD_Ca1-PD_mem1) - k_Cain1 = 0
******
k_Cain2
gCa2*(PD_Ca2-PD_mem2) - k_Cain2 = 0
******
k_NOprod1
Inh_NO*k_NOprod1n*((1 - k_NOpH_frac1 - k_NOq_frac1) + k_NOq_frac1*((P_a - P_1)/(P_an - P_1n))*(r_1/r_1n) + k_NOpH_frac1*_eHy/_eHyn) - k_NOprod1 = 0
******
k_NOprod2
Inh_NO*k_NOprod2n*((1 - k_NOpH_frac2 - k_NOq_frac2) + k_NOq_frac2*((P_2 - P_c)/(P_2n - P_cn))*(r_2/r_2n) + k_NOpH_frac2*_eHy/_eHyn) - k_NOprod2 = 0
******
k_MLC1inac
k_MLC1inac0 + K1_NO1*pow(NO1/NO1n, n_NO)/(pow(K2_NO, n_NO) + pow(NO1/NO1n, n_NO)) - k_MLC1inac = 0
******
k_MLC2inac
k_MLC2inac0 + K1_NO2*pow(NO2/NO2n, n_NO)/(pow(K2_NO, n_NO) + pow(NO2/NO2n, n_NO)) - k_MLC2inac = 0
******
k_MLC1ac
K1_CaM1*pow((Ca_i1/Ca_i1n), n_Ca)/(pow(K2_CaM, n_Ca) + pow((Ca_i1/Ca_i1n), n_Ca)) - k_MLC1ac = 0
******
k_MLC2ac
K1_CaM2*pow((Ca_i2/Ca_i2n), n_Ca)/(pow(K2_CaM, n_Ca) + pow((Ca_i2/Ca_i2n), n_Ca)) - k_MLC2ac = 0
******
k_OO
k_OOn - k_OO = 0
******
k_CC
k_CCn - k_CC = 0
******
PD_K
C_elec*log(_eK/_K_musn) - PD_K = 0
******
PD_Ca1
C_elec/2.0*log(Ca_en/Ca_i1) - PD_Ca1 = 0
******
PD_Ca2
C_elec/2.0*log(Ca_en/Ca_i2) - PD_Ca2 = 0
******
Vol_art
V_1 + V_2 - Vol_art = 0
******
Vol_ven
Vol_ven00 + (C_vi*(P_v-P_ic) + C_ve*P_vs) - Vol_ven = 0
******