C
C There are a total of 34 entries in the algebraic variable array.
C There are a total of 15 entries in each of the rate and state variable arrays.
C There are a total of 59 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C ALGBRC(1) is L_iso in component beta_1_adrenergic_parameters (uM).
C CONSTS(1) is K_H in component beta_1_adrenergic_parameters (uM).
C CONSTS(2) is K_L in component beta_1_adrenergic_parameters (uM).
C CONSTS(3) is K_C in component beta_1_adrenergic_parameters (uM).
C ALGBRC(2) is L_ach in component muscarinic_parameters (uM).
C CONSTS(4) is K_H in component muscarinic_parameters (uM).
C CONSTS(5) is K_L in component muscarinic_parameters (uM).
C CONSTS(6) is K_C in component muscarinic_parameters (uM).
C CONSTS(7) is k_PDE2 in component PDE_parameters (per_sec).
C CONSTS(8) is Km_PDE2 in component PDE_parameters (uM).
C CONSTS(9) is k_PDE3 in component PDE_parameters (per_sec).
C CONSTS(10) is Km_PDE3 in component PDE_parameters (uM).
C CONSTS(11) is k_PDE4 in component PDE_parameters (per_sec).
C CONSTS(12) is Km_PDE4 in component PDE_parameters (uM).
C CONSTS(13) is k_act1 in component G_s_parameters (per_sec).
C CONSTS(14) is k_act2 in component G_s_parameters (per_sec).
C CONSTS(15) is k_hydr in component G_s_parameters (per_sec).
C CONSTS(16) is k_reas in component G_s_parameters (per_uM_per_sec).
C CONSTS(17) is k_act1 in component G_i_parameters (per_sec).
C CONSTS(18) is k_act2 in component G_i_parameters (per_sec).
C CONSTS(19) is k_hydr in component G_i_parameters (per_sec).
C CONSTS(20) is k_reas in component G_i_parameters (per_uM_per_sec).
C ALGBRC(13) is R in component caveolar_beta_1_adrenergic_receptor_module (uM).
C ALGBRC(14) is LR in component caveolar_beta_1_adrenergic_receptor_module (uM).
C ALGBRC(15) is LRG in component caveolar_beta_1_adrenergic_receptor_module (uM).
C ALGBRC(16) is RG in component caveolar_beta_1_adrenergic_receptor_module (uM).
C CONSTS(21) is R_Total in component caveolar_beta_1_adrenergic_receptor_module (uM).
C ALGBRC(6) is Gs_alpha_beta_gamma in component caveolar_G_s_protein_activation_module (uM).
C ALGBRC(17) is R in component caveolar_muscarinic_receptor_module (uM).
C ALGBRC(18) is LR in component caveolar_muscarinic_receptor_module (uM).
C ALGBRC(19) is LRG in component caveolar_muscarinic_receptor_module (uM).
C ALGBRC(20) is RG in component caveolar_muscarinic_receptor_module (uM).
C CONSTS(22) is R_Total in component caveolar_muscarinic_receptor_module (uM).
C ALGBRC(7) is Gi_alpha_beta_gamma in component caveolar_G_i_protein_activation_module (uM).
C STATES(1) is Gs_alpha_GTP in component caveolar_G_s_protein_activation_module (uM).
C STATES(2) is Gs_beta_gamma in component caveolar_G_s_protein_activation_module (uM).
C STATES(3) is Gs_alpha_GDP in component caveolar_G_s_protein_activation_module (uM).
C CONSTS(23) is Gs_Total in component caveolar_G_s_protein_activation_module (uM).
C STATES(4) is Gi_alpha_GTP in component caveolar_G_i_protein_activation_module (uM).
C STATES(5) is Gi_beta_gamma in component caveolar_G_i_protein_activation_module (uM).
C STATES(6) is Gi_alpha_GDP in component caveolar_G_i_protein_activation_module (uM).
C CONSTS(24) is Gi_Total in component caveolar_G_i_protein_activation_module (uM).
C ALGBRC(21) is R in component extracaveolar_beta_1_adrenergic_receptor_module (uM).
C ALGBRC(22) is LR in component extracaveolar_beta_1_adrenergic_receptor_module (uM).
C ALGBRC(23) is LRG in component extracaveolar_beta_1_adrenergic_receptor_module (uM).
C ALGBRC(24) is RG in component extracaveolar_beta_1_adrenergic_receptor_module (uM).
C CONSTS(25) is R_Total in component extracaveolar_beta_1_adrenergic_receptor_module (uM).
C ALGBRC(8) is Gs_alpha_beta_gamma in component extracaveolar_G_s_protein_activation_module (uM).
C ALGBRC(25) is R in component extracaveolar_muscarinic_receptor_module (uM).
C ALGBRC(26) is LR in component extracaveolar_muscarinic_receptor_module (uM).
C ALGBRC(27) is LRG in component extracaveolar_muscarinic_receptor_module (uM).
C ALGBRC(28) is RG in component extracaveolar_muscarinic_receptor_module (uM).
C CONSTS(26) is R_Total in component extracaveolar_muscarinic_receptor_module (uM).
C ALGBRC(9) is Gi_alpha_beta_gamma in component extracaveolar_G_i_protein_activation_module (uM).
C STATES(7) is Gs_alpha_GTP in component extracaveolar_G_s_protein_activation_module (uM).
C STATES(8) is Gs_beta_gamma in component extracaveolar_G_s_protein_activation_module (uM).
C STATES(9) is Gs_alpha_GDP in component extracaveolar_G_s_protein_activation_module (uM).
C CONSTS(27) is Gs_Total in component extracaveolar_G_s_protein_activation_module (uM).
C STATES(10) is Gi_alpha_GTP in component extracaveolar_G_i_protein_activation_module (uM).
C STATES(11) is Gi_beta_gamma in component extracaveolar_G_i_protein_activation_module (uM).
C STATES(12) is Gi_alpha_GDP in component extracaveolar_G_i_protein_activation_module (uM).
C CONSTS(28) is Gi_Total in component extracaveolar_G_i_protein_activation_module (uM).
C ALGBRC(10) is dcAMP_AC_56_dt in component AC56_module (uM_per_sec).
C ALGBRC(3) is k_AC56 in component AC56_module (per_sec).
C CONSTS(29) is AC_56 in component AC56_module (uM).
C CONSTS(30) is AF56 in component AC56_module (dimensionless).
C CONSTS(31) is MW_AC56 in component AC56_module (kDa).
C CONSTS(32) is ATP in component AC56_module (uM).
C CONSTS(33) is Km_ATP in component AC56_module (uM).
C ALGBRC(11) is dcAMP_AC_47_ecav_dt in component AC47_ecav_module (uM_per_sec).
C ALGBRC(4) is k_AC47_ecav in component AC47_ecav_module (per_sec).
C CONSTS(34) is AC_47_ecav in component AC47_ecav_module (uM).
C CONSTS(35) is AF47 in component AC47_ecav_module (dimensionless).
C CONSTS(36) is MW_AC47 in component AC47_ecav_module (kDa).
C CONSTS(37) is ATP in component AC47_ecav_module (uM).
C CONSTS(38) is Km_ATP in component AC47_ecav_module (uM).
C CONSTS(56) is dcAMP_AC_47_cyt_dt in component AC47_cyt_module (uM_per_sec).
C CONSTS(39) is k_AC47_cyt in component AC47_cyt_module (per_sec).
C CONSTS(40) is AC_47_cyt in component AC47_cyt_module (uM).
C CONSTS(41) is AF47 in component AC47_cyt_module (dimensionless).
C CONSTS(42) is ATP in component AC47_cyt_module (uM).
C CONSTS(43) is Km_ATP in component AC47_cyt_module (uM).
C ALGBRC(29) is dcAMP_cav_PDE2_dt in component caveolar_PDE_module (uM_per_sec).
C ALGBRC(32) is dcAMP_cav_PDE3_dt in component caveolar_PDE_module (uM_per_sec).
C ALGBRC(34) is dcAMP_cav_PDE4_dt in component caveolar_PDE_module (uM_per_sec).
C STATES(13) is cAMP_cav in component cAMP_flux_module (uM).
C CONSTS(44) is PDE2 in component caveolar_PDE_module (uM).
C CONSTS(45) is PDE3 in component caveolar_PDE_module (uM).
C CONSTS(46) is PDE4 in component caveolar_PDE_module (uM).
C ALGBRC(30) is dcAMP_ecav_PDE2_dt in component extracaveolar_PDE_module (uM_per_sec).
C ALGBRC(33) is dcAMP_ecav_PDE4_dt in component extracaveolar_PDE_module (uM_per_sec).
C STATES(14) is cAMP_ecav in component cAMP_flux_module (uM).
C CONSTS(47) is PDE2 in component extracaveolar_PDE_module (uM).
C CONSTS(48) is PDE4 in component extracaveolar_PDE_module (uM).
C ALGBRC(5) is dcAMP_cyt_PDE2_dt in component bulk_cytoplasmic_PDE_module (uM_per_sec).
C ALGBRC(12) is dcAMP_cyt_PDE3_dt in component bulk_cytoplasmic_PDE_module (uM_per_sec).
C ALGBRC(31) is dcAMP_cyt_PDE4_dt in component bulk_cytoplasmic_PDE_module (uM_per_sec).
C STATES(15) is cAMP_cyt in component cAMP_flux_module (uM).
C CONSTS(49) is PDE2 in component bulk_cytoplasmic_PDE_module (uM).
C CONSTS(50) is PDE3 in component bulk_cytoplasmic_PDE_module (uM).
C CONSTS(51) is PDE4 in component bulk_cytoplasmic_PDE_module (uM).
C CONSTS(57) is V_cav in component cAMP_flux_module (liter).
C CONSTS(58) is V_ecav in component cAMP_flux_module (liter).
C CONSTS(59) is V_cyt in component cAMP_flux_module (liter).
C CONSTS(52) is V_cell in component cAMP_flux_module (liter).
C CONSTS(53) is J_cav_ecav in component cAMP_flux_module (liters_per_second).
C CONSTS(54) is J_cav_cyt in component cAMP_flux_module (liters_per_second).
C CONSTS(55) is J_ecav_cyt in component cAMP_flux_module (liters_per_second).
C RATES(1) is d/dt Gs_alpha_GTP in component caveolar_G_s_protein_activation_module (uM).
C RATES(2) is d/dt Gs_beta_gamma in component caveolar_G_s_protein_activation_module (uM).
C RATES(3) is d/dt Gs_alpha_GDP in component caveolar_G_s_protein_activation_module (uM).
C RATES(4) is d/dt Gi_alpha_GTP in component caveolar_G_i_protein_activation_module (uM).
C RATES(5) is d/dt Gi_beta_gamma in component caveolar_G_i_protein_activation_module (uM).
C RATES(6) is d/dt Gi_alpha_GDP in component caveolar_G_i_protein_activation_module (uM).
C RATES(7) is d/dt Gs_alpha_GTP in component extracaveolar_G_s_protein_activation_module (uM).
C RATES(8) is d/dt Gs_beta_gamma in component extracaveolar_G_s_protein_activation_module (uM).
C RATES(9) is d/dt Gs_alpha_GDP in component extracaveolar_G_s_protein_activation_module (uM).
C RATES(10) is d/dt Gi_alpha_GTP in component extracaveolar_G_i_protein_activation_module (uM).
C RATES(11) is d/dt Gi_beta_gamma in component extracaveolar_G_i_protein_activation_module (uM).
C RATES(12) is d/dt Gi_alpha_GDP in component extracaveolar_G_i_protein_activation_module (uM).
C RATES(13) is d/dt cAMP_cav in component cAMP_flux_module (uM).
C RATES(14) is d/dt cAMP_ecav in component cAMP_flux_module (uM).
C RATES(15) is d/dt cAMP_cyt in component cAMP_flux_module (uM).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0.035
      CONSTS(2) = 0.386
      CONSTS(3) = 8.809
      CONSTS(4) = 0.16
      CONSTS(5) = 11
      CONSTS(6) = 30
      CONSTS(7) = 20
      CONSTS(8) = 50
      CONSTS(9) = 1.25
      CONSTS(10) = 0.08
      CONSTS(11) = 2.5
      CONSTS(12) = 2.2
      CONSTS(13) = 5
      CONSTS(14) = 0.1
      CONSTS(15) = 0.8
      CONSTS(16) = 1.21e3
      CONSTS(17) = 2.5
      CONSTS(18) = 0.05
      CONSTS(19) = 0.8
      CONSTS(20) = 1.21e3
      CONSTS(21) = 0.633
      CONSTS(22) = 0.633
      STATES(1) = 0.041983438
      STATES(2) = 0.042634499
      STATES(3) = 0.000651061
      CONSTS(23) = 10
      STATES(4) = 0.012644961
      STATES(5) = 0.013274751
      STATES(6) = 0.00062979
      CONSTS(24) = 20
      CONSTS(25) = 0.633
      CONSTS(26) = 0.633
      STATES(7) = 0.083866891
      STATES(8) = 0.084522918
      STATES(9) = 0.000656025
      CONSTS(27) = 10
      STATES(10) = 0.001018705
      STATES(11) = 0.001475253
      STATES(12) = 0.000456548
      CONSTS(28) = 1
      CONSTS(29) = 3.379
      CONSTS(30) = 500
      CONSTS(31) = 130
      CONSTS(32) = 5000
      CONSTS(33) = 315
      CONSTS(34) = 0.2
      CONSTS(35) = 130
      CONSTS(36) = 130
      CONSTS(37) = 5000
      CONSTS(38) = 315
      CONSTS(39) = 1.08e-3
      CONSTS(40) = 0.136
      CONSTS(41) = 130
      CONSTS(42) = 5000
      CONSTS(43) = 315
      STATES(13) = 0.11750433
      CONSTS(44) = 4.5
      CONSTS(45) = 5.6
      CONSTS(46) = 2
      STATES(14) = 1.092200547
      CONSTS(47) = 0.02
      CONSTS(48) = 0.16
      STATES(15) = 0.992583576
      CONSTS(49) = 5e-3
      CONSTS(50) = 7.5e-3
      CONSTS(51) = 5e-3
      CONSTS(52) = 38e-12
      CONSTS(53) = 7.5e-15
      CONSTS(54) = 7.5e-14
      CONSTS(55) = 1.5e-17
      CONSTS(56) = ( CONSTS(39)*CONSTS(40)*CONSTS(41)*CONSTS(42))/(CONSTS(43)+CONSTS(42))
      CONSTS(57) =  0.0100000*CONSTS(52)
      CONSTS(58) =  0.0200000*CONSTS(52)
      CONSTS(59) =  0.500000*CONSTS(52)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(3) =  STATES(1)*CONSTS(15) -  STATES(3)*STATES(2)*CONSTS(16)
      RATES(6) =  STATES(4)*CONSTS(19) -  STATES(6)*STATES(5)*CONSTS(20)
      RATES(9) =  STATES(7)*CONSTS(15) -  STATES(9)*STATES(8)*CONSTS(16)
      RATES(12) =  STATES(10)*CONSTS(19) -  STATES(12)*STATES(11)*CONSTS(20)
      ALGBRC(1) = TERNRY(VOI.GT.120.000.AND.VOI.LE.720.000, 1.00000, 1.00000)
      ALGBRC(6) = (CONSTS(23) - STATES(1)) - STATES(3)
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
      RATES(1) = ( ALGBRC(16)*CONSTS(14)+ ALGBRC(15)*CONSTS(13)) -  STATES(1)*CONSTS(15)
      RATES(2) = ( ALGBRC(16)*CONSTS(14)+ ALGBRC(15)*CONSTS(13)) -  STATES(3)*STATES(2)*CONSTS(16)
      ALGBRC(2) = TERNRY(VOI.GT.240.000.AND.VOI.LE.540.000, 0.00000, 0.00000)
      ALGBRC(7) = (CONSTS(24) - STATES(4)) - STATES(6)
