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 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