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