CALL minimize(minfunc_1, CONSTS, VARIABLES, <VAR>)
      RATES(4) = ( ALGBRC(20)*CONSTS(18)+ ALGBRC(19)*CONSTS(17)) -  STATES(4)*CONSTS(19)
      RATES(5) = ( ALGBRC(20)*CONSTS(18)+ ALGBRC(19)*CONSTS(17)) -  STATES(6)*STATES(5)*CONSTS(20)
      ALGBRC(8) = (CONSTS(27) - STATES(7)) - STATES(9)
CALL minimize(minfunc_2, CONSTS, VARIABLES, <VAR>)
      RATES(7) = ( ALGBRC(24)*CONSTS(14)+ ALGBRC(23)*CONSTS(13)) -  STATES(7)*CONSTS(15)
      RATES(8) = ( ALGBRC(24)*CONSTS(14)+ ALGBRC(23)*CONSTS(13)) -  STATES(9)*STATES(8)*CONSTS(16)
      ALGBRC(9) = (CONSTS(28) - STATES(10)) - STATES(12)
CALL minimize(minfunc_3, CONSTS, VARIABLES, <VAR>)
      RATES(10) = ( ALGBRC(28)*CONSTS(18)+ ALGBRC(27)*CONSTS(17)) -  STATES(10)*CONSTS(19)
      RATES(11) = ( ALGBRC(28)*CONSTS(18)+ ALGBRC(27)*CONSTS(17)) -  STATES(12)*STATES(11)*CONSTS(20)
      ALGBRC(5) = ( CONSTS(7)*CONSTS(49)*STATES(15))/(CONSTS(8)+STATES(15))
      ALGBRC(12) = ( CONSTS(9)*CONSTS(50)*STATES(15))/(CONSTS(10)+STATES(15))
      ALGBRC(31) = ( CONSTS(11)*CONSTS(51)*STATES(15))/(CONSTS(12)+STATES(15))
      RATES(15) = (CONSTS(56) - (ALGBRC(5)+ALGBRC(12)+ALGBRC(31)))+( CONSTS(54)*(STATES(13) - STATES(15)))/CONSTS(59)+( CONSTS(55)*(STATES(14) - STATES(15)))/CONSTS(59)
      ALGBRC(4) =  (( (0.0630000+( 2.01000* STATES(7)*1000.00 ** 1.00430)/(31.5440+ STATES(7)*1000.00 ** 1.00430))*(1.00000+( (1.00000/3.01000)*49.1000* STATES(11)*1000.00 ** 0.892100)/(25.4400+ STATES(11)*1000.00 ** 0.892100))*CONSTS(36))/60.0000)*0.00100000
      ALGBRC(11) = ( ALGBRC(4)*CONSTS(34)*CONSTS(35)*CONSTS(37))/(CONSTS(38)+CONSTS(37))
      ALGBRC(30) = ( CONSTS(7)*CONSTS(47)*STATES(14))/(CONSTS(8)+STATES(14))
      ALGBRC(33) = ( CONSTS(11)*CONSTS(48)*STATES(14))/(CONSTS(12)+STATES(14))
      RATES(14) = ((ALGBRC(11) - (ALGBRC(30)+ALGBRC(33)))+( CONSTS(53)*(STATES(13) - STATES(14)))/CONSTS(58)) - ( CONSTS(55)*(STATES(14) - STATES(15)))/CONSTS(58)
      ALGBRC(3) =  (( (0.700000+( 3.82340*STATES(1)/1.00000 ** 0.978700)/(0.198600+STATES(1)/1.00000 ** 0.978700))*(1.00000+( (1.00000/1.44320)*- 1.00610*STATES(4)/1.00000 ** 0.835600)/(0.191800+STATES(4)/1.00000 ** 0.835600))*CONSTS(31))/60.0000)*0.00100000
      ALGBRC(10) = ( ALGBRC(3)*CONSTS(29)*CONSTS(30)*CONSTS(32))/(CONSTS(33)+CONSTS(32))
      ALGBRC(29) = ( CONSTS(7)*CONSTS(44)*STATES(13))/(CONSTS(8)+STATES(13))
      ALGBRC(32) = ( CONSTS(9)*CONSTS(45)*STATES(13))/(CONSTS(10)+STATES(13))
      ALGBRC(34) = ( CONSTS(11)*CONSTS(46)*STATES(13))/(CONSTS(12)+STATES(13))
      RATES(13) = ((ALGBRC(10) - (ALGBRC(29)+ALGBRC(32)+ALGBRC(34))) - ( CONSTS(53)*(STATES(13) - STATES(14)))/CONSTS(57)) - ( CONSTS(54)*(STATES(13) - STATES(15)))/CONSTS(57)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = TERNRY(VOI.GT.120.000.AND.VOI.LE.720.000, 1.00000, 1.00000)
      ALGBRC(6) = (CONSTS(23) - STATES(1)) - STATES(3)
      ALGBRC(2) = TERNRY(VOI.GT.240.000.AND.VOI.LE.540.000, 0.00000, 0.00000)
      ALGBRC(7) = (CONSTS(24) - STATES(4)) - STATES(6)
      ALGBRC(8) = (CONSTS(27) - STATES(7)) - STATES(9)
      ALGBRC(9) = (CONSTS(28) - STATES(10)) - STATES(12)
      ALGBRC(5) = ( CONSTS(7)*CONSTS(49)*STATES(15))/(CONSTS(8)+STATES(15))
      ALGBRC(12) = ( CONSTS(9)*CONSTS(50)*STATES(15))/(CONSTS(10)+STATES(15))
      ALGBRC(31) = ( CONSTS(11)*CONSTS(51)*STATES(15))/(CONSTS(12)+STATES(15))
      ALGBRC(4) =  (( (0.0630000+( 2.01000* STATES(7)*1000.00 ** 1.00430)/(31.5440+ STATES(7)*1000.00 ** 1.00430))*(1.00000+( (1.00000/3.01000)*49.1000* STATES(11)*1000.00 ** 0.892100)/(25.4400+ STATES(11)*1000.00 ** 0.892100))*CONSTS(36))/60.0000)*0.00100000
      ALGBRC(11) = ( ALGBRC(4)*CONSTS(34)*CONSTS(35)*CONSTS(37))/(CONSTS(38)+CONSTS(37))
      ALGBRC(30) = ( CONSTS(7)*CONSTS(47)*STATES(14))/(CONSTS(8)+STATES(14))
      ALGBRC(33) = ( CONSTS(11)*CONSTS(48)*STATES(14))/(CONSTS(12)+STATES(14))
      ALGBRC(3) =  (( (0.700000+( 3.82340*STATES(1)/1.00000 ** 0.978700)/(0.198600+STATES(1)/1.00000 ** 0.978700))*(1.00000+( (1.00000/1.44320)*- 1.00610*STATES(4)/1.00000 ** 0.835600)/(0.191800+STATES(4)/1.00000 ** 0.835600))*CONSTS(31))/60.0000)*0.00100000
      ALGBRC(10) = ( ALGBRC(3)*CONSTS(29)*CONSTS(30)*CONSTS(32))/(CONSTS(33)+CONSTS(32))
      ALGBRC(29) = ( CONSTS(7)*CONSTS(44)*STATES(13))/(CONSTS(8)+STATES(13))
      ALGBRC(32) = ( CONSTS(9)*CONSTS(45)*STATES(13))/(CONSTS(10)+STATES(13))
      ALGBRC(34) = ( CONSTS(11)*CONSTS(46)*STATES(13))/(CONSTS(12)+STATES(13))
      RETURN
      END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_0[1] = abs(ALGBRC(13) - (((CONSTS(21) - ALGBRC(14)) - ALGBRC(15)) - ALGBRC(16)))
      minfunc_0[2] = abs(ALGBRC(14) - ( ALGBRC(1)*ALGBRC(13))/CONSTS(2))
      minfunc_0[3] = abs(ALGBRC(15) - ( ALGBRC(1)*ALGBRC(13)*ALGBRC(6))/( CONSTS(1)*CONSTS(3)))
      minfunc_0[4] = abs(ALGBRC(16) - ( ALGBRC(13)*ALGBRC(6))/CONSTS(3))
      RETURN
      END
