C
C There are a total of 70 entries in the algebraic variable array.
C There are a total of 7 entries in each of the rate and state variable arrays.
C There are a total of 43 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C STATES(1) is V in component membrane (millivolt).
C ALGBRC(2) is IK in component membrane (femtoA).
C ALGBRC(5) is ICa in component membrane (femtoA).
C ALGBRC(3) is IKCa in component membrane (femtoA).
C ALGBRC(30) is IKATP in component KATP (femtoA).
C CONSTS(1) is Cm in component membrane (femtoF).
C CONSTS(2) is gK in component membrane (picoS).
C CONSTS(3) is VK in component KATP (millivolt).
C STATES(2) is n in component membrane (dimensionless).
C CONSTS(4) is gKCa in component membrane (picoS).
C CONSTS(5) is kd in component membrane (micromolar).
C STATES(3) is c in component calcium_handling (micromolar).
C CONSTS(6) is gCa in component membrane (picoS).
C ALGBRC(4) is minf in component membrane (dimensionless).
C CONSTS(7) is VCa in component membrane (millivolt).
C CONSTS(8) is taun in component membrane (millisecond).
C ALGBRC(1) is ninf in component membrane (dimensionless).
C STATES(4) is cer in component calcium_handling (micromolar).
C CONSTS(9) is fcyt in component calcium_handling (dimensionless).
C ALGBRC(12) is Jmem in component calcium_handling (flux).
C ALGBRC(10) is Jer in component calcium_handling (flux).
C CONSTS(10) is fer in component calcium_handling (dimensionless).
C CONSTS(11) is sigmaV in component calcium_handling (dimensionless).
C CONSTS(12) is pleak in component calcium_handling (first_order_rate_constant).
C CONSTS(13) is Kserca in component calcium_handling (first_order_rate_constant).
C CONSTS(14) is lambdaer in component calcium_handling (dimensionless).
C CONSTS(15) is epser in component calcium_handling (dimensionless).
C CONSTS(16) is alpha in component calcium_handling (micromolar_per_femtoA_millisecond).
C CONSTS(17) is kpmca in component calcium_handling (first_order_rate_constant).
C ALGBRC(6) is Jserca in component calcium_handling (flux).
C ALGBRC(8) is Jleak in component calcium_handling (flux).
C ALGBRC(9) is rgpdh in component glycolysis (flux).
C CONSTS(18) is Rgk in component glycolysis (per_second).
C CONSTS(19) is atot in component glycolysis (micromolar).
C CONSTS(20) is pfkbas in component glycolysis (dimensionless).
C ALGBRC(7) is f6p in component glycolysis (micromolar).
C CONSTS(21) is lambda in component glycolysis (dimensionless).
C ALGBRC(70) is pfk in component pfk (micromolar).
C STATES(5) is g6p in component glycolysis (micromolar).
C STATES(6) is fbp in component glycolysis (micromolar).
C CONSTS(22) is bottom1 in component pfk (dimensionless).
C CONSTS(23) is topa1 in component pfk (dimensionless).
C CONSTS(24) is k1 in component pfk (micromolar).
C CONSTS(25) is k2 in component pfk (micromolar).
C CONSTS(26) is k3 in component pfk (micromolar).
C CONSTS(27) is k4 in component pfk (micromolar).
C CONSTS(28) is cat in component pfk (dimensionless).
C ALGBRC(22) is atp in component nucleotides (micromolar).
C ALGBRC(23) is weight2 in component pfk (dimensionless).
C CONSTS(43) is topa2 in component pfk (dimensionless).
C ALGBRC(26) is bottom2 in component pfk (dimensionless).
C ALGBRC(13) is topa3 in component pfk (dimensionless).
C ALGBRC(11) is weight3 in component pfk (dimensionless).
C ALGBRC(29) is bottom3 in component pfk (dimensionless).
C CONSTS(29) is famp in component pfk (dimensionless).
C CONSTS(30) is fatp in component pfk (dimensionless).
C CONSTS(31) is ffbp in component pfk (dimensionless).
C CONSTS(32) is fbt in component pfk (dimensionless).
C CONSTS(33) is fmt in component pfk (dimensionless).
C ALGBRC(31) is weight4 in component pfk (dimensionless).
C ALGBRC(32) is topa4 in component pfk (dimensionless).
C ALGBRC(33) is bottom4 in component pfk (dimensionless).
C ALGBRC(14) is weight5 in component pfk (dimensionless).
C ALGBRC(34) is topa5 in component pfk (dimensionless).
C ALGBRC(35) is bottom5 in component pfk (dimensionless).
C ALGBRC(36) is weight6 in component pfk (dimensionless).
C ALGBRC(37) is topa6 in component pfk (dimensionless).
C ALGBRC(38) is bottom6 in component pfk (dimensionless).
C ALGBRC(15) is weight7 in component pfk (dimensionless).
C ALGBRC(39) is topa7 in component pfk (dimensionless).
C ALGBRC(40) is bottom7 in component pfk (dimensionless).
C ALGBRC(41) is weight8 in component pfk (dimensionless).
C ALGBRC(42) is topa8 in component pfk (dimensionless).
C ALGBRC(43) is bottom8 in component pfk (dimensionless).
C ALGBRC(47) is weight9 in component pfk (dimensionless).
C ALGBRC(44) is topa9 in component pfk (dimensionless).
C ALGBRC(48) is bottom9 in component pfk (dimensionless).
C ALGBRC(49) is weight10 in component pfk (dimensionless).
C ALGBRC(45) is topa10 in component pfk (dimensionless).
C ALGBRC(50) is bottom10 in component pfk (dimensionless).
C ALGBRC(51) is weight11 in component pfk (dimensionless).
C ALGBRC(52) is topa11 in component pfk (dimensionless).
C ALGBRC(53) is bottom11 in component pfk (dimensionless).
C ALGBRC(54) is weight12 in component pfk (dimensionless).
C ALGBRC(55) is topa12 in component pfk (dimensionless).
C ALGBRC(56) is bottom12 in component pfk (dimensionless).
C ALGBRC(57) is weight13 in component pfk (dimensionless).
C ALGBRC(58) is topa13 in component pfk (dimensionless).
C ALGBRC(59) is bottom13 in component pfk (dimensionless).
C ALGBRC(60) is weight14 in component pfk (dimensionless).
C ALGBRC(61) is topa14 in component pfk (dimensionless).
C ALGBRC(62) is bottom14 in component pfk (dimensionless).
C ALGBRC(63) is weight15 in component pfk (dimensionless).
C ALGBRC(64) is topa15 in component pfk (dimensionless).
C ALGBRC(65) is bottom15 in component pfk (dimensionless).
C ALGBRC(67) is weight16 in component pfk (dimensionless).
C ALGBRC(68) is topa16 in component pfk (dimensionless).
C ALGBRC(69) is bottom16 in component pfk (dimensionless).
C ALGBRC(66) is topb in component pfk (dimensionless).
C ALGBRC(46) is amp in component nucleotides (micromolar).
C ALGBRC(16) is mgadp in component KATP (micromolar).
C ALGBRC(17) is adp3m in component KATP (micromolar).
C ALGBRC(24) is atp4m in component KATP (micromolar).
C ALGBRC(18) is topo in component KATP (dimensionless).
C ALGBRC(27) is bottomo in component KATP (dimensionless).
C ALGBRC(28) is katpo in component KATP (dimensionless).
C CONSTS(34) is gkatpbar in component KATP (picoS).
C STATES(7) is adp in component nucleotides (micromolar).
C CONSTS(35) is kdd in component KATP (dimensionless).
C CONSTS(36) is ktd in component KATP (dimensionless).
C CONSTS(37) is ktt in component KATP (dimensionless).
C ALGBRC(20) is fback in component nucleotides (dimensionless).
C CONSTS(38) is taua in component nucleotides (dimensionless).
C CONSTS(39) is r1 in component nucleotides (micromolar).
C CONSTS(40) is r in component nucleotides (dimensionless).
C ALGBRC(19) is y in component nucleotides (dimensionless).
C CONSTS(41) is vg in component nucleotides (dimensionless).
C CONSTS(42) is kg in component nucleotides (flux).
C ALGBRC(21) is rad in component nucleotides (dimensionless).
C ALGBRC(25) is ratio in component nucleotides (dimensionless).
C RATES(1) is d/dt V in component membrane (millivolt).
C RATES(2) is d/dt n in component membrane (dimensionless).
C RATES(3) is d/dt c in component calcium_handling (micromolar).
C RATES(4) is d/dt cer in component calcium_handling (micromolar).
C RATES(6) is d/dt fbp in component glycolysis (micromolar).
C RATES(5) is d/dt g6p in component glycolysis (micromolar).
C RATES(7) is d/dt adp in component nucleotides (micromolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = -60
      CONSTS(1) = 5300
      CONSTS(2) = 2700
      CONSTS(3) = -75
      STATES(2) = 0
      CONSTS(4) = 600
      CONSTS(5) = 0.5
      STATES(3) = 0.25
      CONSTS(6) = 1000
      CONSTS(7) = 25
      CONSTS(8) = 20
      STATES(4) = 185
      CONSTS(9) = 0.01
      CONSTS(10) = 0.01
      CONSTS(11) = 31
      CONSTS(12) = 0.0002
      CONSTS(13) = 0.4
      CONSTS(14) = 1
      CONSTS(15) = 1
      CONSTS(16) = 0.00000450
      CONSTS(17) = 0.2
      CONSTS(18) = 0.2
      CONSTS(19) = 3000
      CONSTS(20) = 0.06
      CONSTS(21) = 0.005
      STATES(5) = 200
      STATES(6) = 40
      CONSTS(22) = 1
      CONSTS(23) = 0
      CONSTS(24) = 30
      CONSTS(25) = 1
      CONSTS(26) = 50000
      CONSTS(27) = 1000
      CONSTS(28) = 2
      CONSTS(29) = 0.02
      CONSTS(30) = 20
      CONSTS(31) = 0.2
      CONSTS(32) = 20
      CONSTS(33) = 20
      CONSTS(34) = 25000
      STATES(7) = 780
      CONSTS(35) = 17
      CONSTS(36) = 26
      CONSTS(37) = 1
      CONSTS(38) = 300000
      CONSTS(39) = 0.35
      CONSTS(40) = 1
      CONSTS(41) = 2.2
      CONSTS(42) = 10
      CONSTS(43) = CONSTS(23)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = 1.00000/(1.00000+EXP(- (16.0000+STATES(1)/1.00000)/5.00000))
      RATES(2) = (ALGBRC(1) - STATES(2))/CONSTS(8)
      ALGBRC(6) =  CONSTS(13)*STATES(3)
      ALGBRC(8) =  CONSTS(12)*(STATES(4) - STATES(3))
      ALGBRC(10) = ( CONSTS(15)*(ALGBRC(8) - ALGBRC(6)))/CONSTS(14)
      RATES(4) =  - CONSTS(10)*CONSTS(11)*ALGBRC(10)
      ALGBRC(4) = 1.00000/(1.00000+EXP(- (20.0000+STATES(1)/1.00000)/12.0000))
      ALGBRC(5) =  CONSTS(6)*ALGBRC(4)*(STATES(1) - CONSTS(7))
      ALGBRC(12) = - ( CONSTS(16)*ALGBRC(5)+ CONSTS(17)*STATES(3))
      RATES(3) =  CONSTS(9)*(ALGBRC(12)+ALGBRC(10))
      ALGBRC(21) =  ABS(STATES(7) - CONSTS(19) ** 2.00000 -  4.00000*STATES(7) ** 2.00000) ** (1.0 / 2)/1.00000
      ALGBRC(22) =  0.500000*((CONSTS(19) - STATES(7))+ ALGBRC(21)*1.00000)
      ALGBRC(9) =  0.200000* ABS(( STATES(6)*1.00000)/1.00000 ** 2.00000) ** (1.0 / 2)
      ALGBRC(19) =  CONSTS(41)*(ALGBRC(9)/(CONSTS(42)+ALGBRC(9)))
      ALGBRC(20) = CONSTS(40)+ALGBRC(19)
      RATES(7) = (ALGBRC(22) -  STATES(7)*EXP( ALGBRC(20)*(1.00000 - STATES(3)/CONSTS(39))))/( CONSTS(38)*1.00000)
      ALGBRC(2) =  CONSTS(2)*STATES(2)*(STATES(1) - CONSTS(3))
      ALGBRC(3) =  (CONSTS(4)/(1.00000+CONSTS(5)/STATES(3) ** 2.00000))*(STATES(1) - CONSTS(3))
      ALGBRC(16) =  0.165000*STATES(7)
      ALGBRC(18) =  0.0800000*(1.00000+( 2.00000*ALGBRC(16))/( CONSTS(35)*1.00000))+ 0.890000*ALGBRC(16)/( CONSTS(35)*1.00000) ** 2.00000
      ALGBRC(17) =  0.135000*STATES(7)
      ALGBRC(24) =  0.0500000*ALGBRC(22)
      ALGBRC(27) =  1.00000+ALGBRC(16)/( CONSTS(35)*1.00000) ** 2.00000*(1.00000+ALGBRC(17)/( CONSTS(36)*1.00000)+ALGBRC(24)/( CONSTS(37)*1.00000))
      ALGBRC(28) = ALGBRC(18)/ALGBRC(27)
      ALGBRC(30) =  CONSTS(34)*ALGBRC(28)*(STATES(1) - CONSTS(3))
      RATES(1) = - (ALGBRC(2)+ALGBRC(5)+ALGBRC(3)+ALGBRC(30))/CONSTS(1)
      ALGBRC(7) =  0.300000*STATES(5)
      ALGBRC(11) = ALGBRC(7) ** 2.00000/( CONSTS(26)*1.00000)
      ALGBRC(13) = CONSTS(43)+ALGBRC(11)
      ALGBRC(31) =  ALGBRC(7)*ALGBRC(22) ** 2.00000/( CONSTS(30)*CONSTS(26)*CONSTS(27)*1.00000 ** 2.00000)
      ALGBRC(32) = ALGBRC(13)+ALGBRC(31)
      ALGBRC(34) = ALGBRC(32)
      ALGBRC(37) = ALGBRC(34)
      ALGBRC(15) = ( STATES(6)*ALGBRC(7) ** 2.00000)/( CONSTS(25)*CONSTS(26)*CONSTS(31)*1.00000)
      ALGBRC(39) = ALGBRC(37)+ALGBRC(15)
      ALGBRC(41) = ( STATES(6)*ALGBRC(7) ** 2.00000*ALGBRC(22) ** 2.00000)/( CONSTS(25)*CONSTS(26)*CONSTS(27)*CONSTS(31)*CONSTS(32)*CONSTS(30)*1.00000 ** 2.00000)
      ALGBRC(42) = ALGBRC(39)+ALGBRC(41)
      ALGBRC(44) = ALGBRC(42)
      ALGBRC(45) = ALGBRC(44)
      ALGBRC(46) = ( STATES(7)*STATES(7))/ALGBRC(22)
      ALGBRC(51) = ( ALGBRC(46)*ALGBRC(7) ** 2.00000)/( CONSTS(24)*CONSTS(26)*CONSTS(29)*1.00000)
      ALGBRC(52) = ALGBRC(45)+ALGBRC(51)
      ALGBRC(54) = ( ALGBRC(46)*ALGBRC(7) ** 2.00000*ALGBRC(22) ** 2.00000)/( CONSTS(24)*CONSTS(26)*CONSTS(27)*CONSTS(29)*CONSTS(33)*CONSTS(30)*1.00000 ** 2.00000)
      ALGBRC(55) = ALGBRC(52)+ALGBRC(54)
      ALGBRC(58) = ALGBRC(55)
      ALGBRC(61) = ALGBRC(58)
      ALGBRC(64) = ALGBRC(61)
      ALGBRC(67) = ( ALGBRC(46)*STATES(6)*ALGBRC(7) ** 2.00000*ALGBRC(22) ** 2.00000)/( CONSTS(24)*CONSTS(25)*CONSTS(26)*CONSTS(27)*CONSTS(31)*CONSTS(29)*CONSTS(32)*CONSTS(33)*CONSTS(30)*1.00000 ** 2.00000)
      ALGBRC(68) = ALGBRC(64)+ALGBRC(67)
      ALGBRC(23) = ALGBRC(22) ** 2.00000/( CONSTS(27)*1.00000)
      ALGBRC(26) = CONSTS(22)+ALGBRC(23)
      ALGBRC(29) = ALGBRC(26)+ALGBRC(11)
      ALGBRC(33) = ALGBRC(29)+ALGBRC(31)
      ALGBRC(14) = STATES(6)/CONSTS(25)
      ALGBRC(35) = ALGBRC(33)+ALGBRC(14)
      ALGBRC(36) = ( STATES(6)*ALGBRC(22) ** 2.00000)/( CONSTS(25)*CONSTS(27)*CONSTS(32)*1.00000)
      ALGBRC(38) = ALGBRC(35)+ALGBRC(36)
      ALGBRC(40) = ALGBRC(38)+ALGBRC(15)
      ALGBRC(43) = ALGBRC(40)+ALGBRC(41)
      ALGBRC(47) = ALGBRC(46)/CONSTS(24)
      ALGBRC(48) = ALGBRC(43)+ALGBRC(47)
      ALGBRC(49) = ( ALGBRC(46)*ALGBRC(22) ** 2.00000)/( CONSTS(24)*CONSTS(27)*CONSTS(33)*1.00000)
      ALGBRC(50) = ALGBRC(48)+ALGBRC(49)
      ALGBRC(53) = ALGBRC(50)+ALGBRC(51)
      ALGBRC(56) = ALGBRC(53)+ALGBRC(54)
      ALGBRC(57) = ( ALGBRC(46)*STATES(6))/( CONSTS(24)*CONSTS(25))
      ALGBRC(59) = ALGBRC(56)+ALGBRC(57)
      ALGBRC(60) = ( ALGBRC(46)*STATES(6)*ALGBRC(22) ** 2.00000)/( CONSTS(24)*CONSTS(25)*CONSTS(27)*CONSTS(32)*CONSTS(33)*1.00000)
      ALGBRC(62) = ALGBRC(59)+ALGBRC(60)
      ALGBRC(63) = ( ALGBRC(46)*STATES(6)*ALGBRC(7) ** 2.00000)/( CONSTS(24)*CONSTS(25)*CONSTS(26)*CONSTS(31)*CONSTS(29)*1.00000)
      ALGBRC(65) = ALGBRC(62)+ALGBRC(63)
      ALGBRC(69) = ALGBRC(65)+ALGBRC(67)
      ALGBRC(66) = ALGBRC(63)
      ALGBRC(70) =  1.00000*(( CONSTS(20)*CONSTS(28)*ALGBRC(68)+ CONSTS(28)*ALGBRC(66))/ALGBRC(69))
      RATES(6) =  CONSTS(21)*(ALGBRC(70)/1.00000 -  0.500000*ALGBRC(9))
      RATES(5) =  CONSTS(21)*( CONSTS(18)*1.00000 - ALGBRC(70)/1.00000)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = 1.00000/(1.00000+EXP(- (16.0000+STATES(1)/1.00000)/5.00000))
      ALGBRC(6) =  CONSTS(13)*STATES(3)
      ALGBRC(8) =  CONSTS(12)*(STATES(4) - STATES(3))
      ALGBRC(10) = ( CONSTS(15)*(ALGBRC(8) - ALGBRC(6)))/CONSTS(14)
      ALGBRC(4) = 1.00000/(1.00000+EXP(- (20.0000+STATES(1)/1.00000)/12.0000))
      ALGBRC(5) =  CONSTS(6)*ALGBRC(4)*(STATES(1) - CONSTS(7))
      ALGBRC(12) = - ( CONSTS(16)*ALGBRC(5)+ CONSTS(17)*STATES(3))
      ALGBRC(21) =  ABS(STATES(7) - CONSTS(19) ** 2.00000 -  4.00000*STATES(7) ** 2.00000) ** (1.0 / 2)/1.00000
      ALGBRC(22) =  0.500000*((CONSTS(19) - STATES(7))+ ALGBRC(21)*1.00000)
      ALGBRC(9) =  0.200000* ABS(( STATES(6)*1.00000)/1.00000 ** 2.00000) ** (1.0 / 2)
      ALGBRC(19) =  CONSTS(41)*(ALGBRC(9)/(CONSTS(42)+ALGBRC(9)))
      ALGBRC(20) = CONSTS(40)+ALGBRC(19)
      ALGBRC(2) =  CONSTS(2)*STATES(2)*(STATES(1) - CONSTS(3))
      ALGBRC(3) =  (CONSTS(4)/(1.00000+CONSTS(5)/STATES(3) ** 2.00000))*(STATES(1) - CONSTS(3))
      ALGBRC(16) =  0.165000*STATES(7)
      ALGBRC(18) =  0.0800000*(1.00000+( 2.00000*ALGBRC(16))/( CONSTS(35)*1.00000))+ 0.890000*ALGBRC(16)/( CONSTS(35)*1.00000) ** 2.00000
      ALGBRC(17) =  0.135000*STATES(7)
      ALGBRC(24) =  0.0500000*ALGBRC(22)
      ALGBRC(27) =  1.00000+ALGBRC(16)/( CONSTS(35)*1.00000) ** 2.00000*(1.00000+ALGBRC(17)/( CONSTS(36)*1.00000)+ALGBRC(24)/( CONSTS(37)*1.00000))
      ALGBRC(28) = ALGBRC(18)/ALGBRC(27)
      ALGBRC(30) =  CONSTS(34)*ALGBRC(28)*(STATES(1) - CONSTS(3))
      ALGBRC(7) =  0.300000*STATES(5)
      ALGBRC(11) = ALGBRC(7) ** 2.00000/( CONSTS(26)*1.00000)
      ALGBRC(13) = CONSTS(43)+ALGBRC(11)
      ALGBRC(31) =  ALGBRC(7)*ALGBRC(22) ** 2.00000/( CONSTS(30)*CONSTS(26)*CONSTS(27)*1.00000 ** 2.00000)
      ALGBRC(32) = ALGBRC(13)+ALGBRC(31)
      ALGBRC(34) = ALGBRC(32)
      ALGBRC(37) = ALGBRC(34)
      ALGBRC(15) = ( STATES(6)*ALGBRC(7) ** 2.00000)/( CONSTS(25)*CONSTS(26)*CONSTS(31)*1.00000)
      ALGBRC(39) = ALGBRC(37)+ALGBRC(15)
      ALGBRC(41) = ( STATES(6)*ALGBRC(7) ** 2.00000*ALGBRC(22) ** 2.00000)/( CONSTS(25)*CONSTS(26)*CONSTS(27)*CONSTS(31)*CONSTS(32)*CONSTS(30)*1.00000 ** 2.00000)
      ALGBRC(42) = ALGBRC(39)+ALGBRC(41)
      ALGBRC(44) = ALGBRC(42)
      ALGBRC(45) = ALGBRC(44)
      ALGBRC(46) = ( STATES(7)*STATES(7))/ALGBRC(22)
      ALGBRC(51) = ( ALGBRC(46)*ALGBRC(7) ** 2.00000)/( CONSTS(24)*CONSTS(26)*CONSTS(29)*1.00000)
      ALGBRC(52) = ALGBRC(45)+ALGBRC(51)
      ALGBRC(54) = ( ALGBRC(46)*ALGBRC(7) ** 2.00000*ALGBRC(22) ** 2.00000)/( CONSTS(24)*CONSTS(26)*CONSTS(27)*CONSTS(29)*CONSTS(33)*CONSTS(30)*1.00000 ** 2.00000)
      ALGBRC(55) = ALGBRC(52)+ALGBRC(54)
      ALGBRC(58) = ALGBRC(55)
      ALGBRC(61) = ALGBRC(58)
      ALGBRC(64) = ALGBRC(61)
      ALGBRC(67) = ( ALGBRC(46)*STATES(6)*ALGBRC(7) ** 2.00000*ALGBRC(22) ** 2.00000)/( CONSTS(24)*CONSTS(25)*CONSTS(26)*CONSTS(27)*CONSTS(31)*CONSTS(29)*CONSTS(32)*CONSTS(33)*CONSTS(30)*1.00000 ** 2.00000)
      ALGBRC(68) = ALGBRC(64)+ALGBRC(67)
      ALGBRC(23) = ALGBRC(22) ** 2.00000/( CONSTS(27)*1.00000)
      ALGBRC(26) = CONSTS(22)+ALGBRC(23)
      ALGBRC(29) = ALGBRC(26)+ALGBRC(11)
      ALGBRC(33) = ALGBRC(29)+ALGBRC(31)
      ALGBRC(14) = STATES(6)/CONSTS(25)
      ALGBRC(35) = ALGBRC(33)+ALGBRC(14)
      ALGBRC(36) = ( STATES(6)*ALGBRC(22) ** 2.00000)/( CONSTS(25)*CONSTS(27)*CONSTS(32)*1.00000)
      ALGBRC(38) = ALGBRC(35)+ALGBRC(36)
      ALGBRC(40) = ALGBRC(38)+ALGBRC(15)
      ALGBRC(43) = ALGBRC(40)+ALGBRC(41)
      ALGBRC(47) = ALGBRC(46)/CONSTS(24)
      ALGBRC(48) = ALGBRC(43)+ALGBRC(47)
      ALGBRC(49) = ( ALGBRC(46)*ALGBRC(22) ** 2.00000)/( CONSTS(24)*CONSTS(27)*CONSTS(33)*1.00000)
      ALGBRC(50) = ALGBRC(48)+ALGBRC(49)
      ALGBRC(53) = ALGBRC(50)+ALGBRC(51)
      ALGBRC(56) = ALGBRC(53)+ALGBRC(54)
      ALGBRC(57) = ( ALGBRC(46)*STATES(6))/( CONSTS(24)*CONSTS(25))
      ALGBRC(59) = ALGBRC(56)+ALGBRC(57)
      ALGBRC(60) = ( ALGBRC(46)*STATES(6)*ALGBRC(22) ** 2.00000)/( CONSTS(24)*CONSTS(25)*CONSTS(27)*CONSTS(32)*CONSTS(33)*1.00000)
      ALGBRC(62) = ALGBRC(59)+ALGBRC(60)
      ALGBRC(63) = ( ALGBRC(46)*STATES(6)*ALGBRC(7) ** 2.00000)/( CONSTS(24)*CONSTS(25)*CONSTS(26)*CONSTS(31)*CONSTS(29)*1.00000)
      ALGBRC(65) = ALGBRC(62)+ALGBRC(63)
      ALGBRC(69) = ALGBRC(65)+ALGBRC(67)
      ALGBRC(66) = ALGBRC(63)
      ALGBRC(70) =  1.00000*(( CONSTS(20)*CONSTS(28)*ALGBRC(68)+ CONSTS(28)*ALGBRC(66))/ALGBRC(69))
      ALGBRC(25) = ALGBRC(22)/STATES(7)
      RETURN
      END