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 41 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_M2_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_f in component parameters (per_M2_per_second).
C CONSTS(13) is k0_61_f in component parameters (per_second).
C CONSTS(14) is k0_16_f in component parameters (per_second).
C CONSTS(15) is k0_65_r in component parameters (per_M2_per_second).
C CONSTS(16) is k0_61_r in component parameters (per_second).
C CONSTS(17) is k0_16_r in component parameters (per_second).
C CONSTS(33) is k0_65 in component parameters (per_M2_per_second).
C CONSTS(34) is k0_61 in component parameters (per_second).
C CONSTS(35) is k0_16 in component parameters (per_second).
C CONSTS(18) is delta in component parameters (dimensionless).
C CONSTS(36) is alpha_p in component parameters (dimensionless).
C CONSTS(19) is alpha_pp in component parameters (dimensionless).
C CONSTS(20) is N_C in component parameters (dimensionless).
C CONSTS(21) is N_Avo in component parameters (per_mol).
C CONSTS(22) is area in component parameters (um2).
C CONSTS(37) is C_T in component parameters (umol).
C CONSTS(23) is n in component parameters (dimensionless).
C CONSTS(24) is z_c in component parameters (dimensionless).
C CONSTS(25) is z_Na in component parameters (dimensionless).
C CONSTS(26) is F in component parameters (C_per_mol).
C CONSTS(27) is R in component parameters (J_per_K_per_mol).
C CONSTS(28) 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(29) is Na_o in component ion_concentrations (M).
C CONSTS(30) is Na_i in component ion_concentrations (M).
C CONSTS(31) is glucose_i in component ion_concentrations (M).
C CONSTS(32) 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(38) is k_23 in component rate_constants (per_second).
C CONSTS(39) 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(40) 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_M2_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(36) is C_sum in component king_altman_states (per_second5).
C ALGBRC(37) is C1 in component king_altman_states (umol).
C ALGBRC(38) is C2 in component king_altman_states (umol).
C ALGBRC(39) is C3 in component king_altman_states (umol).
C ALGBRC(40) is C4 in component king_altman_states (umol).
C ALGBRC(41) is C5 in component king_altman_states (umol).
C ALGBRC(42) 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(29) is epsilon in component phenomonological_constants (per_second).
C ALGBRC(18) is lambda in component phenomonological_constants (per_M3_per_second5).
C ALGBRC(21) is chi in component phenomonological_constants (M).
C ALGBRC(27) is alpha in component phenomonological_constants (M3).
C ALGBRC(24) is beta in component phenomonological_constants (M2).
C ALGBRC(22) is gamma in component phenomonological_constants (M3_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_sq in component phenomonological_constants (M2).
C ALGBRC(35) 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) = 80000
      CONSTS(2) = 500
      CONSTS(3) = 1e5
      CONSTS(4) = 20
      CONSTS(5) = 50
      CONSTS(6) = 50
      CONSTS(7) = 800
      CONSTS(8) = 1.219e4
      CONSTS(9) = 0.3
      CONSTS(10) = 9.1e-4
      CONSTS(11) = 10
      CONSTS(12) = 50
      CONSTS(13) = 5
      CONSTS(14) = 35
      CONSTS(15) = 4500
      CONSTS(16) = 3
      CONSTS(17) = 350
      CONSTS(18) = 0.7
      CONSTS(19) = 0
      CONSTS(20) = 3e6
      CONSTS(21) = 6.022e23
      CONSTS(22) = 1e6
      CONSTS(23) = 2
      CONSTS(24) = -2
      CONSTS(25) = 1
      CONSTS(26) = 96485.34
      CONSTS(27) = 8.314
      CONSTS(28) = 310
      STATES(1) = -150e-3
      CONSTS(29) = 10e-3
      CONSTS(30) = 500e-3
      CONSTS(31) = 100e-3
      CONSTS(32) = 0e-3
      STATES(2) = 0
      STATES(3) = 0
      STATES(4) = 0
      STATES(5) = 0
      STATES(6) = 0
      STATES(7) = 0
      CONSTS(33) = TERNRY(CONSTS(29).LT.0.100000, CONSTS(15), CONSTS(12))
      CONSTS(40) = 0.200000
      CONSTS(34) = TERNRY(CONSTS(29).LT.0.100000, CONSTS(16), CONSTS(13))
      CONSTS(35) = TERNRY(CONSTS(29).LT.0.100000, CONSTS(17), CONSTS(14))
      CONSTS(36) = (1.00000 - CONSTS(18)) - CONSTS(19)
      CONSTS(37) = ( 1.00000e+06*CONSTS(20))/CONSTS(21)
      CONSTS(38) =  CONSTS(3)*CONSTS(32)
      CONSTS(39) = CONSTS(4)
      CONSTS(40) = CONSTS(7)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = CONSTS(40)
      ALGBRC(1) = ( CONSTS(26)*STATES(1))/( CONSTS(27)*CONSTS(28))
      ALGBRC(6) =  CONSTS(5)*EXP(( - (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(7) =  CONSTS(6)*EXP(( (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      RATES(4) = ( CONSTS(38)*STATES(3)+ ALGBRC(7)*STATES(5)) -  (CONSTS(39)+ALGBRC(6))*STATES(4)
      ALGBRC(9) =  CONSTS(11)*EXP(( - CONSTS(23)*CONSTS(25)*CONSTS(19)*ALGBRC(1))/2.00000)
      ALGBRC(10) =  CONSTS(33)*CONSTS(30) ** CONSTS(23)*EXP(( CONSTS(23)*CONSTS(25)*CONSTS(19)*ALGBRC(1))/2.00000)
      ALGBRC(11) =  CONSTS(34)*EXP(( CONSTS(24)*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(12) =  CONSTS(35)*EXP(( - CONSTS(24)*CONSTS(18)*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(23)*CONSTS(36)*ALGBRC(1))/2.00000)
      ALGBRC(4) =  ALGBRC(2)*CONSTS(29) ** CONSTS(23)
      ALGBRC(5) =  CONSTS(2)*EXP(( CONSTS(23)*CONSTS(25)*CONSTS(36)*ALGBRC(1))/2.00000)
      ALGBRC(8) =  CONSTS(9)*EXP(( - (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(13) = ( CONSTS(1)*ALGBRC(8)*CONSTS(11)*CONSTS(34))/( CONSTS(2)*CONSTS(35)*CONSTS(33))
      RATES(3) = ( ALGBRC(4)*STATES(2)+ CONSTS(39)*STATES(4)+ ALGBRC(13)*STATES(6)) -  (ALGBRC(5)+CONSTS(38)+ALGBRC(8))*STATES(3)
      ALGBRC(14) = ( CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(13))/( ALGBRC(7)*CONSTS(39)*ALGBRC(8))
      ALGBRC(15) =  ALGBRC(14)*CONSTS(31)
      RATES(5) = ( ALGBRC(6)*STATES(4)+ ALGBRC(15)*STATES(6)) -  (CONSTS(40)+ALGBRC(7))*STATES(5)
      ALGBRC(16) = CONSTS(37) - (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(40)*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(26)*STATES(1))/( CONSTS(27)*CONSTS(28))
      ALGBRC(6) =  CONSTS(5)*EXP(( - (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(7) =  CONSTS(6)*EXP(( (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(9) =  CONSTS(11)*EXP(( - CONSTS(23)*CONSTS(25)*CONSTS(19)*ALGBRC(1))/2.00000)
      ALGBRC(10) =  CONSTS(33)*CONSTS(30) ** CONSTS(23)*EXP(( CONSTS(23)*CONSTS(25)*CONSTS(19)*ALGBRC(1))/2.00000)
      ALGBRC(11) =  CONSTS(34)*EXP(( CONSTS(24)*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(12) =  CONSTS(35)*EXP(( - CONSTS(24)*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(2) =  CONSTS(1)*EXP(( - CONSTS(23)*CONSTS(36)*ALGBRC(1))/2.00000)
      ALGBRC(4) =  ALGBRC(2)*CONSTS(29) ** CONSTS(23)
      ALGBRC(5) =  CONSTS(2)*EXP(( CONSTS(23)*CONSTS(25)*CONSTS(36)*ALGBRC(1))/2.00000)
      ALGBRC(8) =  CONSTS(9)*EXP(( - (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(13) = ( CONSTS(1)*ALGBRC(8)*CONSTS(11)*CONSTS(34))/( CONSTS(2)*CONSTS(35)*CONSTS(33))
      ALGBRC(14) = ( CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(13))/( ALGBRC(7)*CONSTS(39)*ALGBRC(8))
      ALGBRC(15) =  ALGBRC(14)*CONSTS(31)
      ALGBRC(16) = CONSTS(37) - (STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6))
      ALGBRC(3) =  CONSTS(10)*EXP(( (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(17) =  ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*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(40)*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(40)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*CONSTS(40)*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(40)*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(40)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)
      ALGBRC(19) =  ALGBRC(12)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(4)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(4)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)
      ALGBRC(20) =  - CONSTS(26)*( CONSTS(23)*CONSTS(25)*CONSTS(36)*( ALGBRC(4)*STATES(2) -  ALGBRC(5)*STATES(3))+ CONSTS(24)*CONSTS(18)*( ALGBRC(12)*STATES(2) -  ALGBRC(11)*ALGBRC(16))+ CONSTS(23)*CONSTS(25)*CONSTS(19)*( ALGBRC(9)*STATES(6) -  ALGBRC(10)*ALGBRC(16)))
      ALGBRC(21) =  (1.00000/ALGBRC(18))*( ALGBRC(2)*CONSTS(39)*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(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(39)*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(39)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9))
      ALGBRC(22) =  (1.00000/ALGBRC(18))*( ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*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(38)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(38)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(38)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*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(40)*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(40)*ALGBRC(10)+ CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ CONSTS(3)*ALGBRC(12)*ALGBRC(6)*CONSTS(40)*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(38)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(38)*ALGBRC(6)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*ALGBRC(13)*ALGBRC(11)
      ALGBRC(26) =  (1.00000/ALGBRC(18))*(( - ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11) -  ALGBRC(2)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)) -  ALGBRC(2)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11))
      ALGBRC(27) =  (1.00000/ALGBRC(18))*( ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(39)*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(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13))
      ALGBRC(28) =  ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(11)
      ALGBRC(29) =  (1.00000/ALGBRC(18))*- ALGBRC(2)*CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)
      ALGBRC(30) = (ALGBRC(27)+ ALGBRC(24)*CONSTS(32))/(ALGBRC(21)+CONSTS(32))
      ALGBRC(31) = (ALGBRC(27)+ ALGBRC(21)*CONSTS(29) ** 2.00000)/(ALGBRC(24)+CONSTS(29) ** 2.00000)
      ALGBRC(32) =  ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)
      ALGBRC(33) = ( 2.00000*CONSTS(26)*CONSTS(37)*(ALGBRC(26)+ ALGBRC(29)*CONSTS(32)))/(ALGBRC(21)+CONSTS(32))
      ALGBRC(34) = ( 2.00000*CONSTS(26)*CONSTS(37)*ALGBRC(29)*CONSTS(29) ** 2.00000)/(ALGBRC(24)+CONSTS(29) ** 2.00000)
      ALGBRC(35) =  ALGBRC(30) ** (1.0 / 2)
      ALGBRC(36) = ALGBRC(17)+ALGBRC(19)+ALGBRC(23)+ALGBRC(25)+ALGBRC(28)+ALGBRC(32)
      ALGBRC(37) = ( CONSTS(37)*ALGBRC(17))/ALGBRC(36)
      ALGBRC(38) = ( CONSTS(37)*ALGBRC(19))/ALGBRC(36)
      ALGBRC(39) = ( CONSTS(37)*ALGBRC(23))/ALGBRC(36)
      ALGBRC(40) = ( CONSTS(37)*ALGBRC(25))/ALGBRC(36)
      ALGBRC(41) = ( CONSTS(37)*ALGBRC(28))/ALGBRC(36)
      ALGBRC(42) = ( CONSTS(37)*ALGBRC(32))/ALGBRC(36)
      ALGBRC(43) =  - CONSTS(26)*( CONSTS(24)*( ALGBRC(12)*ALGBRC(37) -  ALGBRC(11)*ALGBRC(42))+ (CONSTS(24)+ CONSTS(25)*CONSTS(23))*( ALGBRC(8)*ALGBRC(38) -  ALGBRC(13)*ALGBRC(41))+ (CONSTS(24)+ CONSTS(25)*CONSTS(23))*( ALGBRC(6)*ALGBRC(39) -  ALGBRC(7)*ALGBRC(40)))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END