REAL FUNCTION minfunc_1(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_1[1] = abs(ALGBRC(17) - (((CONSTS(22) - ALGBRC(18)) - ALGBRC(19)) - ALGBRC(20)))
      minfunc_1[2] = abs(ALGBRC(18) - ( ALGBRC(2)*ALGBRC(17))/CONSTS(5))
      minfunc_1[3] = abs(ALGBRC(19) - ( ALGBRC(2)*ALGBRC(17)*ALGBRC(7))/( CONSTS(4)*CONSTS(6)))
      minfunc_1[4] = abs(ALGBRC(20) - ( ALGBRC(17)*ALGBRC(7))/CONSTS(6))
      RETURN
      END
REAL FUNCTION minfunc_2(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_2[1] = abs(ALGBRC(21) - (((CONSTS(25) - ALGBRC(22)) - ALGBRC(23)) - ALGBRC(24)))
      minfunc_2[2] = abs(ALGBRC(22) - ( ALGBRC(1)*ALGBRC(21))/CONSTS(2))
      minfunc_2[3] = abs(ALGBRC(23) - ( ALGBRC(1)*ALGBRC(21)*ALGBRC(8))/( CONSTS(1)*CONSTS(3)))
      minfunc_2[4] = abs(ALGBRC(24) - ( ALGBRC(21)*ALGBRC(8))/CONSTS(3))
      RETURN
      END
REAL FUNCTION minfunc_3(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_3[1] = abs(ALGBRC(25) - (((CONSTS(26) - ALGBRC(26)) - ALGBRC(27)) - ALGBRC(28)))
      minfunc_3[2] = abs(ALGBRC(26) - ( ALGBRC(2)*ALGBRC(25))/CONSTS(5))
      minfunc_3[3] = abs(ALGBRC(27) - ( ALGBRC(2)*ALGBRC(25)*ALGBRC(9))/( CONSTS(4)*CONSTS(6)))
      minfunc_3[4] = abs(ALGBRC(28) - ( ALGBRC(25)*ALGBRC(9))/CONSTS(6))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END