C
C There are a total of 32 entries in the algebraic variable array.
C There are a total of 24 entries in each of the rate and state variable arrays.
C There are a total of 142 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is k_rel_TPI in component parameters (dimensionless).
C CONSTS(2) is k_rel_GAPDH in component parameters (dimensionless).
C CONSTS(3) is SUMAXP in component parameters (mM).
C CONSTS(4) is CO2 in component extracellular (mM).
C CONSTS(5) is ETOH in component extracellular (mM).
C CONSTS(6) is SUCC in component extracellular (mM).
C CONSTS(7) is GLY in component extracellular (mM).
C CONSTS(8) is GLCo in component extracellular (mM).
C CONSTS(9) is cytoplasm in component cytoplasm (liter).
C STATES(1) is P in component cytoplasm (mM).
C STATES(2) is G6P in component cytoplasm (mM).
C STATES(3) is F6P in component cytoplasm (mM).
C STATES(4) is F16P in component cytoplasm (mM).
C STATES(5) is NADH in component cytoplasm (mM).
C STATES(6) is NAD in component cytoplasm (mM).
C STATES(7) is BPG in component cytoplasm (mM).
C STATES(8) is P3G in component cytoplasm (mM).
C STATES(9) is P2G in component cytoplasm (mM).
C STATES(10) is PEP in component cytoplasm (mM).
C STATES(11) is PYR in component cytoplasm (mM).
C STATES(12) is ACE in component cytoplasm (mM).
C CONSTS(10) is X in component cytoplasm (mM).
C STATES(13) is GA3P in component cytoplasm (mM).
C STATES(14) is DHAP in component cytoplasm (mM).
C STATES(15) is D6PGluconoLactone in component cytoplasm (mM).
C STATES(16) is D6PGluconate in component cytoplasm (mM).
C STATES(17) is NADP in component cytoplasm (mM).
C STATES(18) is NADPH in component cytoplasm (mM).
C STATES(19) is Ribulose5P in component cytoplasm (mM).
C STATES(20) is Ribose5P in component cytoplasm (mM).
C STATES(21) is Xyl5P in component cytoplasm (mM).
C STATES(22) is Seduhept7P in component cytoplasm (mM).
C STATES(23) is Erythrose4P in component cytoplasm (mM).
C STATES(24) is GLCi in component cytoplasm (mM).
C CONSTS(11) is F26BP in component cytoplasm (mM).
C ALGBRC(1) is vGLK in component vGLK (mmol_per_min).
C ALGBRC(3) is vPGI in component vPGI (mmol_per_min).
C ALGBRC(11) is vPFK in component vPFK (mmol_per_min).
C ALGBRC(12) is vALD in component vALD (mmol_per_min).
C ALGBRC(13) is vG3PDH in component vG3PDH (mmol_per_min).
C ALGBRC(14) is vGAPDH in component vGAPDH (mmol_per_min).
C ALGBRC(15) is vPGK in component vPGK (mmol_per_min).
C ALGBRC(17) is vPGM in component vPGM (mmol_per_min).
C ALGBRC(19) is vENO in component vENO (mmol_per_min).
C ALGBRC(21) is vPYK in component vPYK (mmol_per_min).
C ALGBRC(23) is vPDC in component vPDC (mmol_per_min).
C ALGBRC(27) is vSUC in component vSUC (mmol_per_min).
C ALGBRC(29) is vADH in component vADH (mmol_per_min).
C ALGBRC(24) is vATP in component vATP (mmol_per_min).
C ALGBRC(16) is vTPI in component vTPI (mmol_per_min).
C ALGBRC(18) is vG6PDH in component vG6PDH (mmol_per_min).
C ALGBRC(20) is v6PGL in component v6PGL (mmol_per_min).
C ALGBRC(22) is vGluDH in component vGluDH (mmol_per_min).
C ALGBRC(25) is vPPI in component vPPI (mmol_per_min).
C ALGBRC(28) is vTransk1 in component vTransk1 (mmol_per_min).
C ALGBRC(30) is vR5PI in component vR5PI (mmol_per_min).
C ALGBRC(31) is vTransald in component vTransald (mmol_per_min).
C ALGBRC(32) is vTransk2 in component vTransk2 (mmol_per_min).
C ALGBRC(26) is vNADPH in component vNADPH (mmol_per_min).
C ALGBRC(4) is vGLT in component vGLT (mmol_per_min).
C ALGBRC(2) is ratio_NADPH_NADP in component rules (dimensionless).
C CONSTS(12) is VmGLK in component vGLK (mM_per_min).
C CONSTS(13) is KeqAK in component vGLK (dimensionless).
C CONSTS(14) is KeqGLK in component vGLK (dimensionless).
C CONSTS(15) is KmGLKATP in component vGLK (mM).
C CONSTS(16) is KmGLKGLCi in component vGLK (mM).
C CONSTS(17) is KmGLKG6P in component vGLK (mM).
C CONSTS(18) is KmGLKADP in component vGLK (mM).
C CONSTS(19) is VmPGI in component vPGI (mM_per_min).
C CONSTS(20) is KmPGIG6P in component vPGI (mM).
C CONSTS(21) is KeqPGI in component vPGI (dimensionless).
C CONSTS(22) is KmPGIF6P in component vPGI (mM).
C ALGBRC(5) is numerator1 in component vPFK (mmol_mM2_per_minute).
C ALGBRC(6) is numerator2 in component vPFK (dimensionless).
C ALGBRC(10) is denominator in component vPFK (mM2).
C ALGBRC(7) is denom1 in component vPFK (dimensionless).
C ALGBRC(8) is denom2 in component vPFK (dimensionless).
C ALGBRC(9) is denom3 in component vPFK (dimensionless).
C CONSTS(23) is gR in component vPFK (dimensionless).
C CONSTS(24) is VmPFK in component vPFK (mM_per_min).
C CONSTS(25) is KeqAK in component vPFK (dimensionless).
C CONSTS(26) is KmPFKF6P in component vPFK (mM).
C CONSTS(27) is KmPFKATP in component vPFK (mM).
C CONSTS(28) is L0 in component vPFK (dimensionless).
C CONSTS(29) is CPFKF26BP in component vPFK (dimensionless).
C CONSTS(30) is KPFKF26BP in component vPFK (mM).
C CONSTS(31) is CPFKF16BP in component vPFK (dimensionless).
C CONSTS(32) is KPFKF16BP in component vPFK (mM).
C CONSTS(33) is CPFKAMP in component vPFK (dimensionless).
C CONSTS(34) is KPFKAMP in component vPFK (mM).
C CONSTS(35) is CiPFKATP in component vPFK (dimensionless).
C CONSTS(36) is KiPFKATP in component vPFK (mM).
C CONSTS(37) is CPFKATP in component vPFK (dimensionless).
C CONSTS(38) is VmALD in component vALD (mM_per_min).
C CONSTS(39) is KeqTPI in component vALD (dimensionless).
C CONSTS(40) is KeqALD in component vALD (mM).
C CONSTS(41) is KmALDF16P in component vALD (mM).
C CONSTS(42) is KmALDDHAP in component vALD (mM).
C CONSTS(43) is KmALDGAP in component vALD (mM).
C CONSTS(44) is KmALDGAPi in component vALD (mM).
C CONSTS(45) is VmG3PDH in component vG3PDH (mM_per_min).
C CONSTS(46) is KeqG3PDH in component vG3PDH (dimensionless).
C CONSTS(47) is KeqTPI in component vG3PDH (dimensionless).
C CONSTS(48) is KmG3PDHDHAP in component vG3PDH (mM).
C CONSTS(49) is KmG3PDHNADH in component vG3PDH (mM).
C CONSTS(50) is KmG3PDHNAD in component vG3PDH (mM).
C CONSTS(51) is KmG3PDHGLY in component vG3PDH (mM).
C CONSTS(52) is VmGAPDHr in component vGAPDH (mM_per_min).
C CONSTS(53) is KmGAPDHBPG in component vGAPDH (mM).
C CONSTS(54) is KmGAPDHNADH in component vGAPDH (mM).
C CONSTS(55) is KeqTPI in component vGAPDH (dimensionless).
C CONSTS(56) is VmGAPDHf in component vGAPDH (mM_per_min).
C CONSTS(57) is KmGAPDHGAP in component vGAPDH (mM).
C CONSTS(58) is KmGAPDHNAD in component vGAPDH (mM).
C CONSTS(59) is KeqGAPDH in component vGAPDH (dimensionless).
C CONSTS(60) is VmPGK in component vPGK (mM_per_min).
C CONSTS(61) is KeqPGK in component vPGK (dimensionless).
C CONSTS(62) is KeqAK in component vPGK (dimensionless).
C CONSTS(63) is KmPGKATP in component vPGK (mM).
C CONSTS(64) is KmPGKP3G in component vPGK (mM).
C CONSTS(65) is KmPGKADP in component vPGK (mM).
C CONSTS(66) is KmPGKBPG in component vPGK (mM).
C CONSTS(67) is VmPGM in component vPGM (mM_per_min).
C CONSTS(68) is KmPGMP3G in component vPGM (mM).
C CONSTS(69) is KeqPGM in component vPGM (dimensionless).
C CONSTS(70) is KmPGMP2G in component vPGM (mM).
C CONSTS(71) is VmENO in component vENO (mM_per_min).
C CONSTS(72) is KmENOP2G in component vENO (mM).
C CONSTS(73) is KeqENO in component vENO (dimensionless).
C CONSTS(74) is KmENOPEP in component vENO (mM).
C CONSTS(75) is VmPYK in component vPYK (mM_per_min).
C CONSTS(76) is KmPYKPEP in component vPYK (mM).
C CONSTS(77) is KmPYKADP in component vPYK (mM).
C CONSTS(78) is KeqAK in component vPYK (dimensionless).
C CONSTS(79) is KeqPYK in component vPYK (dimensionless).
C CONSTS(80) is KmPYKPYR in component vPYK (mM).
C CONSTS(81) is KmPYKATP in component vPYK (mM).
C CONSTS(82) is VmPDC in component vPDC (mM_per_min).
C CONSTS(83) is nPDC in component vPDC (dimensionless).
C CONSTS(84) is KmPDCPYR in component vPDC (mM).
C CONSTS(85) is KSUCC in component vSUC (per_min).
C CONSTS(86) is VmADH in component vADH (mM_per_min).
C CONSTS(87) is KiADHNAD in component vADH (mM).
C CONSTS(88) is KmADHETOH in component vADH (mM).
C CONSTS(89) is KeqADH in component vADH (dimensionless).
C CONSTS(90) is KmADHNAD in component vADH (mM).
C CONSTS(91) is KmADHNADH in component vADH (mM).
C CONSTS(92) is KiADHNADH in component vADH (mM).
C CONSTS(93) is KmADHACE in component vADH (mM).
C CONSTS(94) is KiADHACE in component vADH (mM).
C CONSTS(95) is KiADHETOH in component vADH (mM).
C CONSTS(96) is KATPASE in component vATP (per_min).
C CONSTS(97) is KeqAK in component vATP (dimensionless).
C CONSTS(98) is KmGA3P in component vTPI (mM).
C CONSTS(99) is KmDHAP in component vTPI (mM).
C CONSTS(100) is VmDHAP in component vTPI (mM_per_min).
C CONSTS(101) is VmGA3P in component vTPI (mM_per_min).
C CONSTS(102) is VmG6PDH in component vG6PDH (mM_per_min).
C CONSTS(103) is KmG6P in component vG6PDH (mM).
C CONSTS(104) is KmNADP in component vG6PDH (mM).
C CONSTS(105) is KiNADPH in component vG6PDH (mM).
C CONSTS(106) is Vm6PGL in component v6PGL (mM_per_min).
C CONSTS(107) is Km6PGL in component v6PGL (mM).
C CONSTS(108) is VmGluDH in component vGluDH (mM_per_min).
C CONSTS(109) is KmGluconate in component vGluDH (mM).
C CONSTS(110) is KmNADP in component vGluDH (mM).
C CONSTS(111) is KiNADPH in component vGluDH (mM).
C CONSTS(112) is VmPPIf in component vPPI (mM_per_min).
C CONSTS(113) is VmPPIr in component vPPI (mM_per_min).
C CONSTS(114) is KmRibu5P in component vPPI (mM).
C CONSTS(115) is KmRibo5P in component vPPI (mM).
C CONSTS(116) is VmTransk1f in component vTransk1 (mM_per_min).
C CONSTS(117) is VmTransk1r in component vTransk1 (mM_per_min).
C CONSTS(118) is KmRibose5P in component vTransk1 (mM).
C CONSTS(119) is KmXyl5P in component vTransk1 (mM).
C CONSTS(120) is KmGA3P in component vTransk1 (mM).
C CONSTS(121) is KmSeduhept in component vTransk1 (mM).
C CONSTS(122) is VmR5PIf in component vR5PI (mM_per_min).
C CONSTS(123) is VmR5PIr in component vR5PI (mM_per_min).
C CONSTS(124) is KmRibu5P in component vR5PI (mM).
C CONSTS(125) is KmXyl in component vR5PI (mM).
C CONSTS(126) is VmTransaldf in component vTransald (mM_per_min).
C CONSTS(127) is VmTransaldr in component vTransald (mM_per_min).
C CONSTS(128) is KmGA3P in component vTransald (mM).
C CONSTS(129) is KmSeduhept in component vTransald (mM).
C CONSTS(130) is KmF6P in component vTransald (mM).
C CONSTS(131) is KmEry4P in component vTransald (mM).
C CONSTS(132) is VmTransk2f in component vTransk2 (mM_per_min).
C CONSTS(133) is VmTransk2r in component vTransk2 (mM_per_min).
C CONSTS(134) is KmXyl5P in component vTransk2 (mM).
C CONSTS(135) is KmEry4P in component vTransk2 (mM).
C CONSTS(136) is KmF6P in component vTransk2 (mM).
C CONSTS(137) is KmGA3P in component vTransk2 (mM).
C CONSTS(138) is kNADPH in component vNADPH (per_min).
C CONSTS(139) is VmGLT in component vGLT (mM_per_min).
C CONSTS(140) is KeqGLT in component vGLT (dimensionless).
C CONSTS(141) is KmGLTGLCo in component vGLT (mM).
C CONSTS(142) is KmGLTGLCi in component vGLT (mM).
C RATES(1) is d/dt P in component cytoplasm (mM).
C RATES(2) is d/dt G6P in component cytoplasm (mM).
C RATES(3) is d/dt F6P in component cytoplasm (mM).
C RATES(4) is d/dt F16P in component cytoplasm (mM).
C RATES(5) is d/dt NADH in component cytoplasm (mM).
C RATES(6) is d/dt NAD in component cytoplasm (mM).
C RATES(7) is d/dt BPG in component cytoplasm (mM).
C RATES(8) is d/dt P3G in component cytoplasm (mM).
C RATES(9) is d/dt P2G in component cytoplasm (mM).
C RATES(10) is d/dt PEP in component cytoplasm (mM).
C RATES(11) is d/dt PYR in component cytoplasm (mM).
C RATES(12) is d/dt ACE in component cytoplasm (mM).
C RATES(13) is d/dt GA3P in component cytoplasm (mM).
C RATES(14) is d/dt DHAP in component cytoplasm (mM).
C RATES(15) is d/dt D6PGluconoLactone in component cytoplasm (mM).
C RATES(16) is d/dt D6PGluconate in component cytoplasm (mM).
C RATES(17) is d/dt NADP in component cytoplasm (mM).
C RATES(18) is d/dt NADPH in component cytoplasm (mM).
C RATES(19) is d/dt Ribulose5P in component cytoplasm (mM).
C RATES(20) is d/dt Ribose5P in component cytoplasm (mM).
C RATES(21) is d/dt Xyl5P in component cytoplasm (mM).
C RATES(22) is d/dt Seduhept7P in component cytoplasm (mM).
C RATES(23) is d/dt Erythrose4P in component cytoplasm (mM).
C RATES(24) is d/dt GLCi in component cytoplasm (mM).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1
      CONSTS(2) = 1
      CONSTS(3) = 4.1
      CONSTS(4) = 1
      CONSTS(5) = 50
      CONSTS(6) = 0.1
      CONSTS(7) = 0.15
      CONSTS(8) = 50
      CONSTS(9) = 1
      STATES(1) = 5
      STATES(2) = 1.39
      STATES(3) = 0.28
      STATES(4) = 0.1
      STATES(5) = 0.39
      STATES(6) = 1.2
      STATES(7) = 0.1
      STATES(8) = 0.1
      STATES(9) = 0.1
      STATES(10) = 0.1
      STATES(11) = 3.36
      STATES(12) = 0.04
      CONSTS(10) = 0.1
      STATES(13) = 0.05
      STATES(14) = 1
      STATES(15) = 0.1
      STATES(16) = 0.1
      STATES(17) = 0.4
      STATES(18) = 1.6
      STATES(19) = 0.1
      STATES(20) = 0.1
      STATES(21) = 0.1
      STATES(22) = 0.1
      STATES(23) = 0
      STATES(24) = 0.087
      CONSTS(11) = 0.02
      CONSTS(12) = 226.452
      CONSTS(13) = 0.45
      CONSTS(14) = 3800
      CONSTS(15) = 0.15
      CONSTS(16) = 0.08
      CONSTS(17) = 30
      CONSTS(18) = 0.23
      CONSTS(19) = 339.677
      CONSTS(20) = 1.4
      CONSTS(21) = 0.314
      CONSTS(22) = 0.3
      CONSTS(23) = 5.12
      CONSTS(24) = 182.903
      CONSTS(25) = 0.45
      CONSTS(26) = 0.1
      CONSTS(27) = 0.71
      CONSTS(28) = 0.66
      CONSTS(29) = 0.0174
      CONSTS(30) = 0.000682
      CONSTS(31) = 0.397
      CONSTS(32) = 0.111
      CONSTS(33) = 0.0845
      CONSTS(34) = 0.0995
      CONSTS(35) = 100
      CONSTS(36) = 0.65
      CONSTS(37) = 3
      CONSTS(38) = 322.258
      CONSTS(39) = 0.045
      CONSTS(40) = 0.069
      CONSTS(41) = 0.3
      CONSTS(42) = 2.4
      CONSTS(43) = 2
      CONSTS(44) = 10
      CONSTS(45) = 70.15
      CONSTS(46) = 4300
      CONSTS(47) = 0.045
      CONSTS(48) = 0.4
      CONSTS(49) = 0.023
      CONSTS(50) = 0.93
      CONSTS(51) = 1
      CONSTS(52) = 6549.68
      CONSTS(53) = 0.0098
      CONSTS(54) = 0.06
      CONSTS(55) = 0.045
      CONSTS(56) = 1184.52
      CONSTS(57) = 0.21
      CONSTS(58) = 0.09
      CONSTS(59) = 0.005
      CONSTS(60) = 1306.45
      CONSTS(61) = 3200
      CONSTS(62) = 0.45
      CONSTS(63) = 0.3
      CONSTS(64) = 0.53
      CONSTS(65) = 0.2
      CONSTS(66) = 0.003
      CONSTS(67) = 2525.81
      CONSTS(68) = 1.2
      CONSTS(69) = 0.19
      CONSTS(70) = 0.08
      CONSTS(71) = 365.806
      CONSTS(72) = 0.04
      CONSTS(73) = 6.7
      CONSTS(74) = 0.5
      CONSTS(75) = 1088.71
      CONSTS(76) = 0.14
      CONSTS(77) = 0.53
      CONSTS(78) = 0.45
      CONSTS(79) = 6500
      CONSTS(80) = 21
      CONSTS(81) = 1.5
      CONSTS(82) = 174.194
      CONSTS(83) = 1.9
      CONSTS(84) = 4.33
      CONSTS(85) = 21.4
      CONSTS(86) = 810
      CONSTS(87) = 0.92
      CONSTS(88) = 17
      CONSTS(89) = 6.9e-5
      CONSTS(90) = 0.17
      CONSTS(91) = 0.11
      CONSTS(92) = 0.031
      CONSTS(93) = 1.11
      CONSTS(94) = 1.1
      CONSTS(95) = 90
      CONSTS(96) = 39.5
      CONSTS(97) = 0.45
      CONSTS(98) = 1.27
      CONSTS(99) = 1.23
      CONSTS(100) = 10900
      CONSTS(101) = 555
      CONSTS(102) = 4
      CONSTS(103) = 0.04
      CONSTS(104) = 0.02
      CONSTS(105) = 0.017
      CONSTS(106) = 4
      CONSTS(107) = 0.8
      CONSTS(108) = 4
      CONSTS(109) = 0.02
      CONSTS(110) = 0.03
      CONSTS(111) = 0.03
      CONSTS(112) = 3458
      CONSTS(113) = 3458
      CONSTS(114) = 1.6
      CONSTS(115) = 1.6
      CONSTS(116) = 4
      CONSTS(117) = 2
      CONSTS(118) = 0.1
      CONSTS(119) = 0.15
      CONSTS(120) = 0.1
      CONSTS(121) = 0.15
      CONSTS(122) = 1039
      CONSTS(123) = 1039
      CONSTS(124) = 1.5
      CONSTS(125) = 1.5
      CONSTS(126) = 55
      CONSTS(127) = 10
      CONSTS(128) = 0.22
      CONSTS(129) = 0.18
      CONSTS(130) = 0.32
      CONSTS(131) = 0.018
      CONSTS(132) = 3.2
      CONSTS(133) = 43
      CONSTS(134) = 0.16
      CONSTS(135) = 0.09
      CONSTS(136) = 1.1
      CONSTS(137) = 2.1
      CONSTS(138) = 2
      CONSTS(139) = 97.264
      CONSTS(140) = 1
      CONSTS(141) = 1.1918
      CONSTS(142) = 1.1918
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = ( CONSTS(9)*CONSTS(12)*(( - STATES(2)*(CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(13)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(13)*STATES(1) ** 2.00000) ** (1.0 / 2)))/( (1.00000 -  4.00000*CONSTS(13))*CONSTS(14))+( STATES(24)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(13)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(13)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(13)*STATES(1) ** 2.00000) ** (1.0 / 2)))/(2.00000 -  8.00000*CONSTS(13))))/( CONSTS(15)*CONSTS(16)*(1.00000+STATES(2)/CONSTS(17)+STATES(24)/CONSTS(16))*(1.00000+(CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(13)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(13)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (1.00000 -  4.00000*CONSTS(13))*CONSTS(18))+(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(13)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(13)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(13)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (2.00000 -  8.00000*CONSTS(13))*CONSTS(15))))
      ALGBRC(4) = ( CONSTS(9)*CONSTS(139)*(CONSTS(8) - STATES(24)/CONSTS(140)))/( CONSTS(141)*(1.00000+CONSTS(8)/CONSTS(141)+STATES(24)/CONSTS(142)+( 0.910000*CONSTS(8)*STATES(24))/( CONSTS(142)*CONSTS(141))))
      RATES(24) = (ALGBRC(4) - ALGBRC(1))/CONSTS(9)
      ALGBRC(5) =  CONSTS(9)*CONSTS(23)*CONSTS(24)*STATES(3)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))
      ALGBRC(6) = 1.00000+STATES(3)/CONSTS(26)+(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(27))+( CONSTS(23)*STATES(3)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2)))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(27)*CONSTS(26))
      ALGBRC(7) =  1.00000+( CONSTS(29)*CONSTS(11))/CONSTS(30)+( CONSTS(31)*STATES(4))/CONSTS(32) ** 2.00000*1.00000+( 2.00000*CONSTS(33)*CONSTS(25)*CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2) ** 2.00000)/( (- 1.00000+ 4.00000*CONSTS(25))*CONSTS(34)*(((CONSTS(3) - STATES(1))+ 4.00000*CONSTS(25)*STATES(1)) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))) ** 2.00000
      ALGBRC(8) =  1.00000+( CONSTS(35)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2)))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(36)) ** 2.00000*1.00000+( CONSTS(37)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2)))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(27)) ** 2.00000
      ALGBRC(9) =  1.00000+CONSTS(11)/CONSTS(30)+STATES(4)/CONSTS(32) ** 2.00000*1.00000+( 2.00000*CONSTS(25)*CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2) ** 2.00000)/( (- 1.00000+ 4.00000*CONSTS(25))*CONSTS(34)*(((CONSTS(3) - STATES(1))+ 4.00000*CONSTS(25)*STATES(1)) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))) ** 2.00000*1.00000+(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(36)) ** 2.00000
      ALGBRC(10) =  (2.00000 -  8.00000*CONSTS(25))*CONSTS(27)*CONSTS(26)*(( CONSTS(28)*ALGBRC(7)*ALGBRC(8))/ALGBRC(9)+1.00000+STATES(3)/CONSTS(26)+(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(27))+( CONSTS(23)*STATES(3)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2)))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(27)*CONSTS(26)) ** 2.00000)
      ALGBRC(11) = ( ALGBRC(5)*ALGBRC(6))/ALGBRC(10)
      ALGBRC(12) = ( (( CONSTS(9)*CONSTS(38)*STATES(4))/CONSTS(41))*(1.00000 - ( STATES(14)*STATES(13))/( STATES(4)*CONSTS(40))))/(1.00000+STATES(4)/CONSTS(41)+STATES(14)/CONSTS(42)+STATES(13)/CONSTS(43)+( STATES(4)*STATES(13))/( CONSTS(41)*CONSTS(44))+( STATES(14)*STATES(13))/( CONSTS(42)*CONSTS(43)))
      RATES(4) = (ALGBRC(11) - ALGBRC(12))/CONSTS(9)
      ALGBRC(14) = ( (( CONSTS(9)*CONSTS(2)*CONSTS(56)*STATES(13)*STATES(6))/( CONSTS(57)*CONSTS(58)))*(1.00000 - ( STATES(7)*STATES(5))/( STATES(13)*STATES(6)*CONSTS(59))))/( (1.00000+STATES(13)/CONSTS(57)+STATES(7)/CONSTS(53))*(1.00000+STATES(6)/CONSTS(58)+STATES(5)/CONSTS(54)))
      ALGBRC(15) = ( 1.00000*CONSTS(9)*CONSTS(60)*(( CONSTS(61)*STATES(7)*(CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(62)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(62)*STATES(1) ** 2.00000) ** (1.0 / 2)))/(1.00000 -  4.00000*CONSTS(62)) - ( (((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(62)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(62)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(62)*STATES(1) ** 2.00000) ** (1.0 / 2))*STATES(8))/(2.00000 -  8.00000*CONSTS(62))))/( CONSTS(63)*CONSTS(64)*(1.00000+(CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(62)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(62)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (1.00000 -  4.00000*CONSTS(62))*CONSTS(65))+(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(62)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(62)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(62)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (2.00000 -  8.00000*CONSTS(62))*CONSTS(63)))*(1.00000+STATES(7)/CONSTS(66)+STATES(8)/CONSTS(64)))
      RATES(7) = (ALGBRC(14) - ALGBRC(15))/CONSTS(9)
      ALGBRC(13) = ( CONSTS(9)*CONSTS(45)*(( - CONSTS(7)*STATES(6))/CONSTS(46)+( STATES(5)*STATES(14))/(1.00000+CONSTS(47))))/( CONSTS(48)*CONSTS(49)*(1.00000+STATES(6)/CONSTS(50)+STATES(5)/CONSTS(49))*(1.00000+CONSTS(7)/CONSTS(51)+STATES(14)/( (1.00000+CONSTS(47))*CONSTS(48))))
      ALGBRC(16) = ( CONSTS(9)*CONSTS(1)*(( CONSTS(100)*STATES(13))/CONSTS(98) - ( CONSTS(101)*STATES(14))/CONSTS(99)))/(1.00000+STATES(13)/CONSTS(98)+STATES(14)/CONSTS(99))
      RATES(14) = ((ALGBRC(12) - ALGBRC(13))+ALGBRC(16))/CONSTS(9)
      ALGBRC(3) = ( (( CONSTS(9)*CONSTS(19))/CONSTS(20))*(STATES(2) - STATES(3)/CONSTS(21)))/(1.00000+STATES(2)/CONSTS(20)+STATES(3)/CONSTS(22))
      ALGBRC(18) = (( CONSTS(9)*CONSTS(102)*STATES(2)*STATES(17))/( CONSTS(103)*CONSTS(104)))/( (1.00000+STATES(2)/CONSTS(103)+STATES(18)/CONSTS(105))*(1.00000+STATES(17)/CONSTS(104)))
      RATES(2) = ((ALGBRC(1) - ALGBRC(3)) - ALGBRC(18))/CONSTS(9)
      ALGBRC(17) = ( (( CONSTS(9)*CONSTS(67))/CONSTS(68))*(STATES(8) - STATES(9)/CONSTS(69)))/(1.00000+STATES(8)/CONSTS(68)+STATES(9)/CONSTS(70))
      RATES(8) = (ALGBRC(15) - ALGBRC(17))/CONSTS(9)
      ALGBRC(19) = ( (( CONSTS(9)*CONSTS(71))/CONSTS(72))*(STATES(9) - STATES(10)/CONSTS(73)))/(1.00000+STATES(9)/CONSTS(72)+STATES(10)/CONSTS(74))
      RATES(9) = (ALGBRC(17) - ALGBRC(19))/CONSTS(9)
      ALGBRC(20) = ( CONSTS(9)*CONSTS(106)*STATES(15))/(CONSTS(107)+STATES(15))
      RATES(15) = (ALGBRC(18) - ALGBRC(20))/CONSTS(9)
      ALGBRC(21) = ( (( CONSTS(9)*CONSTS(75))/( CONSTS(76)*CONSTS(77)))*(( STATES(10)*(CONSTS(3) -  (((STATES(1) ** 2.00000 -  4.00000*CONSTS(78)*STATES(1) ** 2.00000) -  2.00000*STATES(1)*CONSTS(3))+ 8.00000*CONSTS(78)*STATES(1)*CONSTS(3)+CONSTS(3) ** 2.00000) ** (1.0 / 2)))/(1.00000 -  4.00000*CONSTS(78)) - (( STATES(11)*(((STATES(1) -  4.00000*CONSTS(78)*STATES(1)) - CONSTS(3))+ (((STATES(1) ** 2.00000 -  4.00000*CONSTS(78)*STATES(1) ** 2.00000) -  2.00000*STATES(1)*CONSTS(3))+ 8.00000*CONSTS(78)*STATES(1)*CONSTS(3)+CONSTS(3) ** 2.00000) ** (1.0 / 2)))/(2.00000 -  8.00000*CONSTS(78)))/CONSTS(79)))/( (1.00000+STATES(10)/CONSTS(76)+STATES(11)/CONSTS(80))*(1.00000+((((STATES(1) -  4.00000*CONSTS(78)*STATES(1)) - CONSTS(3))+ (((STATES(1) ** 2.00000 -  4.00000*CONSTS(78)*STATES(1) ** 2.00000) -  2.00000*STATES(1)*CONSTS(3))+ 8.00000*CONSTS(78)*STATES(1)*CONSTS(3)+CONSTS(3) ** 2.00000) ** (1.0 / 2))/(2.00000 -  8.00000*CONSTS(78)))/CONSTS(81)+((CONSTS(3) -  (((STATES(1) ** 2.00000 -  4.00000*CONSTS(78)*STATES(1) ** 2.00000) -  2.00000*STATES(1)*CONSTS(3))+ 8.00000*CONSTS(78)*STATES(1)*CONSTS(3)+CONSTS(3) ** 2.00000) ** (1.0 / 2))/(1.00000 -  4.00000*CONSTS(78)))/CONSTS(77)))
      RATES(10) = (ALGBRC(19) - ALGBRC(21))/CONSTS(9)
      ALGBRC(22) = (( CONSTS(9)*CONSTS(108)*STATES(16)*STATES(17))/( CONSTS(109)*CONSTS(110)))/( (1.00000+STATES(16)/CONSTS(109)+STATES(18)/CONSTS(111))*(1.00000+STATES(17)/CONSTS(110)))
      RATES(16) = (ALGBRC(20) - ALGBRC(22))/CONSTS(9)
      ALGBRC(24) = ( CONSTS(9)*CONSTS(96)*(((STATES(1) -  4.00000*CONSTS(97)*STATES(1)) - CONSTS(3))+ (((STATES(1) ** 2.00000 -  4.00000*CONSTS(97)*STATES(1) ** 2.00000) -  2.00000*STATES(1)*CONSTS(3))+ 8.00000*CONSTS(97)*STATES(1)*CONSTS(3)+CONSTS(3) ** 2.00000) ** (1.0 / 2)))/(2.00000 -  8.00000*CONSTS(97))
      RATES(1) = (((ALGBRC(15) - (ALGBRC(1)+ALGBRC(11)))+ALGBRC(21)) - ALGBRC(24))/CONSTS(9)
      ALGBRC(23) = (( CONSTS(9)*CONSTS(82)*STATES(11) ** CONSTS(83))/CONSTS(84) ** CONSTS(83))/(1.00000+STATES(11) ** CONSTS(83)/CONSTS(84) ** CONSTS(83))
      RATES(11) = (ALGBRC(21) - ALGBRC(23))/CONSTS(9)
      ALGBRC(26) =  CONSTS(9)*CONSTS(138)*STATES(18)
      RATES(17) = (ALGBRC(26) - (ALGBRC(18)+ALGBRC(22)))/CONSTS(9)
      RATES(18) = ((ALGBRC(18)+ALGBRC(22)) - ALGBRC(26))/CONSTS(9)
      ALGBRC(25) = ( CONSTS(9)*(( CONSTS(112)*STATES(19))/CONSTS(114) - ( CONSTS(113)*STATES(20))/CONSTS(115)))/(1.00000+STATES(19)/CONSTS(114)+STATES(20)/CONSTS(115))
      ALGBRC(28) = ( CONSTS(9)*(( CONSTS(116)*STATES(20)*STATES(21))/( CONSTS(118)*CONSTS(119)) - ( CONSTS(117)*STATES(13)*STATES(22))/( CONSTS(120)*CONSTS(121))))/( (1.00000+STATES(20)/CONSTS(118)+STATES(13)/CONSTS(120))*(1.00000+STATES(21)/CONSTS(119)+STATES(22)/CONSTS(121)))
      RATES(20) = (ALGBRC(25) - ALGBRC(28))/CONSTS(9)
      ALGBRC(27) =  CONSTS(9)*CONSTS(85)*STATES(12)
      ALGBRC(29) = ( (( CONSTS(9)*- CONSTS(86))/( CONSTS(87)*CONSTS(88)))*( STATES(6)*CONSTS(5) - ( STATES(5)*STATES(12))/CONSTS(89)))/(1.00000+STATES(6)/CONSTS(87)+( CONSTS(90)*CONSTS(5))/( CONSTS(87)*CONSTS(88))+( CONSTS(91)*STATES(12))/( CONSTS(92)*CONSTS(93))+STATES(5)/CONSTS(92)+( STATES(6)*CONSTS(5))/( CONSTS(87)*CONSTS(88))+( CONSTS(91)*STATES(6)*STATES(12))/( CONSTS(87)*CONSTS(92)*CONSTS(93))+( CONSTS(90)*CONSTS(5)*STATES(5))/( CONSTS(87)*CONSTS(88)*CONSTS(92))+( STATES(5)*STATES(12))/( CONSTS(92)*CONSTS(93))+( STATES(6)*CONSTS(5)*STATES(12))/( CONSTS(87)*CONSTS(88)*CONSTS(94))+( CONSTS(5)*STATES(5)*STATES(12))/( CONSTS(95)*CONSTS(92)*CONSTS(93)))
      RATES(5) = (((ALGBRC(14) - ALGBRC(13))+ 3.00000*ALGBRC(27)) - ALGBRC(29))/CONSTS(9)
      RATES(6) = (((ALGBRC(13) - ALGBRC(14)) -  3.00000*ALGBRC(27))+ALGBRC(29))/CONSTS(9)
      RATES(12) = ((ALGBRC(23) -  2.00000*ALGBRC(27)) - ALGBRC(29))/CONSTS(9)
      ALGBRC(30) = ( CONSTS(9)*(( CONSTS(122)*STATES(19))/CONSTS(124) - ( CONSTS(123)*STATES(21))/CONSTS(125)))/(1.00000+STATES(19)/CONSTS(124)+STATES(21)/CONSTS(125))
      RATES(19) = ((ALGBRC(22) - ALGBRC(25)) - ALGBRC(30))/CONSTS(9)
      ALGBRC(31) = ( CONSTS(9)*(( CONSTS(126)*STATES(13)*STATES(22))/( CONSTS(128)*CONSTS(129)) - ( CONSTS(127)*STATES(3)*STATES(23))/( CONSTS(130)*CONSTS(131))))/( (1.00000+STATES(13)/CONSTS(128)+STATES(3)/CONSTS(130))*(1.00000+STATES(22)/CONSTS(129)+STATES(23)/CONSTS(131)))
      RATES(22) = (ALGBRC(28) - ALGBRC(31))/CONSTS(9)
      ALGBRC(32) = ( CONSTS(9)*(( CONSTS(132)*STATES(23)*STATES(21))/( CONSTS(135)*CONSTS(134)) - ( CONSTS(133)*STATES(3)*STATES(13))/( CONSTS(136)*CONSTS(137))))/( (1.00000+STATES(21)/CONSTS(134)+STATES(13)/CONSTS(137))*(1.00000+STATES(23)/CONSTS(135)+STATES(3)/CONSTS(136)))
      RATES(3) = ((ALGBRC(3) - ALGBRC(11))+ALGBRC(31)+ALGBRC(32))/CONSTS(9)
      RATES(13) = (((((ALGBRC(12) - ALGBRC(14)) - ALGBRC(16))+ALGBRC(28)) - ALGBRC(31))+ALGBRC(32))/CONSTS(9)
      RATES(21) = ((ALGBRC(30) - ALGBRC(28)) - ALGBRC(32))/CONSTS(9)
      RATES(23) = (ALGBRC(31) - ALGBRC(32))/CONSTS(9)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = ( CONSTS(9)*CONSTS(12)*(( - STATES(2)*(CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(13)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(13)*STATES(1) ** 2.00000) ** (1.0 / 2)))/( (1.00000 -  4.00000*CONSTS(13))*CONSTS(14))+( STATES(24)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(13)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(13)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(13)*STATES(1) ** 2.00000) ** (1.0 / 2)))/(2.00000 -  8.00000*CONSTS(13))))/( CONSTS(15)*CONSTS(16)*(1.00000+STATES(2)/CONSTS(17)+STATES(24)/CONSTS(16))*(1.00000+(CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(13)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(13)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (1.00000 -  4.00000*CONSTS(13))*CONSTS(18))+(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(13)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(13)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(13)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (2.00000 -  8.00000*CONSTS(13))*CONSTS(15))))
      ALGBRC(4) = ( CONSTS(9)*CONSTS(139)*(CONSTS(8) - STATES(24)/CONSTS(140)))/( CONSTS(141)*(1.00000+CONSTS(8)/CONSTS(141)+STATES(24)/CONSTS(142)+( 0.910000*CONSTS(8)*STATES(24))/( CONSTS(142)*CONSTS(141))))
      ALGBRC(5) =  CONSTS(9)*CONSTS(23)*CONSTS(24)*STATES(3)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))
      ALGBRC(6) = 1.00000+STATES(3)/CONSTS(26)+(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(27))+( CONSTS(23)*STATES(3)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2)))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(27)*CONSTS(26))
      ALGBRC(7) =  1.00000+( CONSTS(29)*CONSTS(11))/CONSTS(30)+( CONSTS(31)*STATES(4))/CONSTS(32) ** 2.00000*1.00000+( 2.00000*CONSTS(33)*CONSTS(25)*CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2) ** 2.00000)/( (- 1.00000+ 4.00000*CONSTS(25))*CONSTS(34)*(((CONSTS(3) - STATES(1))+ 4.00000*CONSTS(25)*STATES(1)) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))) ** 2.00000
      ALGBRC(8) =  1.00000+( CONSTS(35)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2)))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(36)) ** 2.00000*1.00000+( CONSTS(37)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2)))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(27)) ** 2.00000
      ALGBRC(9) =  1.00000+CONSTS(11)/CONSTS(30)+STATES(4)/CONSTS(32) ** 2.00000*1.00000+( 2.00000*CONSTS(25)*CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2) ** 2.00000)/( (- 1.00000+ 4.00000*CONSTS(25))*CONSTS(34)*(((CONSTS(3) - STATES(1))+ 4.00000*CONSTS(25)*STATES(1)) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))) ** 2.00000*1.00000+(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(36)) ** 2.00000
      ALGBRC(10) =  (2.00000 -  8.00000*CONSTS(25))*CONSTS(27)*CONSTS(26)*(( CONSTS(28)*ALGBRC(7)*ALGBRC(8))/ALGBRC(9)+1.00000+STATES(3)/CONSTS(26)+(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(27))+( CONSTS(23)*STATES(3)*(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(25)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(25)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(25)*STATES(1) ** 2.00000) ** (1.0 / 2)))/( (2.00000 -  8.00000*CONSTS(25))*CONSTS(27)*CONSTS(26)) ** 2.00000)
      ALGBRC(11) = ( ALGBRC(5)*ALGBRC(6))/ALGBRC(10)
      ALGBRC(12) = ( (( CONSTS(9)*CONSTS(38)*STATES(4))/CONSTS(41))*(1.00000 - ( STATES(14)*STATES(13))/( STATES(4)*CONSTS(40))))/(1.00000+STATES(4)/CONSTS(41)+STATES(14)/CONSTS(42)+STATES(13)/CONSTS(43)+( STATES(4)*STATES(13))/( CONSTS(41)*CONSTS(44))+( STATES(14)*STATES(13))/( CONSTS(42)*CONSTS(43)))
      ALGBRC(14) = ( (( CONSTS(9)*CONSTS(2)*CONSTS(56)*STATES(13)*STATES(6))/( CONSTS(57)*CONSTS(58)))*(1.00000 - ( STATES(7)*STATES(5))/( STATES(13)*STATES(6)*CONSTS(59))))/( (1.00000+STATES(13)/CONSTS(57)+STATES(7)/CONSTS(53))*(1.00000+STATES(6)/CONSTS(58)+STATES(5)/CONSTS(54)))
      ALGBRC(15) = ( 1.00000*CONSTS(9)*CONSTS(60)*(( CONSTS(61)*STATES(7)*(CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(62)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(62)*STATES(1) ** 2.00000) ** (1.0 / 2)))/(1.00000 -  4.00000*CONSTS(62)) - ( (((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(62)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(62)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(62)*STATES(1) ** 2.00000) ** (1.0 / 2))*STATES(8))/(2.00000 -  8.00000*CONSTS(62))))/( CONSTS(63)*CONSTS(64)*(1.00000+(CONSTS(3) -  (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(62)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(62)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (1.00000 -  4.00000*CONSTS(62))*CONSTS(65))+(((- CONSTS(3)+STATES(1)) -  4.00000*CONSTS(62)*STATES(1))+ (((CONSTS(3) ** 2.00000 -  2.00000*CONSTS(3)*STATES(1))+ 8.00000*CONSTS(62)*CONSTS(3)*STATES(1)+STATES(1) ** 2.00000) -  4.00000*CONSTS(62)*STATES(1) ** 2.00000) ** (1.0 / 2))/( (2.00000 -  8.00000*CONSTS(62))*CONSTS(63)))*(1.00000+STATES(7)/CONSTS(66)+STATES(8)/CONSTS(64)))
      ALGBRC(13) = ( CONSTS(9)*CONSTS(45)*(( - CONSTS(7)*STATES(6))/CONSTS(46)+( STATES(5)*STATES(14))/(1.00000+CONSTS(47))))/( CONSTS(48)*CONSTS(49)*(1.00000+STATES(6)/CONSTS(50)+STATES(5)/CONSTS(49))*(1.00000+CONSTS(7)/CONSTS(51)+STATES(14)/( (1.00000+CONSTS(47))*CONSTS(48))))
      ALGBRC(16) = ( CONSTS(9)*CONSTS(1)*(( CONSTS(100)*STATES(13))/CONSTS(98) - ( CONSTS(101)*STATES(14))/CONSTS(99)))/(1.00000+STATES(13)/CONSTS(98)+STATES(14)/CONSTS(99))
      ALGBRC(3) = ( (( CONSTS(9)*CONSTS(19))/CONSTS(20))*(STATES(2) - STATES(3)/CONSTS(21)))/(1.00000+STATES(2)/CONSTS(20)+STATES(3)/CONSTS(22))
      ALGBRC(18) = (( CONSTS(9)*CONSTS(102)*STATES(2)*STATES(17))/( CONSTS(103)*CONSTS(104)))/( (1.00000+STATES(2)/CONSTS(103)+STATES(18)/CONSTS(105))*(1.00000+STATES(17)/CONSTS(104)))
      ALGBRC(17) = ( (( CONSTS(9)*CONSTS(67))/CONSTS(68))*(STATES(8) - STATES(9)/CONSTS(69)))/(1.00000+STATES(8)/CONSTS(68)+STATES(9)/CONSTS(70))
      ALGBRC(19) = ( (( CONSTS(9)*CONSTS(71))/CONSTS(72))*(STATES(9) - STATES(10)/CONSTS(73)))/(1.00000+STATES(9)/CONSTS(72)+STATES(10)/CONSTS(74))
      ALGBRC(20) = ( CONSTS(9)*CONSTS(106)*STATES(15))/(CONSTS(107)+STATES(15))
      ALGBRC(21) = ( (( CONSTS(9)*CONSTS(75))/( CONSTS(76)*CONSTS(77)))*(( STATES(10)*(CONSTS(3) -  (((STATES(1) ** 2.00000 -  4.00000*CONSTS(78)*STATES(1) ** 2.00000) -  2.00000*STATES(1)*CONSTS(3))+ 8.00000*CONSTS(78)*STATES(1)*CONSTS(3)+CONSTS(3) ** 2.00000) ** (1.0 / 2)))/(1.00000 -  4.00000*CONSTS(78)) - (( STATES(11)*(((STATES(1) -  4.00000*CONSTS(78)*STATES(1)) - CONSTS(3))+ (((STATES(1) ** 2.00000 -  4.00000*CONSTS(78)*STATES(1) ** 2.00000) -  2.00000*STATES(1)*CONSTS(3))+ 8.00000*CONSTS(78)*STATES(1)*CONSTS(3)+CONSTS(3) ** 2.00000) ** (1.0 / 2)))/(2.00000 -  8.00000*CONSTS(78)))/CONSTS(79)))/( (1.00000+STATES(10)/CONSTS(76)+STATES(11)/CONSTS(80))*(1.00000+((((STATES(1) -  4.00000*CONSTS(78)*STATES(1)) - CONSTS(3))+ (((STATES(1) ** 2.00000 -  4.00000*CONSTS(78)*STATES(1) ** 2.00000) -  2.00000*STATES(1)*CONSTS(3))+ 8.00000*CONSTS(78)*STATES(1)*CONSTS(3)+CONSTS(3) ** 2.00000) ** (1.0 / 2))/(2.00000 -  8.00000*CONSTS(78)))/CONSTS(81)+((CONSTS(3) -  (((STATES(1) ** 2.00000 -  4.00000*CONSTS(78)*STATES(1) ** 2.00000) -  2.00000*STATES(1)*CONSTS(3))+ 8.00000*CONSTS(78)*STATES(1)*CONSTS(3)+CONSTS(3) ** 2.00000) ** (1.0 / 2))/(1.00000 -  4.00000*CONSTS(78)))/CONSTS(77)))
      ALGBRC(22) = (( CONSTS(9)*CONSTS(108)*STATES(16)*STATES(17))/( CONSTS(109)*CONSTS(110)))/( (1.00000+STATES(16)/CONSTS(109)+STATES(18)/CONSTS(111))*(1.00000+STATES(17)/CONSTS(110)))
      ALGBRC(24) = ( CONSTS(9)*CONSTS(96)*(((STATES(1) -  4.00000*CONSTS(97)*STATES(1)) - CONSTS(3))+ (((STATES(1) ** 2.00000 -  4.00000*CONSTS(97)*STATES(1) ** 2.00000) -  2.00000*STATES(1)*CONSTS(3))+ 8.00000*CONSTS(97)*STATES(1)*CONSTS(3)+CONSTS(3) ** 2.00000) ** (1.0 / 2)))/(2.00000 -  8.00000*CONSTS(97))
      ALGBRC(23) = (( CONSTS(9)*CONSTS(82)*STATES(11) ** CONSTS(83))/CONSTS(84) ** CONSTS(83))/(1.00000+STATES(11) ** CONSTS(83)/CONSTS(84) ** CONSTS(83))
      ALGBRC(26) =  CONSTS(9)*CONSTS(138)*STATES(18)
      ALGBRC(25) = ( CONSTS(9)*(( CONSTS(112)*STATES(19))/CONSTS(114) - ( CONSTS(113)*STATES(20))/CONSTS(115)))/(1.00000+STATES(19)/CONSTS(114)+STATES(20)/CONSTS(115))
      ALGBRC(28) = ( CONSTS(9)*(( CONSTS(116)*STATES(20)*STATES(21))/( CONSTS(118)*CONSTS(119)) - ( CONSTS(117)*STATES(13)*STATES(22))/( CONSTS(120)*CONSTS(121))))/( (1.00000+STATES(20)/CONSTS(118)+STATES(13)/CONSTS(120))*(1.00000+STATES(21)/CONSTS(119)+STATES(22)/CONSTS(121)))
      ALGBRC(27) =  CONSTS(9)*CONSTS(85)*STATES(12)
      ALGBRC(29) = ( (( CONSTS(9)*- CONSTS(86))/( CONSTS(87)*CONSTS(88)))*( STATES(6)*CONSTS(5) - ( STATES(5)*STATES(12))/CONSTS(89)))/(1.00000+STATES(6)/CONSTS(87)+( CONSTS(90)*CONSTS(5))/( CONSTS(87)*CONSTS(88))+( CONSTS(91)*STATES(12))/( CONSTS(92)*CONSTS(93))+STATES(5)/CONSTS(92)+( STATES(6)*CONSTS(5))/( CONSTS(87)*CONSTS(88))+( CONSTS(91)*STATES(6)*STATES(12))/( CONSTS(87)*CONSTS(92)*CONSTS(93))+( CONSTS(90)*CONSTS(5)*STATES(5))/( CONSTS(87)*CONSTS(88)*CONSTS(92))+( STATES(5)*STATES(12))/( CONSTS(92)*CONSTS(93))+( STATES(6)*CONSTS(5)*STATES(12))/( CONSTS(87)*CONSTS(88)*CONSTS(94))+( CONSTS(5)*STATES(5)*STATES(12))/( CONSTS(95)*CONSTS(92)*CONSTS(93)))
      ALGBRC(30) = ( CONSTS(9)*(( CONSTS(122)*STATES(19))/CONSTS(124) - ( CONSTS(123)*STATES(21))/CONSTS(125)))/(1.00000+STATES(19)/CONSTS(124)+STATES(21)/CONSTS(125))
      ALGBRC(31) = ( CONSTS(9)*(( CONSTS(126)*STATES(13)*STATES(22))/( CONSTS(128)*CONSTS(129)) - ( CONSTS(127)*STATES(3)*STATES(23))/( CONSTS(130)*CONSTS(131))))/( (1.00000+STATES(13)/CONSTS(128)+STATES(3)/CONSTS(130))*(1.00000+STATES(22)/CONSTS(129)+STATES(23)/CONSTS(131)))
      ALGBRC(32) = ( CONSTS(9)*(( CONSTS(132)*STATES(23)*STATES(21))/( CONSTS(135)*CONSTS(134)) - ( CONSTS(133)*STATES(3)*STATES(13))/( CONSTS(136)*CONSTS(137))))/( (1.00000+STATES(21)/CONSTS(134)+STATES(13)/CONSTS(137))*(1.00000+STATES(23)/CONSTS(135)+STATES(3)/CONSTS(136)))
      ALGBRC(2) = STATES(18)/STATES(17)
      RETURN
      END