C
C There are a total of 123 entries in the algebraic variable array.
C There are a total of 17 entries in each of the rate and state variable arrays.
C There are a total of 59 entries in the constant variable array.
C
C
C VOI is t in component environment (second).
C CONSTS(1) is C_m in component environment (fF).
C STATES(1) is q_Ca_o in component environment (fmol).
C STATES(2) is q_Ca_i in component environment (fmol).
C STATES(3) is q_K_o in component environment (fmol).
C STATES(4) is q_K_i in component environment (fmol).
C STATES(5) is q_S000_LCC in component environment (fmol).
C STATES(6) is q_S010_LCC in component environment (fmol).
C STATES(7) is q_S020_LCC in component environment (fmol).
C STATES(8) is q_S100_LCC in component environment (fmol).
C STATES(9) is q_S110_LCC in component environment (fmol).
C STATES(10) is q_S120_LCC in component environment (fmol).
C STATES(11) is q_S001_LCC in component environment (fmol).
C STATES(12) is q_S011_LCC in component environment (fmol).
C STATES(13) is q_S021_LCC in component environment (fmol).
C STATES(14) is q_S101_LCC in component environment (fmol).
C STATES(15) is q_S111_LCC in component environment (fmol).
C STATES(16) is q_S121_LCC in component environment (fmol).
C STATES(17) is q_mem in component environment (fC).
C ALGBRC(120) is I_mem_LCC in component LCC (fA).
C ALGBRC(116) is v_Ca_o_LCC in component LCC (fmol_per_sec).
C ALGBRC(117) is v_Ca_i_LCC in component LCC (fmol_per_sec).
C ALGBRC(121) is v_K_o_LCC in component LCC (fmol_per_sec).
C ALGBRC(122) is v_K_i_LCC in component LCC (fmol_per_sec).
C CONSTS(2) is kappa_LCC_Ca1 in component LCC_parameters (fmol_per_sec).
C CONSTS(3) is kappa_LCC_Ca2 in component LCC_parameters (fmol_per_sec).
C CONSTS(4) is kappa_LCC_K1 in component LCC_parameters (fmol_per_sec).
C CONSTS(5) is kappa_LCC_K2 in component LCC_parameters (fmol_per_sec).
C CONSTS(6) is kappa_d000 in component LCC_parameters (fmol_per_sec).
C CONSTS(7) is kappa_d010 in component LCC_parameters (fmol_per_sec).
C CONSTS(8) is kappa_d020 in component LCC_parameters (fmol_per_sec).
C CONSTS(9) is kappa_d001 in component LCC_parameters (fmol_per_sec).
C CONSTS(10) is kappa_d011 in component LCC_parameters (fmol_per_sec).
C CONSTS(11) is kappa_d021 in component LCC_parameters (fmol_per_sec).
C CONSTS(12) is kappa_f1_000 in component LCC_parameters (fmol_per_sec).
C CONSTS(13) is kappa_f1_100 in component LCC_parameters (fmol_per_sec).
C CONSTS(14) is kappa_f1_001 in component LCC_parameters (fmol_per_sec).
C CONSTS(15) is kappa_f1_101 in component LCC_parameters (fmol_per_sec).
C CONSTS(16) is kappa_f2_000 in component LCC_parameters (fmol_per_sec).
C CONSTS(17) is kappa_f2_100 in component LCC_parameters (fmol_per_sec).
C CONSTS(18) is kappa_f2_001 in component LCC_parameters (fmol_per_sec).
C CONSTS(19) is kappa_f2_101 in component LCC_parameters (fmol_per_sec).
C CONSTS(20) is kappa_f3_010 in component LCC_parameters (fmol_per_sec).
C CONSTS(21) is kappa_f3_110 in component LCC_parameters (fmol_per_sec).
C CONSTS(22) is kappa_f3_011 in component LCC_parameters (fmol_per_sec).
C CONSTS(23) is kappa_f3_111 in component LCC_parameters (fmol_per_sec).
C CONSTS(24) is kappa_fCa000 in component LCC_parameters (fmol_per_sec).
C CONSTS(25) is kappa_fCa100 in component LCC_parameters (fmol_per_sec).
C CONSTS(26) is kappa_fCa010 in component LCC_parameters (fmol_per_sec).
C CONSTS(27) is kappa_fCa110 in component LCC_parameters (fmol_per_sec).
C CONSTS(28) is kappa_fCa020 in component LCC_parameters (fmol_per_sec).
C CONSTS(29) is kappa_fCa120 in component LCC_parameters (fmol_per_sec).
C CONSTS(30) is K_Ca_i in component LCC_parameters (per_fmol).
C CONSTS(31) is K_Ca_o in component LCC_parameters (per_fmol).
C CONSTS(32) is K_K_i in component LCC_parameters (per_fmol).
C CONSTS(33) is K_K_o in component LCC_parameters (per_fmol).
C CONSTS(34) is K_000_LCC in component LCC_parameters (per_fmol).
C CONSTS(35) is K_010_LCC in component LCC_parameters (per_fmol).
C CONSTS(36) is K_020_LCC in component LCC_parameters (per_fmol).
C CONSTS(37) is K_100_LCC in component LCC_parameters (per_fmol).
C CONSTS(38) is K_110_LCC in component LCC_parameters (per_fmol).
C CONSTS(39) is K_120_LCC in component LCC_parameters (per_fmol).
C CONSTS(40) is K_001_LCC in component LCC_parameters (per_fmol).
C CONSTS(41) is K_011_LCC in component LCC_parameters (per_fmol).
C CONSTS(42) is K_021_LCC in component LCC_parameters (per_fmol).
C CONSTS(43) is K_101_LCC in component LCC_parameters (per_fmol).
C CONSTS(44) is K_111_LCC in component LCC_parameters (per_fmol).
C CONSTS(45) is K_121_LCC in component LCC_parameters (per_fmol).
C CONSTS(46) is zCa in component LCC_parameters (dimensionless).
C CONSTS(47) is zK in component LCC_parameters (dimensionless).
C CONSTS(48) is z_fd in component LCC_parameters (dimensionless).
C CONSTS(49) is z_ff1 in component LCC_parameters (dimensionless).
C CONSTS(50) is z_ff2 in component LCC_parameters (dimensionless).
C CONSTS(51) is z_ff3 in component LCC_parameters (dimensionless).
C CONSTS(52) is z_rCa in component LCC_parameters (dimensionless).
C CONSTS(53) is z_rd in component LCC_parameters (dimensionless).
C CONSTS(54) is z_rf1 in component LCC_parameters (dimensionless).
C CONSTS(55) is z_rf2 in component LCC_parameters (dimensionless).
C CONSTS(56) is z_rf3 in component LCC_parameters (dimensionless).
C CONSTS(57) is R in component constants (J_per_K_per_mol).
C CONSTS(58) is T in component constants (kelvin).
C CONSTS(59) is F in component constants (C_per_mol).
C ALGBRC(7) is mu_Ca_o in component LCC (J_per_mol).
C ALGBRC(8) is mu_Ca_i in component LCC (J_per_mol).
C ALGBRC(9) is mu_K_o in component LCC (J_per_mol).
C ALGBRC(10) is mu_K_i in component LCC (J_per_mol).
C ALGBRC(1) is V_mem in component LCC (volt).
C ALGBRC(19) is Af_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(51) is Ar_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(47) is Am_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(112) is v_LCC_Ca1 in component LCC (fmol_per_sec).
C ALGBRC(20) is Af_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(52) is Ar_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(48) is Am_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(114) is v_LCC_Ca2 in component LCC (fmol_per_sec).
C ALGBRC(21) is Af_LCC_K1 in component LCC (J_per_mol).
C ALGBRC(53) is Ar_LCC_K1 in component LCC (J_per_mol).
C ALGBRC(49) is Am_LCC_K1 in component LCC (J_per_mol).
C ALGBRC(115) is v_LCC_K1 in component LCC (fmol_per_sec).
C ALGBRC(22) is Af_LCC_K2 in component LCC (J_per_mol).
C ALGBRC(54) is Ar_LCC_K2 in component LCC (J_per_mol).
C ALGBRC(50) is Am_LCC_K2 in component LCC (J_per_mol).
C ALGBRC(118) is v_LCC_K2 in component LCC (fmol_per_sec).
C ALGBRC(41) is Af_fCa00 in component LCC (J_per_mol).
C ALGBRC(91) is Ar_fCa00 in component LCC (J_per_mol).
C ALGBRC(92) is v_fCa00 in component LCC (fmol_per_sec).
C ALGBRC(42) is Af_fCa01 in component LCC (J_per_mol).
C ALGBRC(93) is Ar_fCa01 in component LCC (J_per_mol).
C ALGBRC(96) is v_fCa01 in component LCC (fmol_per_sec).
C ALGBRC(43) is Af_fCa02 in component LCC (J_per_mol).
C ALGBRC(97) is Ar_fCa02 in component LCC (J_per_mol).
C ALGBRC(100) is v_fCa02 in component LCC (fmol_per_sec).
C ALGBRC(44) is Af_fCa10 in component LCC (J_per_mol).
C ALGBRC(101) is Ar_fCa10 in component LCC (J_per_mol).
C ALGBRC(104) is v_fCa10 in component LCC (fmol_per_sec).
C ALGBRC(45) is Af_fCa11 in component LCC (J_per_mol).
C ALGBRC(105) is Ar_fCa11 in component LCC (J_per_mol).
C ALGBRC(108) is v_fCa11 in component LCC (fmol_per_sec).
C ALGBRC(46) is Af_fCa12 in component LCC (J_per_mol).
C ALGBRC(109) is Ar_fCa12 in component LCC (J_per_mol).
C ALGBRC(111) is v_fCa12 in component LCC (fmol_per_sec).
C ALGBRC(11) is mu_S000_LCC in component LCC (J_per_mol).
C ALGBRC(94) is v_S000_LCC in component LCC (fmol_per_sec).
C ALGBRC(12) is mu_S010_LCC in component LCC (J_per_mol).
C ALGBRC(98) is v_S010_LCC in component LCC (fmol_per_sec).
C ALGBRC(13) is mu_S020_LCC in component LCC (J_per_mol).
C ALGBRC(102) is v_S020_LCC in component LCC (fmol_per_sec).
C ALGBRC(14) is mu_S100_LCC in component LCC (J_per_mol).
C ALGBRC(106) is v_S100_LCC in component LCC (fmol_per_sec).
C ALGBRC(15) is mu_S110_LCC in component LCC (J_per_mol).
C ALGBRC(110) is v_S110_LCC in component LCC (fmol_per_sec).
C ALGBRC(17) is mu_S120_LCC in component LCC (J_per_mol).
C ALGBRC(113) is v_S120_LCC in component LCC (fmol_per_sec).
C ALGBRC(2) is mu_S001_LCC in component LCC (J_per_mol).
C ALGBRC(95) is v_S001_LCC in component LCC (fmol_per_sec).
C ALGBRC(3) is mu_S011_LCC in component LCC (J_per_mol).
C ALGBRC(99) is v_S011_LCC in component LCC (fmol_per_sec).
C ALGBRC(4) is mu_S021_LCC in component LCC (J_per_mol).
C ALGBRC(103) is v_S021_LCC in component LCC (fmol_per_sec).
C ALGBRC(5) is mu_S101_LCC in component LCC (J_per_mol).
C ALGBRC(107) is v_S101_LCC in component LCC (fmol_per_sec).
C ALGBRC(16) is mu_S111_LCC in component LCC (J_per_mol).
C ALGBRC(119) is v_S111_LCC in component LCC (fmol_per_sec).
C ALGBRC(18) is mu_S121_LCC in component LCC (J_per_mol).
C ALGBRC(123) is v_S121_LCC in component LCC (fmol_per_sec).
C ALGBRC(23) is Af_d000 in component LCC (J_per_mol).
C ALGBRC(55) is Ar_d000 in component LCC (J_per_mol).
C ALGBRC(56) is v_d000 in component LCC (fmol_per_sec).
C ALGBRC(25) is Af_d010 in component LCC (J_per_mol).
C ALGBRC(59) is Ar_d010 in component LCC (J_per_mol).
C ALGBRC(60) is v_d010 in component LCC (fmol_per_sec).
C ALGBRC(27) is Af_d020 in component LCC (J_per_mol).
C ALGBRC(63) is Ar_d020 in component LCC (J_per_mol).
C ALGBRC(64) is v_d020 in component LCC (fmol_per_sec).
C ALGBRC(29) is Af_f1_000 in component LCC (J_per_mol).
C ALGBRC(67) is Ar_f1_000 in component LCC (J_per_mol).
C ALGBRC(68) is v_f1_000 in component LCC (fmol_per_sec).
C ALGBRC(33) is Af_f2_000 in component LCC (J_per_mol).
C ALGBRC(75) is Ar_f2_000 in component LCC (J_per_mol).
C ALGBRC(76) is v_f2_000 in component LCC (fmol_per_sec).
C ALGBRC(37) is Af_f3_010 in component LCC (J_per_mol).
C ALGBRC(83) is Ar_f3_010 in component LCC (J_per_mol).
C ALGBRC(84) is v_f3_010 in component LCC (fmol_per_sec).
C ALGBRC(31) is Af_f1_100 in component LCC (J_per_mol).
C ALGBRC(71) is Ar_f1_100 in component LCC (J_per_mol).
C ALGBRC(72) is v_f1_100 in component LCC (fmol_per_sec).
C ALGBRC(35) is Af_f2_100 in component LCC (J_per_mol).
C ALGBRC(79) is Ar_f2_100 in component LCC (J_per_mol).
C ALGBRC(80) is v_f2_100 in component LCC (fmol_per_sec).
C ALGBRC(39) is Af_f3_110 in component LCC (J_per_mol).
C ALGBRC(87) is Ar_f3_110 in component LCC (J_per_mol).
C ALGBRC(88) is v_f3_110 in component LCC (fmol_per_sec).
C ALGBRC(24) is Af_d001 in component LCC (J_per_mol).
C ALGBRC(57) is Ar_d001 in component LCC (J_per_mol).
C ALGBRC(58) is v_d001 in component LCC (fmol_per_sec).
C ALGBRC(26) is Af_d011 in component LCC (J_per_mol).
C ALGBRC(61) is Ar_d011 in component LCC (J_per_mol).
C ALGBRC(62) is v_d011 in component LCC (fmol_per_sec).
C ALGBRC(28) is Af_d021 in component LCC (J_per_mol).
C ALGBRC(65) is Ar_d021 in component LCC (J_per_mol).
C ALGBRC(66) is v_d021 in component LCC (fmol_per_sec).
C ALGBRC(30) is Af_f1_001 in component LCC (J_per_mol).
C ALGBRC(69) is Ar_f1_001 in component LCC (J_per_mol).
C ALGBRC(70) is v_f1_001 in component LCC (fmol_per_sec).
C ALGBRC(34) is Af_f2_001 in component LCC (J_per_mol).
C ALGBRC(77) is Ar_f2_001 in component LCC (J_per_mol).
C ALGBRC(78) is v_f2_001 in component LCC (fmol_per_sec).
C ALGBRC(38) is Af_f3_011 in component LCC (J_per_mol).
C ALGBRC(85) is Ar_f3_011 in component LCC (J_per_mol).
C ALGBRC(86) is v_f3_011 in component LCC (fmol_per_sec).
C ALGBRC(32) is Af_f1_101 in component LCC (J_per_mol).
C ALGBRC(73) is Ar_f1_101 in component LCC (J_per_mol).
C ALGBRC(74) is v_f1_101 in component LCC (fmol_per_sec).
C ALGBRC(36) is Af_f2_101 in component LCC (J_per_mol).
C ALGBRC(81) is Ar_f2_101 in component LCC (J_per_mol).
C ALGBRC(82) is v_f2_101 in component LCC (fmol_per_sec).
C ALGBRC(40) is Af_f3_111 in component LCC (J_per_mol).
C ALGBRC(89) is Ar_f3_111 in component LCC (J_per_mol).
C ALGBRC(90) is v_f3_111 in component LCC (fmol_per_sec).
C ALGBRC(6) is Ca_tot in component LCC (fmol).
C RATES(1) is d/dt q_Ca_o in component environment (fmol).
C RATES(2) is d/dt q_Ca_i in component environment (fmol).
C RATES(3) is d/dt q_K_o in component environment (fmol).
C RATES(4) is d/dt q_K_i in component environment (fmol).
C RATES(17) is d/dt q_mem in component environment (fC).
C RATES(5) is d/dt q_S000_LCC in component environment (fmol).
C RATES(6) is d/dt q_S010_LCC in component environment (fmol).
C RATES(7) is d/dt q_S020_LCC in component environment (fmol).
C RATES(8) is d/dt q_S100_LCC in component environment (fmol).
C RATES(9) is d/dt q_S110_LCC in component environment (fmol).
C RATES(10) is d/dt q_S120_LCC in component environment (fmol).
C RATES(11) is d/dt q_S001_LCC in component environment (fmol).
C RATES(12) is d/dt q_S011_LCC in component environment (fmol).
C RATES(13) is d/dt q_S021_LCC in component environment (fmol).
C RATES(14) is d/dt q_S101_LCC in component environment (fmol).
C RATES(15) is d/dt q_S111_LCC in component environment (fmol).
C RATES(16) is d/dt q_S121_LCC in component environment (fmol).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 60000
      STATES(1) = 0.245463158
      STATES(2) = 6e-5
      STATES(3) = 27.9828
      STATES(4) = 5510
      STATES(5) = 5.30481E-08
      STATES(6) = 1.05036E-05
      STATES(7) = 5.30481E-08
      STATES(8) = 5.35845E-10
      STATES(9) = 1.06098E-07
      STATES(10) = 5.35845E-10
      STATES(11) = 5.89425E-09
      STATES(12) = 1.16708E-06
      STATES(13) = 5.89425E-09
      STATES(14) = 5.95377E-11
      STATES(15) = 1.17886E-08
      STATES(16) = 5.95377E-11
      STATES(17) = -13039
      CONSTS(2) = 74.6745
      CONSTS(3) = 98.3196
      CONSTS(4) = 0.131085
      CONSTS(5) = 0.172592
      CONSTS(6) = 53.6616
      CONSTS(7) = 2.7099
      CONSTS(8) = 3.56797
      CONSTS(9) = 0.28538
      CONSTS(10) = 0.0144117
      CONSTS(11) = 0.018975
      CONSTS(12) = 5.16892
      CONSTS(13) = 25.6676
      CONSTS(14) = 0.0274891
      CONSTS(15) = 0.136504
      CONSTS(16) = 0.335961
      CONSTS(17) = 1.6683
      CONSTS(18) = 0.00178669
      CONSTS(19) = 0.00887227
      CONSTS(20) = 733.001
      CONSTS(21) = 3639.9
      CONSTS(22) = 3.89821
      CONSTS(23) = 19.3575
      CONSTS(24) = 97713.8
      CONSTS(25) = 485222
      CONSTS(26) = 4934.53
      CONSTS(27) = 24503.7
      CONSTS(28) = 6497.01
      CONSTS(29) = 32262.5
      CONSTS(30) = 0.090536
      CONSTS(31) = 0.663931
      CONSTS(32) = 0.0405426
      CONSTS(33) = 0.297313
      CONSTS(34) = 9.07096
      CONSTS(35) = 179.624
      CONSTS(36) = 136.425
      CONSTS(37) = 1.8267
      CONSTS(38) = 36.1725
      CONSTS(39) = 27.4733
      CONSTS(40) = 1705.66
      CONSTS(41) = 33775.6
      CONSTS(42) = 25652.8
      CONSTS(43) = 343.485
      CONSTS(44) = 6801.71
      CONSTS(45) = 5165.95
      CONSTS(46) = 2
      CONSTS(47) = 1
      CONSTS(48) = 2.1404
      CONSTS(49) = -1.1495
      CONSTS(50) = 0.72162
      CONSTS(51) = 4.2933
      CONSTS(52) = 2
      CONSTS(53) = -2.1404
      CONSTS(54) = 1.8993
      CONSTS(55) = -0.52288
      CONSTS(56) = 0
      CONSTS(57) = 8.31
      CONSTS(58) = 310
      CONSTS(59) = 96485
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(11) =  CONSTS(57)*CONSTS(58)*log( CONSTS(34)*STATES(5))
      ALGBRC(41) = ALGBRC(11)
      ALGBRC(8) =  CONSTS(57)*CONSTS(58)*log( CONSTS(30)*STATES(2))
      ALGBRC(91) =  CONSTS(57)*CONSTS(58)*log( CONSTS(40)*STATES(11))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(92) =  CONSTS(24)*(EXP(ALGBRC(41)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(91)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(1) = STATES(17)/CONSTS(1)
      ALGBRC(23) =  CONSTS(57)*CONSTS(58)*log( CONSTS(34)*STATES(5))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(55) =  CONSTS(57)*CONSTS(58)*log( CONSTS(37)*STATES(8))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(56) =  CONSTS(6)*(EXP(ALGBRC(23)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(55)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(29) =  CONSTS(57)*CONSTS(58)*log( CONSTS(34)*STATES(5))+ CONSTS(49)*CONSTS(59)*ALGBRC(1)
      ALGBRC(67) =  CONSTS(57)*CONSTS(58)*log( CONSTS(35)*STATES(6))+ CONSTS(54)*CONSTS(59)*ALGBRC(1)
      ALGBRC(68) =  CONSTS(12)*(EXP(ALGBRC(29)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(67)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(33) =  CONSTS(57)*CONSTS(58)*log( CONSTS(34)*STATES(5))+ CONSTS(50)*CONSTS(59)*ALGBRC(1)
      ALGBRC(75) =  CONSTS(57)*CONSTS(58)*log( CONSTS(36)*STATES(7))+ CONSTS(55)*CONSTS(59)*ALGBRC(1)
      ALGBRC(76) =  CONSTS(16)*(EXP(ALGBRC(33)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(75)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(94) = ((- ALGBRC(56) - ALGBRC(68)) - ALGBRC(76)) - ALGBRC(92)
      RATES(5) = ALGBRC(94)
      ALGBRC(24) =  CONSTS(57)*CONSTS(58)*log( CONSTS(40)*STATES(11))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(57) =  CONSTS(57)*CONSTS(58)*log( CONSTS(43)*STATES(14))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(58) =  CONSTS(9)*(EXP(ALGBRC(24)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(57)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(30) =  CONSTS(57)*CONSTS(58)*log( CONSTS(40)*STATES(11))+ CONSTS(49)*CONSTS(59)*ALGBRC(1)
      ALGBRC(69) =  CONSTS(57)*CONSTS(58)*log( CONSTS(41)*STATES(12))+ CONSTS(54)*CONSTS(59)*ALGBRC(1)
      ALGBRC(70) =  CONSTS(14)*(EXP(ALGBRC(30)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(69)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(34) =  CONSTS(57)*CONSTS(58)*log( CONSTS(40)*STATES(11))+ CONSTS(50)*CONSTS(59)*ALGBRC(1)
      ALGBRC(77) =  CONSTS(57)*CONSTS(58)*log( CONSTS(42)*STATES(13))+ CONSTS(55)*CONSTS(59)*ALGBRC(1)
      ALGBRC(78) =  CONSTS(18)*(EXP(ALGBRC(34)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(77)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(95) = ((- ALGBRC(58) - ALGBRC(70)) - ALGBRC(78)) - - ALGBRC(92)
      RATES(11) = ALGBRC(95)
      ALGBRC(12) =  CONSTS(57)*CONSTS(58)*log( CONSTS(35)*STATES(6))
      ALGBRC(42) = ALGBRC(12)
      ALGBRC(93) =  CONSTS(57)*CONSTS(58)*log( CONSTS(41)*STATES(12))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(96) =  CONSTS(26)*(EXP(ALGBRC(42)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(93)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(25) =  CONSTS(57)*CONSTS(58)*log( CONSTS(35)*STATES(6))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(59) =  CONSTS(57)*CONSTS(58)*log( CONSTS(38)*STATES(9))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(60) =  CONSTS(7)*(EXP(ALGBRC(25)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(59)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(37) =  CONSTS(57)*CONSTS(58)*log( CONSTS(35)*STATES(6))+ CONSTS(51)*CONSTS(59)*ALGBRC(1)
      ALGBRC(83) =  CONSTS(57)*CONSTS(58)*log( CONSTS(36)*STATES(7))+ CONSTS(56)*CONSTS(59)*ALGBRC(1)
      ALGBRC(84) =  CONSTS(20)*(EXP(ALGBRC(37)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(83)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(98) = ((- ALGBRC(60)+ALGBRC(68)) - ALGBRC(96)) - ALGBRC(84)
      RATES(6) = ALGBRC(98)
      ALGBRC(26) =  CONSTS(57)*CONSTS(58)*log( CONSTS(41)*STATES(12))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(61) =  CONSTS(57)*CONSTS(58)*log( CONSTS(44)*STATES(15))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(62) =  CONSTS(10)*(EXP(ALGBRC(26)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(61)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(38) =  CONSTS(57)*CONSTS(58)*log( CONSTS(41)*STATES(12))+ CONSTS(51)*CONSTS(59)*ALGBRC(1)
      ALGBRC(85) =  CONSTS(57)*CONSTS(58)*log( CONSTS(42)*STATES(13))+ CONSTS(56)*CONSTS(59)*ALGBRC(1)
      ALGBRC(86) =  CONSTS(22)*(EXP(ALGBRC(38)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(85)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(99) = ((- ALGBRC(62)+ALGBRC(70)) - - ALGBRC(96)) - ALGBRC(86)
      RATES(12) = ALGBRC(99)
      ALGBRC(13) =  CONSTS(57)*CONSTS(58)*log( CONSTS(36)*STATES(7))
      ALGBRC(43) = ALGBRC(13)
      ALGBRC(97) =  CONSTS(57)*CONSTS(58)*log( CONSTS(42)*STATES(13))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(100) =  CONSTS(28)*(EXP(ALGBRC(43)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(97)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(27) =  CONSTS(57)*CONSTS(58)*log( CONSTS(36)*STATES(7))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(63) =  CONSTS(57)*CONSTS(58)*log( CONSTS(39)*STATES(10))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(64) =  CONSTS(8)*(EXP(ALGBRC(27)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(63)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(102) = ((- ALGBRC(64)+ALGBRC(76)) - ALGBRC(100))+ALGBRC(84)
      RATES(7) = ALGBRC(102)
      ALGBRC(28) =  CONSTS(57)*CONSTS(58)*log( CONSTS(42)*STATES(13))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(65) =  CONSTS(57)*CONSTS(58)*log( CONSTS(45)*STATES(16))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(66) =  CONSTS(11)*(EXP(ALGBRC(28)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(65)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(103) = ((- ALGBRC(66)+ALGBRC(78)) - - ALGBRC(100))+ALGBRC(86)
      RATES(13) = ALGBRC(103)
      ALGBRC(14) =  CONSTS(57)*CONSTS(58)*log( CONSTS(37)*STATES(8))
      ALGBRC(44) = ALGBRC(14)
      ALGBRC(101) =  CONSTS(57)*CONSTS(58)*log( CONSTS(43)*STATES(14))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(104) =  CONSTS(25)*(EXP(ALGBRC(44)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(101)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(31) =  CONSTS(57)*CONSTS(58)*log( CONSTS(37)*STATES(8))+ CONSTS(49)*CONSTS(59)*ALGBRC(1)
      ALGBRC(71) =  CONSTS(57)*CONSTS(58)*log( CONSTS(38)*STATES(9))+ CONSTS(54)*CONSTS(59)*ALGBRC(1)
      ALGBRC(72) =  CONSTS(13)*(EXP(ALGBRC(31)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(71)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(35) =  CONSTS(57)*CONSTS(58)*log( CONSTS(37)*STATES(8))+ CONSTS(50)*CONSTS(59)*ALGBRC(1)
      ALGBRC(79) =  CONSTS(57)*CONSTS(58)*log( CONSTS(39)*STATES(10))+ CONSTS(55)*CONSTS(59)*ALGBRC(1)
      ALGBRC(80) =  CONSTS(17)*(EXP(ALGBRC(35)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(79)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(106) = ((- ALGBRC(80)+ALGBRC(56)) - ALGBRC(72)) - ALGBRC(104)
      RATES(8) = ALGBRC(106)
      ALGBRC(32) =  CONSTS(57)*CONSTS(58)*log( CONSTS(43)*STATES(14))+ CONSTS(49)*CONSTS(59)*ALGBRC(1)
      ALGBRC(73) =  CONSTS(57)*CONSTS(58)*log( CONSTS(44)*STATES(15))+ CONSTS(54)*CONSTS(59)*ALGBRC(1)
      ALGBRC(74) =  CONSTS(15)*(EXP(ALGBRC(32)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(73)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(36) =  CONSTS(57)*CONSTS(58)*log( CONSTS(43)*STATES(14))+ CONSTS(50)*CONSTS(59)*ALGBRC(1)
      ALGBRC(81) =  CONSTS(57)*CONSTS(58)*log( CONSTS(45)*STATES(16))+ CONSTS(55)*CONSTS(59)*ALGBRC(1)
      ALGBRC(82) =  CONSTS(19)*(EXP(ALGBRC(36)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(81)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(107) = ((- ALGBRC(82)+ALGBRC(58)) - ALGBRC(74)) - - ALGBRC(104)
      RATES(14) = ALGBRC(107)
      ALGBRC(15) =  CONSTS(57)*CONSTS(58)*log( CONSTS(38)*STATES(9))
      ALGBRC(45) = ALGBRC(15)
      ALGBRC(105) =  CONSTS(57)*CONSTS(58)*log( CONSTS(44)*STATES(15))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(108) =  CONSTS(27)*(EXP(ALGBRC(45)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(105)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(39) =  CONSTS(57)*CONSTS(58)*log( CONSTS(38)*STATES(9))+ CONSTS(51)*CONSTS(59)*ALGBRC(1)
      ALGBRC(87) =  CONSTS(57)*CONSTS(58)*log( CONSTS(39)*STATES(10))+ CONSTS(56)*CONSTS(59)*ALGBRC(1)
      ALGBRC(88) =  CONSTS(21)*(EXP(ALGBRC(39)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(87)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(110) = ((ALGBRC(72)+ALGBRC(60)) - ALGBRC(108)) - ALGBRC(88)
      RATES(9) = ALGBRC(110)
      ALGBRC(17) =  CONSTS(57)*CONSTS(58)*log( CONSTS(39)*STATES(10))
      ALGBRC(46) = ALGBRC(17)
      ALGBRC(109) =  CONSTS(57)*CONSTS(58)*log( CONSTS(45)*STATES(16))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(111) =  CONSTS(29)*(EXP(ALGBRC(46)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(109)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(113) = ((ALGBRC(80)+ALGBRC(64)) - ALGBRC(111))+ALGBRC(88)
      RATES(10) = ALGBRC(113)
      ALGBRC(16) =  CONSTS(57)*CONSTS(58)*log( CONSTS(44)*STATES(15))
      ALGBRC(19) = ALGBRC(8)+ CONSTS(46)*CONSTS(59)*ALGBRC(1)+ALGBRC(16)
      ALGBRC(7) =  CONSTS(57)*CONSTS(58)*log( CONSTS(31)*STATES(1))
      ALGBRC(51) = ALGBRC(7)+ALGBRC(16)
      ALGBRC(47) =  CONSTS(46)*CONSTS(59)*ALGBRC(1)
      ALGBRC(112) = TERNRY(ALGBRC(47).EQ.0.00000,  CONSTS(2)*(EXP(ALGBRC(19)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(51)/( CONSTS(57)*CONSTS(58)))),  ((( CONSTS(2)*ALGBRC(47))/( CONSTS(57)*CONSTS(58)))/(EXP(ALGBRC(47)/( CONSTS(57)*CONSTS(58))) - 1.00000))*(EXP(ALGBRC(19)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(51)/( CONSTS(57)*CONSTS(58)))))
      ALGBRC(18) =  CONSTS(57)*CONSTS(58)*log( CONSTS(45)*STATES(16))
      ALGBRC(20) = ALGBRC(8)+ CONSTS(46)*CONSTS(59)*ALGBRC(1)+ALGBRC(18)
      ALGBRC(52) = ALGBRC(7)+ALGBRC(18)
      ALGBRC(48) =  CONSTS(46)*CONSTS(59)*ALGBRC(1)
      ALGBRC(114) = TERNRY(ALGBRC(48).EQ.0.00000,  CONSTS(3)*(EXP(ALGBRC(20)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(52)/( CONSTS(57)*CONSTS(58)))),  ((( CONSTS(3)*ALGBRC(48))/( CONSTS(57)*CONSTS(58)))/(EXP(ALGBRC(48)/( CONSTS(57)*CONSTS(58))) - 1.00000))*(EXP(ALGBRC(20)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(52)/( CONSTS(57)*CONSTS(58)))))
      ALGBRC(116) = ALGBRC(114)+ALGBRC(112)
      RATES(1) = ALGBRC(116)
      ALGBRC(117) = (( (CONSTS(52))*ALGBRC(104)+ CONSTS(52)*ALGBRC(108)+ CONSTS(52)*ALGBRC(100)+ CONSTS(52)*ALGBRC(96)+ CONSTS(52)*ALGBRC(92)+ CONSTS(52)*ALGBRC(111)) - ALGBRC(112)) - ALGBRC(114)
      RATES(2) = ALGBRC(117)
      ALGBRC(10) =  CONSTS(57)*CONSTS(58)*log( CONSTS(32)*STATES(4))
      ALGBRC(21) = ALGBRC(10)+ CONSTS(47)*CONSTS(59)*ALGBRC(1)+ALGBRC(16)
      ALGBRC(9) =  CONSTS(57)*CONSTS(58)*log( CONSTS(33)*STATES(3))
      ALGBRC(53) = ALGBRC(9)+ALGBRC(16)
      ALGBRC(49) =  CONSTS(47)*CONSTS(59)*ALGBRC(1)
      ALGBRC(115) = TERNRY(ALGBRC(49).EQ.0.00000,  CONSTS(4)*(EXP(ALGBRC(21)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(53)/( CONSTS(57)*CONSTS(58)))),  ((( CONSTS(4)*ALGBRC(49))/( CONSTS(57)*CONSTS(58)))/(EXP(ALGBRC(49)/( CONSTS(57)*CONSTS(58))) - 1.00000))*(EXP(ALGBRC(21)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(53)/( CONSTS(57)*CONSTS(58)))))
      ALGBRC(40) =  CONSTS(57)*CONSTS(58)*log( CONSTS(44)*STATES(15))+ CONSTS(51)*CONSTS(59)*ALGBRC(1)
      ALGBRC(89) =  CONSTS(57)*CONSTS(58)*log( CONSTS(45)*STATES(16))+ CONSTS(56)*CONSTS(59)*ALGBRC(1)
      ALGBRC(90) =  CONSTS(23)*(EXP(ALGBRC(40)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(89)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(119) = ((ALGBRC(74)+ALGBRC(62)) - ((((- ALGBRC(108)+ALGBRC(112)) - ALGBRC(112))+ALGBRC(115)) - ALGBRC(115))) - ALGBRC(90)
      RATES(15) = ALGBRC(119)
      ALGBRC(22) = ALGBRC(10)+ CONSTS(47)*CONSTS(59)*ALGBRC(1)+ALGBRC(18)
      ALGBRC(54) = ALGBRC(9)+ALGBRC(18)
      ALGBRC(50) =  CONSTS(47)*CONSTS(59)*ALGBRC(1)
      ALGBRC(118) = TERNRY(ALGBRC(50).EQ.0.00000,  CONSTS(5)*(EXP(ALGBRC(22)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(54)/( CONSTS(57)*CONSTS(58)))),  ((( CONSTS(5)*ALGBRC(50))/( CONSTS(57)*CONSTS(58)))/(EXP(ALGBRC(50)/( CONSTS(57)*CONSTS(58))) - 1.00000))*(EXP(ALGBRC(22)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(54)/( CONSTS(57)*CONSTS(58)))))
      ALGBRC(121) = ALGBRC(118)+ALGBRC(115)
      RATES(3) = ALGBRC(121)
      ALGBRC(122) = - ALGBRC(115) - ALGBRC(118)
      RATES(4) = ALGBRC(122)
      ALGBRC(120) =  CONSTS(59)*((((((((((((((((((((((((((((((((((((( CONSTS(55)*ALGBRC(78) -  CONSTS(50)*ALGBRC(78))+ CONSTS(56)*ALGBRC(86)) -  CONSTS(51)*ALGBRC(86))+ CONSTS(54)*ALGBRC(74)) -  CONSTS(49)*ALGBRC(74))+ CONSTS(55)*ALGBRC(82)) -  CONSTS(50)*ALGBRC(82))+ CONSTS(56)*ALGBRC(90)) -  CONSTS(51)*ALGBRC(90))+ CONSTS(54)*ALGBRC(70)) -  CONSTS(49)*ALGBRC(70))+ CONSTS(53)*ALGBRC(66)) -  CONSTS(48)*ALGBRC(66))+ CONSTS(53)*ALGBRC(62)) -  CONSTS(48)*ALGBRC(62))+ CONSTS(53)*ALGBRC(58)) -  CONSTS(48)*ALGBRC(58))+ CONSTS(56)*ALGBRC(88)) -  CONSTS(51)*ALGBRC(88))+ CONSTS(55)*ALGBRC(80)) -  CONSTS(50)*ALGBRC(80))+ CONSTS(54)*ALGBRC(72)) -  CONSTS(49)*ALGBRC(72))+ CONSTS(56)*ALGBRC(84)) -  CONSTS(51)*ALGBRC(84))+ CONSTS(55)*ALGBRC(76)) -  CONSTS(50)*ALGBRC(76))+ CONSTS(54)*ALGBRC(68)) -  CONSTS(49)*ALGBRC(68))+ CONSTS(53)*ALGBRC(64)) -  CONSTS(48)*ALGBRC(64))+ CONSTS(53)*ALGBRC(60)) -  CONSTS(48)*ALGBRC(60))+ CONSTS(53)*ALGBRC(56)) -  CONSTS(48)*ALGBRC(56)) - ( CONSTS(46)*ALGBRC(112)+ CONSTS(46)*ALGBRC(114))) - ( CONSTS(47)*ALGBRC(115)+ CONSTS(47)*ALGBRC(118)))
      RATES(17) = ALGBRC(120)
      ALGBRC(123) = ((ALGBRC(82)+ALGBRC(66)) - ((((- ALGBRC(111)+ALGBRC(114)) - ALGBRC(114))+ALGBRC(118)) - ALGBRC(118)))+ALGBRC(90)
      RATES(16) = ALGBRC(123)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(11) =  CONSTS(57)*CONSTS(58)*log( CONSTS(34)*STATES(5))
      ALGBRC(41) = ALGBRC(11)
      ALGBRC(8) =  CONSTS(57)*CONSTS(58)*log( CONSTS(30)*STATES(2))
      ALGBRC(91) =  CONSTS(57)*CONSTS(58)*log( CONSTS(40)*STATES(11))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(92) =  CONSTS(24)*(EXP(ALGBRC(41)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(91)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(1) = STATES(17)/CONSTS(1)
      ALGBRC(23) =  CONSTS(57)*CONSTS(58)*log( CONSTS(34)*STATES(5))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(55) =  CONSTS(57)*CONSTS(58)*log( CONSTS(37)*STATES(8))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(56) =  CONSTS(6)*(EXP(ALGBRC(23)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(55)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(29) =  CONSTS(57)*CONSTS(58)*log( CONSTS(34)*STATES(5))+ CONSTS(49)*CONSTS(59)*ALGBRC(1)
      ALGBRC(67) =  CONSTS(57)*CONSTS(58)*log( CONSTS(35)*STATES(6))+ CONSTS(54)*CONSTS(59)*ALGBRC(1)
      ALGBRC(68) =  CONSTS(12)*(EXP(ALGBRC(29)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(67)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(33) =  CONSTS(57)*CONSTS(58)*log( CONSTS(34)*STATES(5))+ CONSTS(50)*CONSTS(59)*ALGBRC(1)
      ALGBRC(75) =  CONSTS(57)*CONSTS(58)*log( CONSTS(36)*STATES(7))+ CONSTS(55)*CONSTS(59)*ALGBRC(1)
      ALGBRC(76) =  CONSTS(16)*(EXP(ALGBRC(33)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(75)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(94) = ((- ALGBRC(56) - ALGBRC(68)) - ALGBRC(76)) - ALGBRC(92)
      ALGBRC(24) =  CONSTS(57)*CONSTS(58)*log( CONSTS(40)*STATES(11))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(57) =  CONSTS(57)*CONSTS(58)*log( CONSTS(43)*STATES(14))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(58) =  CONSTS(9)*(EXP(ALGBRC(24)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(57)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(30) =  CONSTS(57)*CONSTS(58)*log( CONSTS(40)*STATES(11))+ CONSTS(49)*CONSTS(59)*ALGBRC(1)
      ALGBRC(69) =  CONSTS(57)*CONSTS(58)*log( CONSTS(41)*STATES(12))+ CONSTS(54)*CONSTS(59)*ALGBRC(1)
      ALGBRC(70) =  CONSTS(14)*(EXP(ALGBRC(30)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(69)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(34) =  CONSTS(57)*CONSTS(58)*log( CONSTS(40)*STATES(11))+ CONSTS(50)*CONSTS(59)*ALGBRC(1)
      ALGBRC(77) =  CONSTS(57)*CONSTS(58)*log( CONSTS(42)*STATES(13))+ CONSTS(55)*CONSTS(59)*ALGBRC(1)
      ALGBRC(78) =  CONSTS(18)*(EXP(ALGBRC(34)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(77)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(95) = ((- ALGBRC(58) - ALGBRC(70)) - ALGBRC(78)) - - ALGBRC(92)
      ALGBRC(12) =  CONSTS(57)*CONSTS(58)*log( CONSTS(35)*STATES(6))
      ALGBRC(42) = ALGBRC(12)
      ALGBRC(93) =  CONSTS(57)*CONSTS(58)*log( CONSTS(41)*STATES(12))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(96) =  CONSTS(26)*(EXP(ALGBRC(42)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(93)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(25) =  CONSTS(57)*CONSTS(58)*log( CONSTS(35)*STATES(6))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(59) =  CONSTS(57)*CONSTS(58)*log( CONSTS(38)*STATES(9))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(60) =  CONSTS(7)*(EXP(ALGBRC(25)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(59)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(37) =  CONSTS(57)*CONSTS(58)*log( CONSTS(35)*STATES(6))+ CONSTS(51)*CONSTS(59)*ALGBRC(1)
      ALGBRC(83) =  CONSTS(57)*CONSTS(58)*log( CONSTS(36)*STATES(7))+ CONSTS(56)*CONSTS(59)*ALGBRC(1)
      ALGBRC(84) =  CONSTS(20)*(EXP(ALGBRC(37)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(83)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(98) = ((- ALGBRC(60)+ALGBRC(68)) - ALGBRC(96)) - ALGBRC(84)
      ALGBRC(26) =  CONSTS(57)*CONSTS(58)*log( CONSTS(41)*STATES(12))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(61) =  CONSTS(57)*CONSTS(58)*log( CONSTS(44)*STATES(15))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(62) =  CONSTS(10)*(EXP(ALGBRC(26)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(61)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(38) =  CONSTS(57)*CONSTS(58)*log( CONSTS(41)*STATES(12))+ CONSTS(51)*CONSTS(59)*ALGBRC(1)
      ALGBRC(85) =  CONSTS(57)*CONSTS(58)*log( CONSTS(42)*STATES(13))+ CONSTS(56)*CONSTS(59)*ALGBRC(1)
      ALGBRC(86) =  CONSTS(22)*(EXP(ALGBRC(38)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(85)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(99) = ((- ALGBRC(62)+ALGBRC(70)) - - ALGBRC(96)) - ALGBRC(86)
      ALGBRC(13) =  CONSTS(57)*CONSTS(58)*log( CONSTS(36)*STATES(7))
      ALGBRC(43) = ALGBRC(13)
      ALGBRC(97) =  CONSTS(57)*CONSTS(58)*log( CONSTS(42)*STATES(13))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(100) =  CONSTS(28)*(EXP(ALGBRC(43)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(97)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(27) =  CONSTS(57)*CONSTS(58)*log( CONSTS(36)*STATES(7))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(63) =  CONSTS(57)*CONSTS(58)*log( CONSTS(39)*STATES(10))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(64) =  CONSTS(8)*(EXP(ALGBRC(27)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(63)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(102) = ((- ALGBRC(64)+ALGBRC(76)) - ALGBRC(100))+ALGBRC(84)
      ALGBRC(28) =  CONSTS(57)*CONSTS(58)*log( CONSTS(42)*STATES(13))+ CONSTS(48)*CONSTS(59)*ALGBRC(1)
      ALGBRC(65) =  CONSTS(57)*CONSTS(58)*log( CONSTS(45)*STATES(16))+ CONSTS(53)*CONSTS(59)*ALGBRC(1)
      ALGBRC(66) =  CONSTS(11)*(EXP(ALGBRC(28)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(65)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(103) = ((- ALGBRC(66)+ALGBRC(78)) - - ALGBRC(100))+ALGBRC(86)
      ALGBRC(14) =  CONSTS(57)*CONSTS(58)*log( CONSTS(37)*STATES(8))
      ALGBRC(44) = ALGBRC(14)
      ALGBRC(101) =  CONSTS(57)*CONSTS(58)*log( CONSTS(43)*STATES(14))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(104) =  CONSTS(25)*(EXP(ALGBRC(44)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(101)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(31) =  CONSTS(57)*CONSTS(58)*log( CONSTS(37)*STATES(8))+ CONSTS(49)*CONSTS(59)*ALGBRC(1)
      ALGBRC(71) =  CONSTS(57)*CONSTS(58)*log( CONSTS(38)*STATES(9))+ CONSTS(54)*CONSTS(59)*ALGBRC(1)
      ALGBRC(72) =  CONSTS(13)*(EXP(ALGBRC(31)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(71)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(35) =  CONSTS(57)*CONSTS(58)*log( CONSTS(37)*STATES(8))+ CONSTS(50)*CONSTS(59)*ALGBRC(1)
      ALGBRC(79) =  CONSTS(57)*CONSTS(58)*log( CONSTS(39)*STATES(10))+ CONSTS(55)*CONSTS(59)*ALGBRC(1)
      ALGBRC(80) =  CONSTS(17)*(EXP(ALGBRC(35)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(79)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(106) = ((- ALGBRC(80)+ALGBRC(56)) - ALGBRC(72)) - ALGBRC(104)
      ALGBRC(32) =  CONSTS(57)*CONSTS(58)*log( CONSTS(43)*STATES(14))+ CONSTS(49)*CONSTS(59)*ALGBRC(1)
      ALGBRC(73) =  CONSTS(57)*CONSTS(58)*log( CONSTS(44)*STATES(15))+ CONSTS(54)*CONSTS(59)*ALGBRC(1)
      ALGBRC(74) =  CONSTS(15)*(EXP(ALGBRC(32)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(73)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(36) =  CONSTS(57)*CONSTS(58)*log( CONSTS(43)*STATES(14))+ CONSTS(50)*CONSTS(59)*ALGBRC(1)
      ALGBRC(81) =  CONSTS(57)*CONSTS(58)*log( CONSTS(45)*STATES(16))+ CONSTS(55)*CONSTS(59)*ALGBRC(1)
      ALGBRC(82) =  CONSTS(19)*(EXP(ALGBRC(36)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(81)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(107) = ((- ALGBRC(82)+ALGBRC(58)) - ALGBRC(74)) - - ALGBRC(104)
      ALGBRC(15) =  CONSTS(57)*CONSTS(58)*log( CONSTS(38)*STATES(9))
      ALGBRC(45) = ALGBRC(15)
      ALGBRC(105) =  CONSTS(57)*CONSTS(58)*log( CONSTS(44)*STATES(15))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(108) =  CONSTS(27)*(EXP(ALGBRC(45)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(105)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(39) =  CONSTS(57)*CONSTS(58)*log( CONSTS(38)*STATES(9))+ CONSTS(51)*CONSTS(59)*ALGBRC(1)
      ALGBRC(87) =  CONSTS(57)*CONSTS(58)*log( CONSTS(39)*STATES(10))+ CONSTS(56)*CONSTS(59)*ALGBRC(1)
      ALGBRC(88) =  CONSTS(21)*(EXP(ALGBRC(39)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(87)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(110) = ((ALGBRC(72)+ALGBRC(60)) - ALGBRC(108)) - ALGBRC(88)
      ALGBRC(17) =  CONSTS(57)*CONSTS(58)*log( CONSTS(39)*STATES(10))
      ALGBRC(46) = ALGBRC(17)
      ALGBRC(109) =  CONSTS(57)*CONSTS(58)*log( CONSTS(45)*STATES(16))+ CONSTS(52)*ALGBRC(8)
      ALGBRC(111) =  CONSTS(29)*(EXP(ALGBRC(46)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(109)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(113) = ((ALGBRC(80)+ALGBRC(64)) - ALGBRC(111))+ALGBRC(88)
      ALGBRC(16) =  CONSTS(57)*CONSTS(58)*log( CONSTS(44)*STATES(15))
      ALGBRC(19) = ALGBRC(8)+ CONSTS(46)*CONSTS(59)*ALGBRC(1)+ALGBRC(16)
      ALGBRC(7) =  CONSTS(57)*CONSTS(58)*log( CONSTS(31)*STATES(1))
      ALGBRC(51) = ALGBRC(7)+ALGBRC(16)
      ALGBRC(47) =  CONSTS(46)*CONSTS(59)*ALGBRC(1)
      ALGBRC(112) = TERNRY(ALGBRC(47).EQ.0.00000,  CONSTS(2)*(EXP(ALGBRC(19)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(51)/( CONSTS(57)*CONSTS(58)))),  ((( CONSTS(2)*ALGBRC(47))/( CONSTS(57)*CONSTS(58)))/(EXP(ALGBRC(47)/( CONSTS(57)*CONSTS(58))) - 1.00000))*(EXP(ALGBRC(19)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(51)/( CONSTS(57)*CONSTS(58)))))
      ALGBRC(18) =  CONSTS(57)*CONSTS(58)*log( CONSTS(45)*STATES(16))
      ALGBRC(20) = ALGBRC(8)+ CONSTS(46)*CONSTS(59)*ALGBRC(1)+ALGBRC(18)
      ALGBRC(52) = ALGBRC(7)+ALGBRC(18)
      ALGBRC(48) =  CONSTS(46)*CONSTS(59)*ALGBRC(1)
      ALGBRC(114) = TERNRY(ALGBRC(48).EQ.0.00000,  CONSTS(3)*(EXP(ALGBRC(20)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(52)/( CONSTS(57)*CONSTS(58)))),  ((( CONSTS(3)*ALGBRC(48))/( CONSTS(57)*CONSTS(58)))/(EXP(ALGBRC(48)/( CONSTS(57)*CONSTS(58))) - 1.00000))*(EXP(ALGBRC(20)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(52)/( CONSTS(57)*CONSTS(58)))))
      ALGBRC(116) = ALGBRC(114)+ALGBRC(112)
      ALGBRC(117) = (( (CONSTS(52))*ALGBRC(104)+ CONSTS(52)*ALGBRC(108)+ CONSTS(52)*ALGBRC(100)+ CONSTS(52)*ALGBRC(96)+ CONSTS(52)*ALGBRC(92)+ CONSTS(52)*ALGBRC(111)) - ALGBRC(112)) - ALGBRC(114)
      ALGBRC(10) =  CONSTS(57)*CONSTS(58)*log( CONSTS(32)*STATES(4))
      ALGBRC(21) = ALGBRC(10)+ CONSTS(47)*CONSTS(59)*ALGBRC(1)+ALGBRC(16)
      ALGBRC(9) =  CONSTS(57)*CONSTS(58)*log( CONSTS(33)*STATES(3))
      ALGBRC(53) = ALGBRC(9)+ALGBRC(16)
      ALGBRC(49) =  CONSTS(47)*CONSTS(59)*ALGBRC(1)
      ALGBRC(115) = TERNRY(ALGBRC(49).EQ.0.00000,  CONSTS(4)*(EXP(ALGBRC(21)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(53)/( CONSTS(57)*CONSTS(58)))),  ((( CONSTS(4)*ALGBRC(49))/( CONSTS(57)*CONSTS(58)))/(EXP(ALGBRC(49)/( CONSTS(57)*CONSTS(58))) - 1.00000))*(EXP(ALGBRC(21)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(53)/( CONSTS(57)*CONSTS(58)))))
      ALGBRC(40) =  CONSTS(57)*CONSTS(58)*log( CONSTS(44)*STATES(15))+ CONSTS(51)*CONSTS(59)*ALGBRC(1)
      ALGBRC(89) =  CONSTS(57)*CONSTS(58)*log( CONSTS(45)*STATES(16))+ CONSTS(56)*CONSTS(59)*ALGBRC(1)
      ALGBRC(90) =  CONSTS(23)*(EXP(ALGBRC(40)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(89)/( CONSTS(57)*CONSTS(58))))
      ALGBRC(119) = ((ALGBRC(74)+ALGBRC(62)) - ((((- ALGBRC(108)+ALGBRC(112)) - ALGBRC(112))+ALGBRC(115)) - ALGBRC(115))) - ALGBRC(90)
      ALGBRC(22) = ALGBRC(10)+ CONSTS(47)*CONSTS(59)*ALGBRC(1)+ALGBRC(18)
      ALGBRC(54) = ALGBRC(9)+ALGBRC(18)
      ALGBRC(50) =  CONSTS(47)*CONSTS(59)*ALGBRC(1)
      ALGBRC(118) = TERNRY(ALGBRC(50).EQ.0.00000,  CONSTS(5)*(EXP(ALGBRC(22)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(54)/( CONSTS(57)*CONSTS(58)))),  ((( CONSTS(5)*ALGBRC(50))/( CONSTS(57)*CONSTS(58)))/(EXP(ALGBRC(50)/( CONSTS(57)*CONSTS(58))) - 1.00000))*(EXP(ALGBRC(22)/( CONSTS(57)*CONSTS(58))) - EXP(ALGBRC(54)/( CONSTS(57)*CONSTS(58)))))
      ALGBRC(121) = ALGBRC(118)+ALGBRC(115)
      ALGBRC(122) = - ALGBRC(115) - ALGBRC(118)
      ALGBRC(120) =  CONSTS(59)*((((((((((((((((((((((((((((((((((((( CONSTS(55)*ALGBRC(78) -  CONSTS(50)*ALGBRC(78))+ CONSTS(56)*ALGBRC(86)) -  CONSTS(51)*ALGBRC(86))+ CONSTS(54)*ALGBRC(74)) -  CONSTS(49)*ALGBRC(74))+ CONSTS(55)*ALGBRC(82)) -  CONSTS(50)*ALGBRC(82))+ CONSTS(56)*ALGBRC(90)) -  CONSTS(51)*ALGBRC(90))+ CONSTS(54)*ALGBRC(70)) -  CONSTS(49)*ALGBRC(70))+ CONSTS(53)*ALGBRC(66)) -  CONSTS(48)*ALGBRC(66))+ CONSTS(53)*ALGBRC(62)) -  CONSTS(48)*ALGBRC(62))+ CONSTS(53)*ALGBRC(58)) -  CONSTS(48)*ALGBRC(58))+ CONSTS(56)*ALGBRC(88)) -  CONSTS(51)*ALGBRC(88))+ CONSTS(55)*ALGBRC(80)) -  CONSTS(50)*ALGBRC(80))+ CONSTS(54)*ALGBRC(72)) -  CONSTS(49)*ALGBRC(72))+ CONSTS(56)*ALGBRC(84)) -  CONSTS(51)*ALGBRC(84))+ CONSTS(55)*ALGBRC(76)) -  CONSTS(50)*ALGBRC(76))+ CONSTS(54)*ALGBRC(68)) -  CONSTS(49)*ALGBRC(68))+ CONSTS(53)*ALGBRC(64)) -  CONSTS(48)*ALGBRC(64))+ CONSTS(53)*ALGBRC(60)) -  CONSTS(48)*ALGBRC(60))+ CONSTS(53)*ALGBRC(56)) -  CONSTS(48)*ALGBRC(56)) - ( CONSTS(46)*ALGBRC(112)+ CONSTS(46)*ALGBRC(114))) - ( CONSTS(47)*ALGBRC(115)+ CONSTS(47)*ALGBRC(118)))
      ALGBRC(123) = ((ALGBRC(82)+ALGBRC(66)) - ((((- ALGBRC(111)+ALGBRC(114)) - ALGBRC(114))+ALGBRC(118)) - ALGBRC(118)))+ALGBRC(90)
      ALGBRC(2) =  CONSTS(57)*CONSTS(58)*log( CONSTS(40)*STATES(11))
      ALGBRC(3) =  CONSTS(57)*CONSTS(58)*log( CONSTS(41)*STATES(12))
      ALGBRC(4) =  CONSTS(57)*CONSTS(58)*log( CONSTS(42)*STATES(13))
      ALGBRC(5) =  CONSTS(57)*CONSTS(58)*log( CONSTS(43)*STATES(14))
      ALGBRC(6) = STATES(1)+STATES(2)+STATES(11)+STATES(12)+STATES(13)+STATES(14)+STATES(15)+STATES(16)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END