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 43 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 35 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is k0_12 in component parameters (per_M_per_second).
C CONSTS(2) is k0_21 in component parameters (per_second).
C CONSTS(3) is k0_23 in component parameters (per_M_per_second).
C CONSTS(4) is k0_32 in component parameters (per_second).
C CONSTS(5) is k0_34 in component parameters (per_second).
C CONSTS(6) is k0_43 in component parameters (per_second).
C CONSTS(7) is k0_45 in component parameters (per_second).
C CONSTS(8) is k0_54 in component parameters (per_M_per_second).
C CONSTS(9) is k0_25 in component parameters (per_second).
C CONSTS(10) is k0_52 in component parameters (per_second).
C CONSTS(11) is k0_56 in component parameters (per_second).
C CONSTS(12) is k0_65 in component parameters (per_M_per_second).
C CONSTS(13) is k0_61 in component parameters (per_second).
C CONSTS(14) is k0_16 in component parameters (per_second).
C CONSTS(15) is delta in component parameters (dimensionless).
C CONSTS(30) is alpha_p in component parameters (dimensionless).
C CONSTS(16) is alpha_pp in component parameters (dimensionless).
C CONSTS(17) is N_C in component parameters (dimensionless).
C CONSTS(18) is N_Avo in component parameters (per_mol).
C CONSTS(19) is area in component parameters (um2).
C CONSTS(31) is C_T in component parameters (umol).
C CONSTS(20) is n in component parameters (dimensionless).
C CONSTS(21) is z_c in component parameters (dimensionless).
C CONSTS(22) is z_Na in component parameters (dimensionless).
C CONSTS(23) is F in component parameters (C_per_mol).
C CONSTS(24) is R in component parameters (J_per_K_per_mol).
C CONSTS(25) is T in component parameters (kelvin).
C STATES(1) is V in component ion_concentrations (volt).
C ALGBRC(1) is mu in component parameters (dimensionless).
C CONSTS(26) is Na_i in component ion_concentrations (M).
C CONSTS(27) is Na_o in component ion_concentrations (M).
C CONSTS(28) is glucose_i in component ion_concentrations (M).
C CONSTS(29) is glucose_o in component ion_concentrations (M).
C ALGBRC(4) is k_12 in component rate_constants (per_second).
C ALGBRC(5) is k_21 in component rate_constants (per_second).
C CONSTS(32) is k_23 in component rate_constants (per_second).
C CONSTS(33) is k_32 in component rate_constants (per_second).
C ALGBRC(6) is k_34 in component rate_constants (per_second).
C ALGBRC(7) is k_43 in component rate_constants (per_second).
C CONSTS(34) is k_45 in component rate_constants (per_second).
C ALGBRC(15) is k_54 in component rate_constants (per_second).
C ALGBRC(8) is k_25 in component rate_constants (per_second).
C ALGBRC(13) is k_52 in component rate_constants (per_second).
C ALGBRC(9) is k_56 in component rate_constants (per_second).
C ALGBRC(10) is k_65 in component rate_constants (per_second).
C ALGBRC(11) is k_61 in component rate_constants (per_second).
C ALGBRC(12) is k_16 in component rate_constants (per_second).
C ALGBRC(2) is ks_12 in component rate_constants (per_M_per_second).
C ALGBRC(14) is k0_54_temp in component rate_constants (per_M_per_second).
C ALGBRC(3) is k_52_temp in component rate_constants (per_second).
C STATES(2) is C_1 in component kinetic_equations (umol).
C STATES(3) is C_2 in component kinetic_equations (umol).
C STATES(4) is C_3 in component kinetic_equations (umol).
C STATES(5) is C_4 in component kinetic_equations (umol).
C STATES(6) is C_5 in component kinetic_equations (umol).
C ALGBRC(16) is C_6 in component kinetic_equations (umol).
C STATES(7) is C_6_temp in component kinetic_equations (umol).
C ALGBRC(17) is C1_sum in component king_altman_states (per_second5).
C ALGBRC(19) is C2_sum in component king_altman_states (per_second5).
C ALGBRC(23) is C3_sum in component king_altman_states (per_second5).
C ALGBRC(25) is C4_sum in component king_altman_states (per_second5).
C ALGBRC(28) is C5_sum in component king_altman_states (per_second5).
C ALGBRC(32) is C6_sum in component king_altman_states (per_second5).
C ALGBRC(35) is C_sum in component king_altman_states (per_second5).
C ALGBRC(36) is C1 in component king_altman_states (umol).
C ALGBRC(37) is C2 in component king_altman_states (umol).
C ALGBRC(38) is C3 in component king_altman_states (umol).
C ALGBRC(39) is C4 in component king_altman_states (umol).
C ALGBRC(40) is C5 in component king_altman_states (umol).
C ALGBRC(41) is C6 in component king_altman_states (umol).
C ALGBRC(20) is I_NaGl_pSS in component NBC_current (uA).
C ALGBRC(43) is I_NaGl_SS in component NBC_current (uA).
C ALGBRC(42) is J_Na in component NBC_current (umol_per_second).
C ALGBRC(29) is epsilon in component phenomonological_constants (per_second).
C ALGBRC(18) is lambda in component phenomonological_constants (per_M2_per_second5).
C ALGBRC(21) is chi in component phenomonological_constants (M).
C ALGBRC(27) is alpha in component phenomonological_constants (M2).
C ALGBRC(24) is beta in component phenomonological_constants (M).
C ALGBRC(22) is gamma in component phenomonological_constants (M2_per_second).
C ALGBRC(26) is phi in component phenomonological_constants (M_per_second).
C ALGBRC(33) is Imax_Na in component phenomonological_constants (uA).
C ALGBRC(34) is Imax_gluc in component phenomonological_constants (uA).
C ALGBRC(30) is Khalf_Na in component phenomonological_constants (M).
C ALGBRC(31) is Khalf_gluc in component phenomonological_constants (M).
C RATES(1) is d/dt V in component ion_concentrations (volt).
C RATES(2) is d/dt C_1 in component kinetic_equations (umol).
C RATES(3) is d/dt C_2 in component kinetic_equations (umol).
C RATES(4) is d/dt C_3 in component kinetic_equations (umol).
C RATES(5) is d/dt C_4 in component kinetic_equations (umol).
C RATES(6) is d/dt C_5 in component kinetic_equations (umol).
C RATES(7) is d/dt C_6_temp in component kinetic_equations (umol).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 20000
      CONSTS(2) = 400
      CONSTS(3) = 1e4
      CONSTS(4) = 20
      CONSTS(5) = 50
      CONSTS(6) = 50
      CONSTS(7) = 800
      CONSTS(8) = 6.7e6
      CONSTS(9) = 0.01
      CONSTS(10) = 0.17
      CONSTS(11) = 48
      CONSTS(12) = 50
      CONSTS(13) = 35
      CONSTS(14) = 100
      CONSTS(15) = 0.7
      CONSTS(16) = 0
      CONSTS(17) = 8.4e10
      CONSTS(18) = 6.022e23
      CONSTS(19) = 1e6
      CONSTS(20) = 1
      CONSTS(21) = -1
      CONSTS(22) = 1
      CONSTS(23) = 96485.34
      CONSTS(24) = 8.314
      CONSTS(25) = 310
      STATES(1) = -150e-3
      CONSTS(26) = 10e-3
      CONSTS(27) = 100e-3
      CONSTS(28) = 0e-3
      CONSTS(29) = 1e-3
      STATES(2) = 3.923e-9
      STATES(3) = 1.016e-7
      STATES(4) = 1.515e-8
      STATES(5) = 8.912e-10
      STATES(6) = 1.485e-8
      STATES(7) = 3.06798e-9
      CONSTS(30) = (1.00000 - CONSTS(15)) - CONSTS(16)
      CONSTS(34) = 0.000100000
      CONSTS(31) = ( 1.00000e+06*CONSTS(17))/CONSTS(18)
      CONSTS(32) =  CONSTS(3)*CONSTS(29)
      CONSTS(33) = CONSTS(4)
      CONSTS(34) = CONSTS(7)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = CONSTS(34)
      ALGBRC(1) = ( CONSTS(23)*STATES(1))/( CONSTS(24)*CONSTS(25))
      ALGBRC(6) =  CONSTS(5)*EXP(( - (CONSTS(21)+CONSTS(20))*CONSTS(15)*ALGBRC(1))/2.00000)
      ALGBRC(7) =  CONSTS(6)*EXP(( (CONSTS(21)+CONSTS(20))*CONSTS(15)*ALGBRC(1))/2.00000)
      RATES(4) = ( CONSTS(32)*STATES(3)+ ALGBRC(7)*STATES(5)) -  (CONSTS(33)+ALGBRC(6))*STATES(4)
      ALGBRC(9) =  CONSTS(11)*EXP(( - CONSTS(20)*CONSTS(22)*CONSTS(16)*ALGBRC(1))/2.00000)
      ALGBRC(10) =  CONSTS(12)*CONSTS(26) ** CONSTS(20)*EXP(( CONSTS(20)*CONSTS(22)*CONSTS(16)*ALGBRC(1))/2.00000)
      ALGBRC(11) =  CONSTS(13)*EXP(( CONSTS(21)*CONSTS(15)*ALGBRC(1))/2.00000)
      ALGBRC(12) =  CONSTS(14)*EXP(( - CONSTS(21)*CONSTS(15)*ALGBRC(1))/2.00000)
      RATES(7) = ( ALGBRC(12)*STATES(2)+ ALGBRC(9)*STATES(6)) -  (ALGBRC(11)+ALGBRC(10))*STATES(7)
      ALGBRC(2) =  CONSTS(1)*EXP(( - CONSTS(20)*CONSTS(30)*ALGBRC(1))/2.00000)
      ALGBRC(4) =  ALGBRC(2)*CONSTS(27) ** CONSTS(20)
      ALGBRC(5) =  CONSTS(2)*EXP(( CONSTS(20)*CONSTS(22)*CONSTS(30)*ALGBRC(1))/2.00000)
      ALGBRC(8) =  CONSTS(9)*EXP(( - (CONSTS(21)+CONSTS(20))*CONSTS(15)*ALGBRC(1))/2.00000)
      ALGBRC(13) = ( CONSTS(1)*ALGBRC(8)*CONSTS(11)*CONSTS(13))/( CONSTS(2)*CONSTS(14)*CONSTS(12))
      RATES(3) = ( ALGBRC(4)*STATES(2)+ CONSTS(33)*STATES(4)+ ALGBRC(13)*STATES(6)) -  (ALGBRC(5)+CONSTS(32)+ALGBRC(8))*STATES(3)
      ALGBRC(14) = ( CONSTS(3)*ALGBRC(6)*CONSTS(34)*ALGBRC(13))/( ALGBRC(7)*CONSTS(33)*ALGBRC(8))
      ALGBRC(15) =  ALGBRC(14)*CONSTS(28)
      RATES(5) = ( ALGBRC(6)*STATES(4)+ ALGBRC(15)*STATES(6)) -  (CONSTS(34)+ALGBRC(7))*STATES(5)
      ALGBRC(16) = CONSTS(31) - (STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6))
      RATES(2) = ( ALGBRC(5)*STATES(3)+ ALGBRC(11)*ALGBRC(16)) -  (ALGBRC(4)+ALGBRC(12))*STATES(2)
      RATES(6) = ( CONSTS(34)*STATES(5)+ ALGBRC(10)*ALGBRC(16)+ ALGBRC(8)*STATES(3)) -  (ALGBRC(15)+ALGBRC(13)+ALGBRC(9))*STATES(6)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = ( CONSTS(23)*STATES(1))/( CONSTS(24)*CONSTS(25))
      ALGBRC(6) =  CONSTS(5)*EXP(( - (CONSTS(21)+CONSTS(20))*CONSTS(15)*ALGBRC(1))/2.00000)
      ALGBRC(7) =  CONSTS(6)*EXP(( (CONSTS(21)+CONSTS(20))*CONSTS(15)*ALGBRC(1))/2.00000)
      ALGBRC(9) =  CONSTS(11)*EXP(( - CONSTS(20)*CONSTS(22)*CONSTS(16)*ALGBRC(1))/2.00000)
      ALGBRC(10) =  CONSTS(12)*CONSTS(26) ** CONSTS(20)*EXP(( CONSTS(20)*CONSTS(22)*CONSTS(16)*ALGBRC(1))/2.00000)
      ALGBRC(11) =  CONSTS(13)*EXP(( CONSTS(21)*CONSTS(15)*ALGBRC(1))/2.00000)
      ALGBRC(12) =  CONSTS(14)*EXP(( - CONSTS(21)*CONSTS(15)*ALGBRC(1))/2.00000)
      ALGBRC(2) =  CONSTS(1)*EXP(( - CONSTS(20)*CONSTS(30)*ALGBRC(1))/2.00000)
      ALGBRC(4) =  ALGBRC(2)*CONSTS(27) ** CONSTS(20)
      ALGBRC(5) =  CONSTS(2)*EXP(( CONSTS(20)*CONSTS(22)*CONSTS(30)*ALGBRC(1))/2.00000)
      ALGBRC(8) =  CONSTS(9)*EXP(( - (CONSTS(21)+CONSTS(20))*CONSTS(15)*ALGBRC(1))/2.00000)
      ALGBRC(13) = ( CONSTS(1)*ALGBRC(8)*CONSTS(11)*CONSTS(13))/( CONSTS(2)*CONSTS(14)*CONSTS(12))
      ALGBRC(14) = ( CONSTS(3)*ALGBRC(6)*CONSTS(34)*ALGBRC(13))/( ALGBRC(7)*CONSTS(33)*ALGBRC(8))
      ALGBRC(15) =  ALGBRC(14)*CONSTS(28)
      ALGBRC(16) = CONSTS(31) - (STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6))
      ALGBRC(3) =  CONSTS(10)*EXP(( (CONSTS(21)+CONSTS(20))*CONSTS(15)*ALGBRC(1))/2.00000)
      ALGBRC(17) =  ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ CONSTS(32)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)
      ALGBRC(18) =  ALGBRC(2)*CONSTS(3)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*CONSTS(34)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*CONSTS(34)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*CONSTS(34)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)
      ALGBRC(19) =  ALGBRC(12)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(33)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(33)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(33)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(4)*CONSTS(33)*CONSTS(34)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(4)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)
      ALGBRC(20) =  - CONSTS(23)*( CONSTS(20)*CONSTS(22)*CONSTS(30)*( ALGBRC(4)*STATES(2) -  ALGBRC(5)*STATES(3))+ CONSTS(21)*CONSTS(15)*( ALGBRC(12)*STATES(2) -  ALGBRC(11)*ALGBRC(16))+ CONSTS(20)*CONSTS(22)*CONSTS(16)*( ALGBRC(9)*STATES(6) -  ALGBRC(10)*ALGBRC(16)))
      ALGBRC(21) =  (1.00000/ALGBRC(18))*( ALGBRC(2)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*CONSTS(33)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(33)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(33)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*CONSTS(33)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(33)*CONSTS(34)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)+ ALGBRC(2)*ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(9)+ ALGBRC(2)*ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(9))
      ALGBRC(22) =  (1.00000/ALGBRC(18))*( ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10))
      ALGBRC(23) =  ALGBRC(12)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(32)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*CONSTS(32)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(32)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(32)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(32)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(32)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(32)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(32)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(32)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(32)*CONSTS(34)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(4)*CONSTS(32)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)
      ALGBRC(24) =  (1.00000/ALGBRC(18))*( CONSTS(3)*ALGBRC(12)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ CONSTS(3)*ALGBRC(12)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ CONSTS(3)*ALGBRC(12)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ CONSTS(3)*ALGBRC(12)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ CONSTS(3)*ALGBRC(12)*ALGBRC(6)*ALGBRC(13)*ALGBRC(10)+ CONSTS(3)*ALGBRC(12)*ALGBRC(6)*CONSTS(34)*ALGBRC(10)+ CONSTS(3)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ CONSTS(3)*ALGBRC(12)*ALGBRC(6)*CONSTS(34)*ALGBRC(9))
      ALGBRC(25) =  ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(32)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*CONSTS(32)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(33)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*CONSTS(33)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(32)*ALGBRC(6)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(32)*ALGBRC(6)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(32)*ALGBRC(6)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(32)*ALGBRC(6)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*CONSTS(33)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(32)*ALGBRC(6)*ALGBRC(13)*ALGBRC(11)
      ALGBRC(26) =  (1.00000/ALGBRC(18))*(( - ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)*ALGBRC(11) -  ALGBRC(2)*ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)) -  ALGBRC(2)*ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11))
      ALGBRC(27) =  (1.00000/ALGBRC(18))*( ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(33)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(33)*CONSTS(34)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(13))
      ALGBRC(28) =  ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(10)+ ALGBRC(12)*CONSTS(32)*ALGBRC(6)*CONSTS(34)*ALGBRC(10)+ ALGBRC(4)*CONSTS(32)*ALGBRC(6)*CONSTS(34)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(11)+ ALGBRC(4)*CONSTS(32)*ALGBRC(6)*CONSTS(34)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(11)
      ALGBRC(29) =  (1.00000/ALGBRC(18))*- ALGBRC(2)*CONSTS(3)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)*ALGBRC(11)
      ALGBRC(30) = (ALGBRC(27)+ ALGBRC(24)*CONSTS(29))/(ALGBRC(21)+CONSTS(29))
      ALGBRC(31) = (ALGBRC(27)+ ALGBRC(21)*CONSTS(27))/(ALGBRC(24)+CONSTS(27))
      ALGBRC(32) =  ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)+ ALGBRC(12)*CONSTS(32)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)+ ALGBRC(4)*CONSTS(32)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(9)+ ALGBRC(4)*ALGBRC(8)*CONSTS(33)*CONSTS(34)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(4)*ALGBRC(8)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(15)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(34)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*CONSTS(34)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(33)*ALGBRC(7)*ALGBRC(13)
      ALGBRC(33) = ( 2.00000*CONSTS(23)*CONSTS(31)*(ALGBRC(26)+ ALGBRC(29)*CONSTS(29)))/(ALGBRC(21)+CONSTS(29))
      ALGBRC(34) = ( 2.00000*CONSTS(23)*CONSTS(31)*ALGBRC(29)*CONSTS(27))/(ALGBRC(24)+CONSTS(27))
      ALGBRC(35) = ALGBRC(17)+ALGBRC(19)+ALGBRC(23)+ALGBRC(25)+ALGBRC(28)+ALGBRC(32)
      ALGBRC(36) = ( CONSTS(31)*ALGBRC(17))/ALGBRC(35)
      ALGBRC(37) = ( CONSTS(31)*ALGBRC(19))/ALGBRC(35)
      ALGBRC(38) = ( CONSTS(31)*ALGBRC(23))/ALGBRC(35)
      ALGBRC(39) = ( CONSTS(31)*ALGBRC(25))/ALGBRC(35)
      ALGBRC(40) = ( CONSTS(31)*ALGBRC(28))/ALGBRC(35)
      ALGBRC(41) = ( CONSTS(31)*ALGBRC(32))/ALGBRC(35)
      ALGBRC(42) = ( ALGBRC(8)*ALGBRC(37) -  ALGBRC(13)*ALGBRC(40))+( ALGBRC(6)*ALGBRC(38) -  ALGBRC(7)*ALGBRC(39))
      ALGBRC(43) =  - CONSTS(23)*( CONSTS(21)*( ALGBRC(12)*ALGBRC(36) -  ALGBRC(11)*ALGBRC(41))+ (CONSTS(21)+ CONSTS(22)*CONSTS(20))*( ALGBRC(8)*ALGBRC(37) -  ALGBRC(13)*ALGBRC(40))+ (CONSTS(21)+ CONSTS(22)*CONSTS(20))*( ALGBRC(6)*ALGBRC(38) -  ALGBRC(7)*ALGBRC(39)))
      RETURN
      END