Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
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