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 124 entries in the algebraic variable array.
C There are a total of 19 entries in each of the rate and state variable arrays.
C There are a total of 58 entries in the constant variable array.
C
C
C VOI is t in component environment (second).
C CONSTS(1) is R in component environment (J_per_K_per_mol).
C CONSTS(2) is T in component environment (kelvin).
C CONSTS(3) is F in component environment (C_per_mol).
C CONSTS(4) is C_m in component environment (fF).
C STATES(1) is q_Na_o in component environment (fmol).
C STATES(2) is q_Na_i in component environment (fmol).
C STATES(3) is q_S000_Na in component environment (fmol).
C STATES(4) is q_S010_Na in component environment (fmol).
C STATES(5) is q_S100_Na in component environment (fmol).
C STATES(6) is q_S110_Na in component environment (fmol).
C STATES(7) is q_S200_Na in component environment (fmol).
C STATES(8) is q_S210_Na in component environment (fmol).
C STATES(9) is q_S300_Na in component environment (fmol).
C STATES(10) is q_S310_Na in component environment (fmol).
C STATES(11) is q_S001_Na in component environment (fmol).
C STATES(12) is q_S011_Na in component environment (fmol).
C STATES(13) is q_S101_Na in component environment (fmol).
C STATES(14) is q_S111_Na in component environment (fmol).
C STATES(15) is q_S201_Na in component environment (fmol).
C STATES(16) is q_S211_Na in component environment (fmol).
C STATES(17) is q_S301_Na in component environment (fmol).
C STATES(18) is q_S311_Na in component environment (fmol).
C ALGBRC(121) is v_Na in component fast_Na (fmol_per_sec).
C STATES(19) is q_mem in component environment (fC).
C ALGBRC(123) is I_mem_Na in component fast_Na (fA).
C CONSTS(5) is K_000_Na in component Na_parameters (per_fmol).
C CONSTS(6) is K_001_Na in component Na_parameters (per_fmol).
C CONSTS(7) is K_010_Na in component Na_parameters (per_fmol).
C CONSTS(8) is K_011_Na in component Na_parameters (per_fmol).
C CONSTS(9) is K_100_Na in component Na_parameters (per_fmol).
C CONSTS(10) is K_101_Na in component Na_parameters (per_fmol).
C CONSTS(11) is K_110_Na in component Na_parameters (per_fmol).
C CONSTS(12) is K_111_Na in component Na_parameters (per_fmol).
C CONSTS(13) is K_200_Na in component Na_parameters (per_fmol).
C CONSTS(14) is K_201_Na in component Na_parameters (per_fmol).
C CONSTS(15) is K_210_Na in component Na_parameters (per_fmol).
C CONSTS(16) is K_211_Na in component Na_parameters (per_fmol).
C CONSTS(17) is K_300_Na in component Na_parameters (per_fmol).
C CONSTS(18) is K_301_Na in component Na_parameters (per_fmol).
C CONSTS(19) is K_310_Na in component Na_parameters (per_fmol).
C CONSTS(20) is K_311_Na in component Na_parameters (per_fmol).
C CONSTS(21) is K_Na_o in component Na_parameters (per_fmol).
C CONSTS(22) is K_Na_i in component Na_parameters (per_fmol).
C CONSTS(23) is kappa_Na in component Na_parameters (fmol_per_sec).
C CONSTS(24) is kappa_h000 in component Na_parameters (fmol_per_sec).
C CONSTS(25) is kappa_h001 in component Na_parameters (fmol_per_sec).
C CONSTS(26) is kappa_h100 in component Na_parameters (fmol_per_sec).
C CONSTS(27) is kappa_h101 in component Na_parameters (fmol_per_sec).
C CONSTS(28) is kappa_h200 in component Na_parameters (fmol_per_sec).
C CONSTS(29) is kappa_h201 in component Na_parameters (fmol_per_sec).
C CONSTS(30) is kappa_h300 in component Na_parameters (fmol_per_sec).
C CONSTS(31) is kappa_h301 in component Na_parameters (fmol_per_sec).
C CONSTS(32) is kappa_j000 in component Na_parameters (fmol_per_sec).
C CONSTS(33) is kappa_j010 in component Na_parameters (fmol_per_sec).
C CONSTS(34) is kappa_j100 in component Na_parameters (fmol_per_sec).
C CONSTS(35) is kappa_j110 in component Na_parameters (fmol_per_sec).
C CONSTS(36) is kappa_j200 in component Na_parameters (fmol_per_sec).
C CONSTS(37) is kappa_j210 in component Na_parameters (fmol_per_sec).
C CONSTS(38) is kappa_j300 in component Na_parameters (fmol_per_sec).
C CONSTS(39) is kappa_j310 in component Na_parameters (fmol_per_sec).
C CONSTS(40) is kappa_m000 in component Na_parameters (fmol_per_sec).
C CONSTS(41) is kappa_m001 in component Na_parameters (fmol_per_sec).
C CONSTS(42) is kappa_m010 in component Na_parameters (fmol_per_sec).
C CONSTS(43) is kappa_m011 in component Na_parameters (fmol_per_sec).
C CONSTS(44) is kappa_m100 in component Na_parameters (fmol_per_sec).
C CONSTS(45) is kappa_m101 in component Na_parameters (fmol_per_sec).
C CONSTS(46) is kappa_m110 in component Na_parameters (fmol_per_sec).
C CONSTS(47) is kappa_m111 in component Na_parameters (fmol_per_sec).
C CONSTS(48) is kappa_m200 in component Na_parameters (fmol_per_sec).
C CONSTS(49) is kappa_m201 in component Na_parameters (fmol_per_sec).
C CONSTS(50) is kappa_m210 in component Na_parameters (fmol_per_sec).
C CONSTS(51) is kappa_m211 in component Na_parameters (fmol_per_sec).
C CONSTS(52) is zNa in component Na_parameters (dimensionless).
C CONSTS(53) is z_fh in component Na_parameters (dimensionless).
C CONSTS(54) is z_fj in component Na_parameters (dimensionless).
C CONSTS(55) is z_fm in component Na_parameters (dimensionless).
C CONSTS(56) is z_rh in component Na_parameters (dimensionless).
C CONSTS(57) is z_rj in component Na_parameters (dimensionless).
C CONSTS(58) is z_rm in component Na_parameters (dimensionless).
C ALGBRC(17) is mu_Na_o in component fast_Na (J_per_mol).
C ALGBRC(18) is mu_Na_i in component fast_Na (J_per_mol).
C ALGBRC(2) is mu_S000_Na in component fast_Na (J_per_mol).
C ALGBRC(3) is mu_S010_Na in component fast_Na (J_per_mol).
C ALGBRC(4) is mu_S100_Na in component fast_Na (J_per_mol).
C ALGBRC(5) is mu_S110_Na in component fast_Na (J_per_mol).
C ALGBRC(6) is mu_S200_Na in component fast_Na (J_per_mol).
C ALGBRC(7) is mu_S210_Na in component fast_Na (J_per_mol).
C ALGBRC(8) is mu_S300_Na in component fast_Na (J_per_mol).
C ALGBRC(9) is mu_S310_Na in component fast_Na (J_per_mol).
C ALGBRC(10) is mu_S001_Na in component fast_Na (J_per_mol).
C ALGBRC(11) is mu_S011_Na in component fast_Na (J_per_mol).
C ALGBRC(12) is mu_S101_Na in component fast_Na (J_per_mol).
C ALGBRC(13) is mu_S111_Na in component fast_Na (J_per_mol).
C ALGBRC(14) is mu_S201_Na in component fast_Na (J_per_mol).
C ALGBRC(15) is mu_S211_Na in component fast_Na (J_per_mol).
C ALGBRC(16) is mu_S301_Na in component fast_Na (J_per_mol).
C ALGBRC(19) is mu_S311_Na in component fast_Na (J_per_mol).
C ALGBRC(49) is Am_Na in component fast_Na (J_per_mol).
C ALGBRC(20) is Af_Na in component fast_Na (J_per_mol).
C ALGBRC(50) is Ar_Na in component fast_Na (J_per_mol).
C ALGBRC(21) is Af_h000 in component fast_Na (J_per_mol).
C ALGBRC(51) is Ar_h000 in component fast_Na (J_per_mol).
C ALGBRC(52) is v_h000 in component fast_Na (fmol_per_sec).
C ALGBRC(22) is Af_h001 in component fast_Na (J_per_mol).
C ALGBRC(53) is Ar_h001 in component fast_Na (J_per_mol).
C ALGBRC(54) is v_h001 in component fast_Na (fmol_per_sec).
C ALGBRC(23) is Af_h100 in component fast_Na (J_per_mol).
C ALGBRC(55) is Ar_h100 in component fast_Na (J_per_mol).
C ALGBRC(56) is v_h100 in component fast_Na (fmol_per_sec).
C ALGBRC(24) is Af_h101 in component fast_Na (J_per_mol).
C ALGBRC(57) is Ar_h101 in component fast_Na (J_per_mol).
C ALGBRC(58) is v_h101 in component fast_Na (fmol_per_sec).
C ALGBRC(25) is Af_h200 in component fast_Na (J_per_mol).
C ALGBRC(59) is Ar_h200 in component fast_Na (J_per_mol).
C ALGBRC(60) is v_h200 in component fast_Na (fmol_per_sec).
C ALGBRC(26) is Af_h201 in component fast_Na (J_per_mol).
C ALGBRC(61) is Ar_h201 in component fast_Na (J_per_mol).
C ALGBRC(62) is v_h201 in component fast_Na (fmol_per_sec).
C ALGBRC(27) is Af_h300 in component fast_Na (J_per_mol).
C ALGBRC(63) is Ar_h300 in component fast_Na (J_per_mol).
C ALGBRC(64) is v_h300 in component fast_Na (fmol_per_sec).
C ALGBRC(28) is Af_h301 in component fast_Na (J_per_mol).
C ALGBRC(65) is Ar_h301 in component fast_Na (J_per_mol).
C ALGBRC(66) is v_h301 in component fast_Na (fmol_per_sec).
C ALGBRC(29) is Af_j000 in component fast_Na (J_per_mol).
C ALGBRC(67) is Ar_j000 in component fast_Na (J_per_mol).
C ALGBRC(68) is v_j000 in component fast_Na (fmol_per_sec).
C ALGBRC(30) is Af_j010 in component fast_Na (J_per_mol).
C ALGBRC(69) is Ar_j010 in component fast_Na (J_per_mol).
C ALGBRC(70) is v_j010 in component fast_Na (fmol_per_sec).
C ALGBRC(31) is Af_j100 in component fast_Na (J_per_mol).
C ALGBRC(71) is Ar_j100 in component fast_Na (J_per_mol).
C ALGBRC(72) is v_j100 in component fast_Na (fmol_per_sec).
C ALGBRC(32) is Af_j110 in component fast_Na (J_per_mol).
C ALGBRC(73) is Ar_j110 in component fast_Na (J_per_mol).
C ALGBRC(74) is v_j110 in component fast_Na (fmol_per_sec).
C ALGBRC(33) is Af_j200 in component fast_Na (J_per_mol).
C ALGBRC(75) is Ar_j200 in component fast_Na (J_per_mol).
C ALGBRC(76) is v_j200 in component fast_Na (fmol_per_sec).
C ALGBRC(34) is Af_j210 in component fast_Na (J_per_mol).
C ALGBRC(77) is Ar_j210 in component fast_Na (J_per_mol).
C ALGBRC(78) is v_j210 in component fast_Na (fmol_per_sec).
C ALGBRC(35) is Af_j300 in component fast_Na (J_per_mol).
C ALGBRC(79) is Ar_j300 in component fast_Na (J_per_mol).
C ALGBRC(80) is v_j300 in component fast_Na (fmol_per_sec).
C ALGBRC(36) is Af_j310 in component fast_Na (J_per_mol).
C ALGBRC(81) is Ar_j310 in component fast_Na (J_per_mol).
C ALGBRC(82) is v_j310 in component fast_Na (fmol_per_sec).
C ALGBRC(37) is Af_m000 in component fast_Na (J_per_mol).
C ALGBRC(83) is Ar_m000 in component fast_Na (J_per_mol).
C ALGBRC(84) is v_m000 in component fast_Na (fmol_per_sec).
C ALGBRC(38) is Af_m001 in component fast_Na (J_per_mol).
C ALGBRC(85) is Ar_m001 in component fast_Na (J_per_mol).
C ALGBRC(87) is v_m001 in component fast_Na (fmol_per_sec).
C ALGBRC(39) is Af_m010 in component fast_Na (J_per_mol).
C ALGBRC(88) is Ar_m010 in component fast_Na (J_per_mol).
C ALGBRC(90) is v_m010 in component fast_Na (fmol_per_sec).
C ALGBRC(40) is Af_m011 in component fast_Na (J_per_mol).
C ALGBRC(91) is Ar_m011 in component fast_Na (J_per_mol).
C ALGBRC(93) is v_m011 in component fast_Na (fmol_per_sec).
C ALGBRC(41) is Af_m100 in component fast_Na (J_per_mol).
C ALGBRC(94) is Ar_m100 in component fast_Na (J_per_mol).
C ALGBRC(96) is v_m100 in component fast_Na (fmol_per_sec).
C ALGBRC(42) is Af_m101 in component fast_Na (J_per_mol).
C ALGBRC(97) is Ar_m101 in component fast_Na (J_per_mol).
C ALGBRC(99) is v_m101 in component fast_Na (fmol_per_sec).
C ALGBRC(43) is Af_m110 in component fast_Na (J_per_mol).
C ALGBRC(100) is Ar_m110 in component fast_Na (J_per_mol).
C ALGBRC(102) is v_m110 in component fast_Na (fmol_per_sec).
C ALGBRC(44) is Af_m111 in component fast_Na (J_per_mol).
C ALGBRC(103) is Ar_m111 in component fast_Na (J_per_mol).
C ALGBRC(105) is v_m111 in component fast_Na (fmol_per_sec).
C ALGBRC(45) is Af_m200 in component fast_Na (J_per_mol).
C ALGBRC(106) is Ar_m200 in component fast_Na (J_per_mol).
C ALGBRC(108) is v_m200 in component fast_Na (fmol_per_sec).
C ALGBRC(46) is Af_m201 in component fast_Na (J_per_mol).
C ALGBRC(109) is Ar_m201 in component fast_Na (J_per_mol).
C ALGBRC(112) is v_m201 in component fast_Na (fmol_per_sec).
C ALGBRC(47) is Af_m210 in component fast_Na (J_per_mol).
C ALGBRC(113) is Ar_m210 in component fast_Na (J_per_mol).
C ALGBRC(116) is v_m210 in component fast_Na (fmol_per_sec).
C ALGBRC(48) is Af_m211 in component fast_Na (J_per_mol).
C ALGBRC(117) is Ar_m211 in component fast_Na (J_per_mol).
C ALGBRC(120) is v_m211 in component fast_Na (fmol_per_sec).
C ALGBRC(86) is v_S000_Na in component fast_Na (fmol_per_sec).
C ALGBRC(92) is v_S010_Na in component fast_Na (fmol_per_sec).
C ALGBRC(98) is v_S100_Na in component fast_Na (fmol_per_sec).
C ALGBRC(104) is v_S110_Na in component fast_Na (fmol_per_sec).
C ALGBRC(110) is v_S200_Na in component fast_Na (fmol_per_sec).
C ALGBRC(118) is v_S210_Na in component fast_Na (fmol_per_sec).
C ALGBRC(111) is v_S300_Na in component fast_Na (fmol_per_sec).
C ALGBRC(119) is v_S310_Na in component fast_Na (fmol_per_sec).
C ALGBRC(89) is v_S001_Na in component fast_Na (fmol_per_sec).
C ALGBRC(95) is v_S011_Na in component fast_Na (fmol_per_sec).
C ALGBRC(101) is v_S101_Na in component fast_Na (fmol_per_sec).
C ALGBRC(107) is v_S111_Na in component fast_Na (fmol_per_sec).
C ALGBRC(114) is v_S201_Na in component fast_Na (fmol_per_sec).
C ALGBRC(122) is v_S211_Na in component fast_Na (fmol_per_sec).
C ALGBRC(115) is v_S301_Na in component fast_Na (fmol_per_sec).
C ALGBRC(124) is v_S311_Na in component fast_Na (fmol_per_sec).
C ALGBRC(1) is V_mem in component fast_Na (volt).
C RATES(1) is d/dt q_Na_o in component environment (fmol).
C RATES(2) is d/dt q_Na_i in component environment (fmol).
C RATES(19) is d/dt q_mem in component environment (fC).
C RATES(3) is d/dt q_S000_Na in component environment (fmol).
C RATES(4) is d/dt q_S010_Na in component environment (fmol).
C RATES(5) is d/dt q_S100_Na in component environment (fmol).
C RATES(6) is d/dt q_S110_Na in component environment (fmol).
C RATES(7) is d/dt q_S200_Na in component environment (fmol).
C RATES(8) is d/dt q_S210_Na in component environment (fmol).
C RATES(9) is d/dt q_S300_Na in component environment (fmol).
C RATES(10) is d/dt q_S310_Na in component environment (fmol).
C RATES(11) is d/dt q_S001_Na in component environment (fmol).
C RATES(12) is d/dt q_S011_Na in component environment (fmol).
C RATES(13) is d/dt q_S101_Na in component environment (fmol).
C RATES(14) is d/dt q_S111_Na in component environment (fmol).
C RATES(15) is d/dt q_S201_Na in component environment (fmol).
C RATES(16) is d/dt q_S211_Na in component environment (fmol).
C RATES(17) is d/dt q_S301_Na in component environment (fmol).
C RATES(18) is d/dt q_S311_Na in component environment (fmol).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 8.314
CONSTS(2) = 310
CONSTS(3) = 96485
CONSTS(4) = 153400
STATES(1) = 725.48
STATES(2) = 380
STATES(3) = 1.72572E-05
STATES(4) = 1.91747E-06
STATES(5) = 5.75256E-06
STATES(6) = 6.39176E-07
STATES(7) = 6.39176E-07
STATES(8) = 7.10196E-08
STATES(9) = 2.36737E-08
STATES(10) = 2.63039E-09
STATES(11) = 1.91747E-06
STATES(12) = 2.13059E-07
STATES(13) = 6.39176E-07
STATES(14) = 7.10196E-08
STATES(15) = 7.10196E-08
STATES(16) = 7.89103E-09
STATES(17) = 2.63039E-09
STATES(18) = 2.92267E-10
STATES(19) = -13039
CONSTS(5) = 0.0073376
CONSTS(6) = 23718.0222
CONSTS(7) = 17499.8767
CONSTS(8) = 56566767727.4036
CONSTS(9) = 1.5633e-05
CONSTS(10) = 50.5317
CONSTS(11) = 37.2839
CONSTS(12) = 120516694.0958
CONSTS(13) = 9.9918e-08
CONSTS(14) = 0.32298
CONSTS(15) = 0.2383
CONSTS(16) = 770290.0911
CONSTS(17) = 1.9159e-09
CONSTS(18) = 0.006193
CONSTS(19) = 0.0045694
CONSTS(20) = 14770.0739
CONSTS(21) = 0.017678
CONSTS(22) = 0.0024108
CONSTS(23) = 32.858
CONSTS(24) = 0.045709
CONSTS(25) = 1.4141e-08
CONSTS(26) = 21.4544
CONSTS(27) = 6.6373e-06
CONSTS(28) = 3356.675
CONSTS(29) = 0.0010384
CONSTS(30) = 175057.5863
CONSTS(31) = 0.054157
CONSTS(32) = 0.017824
CONSTS(33) = 7.4736e-09
CONSTS(34) = 8.3662
CONSTS(35) = 3.5079e-06
CONSTS(36) = 1308.9427
CONSTS(37) = 0.00054883
CONSTS(38) = 68264.0876
CONSTS(39) = 0.028623
CONSTS(40) = 5117405.6443
CONSTS(41) = 1.5832
CONSTS(42) = 2.1457
CONSTS(43) = 6.6381e-07
CONSTS(44) = 1601300127.6001
CONSTS(45) = 495.389
CONSTS(46) = 671.4132
CONSTS(47) = 0.00020771
CONSTS(48) = 125266701374.7363
CONSTS(49) = 38753.3515
CONSTS(50) = 52523.3901
CONSTS(51) = 0.016249
CONSTS(52) = 1
CONSTS(53) = -4.1892
CONSTS(54) = -4.0381
CONSTS(55) = 0.49541
CONSTS(56) = 1.2995
CONSTS(57) = 1.4281
CONSTS(58) = -2.4284
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = STATES(19)/CONSTS(4)
ALGBRC(21) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(51) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(52) = CONSTS(24)*(EXP(ALGBRC(21)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(51)/( CONSTS(1)*CONSTS(2))))
ALGBRC(29) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(67) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(68) = CONSTS(32)*(EXP(ALGBRC(29)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(67)/( CONSTS(1)*CONSTS(2))))
ALGBRC(37) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(83) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(84) = CONSTS(40)*(EXP(ALGBRC(37)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(83)/( CONSTS(1)*CONSTS(2))))
ALGBRC(86) = (- ALGBRC(84) - ALGBRC(52)) - ALGBRC(68)
RATES(3) = ALGBRC(86)
ALGBRC(22) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(53) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(54) = CONSTS(25)*(EXP(ALGBRC(22)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(53)/( CONSTS(1)*CONSTS(2))))
ALGBRC(38) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(85) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(87) = CONSTS(41)*(EXP(ALGBRC(38)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(85)/( CONSTS(1)*CONSTS(2))))
ALGBRC(89) = (- ALGBRC(87) - ALGBRC(54)) - - ALGBRC(68)
RATES(11) = ALGBRC(89)
ALGBRC(30) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(69) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(70) = CONSTS(33)*(EXP(ALGBRC(30)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(69)/( CONSTS(1)*CONSTS(2))))
ALGBRC(39) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(88) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(90) = CONSTS(42)*(EXP(ALGBRC(39)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(88)/( CONSTS(1)*CONSTS(2))))
ALGBRC(92) = (- ALGBRC(90)+ALGBRC(52)) - ALGBRC(70)
RATES(4) = ALGBRC(92)
ALGBRC(40) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(91) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(93) = CONSTS(43)*(EXP(ALGBRC(40)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(91)/( CONSTS(1)*CONSTS(2))))
ALGBRC(95) = (- ALGBRC(93)+ALGBRC(54)) - - ALGBRC(70)
RATES(12) = ALGBRC(95)
ALGBRC(23) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(55) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(56) = CONSTS(26)*(EXP(ALGBRC(23)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(55)/( CONSTS(1)*CONSTS(2))))
ALGBRC(31) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(71) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(72) = CONSTS(34)*(EXP(ALGBRC(31)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(71)/( CONSTS(1)*CONSTS(2))))
ALGBRC(41) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(94) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(96) = CONSTS(44)*(EXP(ALGBRC(41)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(94)/( CONSTS(1)*CONSTS(2))))
ALGBRC(98) = ((ALGBRC(84) - ALGBRC(96)) - ALGBRC(56)) - ALGBRC(72)
RATES(5) = ALGBRC(98)
ALGBRC(24) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(57) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(58) = CONSTS(27)*(EXP(ALGBRC(24)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(57)/( CONSTS(1)*CONSTS(2))))
ALGBRC(42) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(97) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(99) = CONSTS(45)*(EXP(ALGBRC(42)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(97)/( CONSTS(1)*CONSTS(2))))
ALGBRC(101) = ((ALGBRC(87) - ALGBRC(99)) - ALGBRC(58)) - - ALGBRC(72)
RATES(13) = ALGBRC(101)
ALGBRC(32) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(73) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(74) = CONSTS(35)*(EXP(ALGBRC(32)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(73)/( CONSTS(1)*CONSTS(2))))
ALGBRC(43) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(100) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(102) = CONSTS(46)*(EXP(ALGBRC(43)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(100)/( CONSTS(1)*CONSTS(2))))
ALGBRC(104) = ((ALGBRC(90) - ALGBRC(102))+ALGBRC(56)) - ALGBRC(74)
RATES(6) = ALGBRC(104)
ALGBRC(44) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(103) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(105) = CONSTS(47)*(EXP(ALGBRC(44)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(103)/( CONSTS(1)*CONSTS(2))))
ALGBRC(107) = ((ALGBRC(93) - ALGBRC(105))+ALGBRC(58)) - - ALGBRC(74)
RATES(14) = ALGBRC(107)
ALGBRC(25) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(59) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(60) = CONSTS(28)*(EXP(ALGBRC(25)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(59)/( CONSTS(1)*CONSTS(2))))
ALGBRC(33) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(75) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(76) = CONSTS(36)*(EXP(ALGBRC(33)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(75)/( CONSTS(1)*CONSTS(2))))
ALGBRC(45) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(106) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(108) = CONSTS(48)*(EXP(ALGBRC(45)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(106)/( CONSTS(1)*CONSTS(2))))
ALGBRC(110) = ((ALGBRC(96) - ALGBRC(108)) - ALGBRC(60)) - ALGBRC(76)
RATES(7) = ALGBRC(110)
ALGBRC(27) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(63) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(64) = CONSTS(30)*(EXP(ALGBRC(27)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(63)/( CONSTS(1)*CONSTS(2))))
ALGBRC(35) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(79) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(80) = CONSTS(38)*(EXP(ALGBRC(35)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(79)/( CONSTS(1)*CONSTS(2))))
ALGBRC(111) = (ALGBRC(108) - ALGBRC(64)) - ALGBRC(80)
RATES(9) = ALGBRC(111)
ALGBRC(26) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(61) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(62) = CONSTS(29)*(EXP(ALGBRC(26)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(61)/( CONSTS(1)*CONSTS(2))))
ALGBRC(46) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(109) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(112) = CONSTS(49)*(EXP(ALGBRC(46)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(109)/( CONSTS(1)*CONSTS(2))))
ALGBRC(114) = ((ALGBRC(99) - ALGBRC(112)) - ALGBRC(62)) - - ALGBRC(76)
RATES(15) = ALGBRC(114)
ALGBRC(28) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(65) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(66) = CONSTS(31)*(EXP(ALGBRC(28)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(65)/( CONSTS(1)*CONSTS(2))))
ALGBRC(115) = (ALGBRC(112) - ALGBRC(66)) - - ALGBRC(80)
RATES(17) = ALGBRC(115)
ALGBRC(34) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(77) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(78) = CONSTS(37)*(EXP(ALGBRC(34)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(77)/( CONSTS(1)*CONSTS(2))))
ALGBRC(47) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(113) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(116) = CONSTS(50)*(EXP(ALGBRC(47)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(113)/( CONSTS(1)*CONSTS(2))))
ALGBRC(118) = ((ALGBRC(102) - ALGBRC(116))+ALGBRC(60)) - ALGBRC(78)
RATES(8) = ALGBRC(118)
ALGBRC(36) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(81) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(82) = CONSTS(39)*(EXP(ALGBRC(36)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(81)/( CONSTS(1)*CONSTS(2))))
ALGBRC(119) = (ALGBRC(116)+ALGBRC(64)) - ALGBRC(82)
RATES(10) = ALGBRC(119)
ALGBRC(49) = CONSTS(52)*CONSTS(3)*ALGBRC(1)
ALGBRC(18) = CONSTS(1)*CONSTS(2)*log( CONSTS(22)*STATES(2))
ALGBRC(19) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))
ALGBRC(20) = ALGBRC(18)+ CONSTS(52)*CONSTS(3)*ALGBRC(1)+ALGBRC(19)
ALGBRC(17) = CONSTS(1)*CONSTS(2)*log( CONSTS(21)*STATES(1))
ALGBRC(50) = ALGBRC(17)+ALGBRC(19)
ALGBRC(121) = TERNRY(ALGBRC(49).EQ.0.00000, CONSTS(23)*(EXP(ALGBRC(20)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(50)/( CONSTS(1)*CONSTS(2)))), ((( CONSTS(23)*ALGBRC(49))/( CONSTS(1)*CONSTS(2)))/(EXP(ALGBRC(49)/( CONSTS(1)*CONSTS(2))) - 1.00000))*(EXP(ALGBRC(20)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(50)/( CONSTS(1)*CONSTS(2)))))
RATES(1) = ALGBRC(121)
RATES(2) = - ALGBRC(121)
ALGBRC(48) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(117) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(120) = CONSTS(51)*(EXP(ALGBRC(48)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(117)/( CONSTS(1)*CONSTS(2))))
ALGBRC(122) = ((ALGBRC(105) - ALGBRC(120))+ALGBRC(62)) - - ALGBRC(78)
RATES(16) = ALGBRC(122)
ALGBRC(123) = CONSTS(3)*(((((((((((((((((((((((((((((((((((((((((((((((((((((((( - CONSTS(52)*ALGBRC(121)+ CONSTS(58)*ALGBRC(84)) - CONSTS(55)*ALGBRC(84))+ CONSTS(58)*ALGBRC(96)) - CONSTS(55)*ALGBRC(96))+ CONSTS(58)*ALGBRC(90)) - CONSTS(55)*ALGBRC(90))+ CONSTS(58)*ALGBRC(102)) - CONSTS(55)*ALGBRC(102))+ CONSTS(58)*ALGBRC(116)) - CONSTS(55)*ALGBRC(116))+ CONSTS(56)*ALGBRC(52)) - CONSTS(53)*ALGBRC(52))+ CONSTS(56)*ALGBRC(60)) - CONSTS(53)*ALGBRC(60))+ CONSTS(57)*ALGBRC(72)) - CONSTS(54)*ALGBRC(72))+ CONSTS(57)*ALGBRC(76)) - CONSTS(54)*ALGBRC(76))+ CONSTS(57)*ALGBRC(80)) - CONSTS(54)*ALGBRC(80))+ CONSTS(58)*ALGBRC(108)) - CONSTS(55)*ALGBRC(108))+ CONSTS(56)*ALGBRC(56)) - CONSTS(53)*ALGBRC(56))+ CONSTS(56)*ALGBRC(64)) - CONSTS(53)*ALGBRC(64))+ CONSTS(57)*ALGBRC(68)) - CONSTS(54)*ALGBRC(68))+ CONSTS(58)*ALGBRC(87)) - CONSTS(55)*ALGBRC(87))+ CONSTS(58)*ALGBRC(99)) - CONSTS(55)*ALGBRC(99))+ CONSTS(58)*ALGBRC(93)) - CONSTS(55)*ALGBRC(93))+ CONSTS(58)*ALGBRC(105)) - CONSTS(55)*ALGBRC(105))+ CONSTS(58)*ALGBRC(120)) - CONSTS(55)*ALGBRC(120))+ CONSTS(56)*ALGBRC(54)) - CONSTS(53)*ALGBRC(54))+ CONSTS(56)*ALGBRC(62)) - CONSTS(53)*ALGBRC(62))+ CONSTS(57)*ALGBRC(70)) - CONSTS(54)*ALGBRC(70))+ CONSTS(57)*ALGBRC(74)) - CONSTS(54)*ALGBRC(74))+ CONSTS(57)*ALGBRC(78)) - CONSTS(54)*ALGBRC(78))+ CONSTS(57)*ALGBRC(82)) - CONSTS(54)*ALGBRC(82))+ CONSTS(58)*ALGBRC(112)) - CONSTS(55)*ALGBRC(112))+ CONSTS(56)*ALGBRC(58)) - CONSTS(53)*ALGBRC(58))+ CONSTS(56)*ALGBRC(66)) - CONSTS(53)*ALGBRC(66))
RATES(19) = ALGBRC(123)
ALGBRC(124) = (ALGBRC(120)+ALGBRC(66)) - ((- ALGBRC(82)+ALGBRC(121)) - ALGBRC(121))
RATES(18) = ALGBRC(124)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = STATES(19)/CONSTS(4)
ALGBRC(21) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(51) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(52) = CONSTS(24)*(EXP(ALGBRC(21)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(51)/( CONSTS(1)*CONSTS(2))))
ALGBRC(29) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(67) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(68) = CONSTS(32)*(EXP(ALGBRC(29)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(67)/( CONSTS(1)*CONSTS(2))))
ALGBRC(37) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(83) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(84) = CONSTS(40)*(EXP(ALGBRC(37)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(83)/( CONSTS(1)*CONSTS(2))))
ALGBRC(86) = (- ALGBRC(84) - ALGBRC(52)) - ALGBRC(68)
ALGBRC(22) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(53) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(54) = CONSTS(25)*(EXP(ALGBRC(22)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(53)/( CONSTS(1)*CONSTS(2))))
ALGBRC(38) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(85) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(87) = CONSTS(41)*(EXP(ALGBRC(38)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(85)/( CONSTS(1)*CONSTS(2))))
ALGBRC(89) = (- ALGBRC(87) - ALGBRC(54)) - - ALGBRC(68)
ALGBRC(30) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(69) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(70) = CONSTS(33)*(EXP(ALGBRC(30)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(69)/( CONSTS(1)*CONSTS(2))))
ALGBRC(39) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(88) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(90) = CONSTS(42)*(EXP(ALGBRC(39)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(88)/( CONSTS(1)*CONSTS(2))))
ALGBRC(92) = (- ALGBRC(90)+ALGBRC(52)) - ALGBRC(70)
ALGBRC(40) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(91) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(93) = CONSTS(43)*(EXP(ALGBRC(40)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(91)/( CONSTS(1)*CONSTS(2))))
ALGBRC(95) = (- ALGBRC(93)+ALGBRC(54)) - - ALGBRC(70)
ALGBRC(23) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(55) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(56) = CONSTS(26)*(EXP(ALGBRC(23)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(55)/( CONSTS(1)*CONSTS(2))))
ALGBRC(31) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(71) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(72) = CONSTS(34)*(EXP(ALGBRC(31)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(71)/( CONSTS(1)*CONSTS(2))))
ALGBRC(41) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(94) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(96) = CONSTS(44)*(EXP(ALGBRC(41)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(94)/( CONSTS(1)*CONSTS(2))))
ALGBRC(98) = ((ALGBRC(84) - ALGBRC(96)) - ALGBRC(56)) - ALGBRC(72)
ALGBRC(24) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(57) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(58) = CONSTS(27)*(EXP(ALGBRC(24)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(57)/( CONSTS(1)*CONSTS(2))))
ALGBRC(42) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(97) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(99) = CONSTS(45)*(EXP(ALGBRC(42)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(97)/( CONSTS(1)*CONSTS(2))))
ALGBRC(101) = ((ALGBRC(87) - ALGBRC(99)) - ALGBRC(58)) - - ALGBRC(72)
ALGBRC(32) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(73) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(74) = CONSTS(35)*(EXP(ALGBRC(32)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(73)/( CONSTS(1)*CONSTS(2))))
ALGBRC(43) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(100) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(102) = CONSTS(46)*(EXP(ALGBRC(43)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(100)/( CONSTS(1)*CONSTS(2))))
ALGBRC(104) = ((ALGBRC(90) - ALGBRC(102))+ALGBRC(56)) - ALGBRC(74)
ALGBRC(44) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(103) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(105) = CONSTS(47)*(EXP(ALGBRC(44)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(103)/( CONSTS(1)*CONSTS(2))))
ALGBRC(107) = ((ALGBRC(93) - ALGBRC(105))+ALGBRC(58)) - - ALGBRC(74)
ALGBRC(25) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(59) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(60) = CONSTS(28)*(EXP(ALGBRC(25)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(59)/( CONSTS(1)*CONSTS(2))))
ALGBRC(33) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(75) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(76) = CONSTS(36)*(EXP(ALGBRC(33)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(75)/( CONSTS(1)*CONSTS(2))))
ALGBRC(45) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(106) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(108) = CONSTS(48)*(EXP(ALGBRC(45)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(106)/( CONSTS(1)*CONSTS(2))))
ALGBRC(110) = ((ALGBRC(96) - ALGBRC(108)) - ALGBRC(60)) - ALGBRC(76)
ALGBRC(27) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(63) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(64) = CONSTS(30)*(EXP(ALGBRC(27)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(63)/( CONSTS(1)*CONSTS(2))))
ALGBRC(35) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(79) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(80) = CONSTS(38)*(EXP(ALGBRC(35)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(79)/( CONSTS(1)*CONSTS(2))))
ALGBRC(111) = (ALGBRC(108) - ALGBRC(64)) - ALGBRC(80)
ALGBRC(26) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(61) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(62) = CONSTS(29)*(EXP(ALGBRC(26)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(61)/( CONSTS(1)*CONSTS(2))))
ALGBRC(46) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(109) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(112) = CONSTS(49)*(EXP(ALGBRC(46)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(109)/( CONSTS(1)*CONSTS(2))))
ALGBRC(114) = ((ALGBRC(99) - ALGBRC(112)) - ALGBRC(62)) - - ALGBRC(76)
ALGBRC(28) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(53)*CONSTS(3)*ALGBRC(1)
ALGBRC(65) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(56)*CONSTS(3)*ALGBRC(1)
ALGBRC(66) = CONSTS(31)*(EXP(ALGBRC(28)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(65)/( CONSTS(1)*CONSTS(2))))
ALGBRC(115) = (ALGBRC(112) - ALGBRC(66)) - - ALGBRC(80)
ALGBRC(34) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(77) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(78) = CONSTS(37)*(EXP(ALGBRC(34)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(77)/( CONSTS(1)*CONSTS(2))))
ALGBRC(47) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(113) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(116) = CONSTS(50)*(EXP(ALGBRC(47)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(113)/( CONSTS(1)*CONSTS(2))))
ALGBRC(118) = ((ALGBRC(102) - ALGBRC(116))+ALGBRC(60)) - ALGBRC(78)
ALGBRC(36) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(54)*CONSTS(3)*ALGBRC(1)
ALGBRC(81) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(57)*CONSTS(3)*ALGBRC(1)
ALGBRC(82) = CONSTS(39)*(EXP(ALGBRC(36)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(81)/( CONSTS(1)*CONSTS(2))))
ALGBRC(119) = (ALGBRC(116)+ALGBRC(64)) - ALGBRC(82)
ALGBRC(49) = CONSTS(52)*CONSTS(3)*ALGBRC(1)
ALGBRC(18) = CONSTS(1)*CONSTS(2)*log( CONSTS(22)*STATES(2))
ALGBRC(19) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))
ALGBRC(20) = ALGBRC(18)+ CONSTS(52)*CONSTS(3)*ALGBRC(1)+ALGBRC(19)
ALGBRC(17) = CONSTS(1)*CONSTS(2)*log( CONSTS(21)*STATES(1))
ALGBRC(50) = ALGBRC(17)+ALGBRC(19)
ALGBRC(121) = TERNRY(ALGBRC(49).EQ.0.00000, CONSTS(23)*(EXP(ALGBRC(20)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(50)/( CONSTS(1)*CONSTS(2)))), ((( CONSTS(23)*ALGBRC(49))/( CONSTS(1)*CONSTS(2)))/(EXP(ALGBRC(49)/( CONSTS(1)*CONSTS(2))) - 1.00000))*(EXP(ALGBRC(20)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(50)/( CONSTS(1)*CONSTS(2)))))
ALGBRC(48) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(55)*CONSTS(3)*ALGBRC(1)
ALGBRC(117) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(58)*CONSTS(3)*ALGBRC(1)
ALGBRC(120) = CONSTS(51)*(EXP(ALGBRC(48)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(117)/( CONSTS(1)*CONSTS(2))))
ALGBRC(122) = ((ALGBRC(105) - ALGBRC(120))+ALGBRC(62)) - - ALGBRC(78)
ALGBRC(123) = CONSTS(3)*(((((((((((((((((((((((((((((((((((((((((((((((((((((((( - CONSTS(52)*ALGBRC(121)+ CONSTS(58)*ALGBRC(84)) - CONSTS(55)*ALGBRC(84))+ CONSTS(58)*ALGBRC(96)) - CONSTS(55)*ALGBRC(96))+ CONSTS(58)*ALGBRC(90)) - CONSTS(55)*ALGBRC(90))+ CONSTS(58)*ALGBRC(102)) - CONSTS(55)*ALGBRC(102))+ CONSTS(58)*ALGBRC(116)) - CONSTS(55)*ALGBRC(116))+ CONSTS(56)*ALGBRC(52)) - CONSTS(53)*ALGBRC(52))+ CONSTS(56)*ALGBRC(60)) - CONSTS(53)*ALGBRC(60))+ CONSTS(57)*ALGBRC(72)) - CONSTS(54)*ALGBRC(72))+ CONSTS(57)*ALGBRC(76)) - CONSTS(54)*ALGBRC(76))+ CONSTS(57)*ALGBRC(80)) - CONSTS(54)*ALGBRC(80))+ CONSTS(58)*ALGBRC(108)) - CONSTS(55)*ALGBRC(108))+ CONSTS(56)*ALGBRC(56)) - CONSTS(53)*ALGBRC(56))+ CONSTS(56)*ALGBRC(64)) - CONSTS(53)*ALGBRC(64))+ CONSTS(57)*ALGBRC(68)) - CONSTS(54)*ALGBRC(68))+ CONSTS(58)*ALGBRC(87)) - CONSTS(55)*ALGBRC(87))+ CONSTS(58)*ALGBRC(99)) - CONSTS(55)*ALGBRC(99))+ CONSTS(58)*ALGBRC(93)) - CONSTS(55)*ALGBRC(93))+ CONSTS(58)*ALGBRC(105)) - CONSTS(55)*ALGBRC(105))+ CONSTS(58)*ALGBRC(120)) - CONSTS(55)*ALGBRC(120))+ CONSTS(56)*ALGBRC(54)) - CONSTS(53)*ALGBRC(54))+ CONSTS(56)*ALGBRC(62)) - CONSTS(53)*ALGBRC(62))+ CONSTS(57)*ALGBRC(70)) - CONSTS(54)*ALGBRC(70))+ CONSTS(57)*ALGBRC(74)) - CONSTS(54)*ALGBRC(74))+ CONSTS(57)*ALGBRC(78)) - CONSTS(54)*ALGBRC(78))+ CONSTS(57)*ALGBRC(82)) - CONSTS(54)*ALGBRC(82))+ CONSTS(58)*ALGBRC(112)) - CONSTS(55)*ALGBRC(112))+ CONSTS(56)*ALGBRC(58)) - CONSTS(53)*ALGBRC(58))+ CONSTS(56)*ALGBRC(66)) - CONSTS(53)*ALGBRC(66))
ALGBRC(124) = (ALGBRC(120)+ALGBRC(66)) - ((- ALGBRC(82)+ALGBRC(121)) - ALGBRC(121))
ALGBRC(2) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))
ALGBRC(3) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))
ALGBRC(4) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))
ALGBRC(5) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))
ALGBRC(6) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))
ALGBRC(7) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))
ALGBRC(8) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))
ALGBRC(9) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))
ALGBRC(10) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))
ALGBRC(11) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))
ALGBRC(12) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))
ALGBRC(13) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))
ALGBRC(14) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))
ALGBRC(15) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))
ALGBRC(16) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END