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