C
C There are a total of 22 entries in the algebraic variable array.
C There are a total of 13 entries in each of the rate and state variable arrays.
C There are a total of 67 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C STATES(1) is GlcI in component GlcI (millimolar).
C ALGBRC(6) is vHK in component vHK (flux).
C ALGBRC(5) is vGlcTr in component vGlcTr (flux).
C STATES(2) is hexose_P in component hexose_P (millimolar).
C ALGBRC(7) is vPFK in component vPFK (flux).
C STATES(3) is Fru16BP in component Fru16BP (millimolar).
C ALGBRC(15) is vALD in component vALD (flux).
C STATES(4) is triose_P in component triose_P (millimolar).
C ALGBRC(17) is vGAPDH in component vGAPDH (flux).
C ALGBRC(18) is vGDH in component vGDH (flux).
C ALGBRC(9) is vGPO in component vGPO (flux).
C STATES(5) is BPGA13 in component BPGA13 (millimolar).
C ALGBRC(21) is vPGK in component vPGK (flux).
C STATES(6) is N in component N (millimolar).
C ALGBRC(22) is vPK in component vPK (flux).
C STATES(7) is Pyr in component Pyr (millimolar).
C ALGBRC(10) is vPyrTr in component vPyrTr (flux).
C STATES(8) is NADH in component NADH (millimolar).
C STATES(9) is NAD in component NAD (millimolar).
C STATES(10) is Gly3P in component Gly3P (millimolar).
C ALGBRC(12) is vGlyK in component vGlyK (flux).
C STATES(11) is Gly in component Gly (millimolar).
C STATES(12) is P in component P (millimolar).
C ALGBRC(11) is vATPase in component vATPase (flux).
C ALGBRC(1) is ATP in component ATP (millimolar).
C CONSTS(1) is sumA in component ATP (millimolar).
C CONSTS(2) is Keq_AK in component ATP (dimensionless).
C ALGBRC(2) is ADP in component ADP (millimolar).
C ALGBRC(13) is DHAP in component DHAP (millimolar).
C ALGBRC(4) is Fru6P in component Fru6P (millimolar).
C ALGBRC(14) is GAP in component GAP (millimolar).
C STATES(13) is Glc6P in component Glc6P (millimolar).
C CONSTS(3) is sumc4 in component DHAP (millimolar).
C CONSTS(4) is sumc5 in component DHAP (millimolar).
C ALGBRC(3) is GlcE in component GlcE (millimolar).
C ALGBRC(8) is vPGI in component vPGI (flux).
C ALGBRC(19) is PGA3 in component PGA3 (millimolar).
C ALGBRC(20) is PEP in component PEP (millimolar).
C CONSTS(5) is Keq_ENO in component PEP (dimensionless).
C CONSTS(6) is Keq_PGM in component PEP (dimensionless).
C CONSTS(7) is K_Glc in component vGlcTr (millimolar).
C CONSTS(8) is alpha in component vGlcTr (dimensionless).
C CONSTS(9) is vGlcTr_max in component vGlcTr (flux).
C CONSTS(10) is K_GlcI in component vHK (millimolar).
C CONSTS(11) is K_Glc6P in component vHK (millimolar).
C CONSTS(12) is K_ATP in component vHK (millimolar).
C CONSTS(13) is K_ADP in component vHK (millimolar).
C CONSTS(14) is vHK_max in component vHK (flux).
C CONSTS(15) is K_Glc6P in component vPGI (millimolar).
C CONSTS(16) is K_Fru6P in component vPGI (millimolar).
C CONSTS(17) is vPGI_max in component vPGI (flux).
C CONSTS(18) is Ki_1 in component vPFK (millimolar).
C CONSTS(19) is Ki_2 in component vPFK (millimolar).
C CONSTS(20) is KM_Fru6P in component vPFK (millimolar).
C CONSTS(21) is KM_ATP in component vPFK (millimolar).
C CONSTS(22) is vPFK_max in component vPFK (flux).
C CONSTS(23) is sumA in component vALD (millimolar).
C CONSTS(24) is KM_GAP in component vALD (millimolar).
C CONSTS(25) is Ki_GAP in component vALD (millimolar).
C CONSTS(26) is KM_DHAP in component vALD (millimolar).
C CONSTS(27) is vALD_max_forward in component vALD (flux).
C CONSTS(28) is vALD_max_reverse in component vALD (flux).
C ALGBRC(16) is vTPI in component vTPI (flux).
C CONSTS(29) is K_DHAP in component vTPI (millimolar).
C CONSTS(30) is K_GAP in component vTPI (millimolar).
C CONSTS(31) is vTPI_max in component vTPI (flux).
C CONSTS(32) is K_NAD in component vGAPDH (millimolar).
C CONSTS(33) is K_GAP in component vGAPDH (millimolar).
C CONSTS(34) is K_BPGA13 in component vGAPDH (millimolar).
C CONSTS(35) is K_NADH in component vGAPDH (millimolar).
C CONSTS(36) is vGAPDH_max_forward in component vGAPDH (flux).
C CONSTS(37) is vGAPDH_max_reverse in component vGAPDH (flux).
C CONSTS(38) is vGAPDH_max in component vGAPDH (dimensionless).
C CONSTS(39) is K_NADH in component vGDH (millimolar).
C CONSTS(40) is K_Gly3P in component vGDH (millimolar).
C CONSTS(41) is K_DHAP in component vGDH (millimolar).
C CONSTS(42) is K_NAD in component vGDH (millimolar).
C CONSTS(43) is vGDH_max_forward in component vGDH (flux).
C CONSTS(44) is vGDH_max_reverse in component vGDH (flux).
C CONSTS(45) is vGDH_max in component vGDH (dimensionless).
C CONSTS(46) is K_Gly3P in component vGPO (millimolar).
C CONSTS(47) is vGPO_max in component vGPO (flux).
C CONSTS(48) is K_pyruvate in component vPyrTr (millimolar).
C CONSTS(49) is vPyrTr_max in component vPyrTr (flux).
C CONSTS(50) is K_ADP in component vPGK (millimolar).
C CONSTS(51) is K_BPGA13 in component vPGK (millimolar).
C CONSTS(52) is K_PGA3 in component vPGK (millimolar).
C CONSTS(53) is K_ATP in component vPGK (millimolar).
C CONSTS(54) is vPGK_max_forward in component vPGK (flux).
C CONSTS(55) is vPGK_max_reverse in component vPGK (flux).
C CONSTS(56) is vPGK_max in component vPGK (dimensionless).
C CONSTS(57) is KM_ADP in component vPK (millimolar).
C CONSTS(58) is n in component vPK (dimensionless).
C CONSTS(59) is vPK_max in component vPK (flux).
C CONSTS(60) is k in component vATPase (flux).
C CONSTS(61) is K_ADP in component vGlyK (millimolar).
C CONSTS(62) is K_Gly3P in component vGlyK (millimolar).
C CONSTS(63) is K_Gly in component vGlyK (millimolar).
C CONSTS(64) is K_ATP in component vGlyK (millimolar).
C CONSTS(65) is vGlyK_max_forward in component vGlyK (flux).
C CONSTS(66) is vGlyK_max_reverse in component vGlyK (flux).
C CONSTS(67) is vGlyK_max in component vGlyK (dimensionless).
C RATES(1) is d/dt GlcI in component GlcI (millimolar).
C RATES(2) is d/dt hexose_P in component hexose_P (millimolar).
C RATES(3) is d/dt Fru16BP in component Fru16BP (millimolar).
C RATES(4) is d/dt triose_P in component triose_P (millimolar).
C RATES(5) is d/dt BPGA13 in component BPGA13 (millimolar).
C RATES(6) is d/dt N in component N (millimolar).
C RATES(7) is d/dt Pyr in component Pyr (millimolar).
C RATES(8) is d/dt NADH in component NADH (millimolar).
C RATES(9) is d/dt NAD in component NAD (millimolar).
C RATES(10) is d/dt Gly3P in component Gly3P (millimolar).
C RATES(11) is d/dt Gly in component Gly (millimolar).
C RATES(12) is d/dt P in component P (millimolar).
C RATES(13) is d/dt Glc6P in component Glc6P (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.0340009
      STATES(2) = 2.583763
      STATES(3) = 16.5371
      STATES(4) = 3.9391429
      STATES(5) = 0.0326745
      STATES(6) = 1.59603
      STATES(7) = 4.77413
      STATES(8) = 0.0448639
      STATES(9) = 0.0448639
      STATES(10) = 0.0
      STATES(11) = 0.0
      STATES(12) = 7.63936
      CONSTS(1) = 3.9
      CONSTS(2) = 0.442
      STATES(13) = 2.07199
      CONSTS(3) = 45.0
      CONSTS(4) = 5.0
      CONSTS(5) = 6.7
      CONSTS(6) = 0.187
      CONSTS(7) = 2.0
      CONSTS(8) = 0.75
      CONSTS(9) = 106.2
      CONSTS(10) = 0.1
      CONSTS(11) = 12.0
      CONSTS(12) = 0.116
      CONSTS(13) = 0.126
      CONSTS(14) = 625.0
      CONSTS(15) = 0.4
      CONSTS(16) = 0.12
      CONSTS(17) = 848.0
      CONSTS(18) = 15.8
      CONSTS(19) = 10.7
      CONSTS(20) = 0.82
      CONSTS(21) = 0.026
      CONSTS(22) = 780.0
      CONSTS(23) = 6.0
      CONSTS(24) = 0.067
      CONSTS(25) = 0.098
      CONSTS(26) = 0.015
      CONSTS(27) = 184.5
      CONSTS(28) = 219.555
      CONSTS(29) = 1.2
      CONSTS(30) = 0.25
      CONSTS(31) = 842.0
      CONSTS(32) = 0.45
      CONSTS(33) = 0.15
      CONSTS(34) = 0.1
      CONSTS(35) = 0.02
      CONSTS(36) = 1470.0
      CONSTS(37) = 984.9
      CONSTS(38) = 1.0
      CONSTS(39) = 0.01
      CONSTS(40) = 2.0
      CONSTS(41) = 0.1
      CONSTS(42) = 0.4
      CONSTS(43) = 533.0
      CONSTS(44) = 149.24
      CONSTS(45) = 1.0
      CONSTS(46) = 1.7
      CONSTS(47) = 368.0
      CONSTS(48) = 1.96
      CONSTS(49) = 200.0
      CONSTS(50) = 0.1
      CONSTS(51) = 0.05
      CONSTS(52) = 1.62
      CONSTS(53) = 0.29
      CONSTS(54) = 640.0
      CONSTS(55) = 18.56
      CONSTS(56) = 1.0
      CONSTS(57) = 0.114
      CONSTS(58) = 2.5
      CONSTS(59) = 2600
      CONSTS(60) = 50
      CONSTS(61) = 0.12
      CONSTS(62) = 5.1
      CONSTS(63) = 0.12
      CONSTS(64) = 0.19
      CONSTS(65) = 220.0
      CONSTS(66) = 334000.0
      CONSTS(67) = 1.0
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = (( STATES(12)*(1.00000 -  4.00000*CONSTS(2)) - CONSTS(1))+CONSTS(1) -  (1.00000 -  4.00000*CONSTS(2))*STATES(12) ** 2.00000+ 4.00000*(1.00000 -  4.00000*CONSTS(2))*( - CONSTS(2)*STATES(12) ** 2.00000) ** 0.500000)/( 2.00000*(1.00000 -  4.00000*CONSTS(2)))
      ALGBRC(2) = STATES(12) -  2.00000*ALGBRC(1)
      ALGBRC(6) = ( CONSTS(14)*STATES(1)*ALGBRC(1))/( CONSTS(12)*CONSTS(10)*(1.00000+STATES(13)/CONSTS(11)+STATES(1)/CONSTS(10))*(1.00000+ALGBRC(1)/CONSTS(12)+ALGBRC(2)/CONSTS(13)))
      ALGBRC(3) = TERNRY(VOI.GE.60.0000.AND.VOI.LT.61.0000, 5.00000, 0.0500000)
      ALGBRC(5) =  CONSTS(9)*((ALGBRC(3) - STATES(1))/(CONSTS(7)+ALGBRC(3)+STATES(1)+ CONSTS(8)*ALGBRC(3)*(STATES(1)/CONSTS(7))))
      RATES(1) = ALGBRC(5) - ALGBRC(6)
      ALGBRC(4) = STATES(2) - STATES(13)
      ALGBRC(7) = ( CONSTS(18)*CONSTS(22)*ALGBRC(4)*ALGBRC(1))/( CONSTS(21)*CONSTS(20)*(STATES(3)+CONSTS(18))*(1.00000+STATES(3)/CONSTS(19)+ALGBRC(4)/CONSTS(20))*(1.00000+ALGBRC(1)/CONSTS(21)))
      RATES(2) = ALGBRC(6) - ALGBRC(7)
      ALGBRC(8) = ( CONSTS(17)*(STATES(13)/CONSTS(15) - ALGBRC(4)/CONSTS(16)))/(1.00000+STATES(13)/CONSTS(15)+ALGBRC(4)/CONSTS(16))
      RATES(13) = ALGBRC(6) - ALGBRC(8)
      ALGBRC(12) = ( CONSTS(67)*(( CONSTS(65)*ALGBRC(2)*STATES(10))/( CONSTS(61)*CONSTS(62)) - ( CONSTS(66)*ALGBRC(1)*STATES(11))/( CONSTS(64)*CONSTS(63))))/( (1.00000+STATES(10)/CONSTS(62)+STATES(11)/CONSTS(63))*(1.00000+ALGBRC(2)/CONSTS(61)+ALGBRC(1)/CONSTS(64)))
      RATES(11) = ALGBRC(12)
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
      ALGBRC(15) = (( CONSTS(27)*STATES(3))/( 0.00900000*(1.00000+ALGBRC(1)/0.680000+ALGBRC(2)/1.51000+(CONSTS(23) - (ALGBRC(1)+ALGBRC(2)))/3.65000)) - ( CONSTS(28)*ALGBRC(14)*ALGBRC(13))/( CONSTS(26)*CONSTS(24)))/(1.00000+ALGBRC(14)/CONSTS(24)+ALGBRC(13)/CONSTS(26)+( ALGBRC(14)*ALGBRC(13))/( CONSTS(26)*CONSTS(24))+STATES(3)/( 0.00900000*(1.00000+ALGBRC(1)/0.680000+ALGBRC(2)/1.51000+(CONSTS(23) - (ALGBRC(1)+ALGBRC(2)))/3.65000))+( STATES(3)*ALGBRC(14))/( CONSTS(25)*0.00900000*(1.00000+ALGBRC(1)/0.680000+ALGBRC(2)/1.51000+(CONSTS(23) - (ALGBRC(1)+ALGBRC(2)))/3.65000)))
      RATES(3) = ALGBRC(7) - ALGBRC(15)
      ALGBRC(17) =  CONSTS(38)*(( CONSTS(36)*( ALGBRC(14)*((STATES(9)/CONSTS(33))/CONSTS(32)) -  (CONSTS(37)/CONSTS(36))*((( STATES(5)*STATES(8))/CONSTS(34))/CONSTS(35))))/( (1.00000+ALGBRC(14)/CONSTS(33)+STATES(5)/CONSTS(34))*(1.00000+STATES(9)/CONSTS(32)+STATES(8)/CONSTS(35))))
      ALGBRC(18) = ( CONSTS(45)*CONSTS(43)*(( STATES(8)*ALGBRC(13))/( CONSTS(39)*CONSTS(41)) - ( CONSTS(44)*STATES(9)*STATES(10))/( CONSTS(40)*CONSTS(42)*CONSTS(43))))/( (1.00000+STATES(9)/CONSTS(42)+STATES(8)/CONSTS(39))*(1.00000+ALGBRC(13)/CONSTS(41)+STATES(10)/CONSTS(40)))
      ALGBRC(9) = ( CONSTS(47)*STATES(10))/(STATES(10)+CONSTS(46))
      RATES(4) = ( 2.00000*ALGBRC(15)+ALGBRC(9)) - (ALGBRC(17)+ALGBRC(18))
      RATES(8) = ALGBRC(17) - ALGBRC(18)
      RATES(9) = ALGBRC(18) - ALGBRC(17)
      RATES(10) = ALGBRC(18) - (ALGBRC(12)+ALGBRC(9))
CALL minimize(minfunc_1, CONSTS, VARIABLES, <VAR>)
      ALGBRC(21) = ( CONSTS(56)*CONSTS(54)*(( - CONSTS(55)*ALGBRC(19)*ALGBRC(1))/( CONSTS(53)*CONSTS(52)*CONSTS(54))+( STATES(5)*ALGBRC(2))/( CONSTS(51)*CONSTS(50))))/( (1.00000+STATES(5)/CONSTS(51)+ALGBRC(19)/CONSTS(52))*(1.00000+ALGBRC(2)/CONSTS(50)+ALGBRC(1)/CONSTS(53)))
      RATES(5) = ALGBRC(17) - ALGBRC(21)
      ALGBRC(22) = (( CONSTS(59)*ALGBRC(20)/( 0.340000*(1.00000+ALGBRC(1)/0.570000+ALGBRC(2)/0.640000)) ** CONSTS(58)*ALGBRC(2))/CONSTS(57))/( (1.00000+ALGBRC(20)/( 0.340000*(1.00000+ALGBRC(1)/0.570000+ALGBRC(2)/0.640000)) ** CONSTS(58))*(1.00000+ALGBRC(2)/CONSTS(57)))
      RATES(6) = ALGBRC(21) - ALGBRC(22)
      ALGBRC(10) = (( CONSTS(49)*STATES(7))/CONSTS(48))/(1.00000+STATES(7)/CONSTS(48))
      RATES(7) = ALGBRC(22) - ALGBRC(10)
      ALGBRC(11) = ( CONSTS(60)*ALGBRC(1))/ALGBRC(2)
      RATES(12) = (ALGBRC(21)+ALGBRC(12)+ALGBRC(22)) - (ALGBRC(6)+ALGBRC(7)+ALGBRC(11))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = (( STATES(12)*(1.00000 -  4.00000*CONSTS(2)) - CONSTS(1))+CONSTS(1) -  (1.00000 -  4.00000*CONSTS(2))*STATES(12) ** 2.00000+ 4.00000*(1.00000 -  4.00000*CONSTS(2))*( - CONSTS(2)*STATES(12) ** 2.00000) ** 0.500000)/( 2.00000*(1.00000 -  4.00000*CONSTS(2)))
      ALGBRC(2) = STATES(12) -  2.00000*ALGBRC(1)
      ALGBRC(6) = ( CONSTS(14)*STATES(1)*ALGBRC(1))/( CONSTS(12)*CONSTS(10)*(1.00000+STATES(13)/CONSTS(11)+STATES(1)/CONSTS(10))*(1.00000+ALGBRC(1)/CONSTS(12)+ALGBRC(2)/CONSTS(13)))
      ALGBRC(3) = TERNRY(VOI.GE.60.0000.AND.VOI.LT.61.0000, 5.00000, 0.0500000)
      ALGBRC(5) =  CONSTS(9)*((ALGBRC(3) - STATES(1))/(CONSTS(7)+ALGBRC(3)+STATES(1)+ CONSTS(8)*ALGBRC(3)*(STATES(1)/CONSTS(7))))
      ALGBRC(4) = STATES(2) - STATES(13)
      ALGBRC(7) = ( CONSTS(18)*CONSTS(22)*ALGBRC(4)*ALGBRC(1))/( CONSTS(21)*CONSTS(20)*(STATES(3)+CONSTS(18))*(1.00000+STATES(3)/CONSTS(19)+ALGBRC(4)/CONSTS(20))*(1.00000+ALGBRC(1)/CONSTS(21)))
      ALGBRC(8) = ( CONSTS(17)*(STATES(13)/CONSTS(15) - ALGBRC(4)/CONSTS(16)))/(1.00000+STATES(13)/CONSTS(15)+ALGBRC(4)/CONSTS(16))
      ALGBRC(12) = ( CONSTS(67)*(( CONSTS(65)*ALGBRC(2)*STATES(10))/( CONSTS(61)*CONSTS(62)) - ( CONSTS(66)*ALGBRC(1)*STATES(11))/( CONSTS(64)*CONSTS(63))))/( (1.00000+STATES(10)/CONSTS(62)+STATES(11)/CONSTS(63))*(1.00000+ALGBRC(2)/CONSTS(61)+ALGBRC(1)/CONSTS(64)))
      ALGBRC(15) = (( CONSTS(27)*STATES(3))/( 0.00900000*(1.00000+ALGBRC(1)/0.680000+ALGBRC(2)/1.51000+(CONSTS(23) - (ALGBRC(1)+ALGBRC(2)))/3.65000)) - ( CONSTS(28)*ALGBRC(14)*ALGBRC(13))/( CONSTS(26)*CONSTS(24)))/(1.00000+ALGBRC(14)/CONSTS(24)+ALGBRC(13)/CONSTS(26)+( ALGBRC(14)*ALGBRC(13))/( CONSTS(26)*CONSTS(24))+STATES(3)/( 0.00900000*(1.00000+ALGBRC(1)/0.680000+ALGBRC(2)/1.51000+(CONSTS(23) - (ALGBRC(1)+ALGBRC(2)))/3.65000))+( STATES(3)*ALGBRC(14))/( CONSTS(25)*0.00900000*(1.00000+ALGBRC(1)/0.680000+ALGBRC(2)/1.51000+(CONSTS(23) - (ALGBRC(1)+ALGBRC(2)))/3.65000)))
      ALGBRC(17) =  CONSTS(38)*(( CONSTS(36)*( ALGBRC(14)*((STATES(9)/CONSTS(33))/CONSTS(32)) -  (CONSTS(37)/CONSTS(36))*((( STATES(5)*STATES(8))/CONSTS(34))/CONSTS(35))))/( (1.00000+ALGBRC(14)/CONSTS(33)+STATES(5)/CONSTS(34))*(1.00000+STATES(9)/CONSTS(32)+STATES(8)/CONSTS(35))))
      ALGBRC(18) = ( CONSTS(45)*CONSTS(43)*(( STATES(8)*ALGBRC(13))/( CONSTS(39)*CONSTS(41)) - ( CONSTS(44)*STATES(9)*STATES(10))/( CONSTS(40)*CONSTS(42)*CONSTS(43))))/( (1.00000+STATES(9)/CONSTS(42)+STATES(8)/CONSTS(39))*(1.00000+ALGBRC(13)/CONSTS(41)+STATES(10)/CONSTS(40)))
      ALGBRC(9) = ( CONSTS(47)*STATES(10))/(STATES(10)+CONSTS(46))
      ALGBRC(21) = ( CONSTS(56)*CONSTS(54)*(( - CONSTS(55)*ALGBRC(19)*ALGBRC(1))/( CONSTS(53)*CONSTS(52)*CONSTS(54))+( STATES(5)*ALGBRC(2))/( CONSTS(51)*CONSTS(50))))/( (1.00000+STATES(5)/CONSTS(51)+ALGBRC(19)/CONSTS(52))*(1.00000+ALGBRC(2)/CONSTS(50)+ALGBRC(1)/CONSTS(53)))
      ALGBRC(22) = (( CONSTS(59)*ALGBRC(20)/( 0.340000*(1.00000+ALGBRC(1)/0.570000+ALGBRC(2)/0.640000)) ** CONSTS(58)*ALGBRC(2))/CONSTS(57))/( (1.00000+ALGBRC(20)/( 0.340000*(1.00000+ALGBRC(1)/0.570000+ALGBRC(2)/0.640000)) ** CONSTS(58))*(1.00000+ALGBRC(2)/CONSTS(57)))
      ALGBRC(10) = (( CONSTS(49)*STATES(7))/CONSTS(48))/(1.00000+STATES(7)/CONSTS(48))
      ALGBRC(11) = ( CONSTS(60)*ALGBRC(1))/ALGBRC(2)
      ALGBRC(16) = ( CONSTS(31)*(ALGBRC(13)/CONSTS(29) - ( 5.70000*ALGBRC(14))/CONSTS(30)))/(1.00000+ALGBRC(14)/CONSTS(30)+ALGBRC(13)/CONSTS(29))
      RETURN
      END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_0[1] = abs(ALGBRC(13) - ( CONSTS(4)*ALGBRC(13))/((CONSTS(3)+CONSTS(4)) - (STATES(5)+ 2.00000*STATES(3)+ALGBRC(4)+ALGBRC(14)+STATES(13)+STATES(12))))
      minfunc_0[2] = abs(ALGBRC(14) - (STATES(4) - ALGBRC(13)))
      RETURN
      END
REAL FUNCTION minfunc_1(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_1[1] = abs(ALGBRC(19) - (STATES(6) - ALGBRC(20)))
      minfunc_1[2] = abs(ALGBRC(20) -  CONSTS(5)*CONSTS(6)*ALGBRC(19))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END