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 56 entries in the algebraic variable array.
C There are a total of 4 entries in each of the rate and state variable arrays.
C There are a total of 72 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is ADHMK in component kidney (dimensionless).
C CONSTS(2) is AMK in component kidney (dimensionless).
C CONSTS(3) is AMNA in component kidney (dimensionless).
C CONSTS(4) is ANM in component kidney (dimensionless).
C CONSTS(5) is ANPX in component kidney (dimensionless).
C CONSTS(6) is AUM in component kidney (dimensionless).
C CONSTS(7) is CKE in component kidney (monovalent_mEq_per_litre).
C CONSTS(8) is CNA in component kidney (monovalent_mEq_per_litre).
C CONSTS(9) is HM1 in component kidney (dimensionless).
C CONSTS(10) is MYOGRS in component kidney (dimensionless).
C CONSTS(11) is PA in component kidney (mmHg).
C CONSTS(12) is PAMKRN in component kidney (dimensionless).
C CONSTS(13) is PPC in component kidney (mmHg).
C CONSTS(14) is VTW in component kidney (litre).
C ALGBRC(1) is PAR in component perfusion_pressure (mmHg).
C CONSTS(15) is GBL in component parameter_values (mmHg).
C CONSTS(16) is RAPRSP in component parameter_values (mmHg).
C CONSTS(17) is RFCDFT in component parameter_values (dimensionless).
C CONSTS(18) is RCDFPC in component parameter_values (dimensionless).
C CONSTS(19) is RCDFDP in component parameter_values (minute).
C STATES(1) is PAR1 in component perfusion_pressure (mmHg).
C ALGBRC(3) is MDFLW in component proximal_tubular_and_macula_densa_flow (L_per_minute).
C ALGBRC(4) is RNAUG2 in component renal_autoregulatory_feedback_factor (dimensionless).
C CONSTS(20) is RNAUGN in component parameter_values (minute_per_L).
C CONSTS(21) is RNAULL in component parameter_values (dimensionless).
C CONSTS(22) is RNAUUL in component parameter_values (dimensionless).
C CONSTS(23) is RNAUAD in component parameter_values (per_minute).
C ALGBRC(5) is RNAUG1 in component renal_autoregulatory_feedback_factor (dimensionless).
C ALGBRC(6) is RNAUG1T in component renal_autoregulatory_feedback_factor (dimensionless).
C STATES(2) is RNAUG3 in component renal_autoregulatory_feedback_factor (dimensionless).
C CONSTS(64) is AUMK in component autonomic_effect_on_AAR (dimensionless).
C CONSTS(24) is ARF in component parameter_values (dimensionless).
C CONSTS(63) is AUMKT in component autonomic_effect_on_AAR (dimensionless).
C CONSTS(66) is ANMAR in component angiotensin_effect_on_AAR (dimensionless).
C CONSTS(25) is ANMAM in component parameter_values (dimensionless).
C CONSTS(26) is ANMARL in component parameter_values (dimensionless).
C CONSTS(65) is ANMAR1 in component angiotensin_effect_on_AAR (dimensionless).
C ALGBRC(7) is AAR1 in component AAR_calculation (mmHg_minute_per_L).
C CONSTS(27) is AARK in component parameter_values (mmHg_minute_per_L).
C ALGBRC(8) is AAR in component atrial_natriuretic_peptide_effect_on_AAR (mmHg_minute_per_L).
C CONSTS(28) is ANPXAF in component parameter_values (mmHg_minute_per_L).
C CONSTS(29) is AARLL in component parameter_values (mmHg_minute_per_L).
C ALGBRC(9) is AART in component atrial_natriuretic_peptide_effect_on_AAR (mmHg_minute_per_L).
C CONSTS(67) is AUMK2 in component autonomic_effect_on_EAR (dimensionless).
C CONSTS(30) is AUMK1 in component parameter_values (dimensionless).
C CONSTS(68) is ANMER in component angiotensin_effect_on_EAR (dimensionless).
C CONSTS(31) is ANMEM in component parameter_values (dimensionless).
C ALGBRC(10) is RNAUG4 in component effect_of_renal_autoregulatory_feedback_on_EAR (dimensionless).
C CONSTS(32) is EFAFR in component parameter_values (dimensionless).
C ALGBRC(11) is EAR in component EAR_calculation (mmHg_minute_per_L).
C CONSTS(33) is EARK in component parameter_values (mmHg_minute_per_L).
C CONSTS(34) is EARLL in component parameter_values (mmHg_minute_per_L).
C ALGBRC(12) is EAR1 in component EAR_calculation (mmHg_minute_per_L).
C ALGBRC(13) is RR in component total_renal_resistance (mmHg_minute_per_L).
C ALGBRC(14) is RFN in component normal_renal_blood_flow (L_per_minute).
C ALGBRC(25) is RBF in component actual_renal_blood_flow (L_per_minute).
C CONSTS(35) is REK in component parameter_values (dimensionless).
C ALGBRC(15) is GFN in component glomerular_filtration_rate (L_per_minute).
C ALGBRC(16) is GLPC in component glomerular_colloid_osmotic_pressure (mmHg).
C CONSTS(36) is GPPD in component parameter_values (dimensionless).
C CONSTS(37) is GLPCA in component parameter_values (mmHg).
C ALGBRC(17) is EFAFPR in component glomerular_colloid_osmotic_pressure (dimensionless).
C ALGBRC(18) is EFAFPR1 in component glomerular_colloid_osmotic_pressure (dimensionless).
C ALGBRC(19) is GLP in component glomerular_pressure (mmHg).
C ALGBRC(20) is APD in component glomerular_pressure (mmHg).
C ALGBRC(26) is GFR in component glomerular_filtration_rate (L_per_minute).
C CONSTS(38) is PXTP in component parameter_values (mmHg).
C CONSTS(39) is GFLC in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(40) is GFNLL in component parameter_values (L_per_minute).
C ALGBRC(21) is PFL in component glomerular_filtration_rate (mmHg).
C ALGBRC(22) is GFN1 in component glomerular_filtration_rate (L_per_minute).
C CONSTS(41) is MDFL1 in component parameter_values (dimensionless).
C ALGBRC(23) is PTFL in component proximal_tubular_and_macula_densa_flow (L_per_minute).
C ALGBRC(24) is MDFLWT in component proximal_tubular_and_macula_densa_flow (L_per_minute).
C ALGBRC(28) is RTSPPC in component renal_tissue_osmotic_pressure (mmHg).
C CONSTS(42) is RTPPR in component parameter_values (dimensionless).
C CONSTS(43) is RTPPRS in component parameter_values (mmHg).
C ALGBRC(27) is RTSPPC1 in component renal_tissue_osmotic_pressure (mmHg).
C ALGBRC(50) is UROD in component actual_urea_excretion_rate (mOsm_per_minute).
C STATES(3) is PLUR in component glomerular_urea_concentration (mOsm).
C CONSTS(44) is URFORM in component parameter_values (mOsm_per_minute).
C ALGBRC(2) is PLURC in component plasma_urea_concentration (mOsm_per_litre).
C ALGBRC(29) is RCPRS in component peritubular_capillary_pressure (mmHg).
C CONSTS(45) is RFABX in component parameter_values (dimensionless).
C CONSTS(46) is RVRS in component parameter_values (mmHg_minute_per_L).
C ALGBRC(34) is RFABD in component peritubular_capillary_reabsorption_factor (dimensionless).
C CONSTS(47) is RTSPRS in component parameter_values (mmHg).
C CONSTS(48) is RABSC in component parameter_values (per_mmHg).
C CONSTS(49) is RFABDP in component parameter_values (dimensionless).
C CONSTS(50) is RFABDM in component parameter_values (dimensionless).
C ALGBRC(30) is RABSPR in component peritubular_capillary_reabsorption_factor (mmHg).
C ALGBRC(31) is RFAB1 in component peritubular_capillary_reabsorption_factor (dimensionless).
C ALGBRC(32) is RFAB in component peritubular_capillary_reabsorption_factor (dimensionless).
C ALGBRC(33) is RFABD1 in component peritubular_capillary_reabsorption_factor (dimensionless).
C ALGBRC(35) is DTNAI in component distal_tubular_Na_delivery (monovalent_mEq_per_minute).
C ALGBRC(37) is DTNARA in component Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(51) is DTNAR in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(52) is DIURET in component parameter_values (dimensionless).
C CONSTS(53) is AHMNAR in component parameter_values (dimensionless).
C CONSTS(54) is DTNARL in component parameter_values (monovalent_mEq_per_minute).
C ALGBRC(36) is DTNARA1 in component Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(70) is DTNANG in component angiotensin_induced_Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(55) is ANMNAM in component parameter_values (dimensionless).
C CONSTS(69) is DTNANG1 in component angiotensin_induced_Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C ALGBRC(38) is DTKI in component distal_tubular_K_delivery (monovalent_mEq_per_minute).
C ALGBRC(39) is RFABK in component effect_of_physical_forces_on_distal_K_reabsorption (monovalent_mEq_per_minute).
C CONSTS(56) is RFABKM in component parameter_values (monovalent_mEq_per_minute).
C ALGBRC(41) is MDFLK in component effect_of_fluid_flow_on_distal_K_reabsorption (monovalent_mEq_per_minute).
C CONSTS(57) is MDFLKM in component parameter_values (monovalent_mEq_per_litre).
C ALGBRC(40) is MDFLK1 in component effect_of_fluid_flow_on_distal_K_reabsorption (monovalent_mEq_per_minute).
C ALGBRC(47) is KODN in component normal_K_excretion (monovalent_mEq_per_minute).
C ALGBRC(55) is VUDN in component normal_urine_volume (L_per_minute).
C STATES(4) is DTKA in component K_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C ALGBRC(42) is DTKSC in component K_secretion_from_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(58) is ANMKEM in component parameter_values (dimensionless).
C CONSTS(59) is ANMKEL in component parameter_values (dimensionless).
C CONSTS(60) is CKEEX in component parameter_values (dimensionless).
C CONSTS(71) is ANMKE1 in component K_secretion_from_distal_tubules (dimensionless).
C CONSTS(72) is ANMKE in component K_secretion_from_distal_tubules (dimensionless).
C ALGBRC(44) is NODN in component normal_Na_excretion (monovalent_mEq_per_minute).
C ALGBRC(43) is NODN1 in component normal_Na_excretion (monovalent_mEq_per_minute).
C ALGBRC(45) is KODN1 in component normal_K_excretion (monovalent_mEq_per_minute).
C ALGBRC(48) is DTURI in component normal_urea_excretion (mOsm_per_minute).
C ALGBRC(51) is OSMOPN1 in component normal_osmolar_and_water_excretion (mOsm_per_minute).
C ALGBRC(52) is OSMOPN in component normal_osmolar_and_water_excretion (mOsm_per_minute).
C ALGBRC(53) is OSMOP1T in component normal_urine_volume (mOsm_per_minute).
C ALGBRC(54) is OSMOP1 in component normal_urine_volume (mOsm_per_minute).
C ALGBRC(46) is NOD in component actual_Na_excretion_rate (monovalent_mEq_per_minute).
C ALGBRC(49) is KOD in component actual_K_excretion_rate (monovalent_mEq_per_minute).
C ALGBRC(56) is VUD in component actual_urine_volume (L_per_minute).
C CONSTS(61) is RNAGTC in component parameter_values (minute).
C CONSTS(62) is GFNDMP in component parameter_values (dimensionless).
C RATES(1) is d/dt PAR1 in component perfusion_pressure (mmHg).
C RATES(2) is d/dt RNAUG3 in component renal_autoregulatory_feedback_factor (dimensionless).
C RATES(3) is d/dt PLUR in component glomerular_urea_concentration (mOsm).
C RATES(4) is d/dt DTKA in component K_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1.0
      CONSTS(2) = 1.037
      CONSTS(3) = 1.0
      CONSTS(4) = 0.987545
      CONSTS(5) = 1.0
      CONSTS(6) = 1.00066
      CONSTS(7) = 4.44092
      CONSTS(8) = 142.035
      CONSTS(9) = 0.39984739
      CONSTS(10) = 1.0
      CONSTS(11) = 103.525
      CONSTS(12) = 1.0
      CONSTS(13) = 29.9941
      CONSTS(14) = 39.8952
      CONSTS(15) = 0
      CONSTS(16) = 0
      CONSTS(17) = 0
      CONSTS(18) = 0
      CONSTS(19) = 2000
      STATES(1) = 103.525
      CONSTS(20) = 0.6
      CONSTS(21) = 0.3
      CONSTS(22) = 10
      CONSTS(23) = 0
      STATES(2) = 0.0
      CONSTS(24) = 0.5
      CONSTS(25) = 1.4
      CONSTS(26) = 0.86
      CONSTS(27) = 1
      CONSTS(28) = 1.5
      CONSTS(29) = 4
      CONSTS(30) = 0.3
      CONSTS(31) = 1.6
      CONSTS(32) = 0
      CONSTS(33) = 1
      CONSTS(34) = 24
      CONSTS(35) = 1
      CONSTS(36) = 1.0
      CONSTS(37) = 1.0
      CONSTS(38) = 8
      CONSTS(39) = 0.0208333
      CONSTS(40) = 0.001
      CONSTS(41) = 10
      CONSTS(42) = 0.9
      CONSTS(43) = 15.2
      STATES(3) = 159.549
      CONSTS(44) = 0.24
      CONSTS(45) = 0.8
      CONSTS(46) = 19.167
      CONSTS(47) = 6
      CONSTS(48) = 0.5
      CONSTS(49) = 1
      CONSTS(50) = 0.3
      CONSTS(51) = 0.675
      CONSTS(52) = 1
      CONSTS(53) = 0.3
      CONSTS(54) = 1e-06
      CONSTS(55) = 1
      CONSTS(56) = 0.03
      CONSTS(57) = 0.667
      STATES(4) = 0.0367573
      CONSTS(58) = 2
      CONSTS(59) = 0.3
      CONSTS(60) = 4
      CONSTS(61) = 15
      CONSTS(62) = 3
      CONSTS(63) =  (CONSTS(6) - 1.00000)*CONSTS(24)+1.00000
      CONSTS(64) = TERNRY(CONSTS(63).LT.0.800000, 0.800000, CONSTS(63))
      CONSTS(65) =  (CONSTS(4) - 1.00000)*CONSTS(25)+1.00000
      CONSTS(66) = TERNRY(CONSTS(65).LT.CONSTS(26), CONSTS(26), CONSTS(65))
      CONSTS(67) =  (CONSTS(64) - 1.00000)*CONSTS(30)+1.00000
      CONSTS(68) =  (CONSTS(4) - 1.00000)*CONSTS(31)+1.00000
      CONSTS(69) =  ( (CONSTS(4) - 1.00000)*CONSTS(55)+1.00000)*0.100000
      CONSTS(70) = TERNRY(CONSTS(69).LT.0.00000, 0.00000, CONSTS(69))
      CONSTS(71) =  (CONSTS(4) - 1.00000)*CONSTS(58)+1.00000
      CONSTS(72) = TERNRY(CONSTS(71).LT.CONSTS(59), CONSTS(59), CONSTS(71))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = ((100.000+ (CONSTS(11) - 100.000)*CONSTS(18)) - STATES(1))/CONSTS(19)
      ALGBRC(1) = TERNRY(CONSTS(16).GT.0.00000.AND.CONSTS(17).LE.0.00000, CONSTS(16), TERNRY(CONSTS(17).GT.0.00000, STATES(1), CONSTS(11) - CONSTS(15))
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
      RATES(2) =  (ALGBRC(4) - 1.00000)*CONSTS(23)
      ALGBRC(2) = STATES(3)/CONSTS(14)
      ALGBRC(48) =  ALGBRC(15) ** 2.00000*ALGBRC(2)*3.84000
      ALGBRC(50) =  ALGBRC(48)*CONSTS(35)
      RATES(3) = CONSTS(44) - ALGBRC(50)
      ALGBRC(35) =  ALGBRC(3)*CONSTS(8)*0.00616190
      ALGBRC(38) = ( ALGBRC(35)*CONSTS(7))/CONSTS(8)
      ALGBRC(27) =  ALGBRC(16)*CONSTS(42) - CONSTS(43)
      ALGBRC(28) = TERNRY(ALGBRC(27).LT.1.00000, 1.00000, ALGBRC(27))
      ALGBRC(29) =  ( (ALGBRC(14) - 1.20000)*CONSTS(45)+1.20000)*CONSTS(46)
      ALGBRC(30) = ((ALGBRC(16)+CONSTS(47)) - ALGBRC(29)) - ALGBRC(28)
      ALGBRC(31) =  ALGBRC(30)*CONSTS(48)
      ALGBRC(32) = ALGBRC(31)
      ALGBRC(33) =  (ALGBRC(32) - 1.00000)*CONSTS(50)+1.00000
      ALGBRC(34) = TERNRY(ALGBRC(33).LT.0.000100000, 0.000100000, ALGBRC(33))
      ALGBRC(39) =  (ALGBRC(34) - 1.00000)*CONSTS(56)
      ALGBRC(40) =  (ALGBRC(3) - 1.00000)*CONSTS(57)+1.00000
      ALGBRC(41) = TERNRY(ALGBRC(40).LT.0.100000, 0.100000, ALGBRC(40))
      ALGBRC(42) = ( CONSTS(7)/4.40000 ** CONSTS(60)*CONSTS(2)*0.0800000*ALGBRC(41))/CONSTS(72)
      ALGBRC(45) = ((ALGBRC(38)+ALGBRC(42)) - STATES(4)) - ALGBRC(39)
      ALGBRC(47) = TERNRY(ALGBRC(45).LT.0.00000, 0.00000, ALGBRC(45))
      ALGBRC(36) =  (( CONSTS(3)*ALGBRC(34)*CONSTS(51))/CONSTS(52))*( (CONSTS(1) - 1.00000)*CONSTS(53)+1.00000)
      ALGBRC(37) = TERNRY(ALGBRC(36).LT.CONSTS(54), CONSTS(54), ALGBRC(36))
      ALGBRC(43) = (ALGBRC(35) - ALGBRC(37)) - CONSTS(70)
      ALGBRC(44) = TERNRY(ALGBRC(43).LT.1.00000e-08, 1.00000e-08, ALGBRC(43))
      ALGBRC(51) = ALGBRC(48)+ 2.00000*(ALGBRC(44)+ALGBRC(47))
      ALGBRC(52) = TERNRY(ALGBRC(51).GT.0.600000, 0.600000, ALGBRC(51))
      ALGBRC(53) = ALGBRC(51) - 0.600000
      ALGBRC(54) = TERNRY(ALGBRC(53).LT.0.00000, 0.00000, ALGBRC(53))
      ALGBRC(55) = ALGBRC(52)/( 600.000*CONSTS(1))+ALGBRC(54)/360.000
      RATES(4) =  ( (ALGBRC(47)/ALGBRC(55))*0.000451800 - STATES(4))*1.00000
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = TERNRY(CONSTS(16).GT.0.00000.AND.CONSTS(17).LE.0.00000, CONSTS(16), TERNRY(CONSTS(17).GT.0.00000, STATES(1), CONSTS(11) - CONSTS(15))
      ALGBRC(2) = STATES(3)/CONSTS(14)
      ALGBRC(48) =  ALGBRC(15) ** 2.00000*ALGBRC(2)*3.84000
      ALGBRC(50) =  ALGBRC(48)*CONSTS(35)
      ALGBRC(35) =  ALGBRC(3)*CONSTS(8)*0.00616190
      ALGBRC(38) = ( ALGBRC(35)*CONSTS(7))/CONSTS(8)
      ALGBRC(27) =  ALGBRC(16)*CONSTS(42) - CONSTS(43)
      ALGBRC(28) = TERNRY(ALGBRC(27).LT.1.00000, 1.00000, ALGBRC(27))
      ALGBRC(29) =  ( (ALGBRC(14) - 1.20000)*CONSTS(45)+1.20000)*CONSTS(46)
      ALGBRC(30) = ((ALGBRC(16)+CONSTS(47)) - ALGBRC(29)) - ALGBRC(28)
      ALGBRC(31) =  ALGBRC(30)*CONSTS(48)
      ALGBRC(32) = ALGBRC(31)
      ALGBRC(33) =  (ALGBRC(32) - 1.00000)*CONSTS(50)+1.00000
      ALGBRC(34) = TERNRY(ALGBRC(33).LT.0.000100000, 0.000100000, ALGBRC(33))
      ALGBRC(39) =  (ALGBRC(34) - 1.00000)*CONSTS(56)
      ALGBRC(40) =  (ALGBRC(3) - 1.00000)*CONSTS(57)+1.00000
      ALGBRC(41) = TERNRY(ALGBRC(40).LT.0.100000, 0.100000, ALGBRC(40))
      ALGBRC(42) = ( CONSTS(7)/4.40000 ** CONSTS(60)*CONSTS(2)*0.0800000*ALGBRC(41))/CONSTS(72)
      ALGBRC(45) = ((ALGBRC(38)+ALGBRC(42)) - STATES(4)) - ALGBRC(39)
      ALGBRC(47) = TERNRY(ALGBRC(45).LT.0.00000, 0.00000, ALGBRC(45))
      ALGBRC(36) =  (( CONSTS(3)*ALGBRC(34)*CONSTS(51))/CONSTS(52))*( (CONSTS(1) - 1.00000)*CONSTS(53)+1.00000)
      ALGBRC(37) = TERNRY(ALGBRC(36).LT.CONSTS(54), CONSTS(54), ALGBRC(36))
      ALGBRC(43) = (ALGBRC(35) - ALGBRC(37)) - CONSTS(70)
      ALGBRC(44) = TERNRY(ALGBRC(43).LT.1.00000e-08, 1.00000e-08, ALGBRC(43))
      ALGBRC(51) = ALGBRC(48)+ 2.00000*(ALGBRC(44)+ALGBRC(47))
      ALGBRC(52) = TERNRY(ALGBRC(51).GT.0.600000, 0.600000, ALGBRC(51))
      ALGBRC(53) = ALGBRC(51) - 0.600000
      ALGBRC(54) = TERNRY(ALGBRC(53).LT.0.00000, 0.00000, ALGBRC(53))
      ALGBRC(55) = ALGBRC(52)/( 600.000*CONSTS(1))+ALGBRC(54)/360.000
      ALGBRC(25) =  CONSTS(35)*ALGBRC(14)
      ALGBRC(26) =  ALGBRC(15)*CONSTS(35)
      ALGBRC(46) =  ALGBRC(44)*CONSTS(35)
      ALGBRC(49) =  ALGBRC(47)*CONSTS(35)
      ALGBRC(56) =  ALGBRC(55)*CONSTS(35)
      RETURN
      END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_0[1] = abs(ALGBRC(6) - ( (ALGBRC(3) - 1.00000)*CONSTS(20)+1.00000))
      minfunc_0[2] = abs(ALGBRC(5) - (TERNRY(ALGBRC(6).LT.CONSTS(21), CONSTS(21), TERNRY(ALGBRC(6).GT.CONSTS(22), CONSTS(22), ALGBRC(6))))
      minfunc_0[3] = abs(ALGBRC(4) - (ALGBRC(5) - STATES(2)))
      minfunc_0[4] = abs(ALGBRC(7) -  CONSTS(27)*CONSTS(12)*CONSTS(64)*ALGBRC(4)*CONSTS(66)*40.0000*CONSTS(10))
      minfunc_0[5] = abs(ALGBRC(9) - ((ALGBRC(7) -  CONSTS(5)*CONSTS(28))+CONSTS(28)))
      minfunc_0[6] = abs(ALGBRC(8) - (TERNRY(ALGBRC(9).LT.CONSTS(29), CONSTS(29), ALGBRC(9))))
      minfunc_0[7] = abs(ALGBRC(10) - ( (ALGBRC(4) - 1.00000)*CONSTS(32)+1.00000))
      minfunc_0[8] = abs(ALGBRC(12) -  43.3330*CONSTS(33)*CONSTS(68)*ALGBRC(10)*CONSTS(10)*CONSTS(67))
      minfunc_0[9] = abs(ALGBRC(11) - (TERNRY(ALGBRC(12).LT.CONSTS(34), CONSTS(34), ALGBRC(12))))
      minfunc_0[10] = abs(ALGBRC(13) - (ALGBRC(8)+ALGBRC(11)))
      minfunc_0[11] = abs(ALGBRC(14) - ALGBRC(1)/ALGBRC(13))
      minfunc_0[12] = abs(ALGBRC(18) - ( ALGBRC(14)*(1.00000 - CONSTS(9)))/( ALGBRC(14)*(1.00000 - CONSTS(9)) - ALGBRC(15)))
      minfunc_0[13] = abs(ALGBRC(17) - (TERNRY(ALGBRC(18).LT.1.00000, 1.00000, ALGBRC(18))))
      minfunc_0[14] = abs(ALGBRC(16) - (TERNRY(CONSTS(37).GT.0.00000,  ALGBRC(17) ** 1.35000*CONSTS(13)*0.980000, CONSTS(13)+4.00000)))
      minfunc_0[15] = abs(ALGBRC(20) -  ALGBRC(8)*ALGBRC(14))
      minfunc_0[16] = abs(ALGBRC(19) - (ALGBRC(1) - ALGBRC(20)))
      minfunc_0[17] = abs(ALGBRC(21) - ((ALGBRC(19) - ALGBRC(16)) - CONSTS(38)))
      minfunc_0[18] = abs(ALGBRC(22) -  ALGBRC(21)*CONSTS(39))
      minfunc_0[19] = abs(ALGBRC(15) - (TERNRY(ALGBRC(22).LT.CONSTS(40), CONSTS(40), ALGBRC(22))))
      minfunc_0[20] = abs(ALGBRC(23) -  ALGBRC(15)*8.00000)
      minfunc_0[21] = abs(ALGBRC(24) - ( (ALGBRC(23) - 1.00000)*CONSTS(41)+1.00000))
      minfunc_0[22] = abs(ALGBRC(3) - (TERNRY(ALGBRC(24).LT.0.00000, 0.00000, ALGBRC(24))))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END