/*
   There are a total of 325 entries in the algebraic variable array.
   There are a total of 43 entries in each of the rate and state variable arrays.
   There are a total of 210 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * CONSTANTS[0] is MYOGRS in component temp_myogrs_and_pamk (dimensionless).
 * CONSTANTS[1] is PAMK in component temp_myogrs_and_pamk (dimensionless).
 * CONSTANTS[2] is PAMKRN in component temp_myogrs_and_pamk (dimensionless).
 * ALGEBRAIC[6] is ANM in component general_angiotensin_multiplier (dimensionless).
 * ALGEBRAIC[84] is CKE in component extracellular_K_concentration (monovalent_mEq_per_litre).
 * ALGEBRAIC[3] is AMK in component aldosterone_effect_on_cell_membrane_K_transport (dimensionless).
 * ALGEBRAIC[5] is AMNA in component aldosterone_effect_on_cell_membrane_Na_transport (dimensionless).
 * ALGEBRAIC[7] is ANMAL in component angiotensin_control_of_aldosterone_secretion (dimensionless).
 * CONSTANTS[3] is ANMALD in component aldosterone_parameter_values (dimensionless).
 * ALGEBRAIC[86] is OSMAL in component osmotic_control_of_aldosterone_secretion (dimensionless).
 * ALGEBRAIC[98] is AMR1 in component aldosterone_secretion (dimensionless).
 * CONSTANTS[4] is AMKMUL in component aldosterone_parameter_values (dimensionless).
 * CONSTANTS[5] is ALDINF in component aldosterone_parameter_values (dimensionless).
 * CONSTANTS[6] is ALDKNS in component aldosterone_parameter_values (dimensionless).
 * ALGEBRAIC[90] is AMRBSC in component aldosterone_secretion (dimensionless).
 * ALGEBRAIC[93] is AMRT in component aldosterone_secretion (dimensionless).
 * ALGEBRAIC[96] is AMR in component aldosterone_secretion (dimensionless).
 * STATES[0] is AMC in component aldosterone_concentration (dimensionless).
 * CONSTANTS[7] is AMT in component aldosterone_parameter_values (minute).
 * ALGEBRAIC[1] is AM in component general_aldosterone_multiplier (dimensionless).
 * CONSTANTS[8] is AM1UL in component aldosterone_parameter_values (dimensionless).
 * CONSTANTS[9] is AM1LL in component aldosterone_parameter_values (dimensionless).
 * CONSTANTS[10] is AMCSNS in component aldosterone_parameter_values (dimensionless).
 * CONSTANTS[11] is ALDMM in component aldosterone_parameter_values (dimensionless).
 * ALGEBRAIC[0] is AM1 in component general_aldosterone_multiplier (dimensionless).
 * CONSTANTS[12] is AMKM in component aldosterone_parameter_values (dimensionless).
 * ALGEBRAIC[2] is AMKT in component aldosterone_effect_on_cell_membrane_K_transport (dimensionless).
 * CONSTANTS[13] is AMNAM in component aldosterone_parameter_values (dimensionless).
 * CONSTANTS[14] is AMNAUL in component aldosterone_parameter_values (dimensionless).
 * CONSTANTS[15] is AMNALL in component aldosterone_parameter_values (dimensionless).
 * ALGEBRAIC[4] is AMNAT in component aldosterone_effect_on_cell_membrane_Na_transport (dimensionless).
 * ALGEBRAIC[232] is MDFLW in component proximal_tubular_and_macula_densa_flow (L_per_minute).
 * ALGEBRAIC[10] is ANUVN in component angiotensin_effect_on_venous_constriction (dimensionless).
 * ALGEBRAIC[9] is ANU in component angiotensin_effect_on_circulation (dimensionless).
 * ALGEBRAIC[268] is ANGSCR in component instantaneous_angiotensin_formation (dimensionless).
 * ALGEBRAIC[262] is MDFLW3 in component instantaneous_angiotensin_formation (L_per_minute).
 * STATES[1] is ANX1 in component time_delayed_angiotensin_formation (dimensionless).
 * CONSTANTS[16] is ANXM in component angiotensin_parameter_values (dimensionless).
 * CONSTANTS[17] is ANV in component angiotensin_parameter_values (minute).
 * ALGEBRAIC[275] is ANX in component time_delayed_angiotensin_formation (dimensionless).
 * ALGEBRAIC[283] is ANPR in component total_angiotensin_formation (dimensionless).
 * CONSTANTS[18] is REK in component angiotensin_parameter_values (dimensionless).
 * ALGEBRAIC[276] is ANPRT in component total_angiotensin_formation (dimensionless).
 * ALGEBRAIC[293] is ANPR1 in component artificial_angiotensin_formation (dimensionless).
 * CONSTANTS[19] is ANGKNS in component angiotensin_parameter_values (dimensionless).
 * CONSTANTS[20] is ANGINF in component angiotensin_parameter_values (dimensionless).
 * STATES[2] is ANC in component angiotensin_concentration (dimensionless).
 * CONSTANTS[21] is ANT in component angiotensin_parameter_values (minute).
 * CONSTANTS[22] is ANMUL in component angiotensin_parameter_values (dimensionless).
 * CONSTANTS[23] is ANMLL in component angiotensin_parameter_values (dimensionless).
 * CONSTANTS[24] is ANCSNS in component angiotensin_parameter_values (dimensionless).
 * CONSTANTS[25] is ANUM in component angiotensin_parameter_values (dimensionless).
 * CONSTANTS[26] is ANULL in component angiotensin_parameter_values (dimensionless).
 * ALGEBRAIC[8] is ANU1 in component angiotensin_effect_on_circulation (dimensionless).
 * CONSTANTS[27] is ANUVM in component angiotensin_parameter_values (dimensionless).
 * CONSTANTS[28] is Z12 in component angiotensin_parameter_values (dimensionless).
 * ALGEBRAIC[82] is CNA in component extracellular_Na_concentration (monovalent_mEq_per_litre).
 * ALGEBRAIC[65] is PA1 in component pressure_driving_autonomic_receptors (mmHg).
 * ALGEBRAIC[14] is ADHMK in component ADH_effect_on_kidney (dimensionless).
 * ALGEBRAIC[12] is ADHMV in component ADH_effect_on_nonrenal_vascular_resistance (dimensionless).
 * STATES[3] is ADHC in component ADH_in_blood (dimensionless).
 * ALGEBRAIC[87] is ADHNA in component osmotic_control_of_ADH_secretion (dimensionless).
 * CONSTANTS[29] is CNR in component antidiuretic_hormone_parameter_values (monovalent_mEq_per_litre).
 * ALGEBRAIC[83] is ADHNA1 in component osmotic_control_of_ADH_secretion (dimensionless).
 * ALGEBRAIC[69] is ADHPR in component pressure_control_of_ADH_secretion (dimensionless).
 * CONSTANTS[30] is ADHPUL in component antidiuretic_hormone_parameter_values (mmHg).
 * CONSTANTS[31] is ADHPAM in component antidiuretic_hormone_parameter_values (per_mmHg2).
 * ALGEBRAIC[67] is ADHPA in component pressure_control_of_ADH_secretion (mmHg).
 * ALGEBRAIC[94] is ADH in component total_ADH_secretion (dimensionless).
 * CONSTANTS[32] is ADHINF in component antidiuretic_hormone_parameter_values (dimensionless).
 * ALGEBRAIC[91] is ADH1 in component total_ADH_secretion (dimensionless).
 * CONSTANTS[33] is ADHTC in component antidiuretic_hormone_parameter_values (minute).
 * CONSTANTS[34] is ADHVUL in component antidiuretic_hormone_parameter_values (dimensionless).
 * CONSTANTS[35] is ADHVLL in component antidiuretic_hormone_parameter_values (dimensionless).
 * ALGEBRAIC[11] is ADHMV1 in component ADH_effect_on_nonrenal_vascular_resistance (dimensionless).
 * CONSTANTS[36] is ADHKLL in component antidiuretic_hormone_parameter_values (dimensionless).
 * CONSTANTS[37] is ADHKUL in component antidiuretic_hormone_parameter_values (dimensionless).
 * ALGEBRAIC[13] is ADHMK1 in component ADH_effect_on_kidney (dimensionless).
 * ALGEBRAIC[47] is PLA in component left_atrial_pressure (mmHg).
 * ALGEBRAIC[34] is PRA in component right_atrial_pressure (mmHg).
 * ALGEBRAIC[16] is ANPX in component ANP_effect_on_renal_afferent_arteriolar_resistance (dimensionless).
 * ALGEBRAIC[50] is ANP in component total_ANP_secreted (dimensionless).
 * ALGEBRAIC[48] is ANPL in component total_ANP_secreted (dimensionless).
 * ALGEBRAIC[35] is ANPR2 in component total_ANP_secreted (dimensionless).
 * ALGEBRAIC[52] is ANP1 in component ANP_into_circulation (dimensionless).
 * CONSTANTS[38] is ANPKNS in component atrial_natriuretic_peptide_parameter_values (dimensionless).
 * CONSTANTS[39] is ANPINF in component atrial_natriuretic_peptide_parameter_values (dimensionless).
 * STATES[4] is ANPC in component ANP_in_plasma (dimensionless).
 * CONSTANTS[40] is ANPTC in component atrial_natriuretic_peptide_parameter_values (minute).
 * CONSTANTS[41] is ANPXUL in component atrial_natriuretic_peptide_parameter_values (dimensionless).
 * ALGEBRAIC[15] is ANPX1 in component ANP_effect_on_renal_afferent_arteriolar_resistance (dimensionless).
 * ALGEBRAIC[64] is PA in component arterial_pressure_and_pressure_gradient (mmHg).
 * ALGEBRAIC[169] is PO2ART in component arterial_PO2 (mmHg).
 * ALGEBRAIC[40] is PPA in component pulmonary_vasculature_pressure (mmHg).
 * ALGEBRAIC[20] is VVR in component autonomic_drive_on_target_organs_and_tissues (litre).
 * ALGEBRAIC[22] is AUH in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * ALGEBRAIC[23] is AUR in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * ALGEBRAIC[25] is AUM in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * ALGEBRAIC[24] is AOM in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * ALGEBRAIC[26] is AVE in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * ALGEBRAIC[19] is AU in component actual_autonomic_stimulation (dimensionless).
 * CONSTANTS[42] is CRRFLX in component autonomics_parameter_values (mmHg).
 * CONSTANTS[43] is EXE in component autonomics_parameter_values (mmHg).
 * ALGEBRAIC[71] is AUC in component chemoreceptors_effect_of_PA (dimensionless).
 * CONSTANTS[44] is AUC1 in component autonomics_parameter_values (dimensionless).
 * ALGEBRAIC[173] is AUC3 in component chemoreceptors_effect_of_art_PO2 (dimensionless).
 * CONSTANTS[45] is O2CHMO in component autonomics_parameter_values (per_mmHg).
 * ALGEBRAIC[171] is AUC2 in component chemoreceptors_effect_of_art_PO2 (dimensionless).
 * ALGEBRAIC[17] is AU6C in component arterial_baroreceptor_reflex (dimensionless).
 * CONSTANTS[46] is AUX in component autonomics_parameter_values (dimensionless).
 * CONSTANTS[47] is AUK in component autonomics_parameter_values (per_minute).
 * CONSTANTS[48] is BAROTC in component autonomics_parameter_values (minute).
 * ALGEBRAIC[66] is AUB in component arterial_baroreceptor_reflex (dimensionless).
 * ALGEBRAIC[68] is A1B in component arterial_baroreceptor_reflex (dimensionless).
 * ALGEBRAIC[70] is AU6A in component arterial_baroreceptor_reflex (dimensionless).
 * CONSTANTS[49] is AU4 in component arterial_baroreceptor_reflex (dimensionless).
 * STATES[5] is AU6 in component arterial_baroreceptor_reflex (dimensionless).
 * ALGEBRAIC[72] is AUN in component CNS_ischemic_reflex (dimensionless).
 * CONSTANTS[50] is AUN1 in component autonomics_parameter_values (dimensionless).
 * ALGEBRAIC[49] is AULP in component autonomic_response_to_vasculature_pressure (dimensionless).
 * CONSTANTS[51] is AULPM in component autonomics_parameter_values (dimensionless).
 * CONSTANTS[206] is AUEX in component autonomic_response_to_exercise (dimensionless).
 * CONSTANTS[52] is EXC in component autonomics_parameter_values (dimensionless).
 * CONSTANTS[53] is EXCXP in component autonomics_parameter_values (dimensionless).
 * ALGEBRAIC[177] is AUTTL in component total_autonomic_stimulation (dimensionless).
 * ALGEBRAIC[175] is AUTTL1 in component total_autonomic_stimulation (dimensionless).
 * CONSTANTS[54] is EXCML in component autonomics_parameter_values (dimensionless).
 * CONSTANTS[55] is AUDMP in component autonomics_parameter_values (minute).
 * CONSTANTS[56] is AUMAX in component autonomics_parameter_values (dimensionless).
 * CONSTANTS[57] is AUMIN in component autonomics_parameter_values (dimensionless).
 * CONSTANTS[58] is AUSLP in component autonomics_parameter_values (dimensionless).
 * ALGEBRAIC[180] is DAU in component actual_autonomic_stimulation (per_minute).
 * STATES[6] is AU1 in component actual_autonomic_stimulation (dimensionless).
 * ALGEBRAIC[18] is AUT in component actual_autonomic_stimulation (dimensionless).
 * CONSTANTS[59] is VV9 in component autonomics_parameter_values (litre).
 * CONSTANTS[60] is AUL in component autonomics_parameter_values (litre).
 * CONSTANTS[61] is AUV in component autonomics_parameter_values (dimensionless).
 * CONSTANTS[62] is AUS in component autonomics_parameter_values (dimensionless).
 * CONSTANTS[63] is O2A in component autonomics_parameter_values (dimensionless).
 * CONSTANTS[64] is AUM1 in component autonomics_parameter_values (dimensionless).
 * CONSTANTS[65] is AUM2 in component autonomics_parameter_values (dimensionless).
 * CONSTANTS[66] is AUY in component autonomics_parameter_values (dimensionless).
 * ALGEBRAIC[21] is AUO in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * CONSTANTS[67] is MDMP in component autonomics_parameter_values (dimensionless).
 * ALGEBRAIC[81] is VEC in component extracellular_fluid_volume (litre).
 * ALGEBRAIC[152] is PPD in component concentration_of_protein_in_pulmonary_interstitium (gram_per_minute).
 * ALGEBRAIC[233] is RVS in component venous_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[149] is DFP in component pulmonary_interstitial_free_fluid_volume (L_per_minute).
 * ALGEBRAIC[123] is VPF in component pulmonary_interstitial_free_fluid_volume (litre).
 * ALGEBRAIC[234] is BFN in component blood_flow_through_NM_NR_tissues (L_per_minute).
 * ALGEBRAIC[209] is PVS in component venous_average_pressure (mmHg).
 * ALGEBRAIC[235] is PC in component capillary_pressure (mmHg).
 * ALGEBRAIC[27] is CPP in component plasma_protein_concentration (gram_per_L).
 * ALGEBRAIC[30] is PPC in component plasma_colloid_osmotic_pressure (mmHg).
 * STATES[7] is VP in component plasma_volume (litre).
 * ALGEBRAIC[134] is PGH in component hydrostatic_pressure_of_tissue_gel (mmHg).
 * ALGEBRAIC[132] is PTC in component total_osmotic_pressure_of_tissue_gel (mmHg).
 * ALGEBRAIC[269] is VTCPL in component plasma_leakage (L_per_minute).
 * ALGEBRAIC[277] is VTC in component rate_of_fluid_out_of_capillaries (L_per_minute).
 * CONSTANTS[68] is CFC in component capillary_dynamics_parameter_values (L_per_minute_per_mmHg).
 * ALGEBRAIC[139] is VTL in component lymph_flow (L_per_minute).
 * CONSTANTS[69] is TRPL in component capillary_dynamics_parameter_values (L_per_minute).
 * ALGEBRAIC[284] is VPD in component plasma_volume (L_per_minute).
 * STATES[8] is PRP in component total_plasma_protein (gram).
 * ALGEBRAIC[29] is DLP in component protein_destruction_and_formation (gram_per_minute).
 * CONSTANTS[70] is CPR in component capillary_dynamics_parameter_values (gram_per_L).
 * CONSTANTS[71] is LPPR in component capillary_dynamics_parameter_values (gram_per_minute).
 * CONSTANTS[72] is LPDE in component capillary_dynamics_parameter_values (dimensionless).
 * CONSTANTS[73] is LPK in component capillary_dynamics_parameter_values (L_per_minute).
 * ALGEBRAIC[28] is CPPD in component protein_destruction_and_formation (gram_per_L).
 * CONSTANTS[74] is PCR in component capillary_dynamics_parameter_values (mmHg).
 * CONSTANTS[75] is CPK in component capillary_dynamics_parameter_values (L_per_minute_per_mmHg).
 * CONSTANTS[76] is PCE in component capillary_dynamics_parameter_values (dimensionless).
 * ALGEBRAIC[263] is PRCD in component plasma_leakage (mmHg).
 * ALGEBRAIC[127] is CPI in component interstitial_protein_concentration (gram_per_L).
 * ALGEBRAIC[278] is DPC in component protein_influx_into_interstitium (gram_per_minute).
 * ALGEBRAIC[140] is DPL in component lymph_protein_flow (gram_per_minute).
 * ALGEBRAIC[285] is DPP in component total_plasma_protein (gram_per_minute).
 * ALGEBRAIC[124] is VTS in component total_systemic_fluid_volume (litre).
 * ALGEBRAIC[125] is VTS1 in component interstitial_fluid_volume (litre).
 * CONSTANTS[77] is TSSLML in component capillary_dynamics_parameter_values (dimensionless).
 * CONSTANTS[78] is TSSLTC in component capillary_dynamics_parameter_values (per_minute).
 * STATES[9] is VTS2 in component interstitial_fluid_volume (litre).
 * STATES[10] is TSP in component total_interstitial_protein (gram).
 * ALGEBRAIC[286] is DPI in component total_interstitial_protein (gram_per_minute).
 * ALGEBRAIC[129] is PTCPR in component interstitial_colloid_osmotic_pressure (mmHg).
 * ALGEBRAIC[133] is PTT in component total_tissue_pressure (mmHg).
 * ALGEBRAIC[130] is CHY in component hydrostatic_pressure_of_tissue_gel (gram_per_L).
 * CONSTANTS[79] is HYL in component capillary_dynamics_parameter_values (gram).
 * CONSTANTS[80] is CMPTSS in component capillary_dynamics_parameter_values (dimensionless).
 * CONSTANTS[81] is PGHF in component capillary_dynamics_parameter_values (L_mmHg_per_gram).
 * ALGEBRAIC[131] is POSHYL in component total_osmotic_pressure_of_tissue_gel (mmHg).
 * CONSTANTS[82] is GCOPF in component capillary_dynamics_parameter_values (per_mmHg).
 * CONSTANTS[83] is VTSF in component capillary_dynamics_parameter_values (litre).
 * ALGEBRAIC[135] is PIF in component interstial_free_fluid_pressure (mmHg).
 * ALGEBRAIC[137] is PTS in component interstitial_solid_tissue_pressure (mmHg).
 * CONSTANTS[84] is PLDF in component capillary_dynamics_parameter_values (mmHg).
 * ALGEBRAIC[136] is PLD1 in component lymph_flow (mmHg).
 * ALGEBRAIC[138] is PLD in component lymph_flow (mmHg).
 * ALGEBRAIC[126] is VG in component interstitial_gel_volume (litre).
 * ALGEBRAIC[128] is VIF in component interstitial_free_fluid_volume (litre).
 * ALGEBRAIC[100] is AMM in component global_M_blood_flow_autoregulation_output (dimensionless).
 * ALGEBRAIC[112] is ARM in component global_NM_blood_flow_autoregulation_output (dimensionless).
 * ALGEBRAIC[198] is ATRRFB in component volume_effect_on_arteries (dimensionless).
 * ALGEBRAIC[205] is ATRVFB in component volume_effect_on_unstressed_venous_volume (litre).
 * ALGEBRAIC[88] is HMD in component heart_deterioration (dimensionless).
 * STATES[11] is HPL in component left_ventricular_hypertrophy (dimensionless).
 * STATES[12] is HPR in component right_ventricular_hypertrophy (dimensionless).
 * ALGEBRAIC[158] is OSA in component arterial_PO2 (dimensionless).
 * ALGEBRAIC[264] is RBF in component actual_renal_blood_flow (L_per_minute).
 * ALGEBRAIC[182] is VIM in component blood_viscosity (dimensionless).
 * STATES[13] is VRC in component RBC_volume (litre).
 * STATES[14] is VV6 in component long_term_stress_relaxation (litre).
 * STATES[15] is VV7 in component short_term_stress_relaxation (litre).
 * ALGEBRAIC[215] is BFM in component blood_flow_through_M_tissues (L_per_minute).
 * ALGEBRAIC[279] is QAO in component systemic_blood_flow (L_per_minute).
 * ALGEBRAIC[168] is QRO in component right_ventricular_output (L_per_minute).
 * ALGEBRAIC[164] is QLO in component left_ventricular_output (L_per_minute).
 * ALGEBRAIC[44] is RPA in component pulmonary_arterial_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[53] is RPV in component pulmonary_venous_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[207] is VVE in component venous_excess_volume (litre).
 * STATES[16] is VVS1 in component venous_blood_volume (litre).
 * STATES[17] is VAS1 in component arterial_blood_volume (litre).
 * STATES[18] is VLA1 in component left_atrial_blood_volume (litre).
 * STATES[19] is VPA1 in component pulmonary_vasculature_blood_volume (litre).
 * STATES[20] is VRA1 in component right_atrial_blood_volume (litre).
 * ALGEBRAIC[31] is VBD in component total_blood_volume_change (litre).
 * ALGEBRAIC[212] is QVO in component rate_of_blood_flow_from_veins_to_right_atrium (L_per_minute).
 * ALGEBRAIC[32] is VRA in component right_atrial_blood_volume (litre).
 * ALGEBRAIC[214] is DRA in component right_atrial_blood_volume (L_per_minute).
 * ALGEBRAIC[33] is VRE in component right_atrial_pressure (litre).
 * ALGEBRAIC[36] is PRA1 in component autonomic_stimulation_effect_on_right_atrial_pressure (mmHg).
 * CONSTANTS[85] is HTAUML in component circulatory_dynamics_parameter_values (dimensionless).
 * ALGEBRAIC[160] is RVM in component pressure_effect_on_right_ventricular_pumping (dimensionless).
 * ALGEBRAIC[159] is PP2 in component pressure_effect_on_right_ventricular_pumping (mmHg).
 * ALGEBRAIC[58] is QLN in component left_ventricular_output (L_per_minute).
 * ALGEBRAIC[165] is HPEF in component pumping_effectiveness_of_right_ventricle (L_per_minute).
 * CONSTANTS[86] is QRF in component circulatory_dynamics_parameter_values (L_per_minute).
 * CONSTANTS[87] is HSR in component circulatory_dynamics_parameter_values (dimensionless).
 * ALGEBRAIC[37] is QRN in component right_ventricular_output (dimensionless).
 * ALGEBRAIC[56] is QPO in component rate_of_blood_flow_from_pulmonary_veins_to_left_atrium (L_per_minute).
 * ALGEBRAIC[38] is VPA in component pulmonary_vasculature_blood_volume (litre).
 * ALGEBRAIC[170] is DPA in component pulmonary_vasculature_blood_volume (L_per_minute).
 * ALGEBRAIC[39] is VPE in component pulmonary_vasculature_pressure (litre).
 * ALGEBRAIC[41] is PP1T in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
 * ALGEBRAIC[42] is PP1 in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
 * ALGEBRAIC[43] is CPA in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
 * ALGEBRAIC[51] is PL1 in component pulmonary_venous_resistance (mmHg).
 * ALGEBRAIC[54] is RPT in component total_pulmonary_vascular_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[55] is PGL in component pressure_gradient_through_the_lungs (mmHg).
 * ALGEBRAIC[45] is VLA in component left_atrial_blood_volume (litre).
 * ALGEBRAIC[166] is DLA in component left_atrial_blood_volume (L_per_minute).
 * ALGEBRAIC[46] is VLE in component left_atrial_pressure (litre).
 * ALGEBRAIC[57] is PLA1 in component autonomic_stimulation_effect_on_left_atrial_pressure (mmHg).
 * ALGEBRAIC[162] is LVM in component pumping_effectiveness_of_left_ventricle (dimensionless).
 * ALGEBRAIC[161] is PA2 in component pumping_effectiveness_of_left_ventricle (mmHg).
 * ALGEBRAIC[163] is QLOT in component left_ventricular_output (L_per_minute).
 * CONSTANTS[88] is HSL in component circulatory_dynamics_parameter_values (dimensionless).
 * ALGEBRAIC[73] is QLO1 in component left_ventricular_output (L_per_minute).
 * ALGEBRAIC[59] is VVS in component venous_blood_volume (litre).
 * ALGEBRAIC[287] is DVS in component venous_blood_volume (L_per_minute).
 * ALGEBRAIC[60] is VVA in component angiotensin_induced_venous_constriction (litre).
 * CONSTANTS[89] is ANY in component circulatory_dynamics_parameter_values (litre).
 * ALGEBRAIC[206] is VVE1 in component venous_excess_volume (litre).
 * CONSTANTS[90] is CV in component circulatory_dynamics_parameter_values (L_per_mmHg).
 * ALGEBRAIC[208] is PVS1 in component venous_average_pressure (mmHg).
 * ALGEBRAIC[61] is PR1 in component venous_outflow_pressure_into_heart (mmHg).
 * CONSTANTS[91] is PR1LL in component circulatory_dynamics_parameter_values (mmHg).
 * ALGEBRAIC[210] is RVG in component resistance_from_veins_to_right_atrium (mmHg_minute_per_L).
 * ALGEBRAIC[211] is PGV in component rate_of_blood_flow_from_veins_to_right_atrium (mmHg).
 * CONSTANTS[92] is CN7 in component circulatory_dynamics_parameter_values (dimensionless).
 * CONSTANTS[93] is CN2 in component circulatory_dynamics_parameter_values (per_mmHg).
 * CONSTANTS[94] is RVSM in component circulatory_dynamics_parameter_values (mmHg_minute_per_L).
 * ALGEBRAIC[236] is CN3 in component venous_resistance (dimensionless).
 * ALGEBRAIC[237] is RV1 in component venous_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[238] is NNRVR in component NM_NR_venous_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[62] is VAS in component arterial_blood_volume (litre).
 * ALGEBRAIC[288] is DAS in component arterial_blood_volume (L_per_minute).
 * ALGEBRAIC[74] is PAG in component arterial_pressure_and_pressure_gradient (mmHg).
 * ALGEBRAIC[63] is VAE in component arterial_pressure_and_pressure_gradient (litre).
 * ALGEBRAIC[75] is PAM in component pressure_effect_on_arterial_distention (dimensionless).
 * CONSTANTS[95] is PAEX in component circulatory_dynamics_parameter_values (dimensionless).
 * ALGEBRAIC[200] is R1 in component non_renal_systemic_arterial_resistance_multiplier (dimensionless).
 * ALGEBRAIC[202] is NNRAR in component NM_NR_arterial_resistance (mmHg_minute_per_L).
 * CONSTANTS[96] is RAR in component circulatory_dynamics_parameter_values (mmHg_minute_per_L).
 * CONSTANTS[97] is RMULT1 in component circulatory_dynamics_parameter_values (dimensionless).
 * ALGEBRAIC[213] is PGS in component pressure_gradient_from_arteries_to_veins (mmHg).
 * ALGEBRAIC[204] is RSM in component M_systemic_resistance (mmHg_minute_per_L).
 * CONSTANTS[98] is RAM in component circulatory_dynamics_parameter_values (mmHg_minute_per_L).
 * ALGEBRAIC[239] is RSN in component total_NM_NR_systemic_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[76] is FISFLO in component blood_flow_through_AV_fistulas (L_per_minute).
 * CONSTANTS[99] is FIS in component circulatory_dynamics_parameter_values (L_per_minute_per_mmHg).
 * ALGEBRAIC[270] is SYSFLO in component systemic_blood_flow (L_per_minute).
 * ALGEBRAIC[289] is RTP in component total_peripheral_resistance (mmHg_minute_per_L).
 * STATES[21] is TVD in component rate_of_fluid_intake (L_per_minute).
 * ALGEBRAIC[311] is NOD in component actual_Na_excretion_rate (monovalent_mEq_per_minute).
 * ALGEBRAIC[187] is STH in component effect_of_salt_appetite_stimulation_on_thirst (dimensionless).
 * ALGEBRAIC[315] is KOD in component actual_K_excretion_rate (monovalent_mEq_per_minute).
 * ALGEBRAIC[324] is VUD in component actual_urine_volume (L_per_minute).
 * STATES[22] is VTW in component total_body_water (litre).
 * CONSTANTS[100] is NID in component electrolytes_parameter_values (monovalent_mEq_per_minute).
 * CONSTANTS[101] is TRPL in component electrolytes_parameter_values (L_per_minute).
 * ALGEBRAIC[313] is NED in component extracellular_Na_concentration (monovalent_mEq_per_minute).
 * STATES[23] is NAE in component extracellular_Na_concentration (monovalent_mEq).
 * ALGEBRAIC[77] is AMK1 in component aldosterone_effect_on_cellular_K_distribution (dimensionless).
 * CONSTANTS[102] is ALCLK in component electrolytes_parameter_values (dimensionless).
 * ALGEBRAIC[78] is KE in component extracellular_K_concentration (monovalent_mEq).
 * STATES[24] is KTOT in component extracellular_K_concentration (monovalent_mEq).
 * CONSTANTS[103] is KID in component electrolytes_parameter_values (monovalent_mEq_per_minute).
 * ALGEBRAIC[317] is KTOTD in component extracellular_K_concentration (monovalent_mEq_per_minute).
 * STATES[25] is VIC in component intracellular_fluid_volume (litre).
 * ALGEBRAIC[80] is CKI in component intracellular_K_concentration (monovalent_mEq_per_litre).
 * ALGEBRAIC[79] is KI in component intracellular_K_concentration (monovalent_mEq).
 * ALGEBRAIC[89] is VID in component intracellular_fluid_volume (L_per_minute).
 * CONSTANTS[104] is VIDML in component electrolytes_parameter_values (litre2_per_monovalent_mEq_per_minute).
 * ALGEBRAIC[85] is CCD in component intracellular_fluid_volume (monovalent_mEq_per_litre).
 * ALGEBRAIC[114] is POT in component pressure_of_O2_in_NM_tissue_cells (mmHg).
 * CONSTANTS[105] is HSL in component heart_hypertrophy_or_deterioration_parameter_values (dimensionless).
 * CONSTANTS[106] is Z13 in component heart_hypertrophy_or_deterioration_parameter_values (dimensionless).
 * CONSTANTS[107] is HSR in component heart_hypertrophy_or_deterioration_parameter_values (dimensionless).
 * CONSTANTS[108] is DHDTR in component heart_hypertrophy_or_deterioration_parameter_values (per_mmHg_per_minute).
 * STATES[26] is HMD1 in component heart_deterioration (dimensionless).
 * ALGEBRAIC[116] is DHM in component heart_deterioration (per_minute).
 * ALGEBRAIC[92] is AUHR in component effect_of_autonomic_stimulation_on_HR (beats_per_minute).
 * CONSTANTS[207] is PRHR in component effect_of_PRA_on_HR (beats_per_minute).
 * CONSTANTS[109] is PR1LL in component HR_and_SV_parameter_values (mmHg).
 * ALGEBRAIC[95] is HDHR in component effect_of_heart_deterioration_on_HR (dimensionless).
 * ALGEBRAIC[97] is HR in component heart_rate (beats_per_minute).
 * ALGEBRAIC[167] is SVO in component stroke_volume_output (litre).
 * ALGEBRAIC[102] is PMO in component pressure_of_O2_in_M_tissue_cells (mmHg).
 * ALGEBRAIC[104] is PDO in component M_autoregulatory_driving_force (mmHg).
 * ALGEBRAIC[106] is POE in component M_ST_sensitivity_control (mmHg).
 * CONSTANTS[110] is POM in component M_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * ALGEBRAIC[99] is AMM1 in component M_ST_time_delay_and_limit (dimensionless).
 * CONSTANTS[111] is A4K in component M_autoregulatory_local_blood_flow_parameter_values (minute).
 * CONSTANTS[112] is AMM4 in component M_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * STATES[27] is AMM1T in component M_ST_time_delay_and_limit (dimensionless).
 * ALGEBRAIC[107] is POF in component M_LT_sensitivity_control (mmHg).
 * CONSTANTS[113] is POM2 in component M_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * STATES[28] is AMM2 in component M_LT_time_delay (dimensionless).
 * CONSTANTS[114] is A4K2 in component M_autoregulatory_local_blood_flow_parameter_values (minute).
 * STATES[29] is OVA in component O2_volume_of_arterial_blood (mL_per_L).
 * ALGEBRAIC[157] is HM in component hematocrit_fraction (dimensionless).
 * ALGEBRAIC[229] is RMO in component delivery_of_O2_to_M_tissues (mL_per_minute).
 * ALGEBRAIC[216] is O2ARTM in component M_O2_blood_supply (mL_per_minute).
 * ALGEBRAIC[230] is PVO in component M_venous_O2_content (mmHg).
 * CONSTANTS[115] is EXC in component M_O2_delivery_parameter_values (dimensionless).
 * CONSTANTS[116] is EXCXP2 in component M_O2_delivery_parameter_values (dimensionless).
 * ALGEBRAIC[231] is OVS in component M_venous_O2_content (dimensionless).
 * ALGEBRAIC[105] is MMO in component metabolic_O2_consumption_by_M_tissue (mL_per_minute).
 * ALGEBRAIC[103] is P2O in component metabolic_O2_consumption_by_M_tissue (mmHg).
 * CONSTANTS[117] is OMM in component M_O2_delivery_parameter_values (mL_per_minute).
 * CONSTANTS[118] is PM5 in component M_O2_delivery_parameter_values (per_mmHg).
 * ALGEBRAIC[101] is QOM in component volume_of_O2_in_M_tissue (mL).
 * ALGEBRAIC[261] is DO2M in component volume_of_O2_in_M_tissue (mL_per_minute).
 * STATES[30] is QOM1 in component volume_of_O2_in_M_tissue (mL).
 * CONSTANTS[119] is PK2 in component M_O2_delivery_parameter_values (mmHg_per_mL).
 * ALGEBRAIC[117] is POD in component NM_autoregulatory_driving_force (mmHg).
 * CONSTANTS[120] is POR in component NM_autoregulatory_local_blood_flow_parameter_values (mmHg).
 * ALGEBRAIC[119] is POB in component NM_ST_sensitivity_control (mmHg).
 * CONSTANTS[121] is POK in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * ALGEBRAIC[108] is AR1 in component NM_ST_time_delay_and_damping (dimensionless).
 * CONSTANTS[122] is A1K in component NM_autoregulatory_local_blood_flow_parameter_values (minute).
 * STATES[31] is AR1T in component NM_ST_time_delay_and_damping (dimensionless).
 * ALGEBRAIC[120] is POA in component NM_I_sensitivity_control (mmHg).
 * CONSTANTS[123] is PON in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * ALGEBRAIC[109] is AR2 in component NM_I_time_delay_and_limit (dimensionless).
 * CONSTANTS[124] is A2K in component NM_autoregulatory_local_blood_flow_parameter_values (minute).
 * STATES[32] is AR2T in component NM_I_time_delay_and_limit (dimensionless).
 * ALGEBRAIC[121] is POC in component NM_LT_sensitivity_control (mmHg).
 * CONSTANTS[125] is POZ in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * ALGEBRAIC[110] is AR3 in component NM_LT_time_delay_and_limit (dimensionless).
 * CONSTANTS[126] is A3K in component NM_autoregulatory_local_blood_flow_parameter_values (minute).
 * STATES[33] is AR3T in component NM_LT_time_delay_and_limit (dimensionless).
 * ALGEBRAIC[111] is ARM1 in component total_NM_autoregulation (dimensionless).
 * CONSTANTS[127] is AUTOSN in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * ALGEBRAIC[271] is DOB in component delivery_of_O2_to_NM_tissues (mL_per_minute).
 * ALGEBRAIC[265] is O2ARTN in component NM_O2_blood_supply (mL_per_minute).
 * ALGEBRAIC[272] is POV in component NM_venous_O2_content (mmHg).
 * ALGEBRAIC[273] is OSV in component NM_venous_O2_content (dimensionless).
 * ALGEBRAIC[118] is MO2 in component O2_consumption_by_NM_tissue (mL_per_minute).
 * CONSTANTS[128] is O2M in component NM_O2_delivery_parameter_values (mL_per_minute).
 * ALGEBRAIC[115] is P1O in component O2_consumption_by_NM_tissue (mmHg).
 * ALGEBRAIC[113] is QO2 in component volume_of_O2_in_NM_tissue (mL).
 * ALGEBRAIC[290] is DO2N in component volume_of_O2_in_NM_tissue (mL_per_minute).
 * ALGEBRAIC[280] is DO2N1 in component volume_of_O2_in_NM_tissue (mL_per_minute).
 * STATES[34] is QO2T in component volume_of_O2_in_NM_tissue (mL).
 * ALGEBRAIC[122] is PCP in component pulmonary_capillary_pressure (mmHg).
 * ALGEBRAIC[144] is POS in component colloid_osmotic_pressure_of_pulmonary_interstitium (mmHg).
 * ALGEBRAIC[141] is PPI in component pulmonary_interstitial_fluid_pressure (mmHg).
 * ALGEBRAIC[145] is PFI in component fluid_filtration_into_pulmonary_interstitium (L_per_minute).
 * CONSTANTS[129] is CPF in component pulmonary_fluid_dynamics_parameter_values (L_per_minute_per_mmHg).
 * ALGEBRAIC[147] is PLF in component lung_lymphatic_protein_flow (L_per_minute).
 * ALGEBRAIC[148] is DFZ in component pulmonary_interstitial_free_fluid_volume (L_per_minute).
 * STATES[35] is VPF1 in component pulmonary_interstitial_free_fluid_volume (litre).
 * ALGEBRAIC[150] is PPO in component lung_lymphatic_protein_flow (gram_per_minute).
 * ALGEBRAIC[146] is PPN in component protein_leakage_into_pulmonary_interstitium (gram_per_minute).
 * ALGEBRAIC[143] is CPN in component concentration_of_protein_in_pulmonary_interstitium (gram_per_L).
 * ALGEBRAIC[151] is PPZ in component concentration_of_protein_in_pulmonary_interstitium (gram_per_minute).
 * STATES[36] is PPR1 in component concentration_of_protein_in_pulmonary_interstitium (gram).
 * ALGEBRAIC[142] is PPR in component concentration_of_protein_in_pulmonary_interstitium (gram).
 * ALGEBRAIC[281] is O2UTIL in component total_O2_utilization (mL_per_minute).
 * ALGEBRAIC[154] is O2VAD2 in component progressive_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
 * ALGEBRAIC[176] is O2VTS2 in component acute_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
 * ALGEBRAIC[291] is ALVENT in component alveolar_ventilation (L_per_minute).
 * CONSTANTS[130] is VNTSTM in component pulmonary_O2_uptake_parameter_values (dimensionless).
 * ALGEBRAIC[294] is PO2ALV in component alveolar_PO2 (mmHg).
 * CONSTANTS[131] is PO2AMB in component pulmonary_O2_uptake_parameter_values (mmHg).
 * ALGEBRAIC[296] is O2DFS in component respiratory_O2_diffusion_into_capillaries (mL_per_minute).
 * CONSTANTS[132] is PL2 in component pulmonary_O2_uptake_parameter_values (L_mL_per_minute_per_mmHg).
 * CONSTANTS[133] is VPTISS in component pulmonary_O2_uptake_parameter_values (litre).
 * ALGEBRAIC[153] is RSPDFC in component respiratory_O2_diffusion_into_capillaries (mL_per_minute_per_mmHg).
 * ALGEBRAIC[298] is DOVA in component O2_volume_of_arterial_blood (mL_per_L_per_minute).
 * ALGEBRAIC[174] is O2VTST in component acute_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
 * ALGEBRAIC[172] is O2VTST1 in component acute_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
 * ALGEBRAIC[179] is DO2VAD in component progressive_chemoreceptor_adaptation_of_alveolar_ventilation (per_minute).
 * STATES[37] is O2VAD1 in component progressive_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
 * ALGEBRAIC[156] is HM1 in component hematocrit_fraction (dimensionless).
 * ALGEBRAIC[155] is VB in component hematocrit_fraction (litre).
 * ALGEBRAIC[178] is VIE in component viscosity_due_to_RBCs (dimensionless).
 * CONSTANTS[134] is HMK in component red_cells_and_viscosity_parameter_values (dimensionless).
 * CONSTANTS[135] is HKM in component red_cells_and_viscosity_parameter_values (dimensionless).
 * ALGEBRAIC[181] is VIB in component blood_viscosity (dimensionless).
 * ALGEBRAIC[188] is HM7 in component oxygen_stimulation (mmHg).
 * CONSTANTS[136] is PO2AMB in component red_cells_and_viscosity_parameter_values (mmHg).
 * CONSTANTS[137] is HM6 in component red_cells_and_viscosity_parameter_values (mmHg).
 * CONSTANTS[208] is PO2AM1 in component oxygen_stimulation (mmHg).
 * ALGEBRAIC[184] is HM3 in component oxygen_stimulation (mmHg).
 * CONSTANTS[209] is HM4 in component oxygen_stimulation (mmHg).
 * ALGEBRAIC[186] is HM5 in component oxygen_stimulation (mmHg).
 * ALGEBRAIC[190] is RC1 in component RBC_production (L_per_minute).
 * CONSTANTS[138] is HM8 in component red_cells_and_viscosity_parameter_values (L_per_minute_per_mmHg).
 * CONSTANTS[139] is REK in component red_cells_and_viscosity_parameter_values (dimensionless).
 * ALGEBRAIC[193] is RC2 in component RBC_destruction (L_per_minute).
 * CONSTANTS[140] is RKC in component red_cells_and_viscosity_parameter_values (per_minute).
 * CONSTANTS[141] is TRRBC in component red_cells_and_viscosity_parameter_values (L_per_minute).
 * ALGEBRAIC[196] is RCD in component RBC_volume (L_per_minute).
 * CONSTANTS[142] is SR in component stress_relaxation_parameter_values (dimensionless).
 * CONSTANTS[143] is SRK in component stress_relaxation_parameter_values (minute).
 * CONSTANTS[144] is SR2 in component stress_relaxation_parameter_values (dimensionless).
 * CONSTANTS[145] is SRK2 in component stress_relaxation_parameter_values (minute).
 * CONSTANTS[146] is ANMSLT in component thirst_drinking_and_salt_appetite_parameter_values (dimensionless).
 * CONSTANTS[147] is Z10 in component thirst_drinking_and_salt_appetite_parameter_values (mmHg).
 * CONSTANTS[148] is Z11 in component thirst_drinking_and_salt_appetite_parameter_values (per_mmHg2).
 * ALGEBRAIC[183] is ANMSML in component effect_of_salt_appetite_stimulation_on_thirst (dimensionless).
 * ALGEBRAIC[185] is STH1 in component effect_of_salt_appetite_stimulation_on_thirst (dimensionless).
 * ALGEBRAIC[191] is AHCM in component effect_of_antidiuretic_hormone_on_thirst (dimensionless).
 * CONSTANTS[149] is AHTHM in component thirst_drinking_and_salt_appetite_parameter_values (dimensionless).
 * ALGEBRAIC[194] is ANMTH in component effect_of_angiotensin_on_thirst (dimensionless).
 * CONSTANTS[150] is ANMTM in component thirst_drinking_and_salt_appetite_parameter_values (dimensionless).
 * CONSTANTS[151] is DR in component thirst_drinking_and_salt_appetite_parameter_values (L_per_minute).
 * CONSTANTS[152] is TVDDL in component thirst_drinking_and_salt_appetite_parameter_values (minute).
 * ALGEBRAIC[199] is AHTH in component rate_of_fluid_intake (dimensionless).
 * ALGEBRAIC[197] is AHTH1 in component rate_of_fluid_intake (dimensionless).
 * ALGEBRAIC[203] is TVZ in component rate_of_fluid_intake (L_per_minute).
 * ALGEBRAIC[201] is TVZ1 in component rate_of_fluid_intake (L_per_minute).
 * ALGEBRAIC[192] is AHZ in component effect_of_pressure_on_volume_receptors (dimensionless).
 * CONSTANTS[153] is AH10 in component volume_receptors_parameter_values (dimensionless).
 * CONSTANTS[154] is AH9 in component volume_receptors_parameter_values (per_mmHg).
 * ALGEBRAIC[189] is AHZ1 in component effect_of_pressure_on_volume_receptors (dimensionless).
 * STATES[38] is AHY in component time_dependent_volume_receptor_adaptation (dimensionless).
 * CONSTANTS[155] is AH11 in component volume_receptors_parameter_values (minute).
 * ALGEBRAIC[195] is AH7 in component total_volume_nervous_feedback (dimensionless).
 * CONSTANTS[156] is ATRFBM in component volume_receptors_parameter_values (dimensionless).
 * CONSTANTS[157] is ATRVM in component volume_receptors_parameter_values (litre).
 * ALGEBRAIC[217] is PAR in component perfusion_pressure (mmHg).
 * CONSTANTS[158] is GBL in component kidney_parameter_values (mmHg).
 * CONSTANTS[159] is RAPRSP in component kidney_parameter_values (mmHg).
 * CONSTANTS[160] is RFCDFT in component kidney_parameter_values (dimensionless).
 * CONSTANTS[161] is RCDFPC in component kidney_parameter_values (dimensionless).
 * CONSTANTS[162] is RCDFDP in component kidney_parameter_values (minute).
 * STATES[39] is PAR1 in component perfusion_pressure (mmHg).
 * ALGEBRAIC[240] is RNAUG2 in component renal_autoregulatory_feedback_factor (dimensionless).
 * CONSTANTS[163] is RNAUGN in component kidney_parameter_values (minute_per_L).
 * CONSTANTS[164] is RNAULL in component kidney_parameter_values (dimensionless).
 * CONSTANTS[165] is RNAUUL in component kidney_parameter_values (dimensionless).
 * CONSTANTS[166] is RNAUAD in component kidney_parameter_values (per_minute).
 * ALGEBRAIC[241] is RNAUG1 in component renal_autoregulatory_feedback_factor (dimensionless).
 * ALGEBRAIC[242] is RNAUG1T in component renal_autoregulatory_feedback_factor (dimensionless).
 * STATES[40] is RNAUG3 in component renal_autoregulatory_feedback_factor (dimensionless).
 * ALGEBRAIC[219] is AUMK in component autonomic_effect_on_AAR (dimensionless).
 * CONSTANTS[167] is ARF in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[218] is AUMKT in component autonomic_effect_on_AAR (dimensionless).
 * ALGEBRAIC[221] is ANMAR in component angiotensin_effect_on_AAR (dimensionless).
 * CONSTANTS[168] is ANMAM in component kidney_parameter_values (dimensionless).
 * CONSTANTS[169] is ANMARL in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[220] is ANMAR1 in component angiotensin_effect_on_AAR (dimensionless).
 * ALGEBRAIC[243] is AAR1 in component AAR_calculation (mmHg_minute_per_L).
 * CONSTANTS[170] is AARK in component kidney_parameter_values (mmHg_minute_per_L).
 * ALGEBRAIC[244] is AAR in component atrial_natriuretic_peptide_effect_on_AAR (mmHg_minute_per_L).
 * CONSTANTS[171] is ANPXAF in component kidney_parameter_values (mmHg_minute_per_L).
 * CONSTANTS[172] is AARLL in component kidney_parameter_values (mmHg_minute_per_L).
 * ALGEBRAIC[245] is AART in component atrial_natriuretic_peptide_effect_on_AAR (mmHg_minute_per_L).
 * ALGEBRAIC[222] is AUMK2 in component autonomic_effect_on_EAR (dimensionless).
 * CONSTANTS[173] is AUMK1 in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[223] is ANMER in component angiotensin_effect_on_EAR (dimensionless).
 * CONSTANTS[174] is ANMEM in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[246] is RNAUG4 in component effect_of_renal_autoregulatory_feedback_on_EAR (dimensionless).
 * CONSTANTS[175] is EFAFR in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[247] is EAR in component EAR_calculation (mmHg_minute_per_L).
 * CONSTANTS[176] is EARK in component kidney_parameter_values (mmHg_minute_per_L).
 * CONSTANTS[177] is EARLL in component kidney_parameter_values (mmHg_minute_per_L).
 * ALGEBRAIC[248] is EAR1 in component EAR_calculation (mmHg_minute_per_L).
 * ALGEBRAIC[249] is RR in component total_renal_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[250] is RFN in component normal_renal_blood_flow (L_per_minute).
 * CONSTANTS[178] is REK in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[251] is GFN in component glomerular_filtration_rate (L_per_minute).
 * ALGEBRAIC[252] is GLPC in component glomerular_colloid_osmotic_pressure (mmHg).
 * CONSTANTS[179] is GPPD in component kidney_parameter_values (dimensionless).
 * CONSTANTS[180] is GLPCA in component kidney_parameter_values (mmHg).
 * ALGEBRAIC[253] is EFAFPR in component glomerular_colloid_osmotic_pressure (dimensionless).
 * ALGEBRAIC[254] is EFAFPR1 in component glomerular_colloid_osmotic_pressure (dimensionless).
 * ALGEBRAIC[255] is GLP in component glomerular_pressure (mmHg).
 * ALGEBRAIC[256] is APD in component glomerular_pressure (mmHg).
 * ALGEBRAIC[266] is GFR in component glomerular_filtration_rate (L_per_minute).
 * CONSTANTS[181] is PXTP in component kidney_parameter_values (mmHg).
 * CONSTANTS[182] is GFLC in component kidney_parameter_values (L_per_minute_per_mmHg).
 * CONSTANTS[183] is GFNLL in component kidney_parameter_values (L_per_minute).
 * ALGEBRAIC[257] is PFL in component glomerular_filtration_rate (mmHg).
 * ALGEBRAIC[258] is GFN1 in component glomerular_filtration_rate (L_per_minute).
 * CONSTANTS[184] is MDFL1 in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[259] is PTFL in component proximal_tubular_and_macula_densa_flow (L_per_minute).
 * ALGEBRAIC[260] is MDFLWT in component proximal_tubular_and_macula_densa_flow (L_per_minute).
 * ALGEBRAIC[274] is RTSPPC in component renal_tissue_osmotic_pressure (mmHg).
 * CONSTANTS[185] is RTPPR in component kidney_parameter_values (dimensionless).
 * CONSTANTS[186] is RTPPRS in component kidney_parameter_values (mmHg).
 * ALGEBRAIC[267] is RTSPPC1 in component renal_tissue_osmotic_pressure (mmHg).
 * ALGEBRAIC[318] is UROD in component actual_urea_excretion_rate (mOsm_per_minute).
 * STATES[41] is PLUR in component glomerular_urea_concentration (mOsm).
 * CONSTANTS[187] is URFORM in component kidney_parameter_values (mOsm_per_minute).
 * ALGEBRAIC[224] is PLURC in component plasma_urea_concentration (mOsm_per_litre).
 * ALGEBRAIC[282] is RCPRS in component peritubular_capillary_pressure (mmHg).
 * CONSTANTS[188] is RFABX in component kidney_parameter_values (dimensionless).
 * CONSTANTS[189] is RVRS in component kidney_parameter_values (mmHg_minute_per_L).
 * ALGEBRAIC[300] is RFABD in component peritubular_capillary_reabsorption_factor (dimensionless).
 * CONSTANTS[190] is RTSPRS in component kidney_parameter_values (mmHg).
 * CONSTANTS[191] is RABSC in component kidney_parameter_values (per_mmHg).
 * CONSTANTS[192] is RFABDP in component kidney_parameter_values (dimensionless).
 * CONSTANTS[193] is RFABDM in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[292] is RABSPR in component peritubular_capillary_reabsorption_factor (mmHg).
 * ALGEBRAIC[295] is RFAB1 in component peritubular_capillary_reabsorption_factor (dimensionless).
 * ALGEBRAIC[297] is RFAB in component peritubular_capillary_reabsorption_factor (dimensionless).
 * ALGEBRAIC[299] is RFABD1 in component peritubular_capillary_reabsorption_factor (dimensionless).
 * ALGEBRAIC[301] is DTNAI in component distal_tubular_Na_delivery (monovalent_mEq_per_minute).
 * ALGEBRAIC[303] is DTNARA in component Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
 * CONSTANTS[194] is DTNAR in component kidney_parameter_values (monovalent_mEq_per_minute).
 * CONSTANTS[195] is DIURET in component kidney_parameter_values (dimensionless).
 * CONSTANTS[196] is AHMNAR in component kidney_parameter_values (dimensionless).
 * CONSTANTS[197] is DTNARL in component kidney_parameter_values (monovalent_mEq_per_minute).
 * ALGEBRAIC[302] is DTNARA1 in component Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
 * ALGEBRAIC[226] is DTNANG in component angiotensin_induced_Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
 * CONSTANTS[198] is ANMNAM in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[225] is DTNANG1 in component angiotensin_induced_Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
 * ALGEBRAIC[304] is DTKI in component distal_tubular_K_delivery (monovalent_mEq_per_minute).
 * ALGEBRAIC[305] is RFABK in component effect_of_physical_forces_on_distal_K_reabsorption (monovalent_mEq_per_minute).
 * CONSTANTS[199] is RFABKM in component kidney_parameter_values (monovalent_mEq_per_minute).
 * ALGEBRAIC[307] is MDFLK in component effect_of_fluid_flow_on_distal_K_reabsorption (monovalent_mEq_per_minute).
 * CONSTANTS[200] is MDFLKM in component kidney_parameter_values (monovalent_mEq_per_litre).
 * ALGEBRAIC[306] is MDFLK1 in component effect_of_fluid_flow_on_distal_K_reabsorption (monovalent_mEq_per_minute).
 * ALGEBRAIC[314] is KODN in component normal_K_excretion (monovalent_mEq_per_minute).
 * ALGEBRAIC[323] is VUDN in component normal_urine_volume (L_per_minute).
 * STATES[42] is DTKA in component K_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
 * ALGEBRAIC[308] is DTKSC in component K_secretion_from_distal_tubules (monovalent_mEq_per_minute).
 * CONSTANTS[201] is ANMKEM in component kidney_parameter_values (dimensionless).
 * CONSTANTS[202] is ANMKEL in component kidney_parameter_values (dimensionless).
 * CONSTANTS[203] is CKEEX in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[227] is ANMKE1 in component K_secretion_from_distal_tubules (dimensionless).
 * ALGEBRAIC[228] is ANMKE in component K_secretion_from_distal_tubules (dimensionless).
 * ALGEBRAIC[310] is NODN in component normal_Na_excretion (monovalent_mEq_per_minute).
 * ALGEBRAIC[309] is NODN1 in component normal_Na_excretion (monovalent_mEq_per_minute).
 * ALGEBRAIC[312] is KODN1 in component normal_K_excretion (monovalent_mEq_per_minute).
 * ALGEBRAIC[316] is DTURI in component normal_urea_excretion (mOsm_per_minute).
 * ALGEBRAIC[319] is OSMOPN1 in component normal_osmolar_and_water_excretion (mOsm_per_minute).
 * ALGEBRAIC[320] is OSMOPN in component normal_osmolar_and_water_excretion (mOsm_per_minute).
 * ALGEBRAIC[321] is OSMOP1T in component normal_urine_volume (mOsm_per_minute).
 * ALGEBRAIC[322] is OSMOP1 in component normal_urine_volume (mOsm_per_minute).
 * CONSTANTS[204] is RNAGTC in component kidney_parameter_values (minute).
 * CONSTANTS[205] is GFNDMP in component kidney_parameter_values (dimensionless).
 * RATES[0] is d/dt AMC in component aldosterone_concentration (dimensionless).
 * RATES[1] is d/dt ANX1 in component time_delayed_angiotensin_formation (dimensionless).
 * RATES[2] is d/dt ANC in component angiotensin_concentration (dimensionless).
 * RATES[3] is d/dt ADHC in component ADH_in_blood (dimensionless).
 * RATES[4] is d/dt ANPC in component ANP_in_plasma (dimensionless).
 * RATES[5] is d/dt AU6 in component arterial_baroreceptor_reflex (dimensionless).
 * RATES[6] is d/dt AU1 in component actual_autonomic_stimulation (dimensionless).
 * RATES[7] is d/dt VP in component plasma_volume (litre).
 * RATES[8] is d/dt PRP in component total_plasma_protein (gram).
 * RATES[9] is d/dt VTS2 in component interstitial_fluid_volume (litre).
 * RATES[10] is d/dt TSP in component total_interstitial_protein (gram).
 * RATES[20] is d/dt VRA1 in component right_atrial_blood_volume (litre).
 * RATES[19] is d/dt VPA1 in component pulmonary_vasculature_blood_volume (litre).
 * RATES[18] is d/dt VLA1 in component left_atrial_blood_volume (litre).
 * RATES[16] is d/dt VVS1 in component venous_blood_volume (litre).
 * RATES[17] is d/dt VAS1 in component arterial_blood_volume (litre).
 * RATES[23] is d/dt NAE in component extracellular_Na_concentration (monovalent_mEq).
 * RATES[24] is d/dt KTOT in component extracellular_K_concentration (monovalent_mEq).
 * RATES[25] is d/dt VIC in component intracellular_fluid_volume (litre).
 * RATES[22] is d/dt VTW in component total_body_water (litre).
 * RATES[11] is d/dt HPL in component left_ventricular_hypertrophy (dimensionless).
 * RATES[12] is d/dt HPR in component right_ventricular_hypertrophy (dimensionless).
 * RATES[26] is d/dt HMD1 in component heart_deterioration (dimensionless).
 * RATES[27] is d/dt AMM1T in component M_ST_time_delay_and_limit (dimensionless).
 * RATES[28] is d/dt AMM2 in component M_LT_time_delay (dimensionless).
 * RATES[30] is d/dt QOM1 in component volume_of_O2_in_M_tissue (mL).
 * RATES[31] is d/dt AR1T in component NM_ST_time_delay_and_damping (dimensionless).
 * RATES[32] is d/dt AR2T in component NM_I_time_delay_and_limit (dimensionless).
 * RATES[33] is d/dt AR3T in component NM_LT_time_delay_and_limit (dimensionless).
 * RATES[34] is d/dt QO2T in component volume_of_O2_in_NM_tissue (mL).
 * RATES[35] is d/dt VPF1 in component pulmonary_interstitial_free_fluid_volume (litre).
 * RATES[36] is d/dt PPR1 in component concentration_of_protein_in_pulmonary_interstitium (gram).
 * RATES[29] is d/dt OVA in component O2_volume_of_arterial_blood (mL_per_L).
 * RATES[37] is d/dt O2VAD1 in component progressive_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
 * RATES[13] is d/dt VRC in component RBC_volume (litre).
 * RATES[15] is d/dt VV7 in component short_term_stress_relaxation (litre).
 * RATES[14] is d/dt VV6 in component long_term_stress_relaxation (litre).
 * RATES[21] is d/dt TVD in component rate_of_fluid_intake (L_per_minute).
 * RATES[38] is d/dt AHY in component time_dependent_volume_receptor_adaptation (dimensionless).
 * RATES[39] is d/dt PAR1 in component perfusion_pressure (mmHg).
 * RATES[40] is d/dt RNAUG3 in component renal_autoregulatory_feedback_factor (dimensionless).
 * RATES[41] is d/dt PLUR in component glomerular_urea_concentration (mOsm).
 * RATES[42] is d/dt DTKA in component K_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1.0;
CONSTANTS[1] = 1.0;
CONSTANTS[2] = 1.0;
CONSTANTS[3] = 2.5;
CONSTANTS[4] = 12;
CONSTANTS[5] = 0;
CONSTANTS[6] = 0;
STATES[0] = 1.0;
CONSTANTS[7] = 60;
CONSTANTS[8] = 5;
CONSTANTS[9] = 0;
CONSTANTS[10] = 0.65;
CONSTANTS[11] = 2.5;
CONSTANTS[12] = 0.5;
CONSTANTS[13] = 0.8;
CONSTANTS[14] = 15;
CONSTANTS[15] = 0.04;
STATES[1] = 0.0;
CONSTANTS[16] = 0;
CONSTANTS[17] = 5000;
CONSTANTS[18] = 1;
CONSTANTS[19] = 0;
CONSTANTS[20] = 0;
STATES[2] = 0.859476;
CONSTANTS[21] = 12;
CONSTANTS[22] = 1.8;
CONSTANTS[23] = 0.7;
CONSTANTS[24] = 0.4;
CONSTANTS[25] = 6;
CONSTANTS[26] = 0.8;
CONSTANTS[27] = 0;
CONSTANTS[28] = 5;
STATES[3] = 1.0;
CONSTANTS[29] = 139;
CONSTANTS[30] = 85;
CONSTANTS[31] = 0.3;
CONSTANTS[32] = 0;
CONSTANTS[33] = 15;
CONSTANTS[34] = 2.5;
CONSTANTS[35] = 0.93617;
CONSTANTS[36] = 0.2;
CONSTANTS[37] = 5;
CONSTANTS[38] = 0;
CONSTANTS[39] = 0;
STATES[4] = 1.0;
CONSTANTS[40] = 4;
CONSTANTS[41] = 10;
CONSTANTS[42] = 0;
CONSTANTS[43] = 0;
CONSTANTS[44] = 0.3;
CONSTANTS[45] = 0.01;
CONSTANTS[46] = 1;
CONSTANTS[47] = 0.004;
CONSTANTS[48] = 0.16;
CONSTANTS[49] = -0.060024;
STATES[5] = 1.00132;
CONSTANTS[50] = 0.5;
CONSTANTS[51] = 0;
CONSTANTS[52] = 1;
CONSTANTS[53] = 1.0;
CONSTANTS[54] = 0.01;
CONSTANTS[55] = 0.3;
CONSTANTS[56] = 5.0;
CONSTANTS[57] = 0.4;
CONSTANTS[58] = 1.5;
STATES[6] = 1.00007;
CONSTANTS[59] = 2.51;
CONSTANTS[60] = 1.5;
CONSTANTS[61] = 0.55;
CONSTANTS[62] = 1;
CONSTANTS[63] = 0.1;
CONSTANTS[64] = 3;
CONSTANTS[65] = 1;
CONSTANTS[66] = 0;
CONSTANTS[67] = 0;
STATES[7] = 3.00449;
CONSTANTS[68] = 0.01167;
CONSTANTS[69] = 0;
STATES[8] = 216.243;
CONSTANTS[70] = 40;
CONSTANTS[71] = 0.03;
CONSTANTS[72] = 8;
CONSTANTS[73] = 2.728e-14;
CONSTANTS[74] = 15;
CONSTANTS[75] = 0.000253;
CONSTANTS[76] = 1;
CONSTANTS[77] = 0.15;
CONSTANTS[78] = 0.005;
STATES[9] = 0.0;
STATES[10] = 279.945;
CONSTANTS[79] = 60;
CONSTANTS[80] = 2;
CONSTANTS[81] = -2;
CONSTANTS[82] = 0.8092;
CONSTANTS[83] = 6;
CONSTANTS[84] = 4.2;
STATES[11] = 1.00163;
STATES[12] = 1.00237;
STATES[13] = 2.00439;
STATES[14] = 0.0101913;
STATES[15] = 0.00366525;
STATES[16] = 3.28246;
STATES[17] = 0.862514;
STATES[18] = 0.379883;
STATES[19] = 0.38131;
STATES[20] = 0.100043;
CONSTANTS[85] = 0.4;
CONSTANTS[86] = 0.15;
CONSTANTS[87] = 1;
CONSTANTS[88] = 1;
CONSTANTS[89] = -0.2;
CONSTANTS[90] = 0.1;
CONSTANTS[91] = 0;
CONSTANTS[92] = 0.2;
CONSTANTS[93] = 0.0212;
CONSTANTS[94] = 1;
CONSTANTS[95] = 2;
CONSTANTS[96] = 30.52;
CONSTANTS[97] = 1;
CONSTANTS[98] = 96.3;
CONSTANTS[99] = 0;
STATES[21] = 0.000980838;
STATES[22] = 39.8952;
CONSTANTS[100] = 0.1;
CONSTANTS[101] = 0;
STATES[23] = 2109.91;
CONSTANTS[102] = 0.3;
STATES[24] = 3622.54;
CONSTANTS[103] = 0.08;
STATES[25] = 25.0404;
CONSTANTS[104] = 0.01;
CONSTANTS[105] = 1;
CONSTANTS[106] = 0.625;
CONSTANTS[107] = 1;
CONSTANTS[108] = 0.05;
STATES[26] = 1.0;
CONSTANTS[109] = 0;
CONSTANTS[110] = 0.04;
CONSTANTS[111] = 0.1;
CONSTANTS[112] = 0.005;
STATES[27] = 1.00269;
CONSTANTS[113] = 2;
STATES[28] = 1.09071;
CONSTANTS[114] = 40000;
STATES[29] = 204.497;
CONSTANTS[115] = 1;
CONSTANTS[116] = 0.17;
CONSTANTS[117] = 57.1;
CONSTANTS[118] = 30;
STATES[30] = 48.0839;
CONSTANTS[119] = 0.79167;
CONSTANTS[120] = 35;
CONSTANTS[121] = 0.1;
CONSTANTS[122] = 0.5;
STATES[31] = 1.02127;
CONSTANTS[123] = 0.1;
CONSTANTS[124] = 60;
STATES[32] = 1.01179;
CONSTANTS[125] = 2;
CONSTANTS[126] = 40000;
STATES[33] = 1.1448;
CONSTANTS[127] = 0.9;
CONSTANTS[128] = 164;
STATES[34] = 72.2362;
CONSTANTS[129] = 0.0003;
STATES[35] = 0.0123238;
STATES[36] = 0.419998;
CONSTANTS[130] = 1;
CONSTANTS[131] = 150;
CONSTANTS[132] = 1.8;
CONSTANTS[133] = 0.0175;
STATES[37] = 2.368e-07;
CONSTANTS[134] = 90;
CONSTANTS[135] = 0.53333;
CONSTANTS[136] = 150;
CONSTANTS[137] = 1850;
CONSTANTS[138] = 4.714e-08;
CONSTANTS[139] = 1;
CONSTANTS[140] = 5.8e-06;
CONSTANTS[141] = 0;
CONSTANTS[142] = 1;
CONSTANTS[143] = 5;
CONSTANTS[144] = 1;
CONSTANTS[145] = 10000;
CONSTANTS[146] = 2;
CONSTANTS[147] = 45;
CONSTANTS[148] = 0.01;
CONSTANTS[149] = 2;
CONSTANTS[150] = 1.5;
CONSTANTS[151] = 0;
CONSTANTS[152] = 30;
CONSTANTS[153] = 0.333;
CONSTANTS[154] = 1;
STATES[38] = 0.301963;
CONSTANTS[155] = 1000;
CONSTANTS[156] = 0;
CONSTANTS[157] = 0;
CONSTANTS[158] = 0;
CONSTANTS[159] = 0;
CONSTANTS[160] = 0;
CONSTANTS[161] = 0;
CONSTANTS[162] = 2000;
STATES[39] = 103.525;
CONSTANTS[163] = 0.6;
CONSTANTS[164] = 0.3;
CONSTANTS[165] = 10;
CONSTANTS[166] = 0;
STATES[40] = 0.0;
CONSTANTS[167] = 0.5;
CONSTANTS[168] = 1.4;
CONSTANTS[169] = 0.86;
CONSTANTS[170] = 1;
CONSTANTS[171] = 1.5;
CONSTANTS[172] = 4;
CONSTANTS[173] = 0.3;
CONSTANTS[174] = 1.6;
CONSTANTS[175] = 0;
CONSTANTS[176] = 1;
CONSTANTS[177] = 24;
CONSTANTS[178] = 1;
CONSTANTS[179] = 1.0;
CONSTANTS[180] = 1.0;
CONSTANTS[181] = 8;
CONSTANTS[182] = 0.0208333;
CONSTANTS[183] = 0.001;
CONSTANTS[184] = 10;
CONSTANTS[185] = 0.9;
CONSTANTS[186] = 15.2;
STATES[41] = 159.549;
CONSTANTS[187] = 0.24;
CONSTANTS[188] = 0.8;
CONSTANTS[189] = 19.167;
CONSTANTS[190] = 6;
CONSTANTS[191] = 0.5;
CONSTANTS[192] = 1;
CONSTANTS[193] = 0.3;
CONSTANTS[194] = 0.675;
CONSTANTS[195] = 1;
CONSTANTS[196] = 0.3;
CONSTANTS[197] = 1e-06;
CONSTANTS[198] = 1;
CONSTANTS[199] = 0.03;
CONSTANTS[200] = 0.667;
STATES[42] = 0.0367573;
CONSTANTS[201] = 2;
CONSTANTS[202] = 0.3;
CONSTANTS[203] = 4;
CONSTANTS[204] = 15;
CONSTANTS[205] = 3;
CONSTANTS[206] = pow(CONSTANTS[52], CONSTANTS[53]);
CONSTANTS[207] =  pow(CONSTANTS[109], 0.500000)*5.00000;
CONSTANTS[208] = (CONSTANTS[136]>80.0000 ? 80.0000 : CONSTANTS[136]);
CONSTANTS[209] = CONSTANTS[136] - 40.0000;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[31] = ((((((STATES[7]+STATES[13]) - STATES[16]) - STATES[17]) - STATES[18]) - STATES[19]) - STATES[20])/2.00000;
ALGEBRAIC[45] = STATES[18]+ ALGEBRAIC[31]*0.128000;
ALGEBRAIC[46] = ALGEBRAIC[45] - 0.380000;
ALGEBRAIC[47] = ALGEBRAIC[46]/0.0100000;
ALGEBRAIC[48] = ( (ALGEBRAIC[47] - 1.00000)*1.00000<0.00000 ? 0.00000 :  (ALGEBRAIC[47] - 1.00000)*1.00000);
ALGEBRAIC[32] = STATES[20]+ ALGEBRAIC[31]*0.0574000;
ALGEBRAIC[33] = ALGEBRAIC[32] - 0.100000;
ALGEBRAIC[34] = ALGEBRAIC[33]/0.00500000;
ALGEBRAIC[35] = ( (ALGEBRAIC[34]+1.00000)*2.00000<0.00000 ? 0.00000 :  (ALGEBRAIC[34]+1.00000)*2.00000);
ALGEBRAIC[50] = (ALGEBRAIC[48]+ALGEBRAIC[35])/3.00000;
ALGEBRAIC[52] = (CONSTANTS[38]>0.00000 ? CONSTANTS[38] : ALGEBRAIC[50]+CONSTANTS[39]);
RATES[4] = (ALGEBRAIC[52] - STATES[4])/CONSTANTS[40];
ALGEBRAIC[62] = STATES[17]+ ALGEBRAIC[31]*0.261000;
ALGEBRAIC[63] = ALGEBRAIC[62] - 0.495000;
ALGEBRAIC[64] = ALGEBRAIC[63]/0.00355000;
RATES[39] = ((100.000+ (ALGEBRAIC[64] - 100.000)*CONSTANTS[161]) - STATES[39])/CONSTANTS[162];
ALGEBRAIC[65] = (CONSTANTS[42]>1.00000e-07 ? CONSTANTS[42] : ALGEBRAIC[64] - CONSTANTS[43]);
ALGEBRAIC[66] = (ALGEBRAIC[65]<160.000&&ALGEBRAIC[65]>=80.0000 ?  0.0166670*(160.000 - ALGEBRAIC[65]) : ALGEBRAIC[65]<80.0000 ? 1.33360 : 0.00000);
ALGEBRAIC[68] =  (ALGEBRAIC[66] - 1.00000)*CONSTANTS[46]+1.00000;
ALGEBRAIC[70] = ALGEBRAIC[68] - CONSTANTS[49];
RATES[5] = (ALGEBRAIC[70] - STATES[5])/CONSTANTS[48];
ALGEBRAIC[81] = STATES[22] - STATES[25];
ALGEBRAIC[82] = STATES[23]/ALGEBRAIC[81];
ALGEBRAIC[0] = CONSTANTS[8] - (CONSTANTS[8] - 1.00000)/( ((CONSTANTS[9] - 1.00000)/(CONSTANTS[9] - CONSTANTS[8]))*(STATES[0] - 1.00000)*CONSTANTS[10]+1.00000);
ALGEBRAIC[1] =  (ALGEBRAIC[0] - 1.00000)*CONSTANTS[11]+1.00000;
ALGEBRAIC[2] =  (ALGEBRAIC[1] - 1.00000)*CONSTANTS[12]+1.00000;
ALGEBRAIC[3] = (ALGEBRAIC[2]<0.200000 ? 0.200000 : ALGEBRAIC[2]);
ALGEBRAIC[77] =  (ALGEBRAIC[3] - 1.00000)*CONSTANTS[102]+1.00000;
ALGEBRAIC[78] = (STATES[24] - 3000.00)/( ALGEBRAIC[77]*9.33330);
ALGEBRAIC[79] = STATES[24] - ALGEBRAIC[78];
ALGEBRAIC[80] = ALGEBRAIC[79]/STATES[25];
ALGEBRAIC[85] = ALGEBRAIC[80] - ALGEBRAIC[82];
ALGEBRAIC[89] =  ALGEBRAIC[85]*CONSTANTS[104];
RATES[25] = ALGEBRAIC[89];
ALGEBRAIC[83] = (ALGEBRAIC[82] - CONSTANTS[29])/(142.000 - CONSTANTS[29]);
ALGEBRAIC[87] = (ALGEBRAIC[83]<0.00000 ? 0.00000 : ALGEBRAIC[83]);
ALGEBRAIC[67] = (ALGEBRAIC[65]>CONSTANTS[30] ? CONSTANTS[30] : ALGEBRAIC[65]);
ALGEBRAIC[69] =  pow(CONSTANTS[30] - ALGEBRAIC[67], 2.00000)*CONSTANTS[31];
ALGEBRAIC[91] = ALGEBRAIC[87]+ALGEBRAIC[69]+CONSTANTS[32];
ALGEBRAIC[94] = (ALGEBRAIC[91]<0.00000 ? 0.00000 : ALGEBRAIC[91]);
RATES[3] = (ALGEBRAIC[94] - STATES[3])/CONSTANTS[33];
ALGEBRAIC[6] = CONSTANTS[22] - (CONSTANTS[22] - 1.00000)/( ((CONSTANTS[23] - 1.00000)/(CONSTANTS[23] - CONSTANTS[22]))*(STATES[2] - 1.00000)*CONSTANTS[24]+1.00000);
ALGEBRAIC[7] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[3]+1.00000;
ALGEBRAIC[84] = ALGEBRAIC[78]/ALGEBRAIC[81];
ALGEBRAIC[86] = (ALGEBRAIC[84] - 3.30000)/1.00000;
ALGEBRAIC[90] =  ALGEBRAIC[7]*0.909000*ALGEBRAIC[86];
ALGEBRAIC[93] =  (ALGEBRAIC[90] - 1.00000)*CONSTANTS[4]+1.00000;
ALGEBRAIC[96] = (ALGEBRAIC[93]<0.00000 ? 0.00000 : ALGEBRAIC[93]);
ALGEBRAIC[98] = (CONSTANTS[6]>0.00000 ? CONSTANTS[6] : ALGEBRAIC[96]+CONSTANTS[5]);
RATES[0] = (ALGEBRAIC[98] - STATES[0])/CONSTANTS[7];
ALGEBRAIC[101] = (STATES[30]<0.000100000 ? 0.000100000 : STATES[30]);
ALGEBRAIC[102] =  CONSTANTS[119]*ALGEBRAIC[101];
ALGEBRAIC[104] = ALGEBRAIC[102] - 38.0000;
ALGEBRAIC[106] =  ALGEBRAIC[104]*CONSTANTS[110]+1.00000;
RATES[27] = ( ALGEBRAIC[106]*1.00000 - STATES[27])/CONSTANTS[111];
ALGEBRAIC[107] =  CONSTANTS[113]*ALGEBRAIC[104]+1.00000;
RATES[28] = ( ALGEBRAIC[107]*1.00000 - STATES[28])/CONSTANTS[114];
ALGEBRAIC[113] = (STATES[34]<0.00000 ? 0.00000 : STATES[34]);
ALGEBRAIC[114] =  ALGEBRAIC[113]*0.486110;
ALGEBRAIC[116] =  (ALGEBRAIC[114] - 10.0000)*CONSTANTS[108];
RATES[26] = ALGEBRAIC[116];
ALGEBRAIC[117] = ALGEBRAIC[114] - CONSTANTS[120];
ALGEBRAIC[119] =  ALGEBRAIC[117]*CONSTANTS[121]+1.00000;
RATES[31] = ( ALGEBRAIC[119]*1.00000 - STATES[31])/CONSTANTS[122];
ALGEBRAIC[120] =  CONSTANTS[123]*ALGEBRAIC[117]+1.00000;
RATES[32] = ( ALGEBRAIC[120]*1.00000 - STATES[32])/CONSTANTS[124];
ALGEBRAIC[121] =  CONSTANTS[125]*ALGEBRAIC[117]+1.00000;
RATES[33] = ( ALGEBRAIC[121]*1.00000 - STATES[33])/CONSTANTS[126];
ALGEBRAIC[123] = (STATES[35]<0.00100000 ? 0.00100000 : STATES[35]);
ALGEBRAIC[124] = (ALGEBRAIC[81] - STATES[7]) - ALGEBRAIC[123];
RATES[9] =  ( (ALGEBRAIC[124] - 12.0000)*CONSTANTS[77] - STATES[9])*CONSTANTS[78];
ALGEBRAIC[27] = STATES[8]/STATES[7];
ALGEBRAIC[30] =  0.280000*ALGEBRAIC[27]+ 0.00190000*pow(ALGEBRAIC[27], 2.00000);
ALGEBRAIC[38] = STATES[19]+ ALGEBRAIC[31]*0.155000;
ALGEBRAIC[39] = ALGEBRAIC[38] - 0.306250;
ALGEBRAIC[40] = ALGEBRAIC[39]/0.00480000;
ALGEBRAIC[41] =  0.0260000*ALGEBRAIC[40];
ALGEBRAIC[42] = (ALGEBRAIC[41]<1.00000e-05 ? 1.00000e-05 : ALGEBRAIC[41]);
ALGEBRAIC[43] = pow(ALGEBRAIC[42], 0.500000);
ALGEBRAIC[44] = 1.00000/ALGEBRAIC[43];
ALGEBRAIC[51] = ALGEBRAIC[47]+18.0000;
ALGEBRAIC[53] = 1.00000/( ALGEBRAIC[51]*0.0357000);
ALGEBRAIC[122] = ( (ALGEBRAIC[40] - ALGEBRAIC[47])*ALGEBRAIC[53])/(ALGEBRAIC[53]+ALGEBRAIC[44])+ALGEBRAIC[47];
ALGEBRAIC[142] = (STATES[36]<0.0250000 ? 0.0250000 : STATES[36]);
ALGEBRAIC[143] = ALGEBRAIC[142]/ALGEBRAIC[123];
ALGEBRAIC[144] =  ALGEBRAIC[143]*0.400000;
ALGEBRAIC[141] = 2.00000 - 0.150000/ALGEBRAIC[123];
ALGEBRAIC[145] =  (((ALGEBRAIC[122] - ALGEBRAIC[141])+ALGEBRAIC[144]) - ALGEBRAIC[30])*CONSTANTS[129];
ALGEBRAIC[147] =  (ALGEBRAIC[141]+11.0000)*0.000300000;
ALGEBRAIC[148] = ALGEBRAIC[145] - ALGEBRAIC[147];
ALGEBRAIC[149] = ALGEBRAIC[148];
RATES[35] = ALGEBRAIC[149];
ALGEBRAIC[150] =  ALGEBRAIC[147]*ALGEBRAIC[143];
ALGEBRAIC[146] =  (ALGEBRAIC[27] - ALGEBRAIC[143])*0.000225000;
ALGEBRAIC[151] = ALGEBRAIC[146] - ALGEBRAIC[150];
ALGEBRAIC[152] = ALGEBRAIC[151];
RATES[36] = ALGEBRAIC[152];
ALGEBRAIC[18] = CONSTANTS[56] - (CONSTANTS[56] - 1.00000)/exp( CONSTANTS[58]*(STATES[6] - 1.00000));
ALGEBRAIC[19] = (ALGEBRAIC[18]<CONSTANTS[57] ? CONSTANTS[57] : ALGEBRAIC[18]);
ALGEBRAIC[21] = ALGEBRAIC[19] - 1.00000;
ALGEBRAIC[22] =  ALGEBRAIC[21]*CONSTANTS[61]+1.00000;
ALGEBRAIC[88] = (STATES[26]>1.00000 ? 1.00000 : STATES[26]);
ALGEBRAIC[57] =  (ALGEBRAIC[47]+4.00000)*( CONSTANTS[85]*(ALGEBRAIC[19] - 1.00000)+1.00000) - 4.00000;
ALGEBRAIC[58] = (ALGEBRAIC[57]<=- 2.00000 ? 0.0100000 : ALGEBRAIC[57]>- 2.00000&&ALGEBRAIC[57]<=1.00000 ? 0.0100000+( (3.60000 - 0.0100000)*(ALGEBRAIC[57] - - 2.00000))/(1.00000 - - 2.00000) : ALGEBRAIC[57]>1.00000&&ALGEBRAIC[57]<=5.00000 ? 3.60000+( (9.40000 - 3.60000)*(ALGEBRAIC[57] - 1.00000))/(5.00000 - 1.00000) : ALGEBRAIC[57]>5.00000&&ALGEBRAIC[57]<=8.00000 ? 9.40000+( (11.6000 - 9.40000)*(ALGEBRAIC[57] - 5.00000))/(8.00000 - 5.00000) : ALGEBRAIC[57]>8.00000&&ALGEBRAIC[57]<=12.0000 ? 11.6000+( (13.5000 - 11.6000)*(ALGEBRAIC[57] - 8.00000))/(12.0000 - 8.00000) : 13.5000);
ALGEBRAIC[155] = STATES[7]+STATES[13];
ALGEBRAIC[156] = STATES[13]/ALGEBRAIC[155];
ALGEBRAIC[157] =  100.000*ALGEBRAIC[156];
ALGEBRAIC[158] = (STATES[29]/ALGEBRAIC[157])/5.25000;
ALGEBRAIC[161] = ALGEBRAIC[64]/( ALGEBRAIC[22]*ALGEBRAIC[158]);
ALGEBRAIC[162] = (ALGEBRAIC[161]<=0.00000 ? 1.04000 : ALGEBRAIC[161]>0.00000&&ALGEBRAIC[161]<=60.0000 ? 1.04000+( (1.02500 - 1.04000)*(ALGEBRAIC[161] - 0.00000))/(60.0000 - 0.00000) : ALGEBRAIC[161]>60.0000&&ALGEBRAIC[161]<=125.000 ? 1.02500+( (0.970000 - 1.02500)*(ALGEBRAIC[161] - 60.0000))/(125.000 - 60.0000) : ALGEBRAIC[161]>125.000&&ALGEBRAIC[161]<=160.000 ? 0.970000+( (0.880000 - 0.970000)*(ALGEBRAIC[161] - 125.000))/(160.000 - 125.000) : ALGEBRAIC[161]>160.000&&ALGEBRAIC[161]<=200.000 ? 0.880000+( (0.590000 - 0.880000)*(ALGEBRAIC[161] - 160.000))/(200.000 - 160.000) : ALGEBRAIC[161]>200.000&&ALGEBRAIC[161]<=240.000 ? 0.590000+( (0.00000 - 0.590000)*(ALGEBRAIC[161] - 200.000))/(240.000 - 200.000) : 0.00000);
ALGEBRAIC[163] =  ALGEBRAIC[162]*ALGEBRAIC[58]*ALGEBRAIC[22]*CONSTANTS[88]*ALGEBRAIC[88]*STATES[11];
ALGEBRAIC[73] = (ALGEBRAIC[47] - ALGEBRAIC[64])/3.00000;
ALGEBRAIC[164] = (ALGEBRAIC[73]>0.00000 ? ALGEBRAIC[163]+ALGEBRAIC[73] : ALGEBRAIC[163]);
ALGEBRAIC[54] = ALGEBRAIC[53]+ALGEBRAIC[44];
ALGEBRAIC[55] = ALGEBRAIC[40] - ALGEBRAIC[47];
ALGEBRAIC[56] = ALGEBRAIC[55]/ALGEBRAIC[54];
ALGEBRAIC[166] = ALGEBRAIC[56] - ALGEBRAIC[164];
RATES[18] = ALGEBRAIC[166];
ALGEBRAIC[159] = (ALGEBRAIC[40]/ALGEBRAIC[22])/ALGEBRAIC[158];
ALGEBRAIC[160] = (ALGEBRAIC[159]<=0.00000 ? 1.06000 : ALGEBRAIC[159]>0.00000&&ALGEBRAIC[159]<=32.0000 ? 1.06000+( (0.970000 - 1.06000)*(ALGEBRAIC[159] - 0.00000))/(32.0000 - 0.00000) : ALGEBRAIC[159]>32.0000&&ALGEBRAIC[159]<=38.4000 ? 0.970000+( (0.930000 - 0.970000)*(ALGEBRAIC[159] - 32.0000))/(38.4000 - 32.0000) : ALGEBRAIC[159]>38.4000&&ALGEBRAIC[159]<=48.0000 ? 0.930000+( (0.800000 - 0.930000)*(ALGEBRAIC[159] - 38.4000))/(48.0000 - 38.4000) : ALGEBRAIC[159]>48.0000&&ALGEBRAIC[159]<=60.8000 ? 0.800000+( (0.460000 - 0.800000)*(ALGEBRAIC[159] - 48.0000))/(60.8000 - 48.0000) : ALGEBRAIC[159]>60.8000&&ALGEBRAIC[159]<=72.0000 ? 0.460000+( (0.00000 - 0.460000)*(ALGEBRAIC[159] - 60.8000))/(72.0000 - 60.8000) : 0.00000);
ALGEBRAIC[165] =  (1.00000 - CONSTANTS[86])*ALGEBRAIC[22]*ALGEBRAIC[160]*CONSTANTS[87]*ALGEBRAIC[88]*STATES[12]+( CONSTANTS[86]*ALGEBRAIC[164])/ALGEBRAIC[58];
ALGEBRAIC[36] =  (ALGEBRAIC[34]+8.00000)*( CONSTANTS[85]*(ALGEBRAIC[19] - 1.00000)+1.00000) - 8.00000;
ALGEBRAIC[37] = (ALGEBRAIC[36]<=- 8.00000 ? 0.00000 : ALGEBRAIC[36]>- 8.00000&&ALGEBRAIC[36]<=- 6.00000 ? 0.00000+( (0.750000 - 0.00000)*(ALGEBRAIC[36] - - 8.00000))/(- 6.00000 - - 8.00000) : ALGEBRAIC[36]>- 6.00000&&ALGEBRAIC[36]<=- 2.00000 ? 0.750000+( (2.60000 - 0.750000)*(ALGEBRAIC[36] - - 6.00000))/(- 2.00000 - - 6.00000) : ALGEBRAIC[36]>- 2.00000&&ALGEBRAIC[36]<=4.00000 ? 2.60000+( (9.80000 - 2.60000)*(ALGEBRAIC[36] - - 2.00000))/(4.00000 - - 2.00000) : ALGEBRAIC[36]>4.00000&&ALGEBRAIC[36]<=12.0000 ? 9.80000+( (13.5000 - 9.80000)*(ALGEBRAIC[36] - 4.00000))/(12.0000 - 4.00000) : 13.5000);
ALGEBRAIC[168] =  ALGEBRAIC[37]*ALGEBRAIC[165];
ALGEBRAIC[170] = ALGEBRAIC[168] - ALGEBRAIC[56];
RATES[19] = ALGEBRAIC[170];
ALGEBRAIC[169] = (ALGEBRAIC[158]>1.00000 ? 114.000+ (ALGEBRAIC[158] - 1.00000)*6667.00 : ALGEBRAIC[158]>0.936000&&ALGEBRAIC[158]<=1.00000 ? 74.0000+ (ALGEBRAIC[158] - 0.936000)*625.000 : ALGEBRAIC[158]>0.800000&&ALGEBRAIC[158]<=0.936000 ? 46.0000+ (ALGEBRAIC[158] - 0.800000)*205.882 :  ALGEBRAIC[158]*57.5000);
ALGEBRAIC[172] = (ALGEBRAIC[169] - 67.0000)/30.0000;
ALGEBRAIC[174] = (ALGEBRAIC[172]>1.00000 ? 1.00000 : ALGEBRAIC[172]<0.600000 ? 0.600000 : ALGEBRAIC[172]);
ALGEBRAIC[176] = 1.00000/ALGEBRAIC[174];
ALGEBRAIC[179] =  ( (ALGEBRAIC[176] - 1.00000)*3.00000 - STATES[37])*0.000500000;
RATES[37] = ALGEBRAIC[179];
ALGEBRAIC[71] = (ALGEBRAIC[65]<80.0000&&ALGEBRAIC[65]>=40.0000 ?  0.00500000*(80.0000 - ALGEBRAIC[65])*CONSTANTS[44] : ALGEBRAIC[65]<40.0000 ?  0.200000*CONSTANTS[44] : 0.00000);
ALGEBRAIC[171] = (ALGEBRAIC[169]<80.0000&&ALGEBRAIC[169]>=40.0000 ?  CONSTANTS[45]*(80.0000 - ALGEBRAIC[169]) : ALGEBRAIC[169]<40.0000 ?  CONSTANTS[45]*40.0000 : 0.00000);
ALGEBRAIC[173] = ALGEBRAIC[71]+ALGEBRAIC[171];
ALGEBRAIC[17] = STATES[5];
ALGEBRAIC[72] = (ALGEBRAIC[65]<40.0000 ?  0.0400000*(40.0000 - ALGEBRAIC[65])*CONSTANTS[50] : 0.00000);
ALGEBRAIC[49] =  (15.0000/(ALGEBRAIC[47]+ALGEBRAIC[34]+ALGEBRAIC[40]) - 1.00000)*CONSTANTS[51]+1.00000;
ALGEBRAIC[175] =  ( CONSTANTS[206]*ALGEBRAIC[49]*(ALGEBRAIC[173]+ALGEBRAIC[17]+ALGEBRAIC[72]) - 1.00000)*CONSTANTS[54]+1.00000;
ALGEBRAIC[177] = (ALGEBRAIC[175]<0.00000 ? 0.00000 : ALGEBRAIC[175]);
ALGEBRAIC[180] = (ALGEBRAIC[177] - STATES[6])/CONSTANTS[55];
RATES[6] = ALGEBRAIC[180];
ALGEBRAIC[189] =  pow(fabs(ALGEBRAIC[34]), CONSTANTS[153])*CONSTANTS[154];
ALGEBRAIC[192] = (ALGEBRAIC[34]<0.00000 ? - ALGEBRAIC[189] : ALGEBRAIC[189]);
RATES[38] = (ALGEBRAIC[192] - STATES[38])/CONSTANTS[155];
ALGEBRAIC[184] =  (CONSTANTS[208] - 40.0000)*ALGEBRAIC[157];
ALGEBRAIC[186] = (ALGEBRAIC[184]+CONSTANTS[209]<0.00000 ? 0.00000 : ALGEBRAIC[184]+CONSTANTS[209]);
ALGEBRAIC[188] = CONSTANTS[137] - ALGEBRAIC[186];
ALGEBRAIC[190] = ( ALGEBRAIC[188]*CONSTANTS[138]*CONSTANTS[139]+5.00000e-06<0.00000 ? 0.00000 :  ALGEBRAIC[188]*CONSTANTS[138]*CONSTANTS[139]+5.00000e-06);
ALGEBRAIC[178] = ALGEBRAIC[157]/( (CONSTANTS[134] - ALGEBRAIC[157])*CONSTANTS[135]);
ALGEBRAIC[181] = ALGEBRAIC[178]+1.50000;
ALGEBRAIC[182] =  0.333300*ALGEBRAIC[181];
ALGEBRAIC[193] =  STATES[13]*CONSTANTS[140]*ALGEBRAIC[182];
ALGEBRAIC[196] = (ALGEBRAIC[190] - ALGEBRAIC[193])+CONSTANTS[141];
RATES[13] = ALGEBRAIC[196];
ALGEBRAIC[194] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[150]*0.00100000;
ALGEBRAIC[183] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[146]+1.00000;
ALGEBRAIC[185] =  pow(CONSTANTS[147] - ALGEBRAIC[114], 2.00000)*CONSTANTS[148]*ALGEBRAIC[183];
ALGEBRAIC[187] = (ALGEBRAIC[185]<0.800000 ? 0.800000 : ALGEBRAIC[185]>8.00000 ? 8.00000 : ALGEBRAIC[185]);
ALGEBRAIC[191] =  (STATES[3] - 1.00000)*CONSTANTS[149]+1.00000;
ALGEBRAIC[197] =  ALGEBRAIC[191]*ALGEBRAIC[187]*0.00100000;
ALGEBRAIC[199] = (ALGEBRAIC[197]<0.00000 ? 0.00000 : ALGEBRAIC[197]);
ALGEBRAIC[201] =  (ALGEBRAIC[194]+ALGEBRAIC[199])*1.00000;
ALGEBRAIC[203] = (ALGEBRAIC[201]<0.00000 ? 0.00000 : ALGEBRAIC[201]);
RATES[21] = ((ALGEBRAIC[203]+CONSTANTS[151]) - STATES[21])/CONSTANTS[152];
ALGEBRAIC[20] = (CONSTANTS[59] -  ALGEBRAIC[19]*CONSTANTS[60])+CONSTANTS[60];
ALGEBRAIC[195] = ALGEBRAIC[192] - STATES[38];
ALGEBRAIC[205] =  ALGEBRAIC[195]*CONSTANTS[157];
ALGEBRAIC[59] = STATES[16]+ ALGEBRAIC[31]*0.398600;
ALGEBRAIC[8] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[25]+1.00000;
ALGEBRAIC[9] = (ALGEBRAIC[8]<CONSTANTS[26] ? CONSTANTS[26] : ALGEBRAIC[8]);
ALGEBRAIC[60] =  (ALGEBRAIC[9] - 1.00000)*CONSTANTS[89];
ALGEBRAIC[206] = ((((ALGEBRAIC[59] - ALGEBRAIC[20]) - ALGEBRAIC[60]) - STATES[15]) - STATES[14]) - ALGEBRAIC[205];
ALGEBRAIC[207] = (ALGEBRAIC[206]<0.000100000 ? 0.000100000 : ALGEBRAIC[206]);
RATES[15] = ( (ALGEBRAIC[207] - 0.740000)*CONSTANTS[142] - STATES[15])/CONSTANTS[143];
RATES[14] = ( (ALGEBRAIC[207] - 0.740000)*CONSTANTS[144] - STATES[14])/CONSTANTS[145];
ALGEBRAIC[208] = 3.70000+(ALGEBRAIC[207] - 0.740000)/CONSTANTS[90];
ALGEBRAIC[209] = (ALGEBRAIC[208]<0.000100000 ? 0.000100000 : ALGEBRAIC[208]);
ALGEBRAIC[210] = 0.740000/pow(ALGEBRAIC[209]/( ALGEBRAIC[182]*3.70000), 0.500000);
ALGEBRAIC[61] = (ALGEBRAIC[34]<CONSTANTS[91] ? CONSTANTS[91] : ALGEBRAIC[34]);
ALGEBRAIC[211] = ALGEBRAIC[209] - ALGEBRAIC[61];
ALGEBRAIC[212] = ALGEBRAIC[211]/ALGEBRAIC[210];
ALGEBRAIC[214] = ALGEBRAIC[212] - ALGEBRAIC[168];
RATES[20] = ALGEBRAIC[214];
ALGEBRAIC[213] = ALGEBRAIC[64] - ALGEBRAIC[209];
ALGEBRAIC[99] = (STATES[27]<CONSTANTS[112] ? CONSTANTS[112] : STATES[27]);
ALGEBRAIC[100] =  ALGEBRAIC[99]*STATES[28];
ALGEBRAIC[11] = CONSTANTS[34] - (CONSTANTS[34] - 1.00000)/( ((CONSTANTS[35] - 1.00000)/(CONSTANTS[35] - CONSTANTS[34]))*(STATES[3] - 1.00000)+1.00000);
ALGEBRAIC[12] = (ALGEBRAIC[11]<CONSTANTS[35] ? CONSTANTS[35] : ALGEBRAIC[11]);
ALGEBRAIC[25] = pow( ALGEBRAIC[21]*CONSTANTS[64]+1.00000, CONSTANTS[65]);
ALGEBRAIC[198] =  ALGEBRAIC[195]*CONSTANTS[156]+1.00000;
ALGEBRAIC[75] = pow(ALGEBRAIC[64]/100.000, CONSTANTS[95]);
ALGEBRAIC[200] = (( ALGEBRAIC[9]*ALGEBRAIC[12]*ALGEBRAIC[25]*ALGEBRAIC[182]*CONSTANTS[1])/ALGEBRAIC[75])/ALGEBRAIC[198];
ALGEBRAIC[204] =  CONSTANTS[98]*ALGEBRAIC[100]*ALGEBRAIC[200]*CONSTANTS[0]*CONSTANTS[97];
ALGEBRAIC[215] = ALGEBRAIC[213]/ALGEBRAIC[204];
ALGEBRAIC[216] =  STATES[29]*ALGEBRAIC[215];
rootfind_0(VOI, CONSTANTS, RATES, STATES, ALGEBRAIC, pret);
ALGEBRAIC[24] =  ALGEBRAIC[21]*CONSTANTS[63]+1.00000;
ALGEBRAIC[103] = (ALGEBRAIC[102]>38.0000 ? 38.0000 : ALGEBRAIC[102]);
ALGEBRAIC[105] =  ALGEBRAIC[24]*CONSTANTS[117]*CONSTANTS[115]*(1.00000 - pow(38.0001 - ALGEBRAIC[103], 3.00000)/54872.0);
ALGEBRAIC[261] = ALGEBRAIC[229] - ALGEBRAIC[105];
RATES[30] = ALGEBRAIC[261];
ALGEBRAIC[10] =  (ALGEBRAIC[9] - 1.00000)*CONSTANTS[27]+1.00000;
ALGEBRAIC[15] = CONSTANTS[41] - CONSTANTS[41]/( 0.555556*(1.00000+STATES[4]));
ALGEBRAIC[16] = (ALGEBRAIC[15]<- 1.00000 ? - 1.00000 : ALGEBRAIC[15]);
ALGEBRAIC[26] =  ALGEBRAIC[21]*CONSTANTS[66]+1.00000;
ALGEBRAIC[108] = (STATES[31]<0.500000 ? 0.500000 : STATES[31]);
ALGEBRAIC[109] = (STATES[32]<0.500000 ? 0.500000 : STATES[32]);
ALGEBRAIC[110] = (STATES[33]<0.300000 ? 0.300000 : STATES[33]);
ALGEBRAIC[111] =  ALGEBRAIC[108]*ALGEBRAIC[109]*ALGEBRAIC[110];
ALGEBRAIC[112] =  (ALGEBRAIC[111] - 1.00000)*CONSTANTS[127]+1.00000;
ALGEBRAIC[202] =  CONSTANTS[96]*ALGEBRAIC[112]*ALGEBRAIC[200]*CONSTANTS[0]*CONSTANTS[97];
ALGEBRAIC[217] = (CONSTANTS[159]>0.00000&&CONSTANTS[160]<=0.00000 ? CONSTANTS[159] : CONSTANTS[160]>0.00000 ? STATES[39] : ALGEBRAIC[64] - CONSTANTS[158]);
ALGEBRAIC[218] =  (ALGEBRAIC[25] - 1.00000)*CONSTANTS[167]+1.00000;
ALGEBRAIC[219] = (ALGEBRAIC[218]<0.800000 ? 0.800000 : ALGEBRAIC[218]);
ALGEBRAIC[220] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[168]+1.00000;
ALGEBRAIC[221] = (ALGEBRAIC[220]<CONSTANTS[169] ? CONSTANTS[169] : ALGEBRAIC[220]);
ALGEBRAIC[222] =  (ALGEBRAIC[219] - 1.00000)*CONSTANTS[173]+1.00000;
ALGEBRAIC[223] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[174]+1.00000;
rootfind_1(VOI, CONSTANTS, RATES, STATES, ALGEBRAIC, pret);
RATES[40] =  (ALGEBRAIC[240] - 1.00000)*CONSTANTS[166];
ALGEBRAIC[262] = ALGEBRAIC[232];
ALGEBRAIC[268] = (ALGEBRAIC[262]>1.00000 ? 1.00000/(1.00000+ (ALGEBRAIC[262] - 1.00000)*72.0000) : 10.0000 - 9.00000/(1.00000+ (1.00000 - ALGEBRAIC[262])*8.00000));
ALGEBRAIC[275] =  (ALGEBRAIC[268] - 1.00000)*CONSTANTS[16];
RATES[1] = (ALGEBRAIC[275] - STATES[1])/CONSTANTS[17];
ALGEBRAIC[74] = ALGEBRAIC[64] - ALGEBRAIC[34];
ALGEBRAIC[76] =  ALGEBRAIC[74]*CONSTANTS[99];
ALGEBRAIC[264] =  CONSTANTS[178]*ALGEBRAIC[250];
ALGEBRAIC[270] = ALGEBRAIC[215]+ALGEBRAIC[234]+ALGEBRAIC[264];
ALGEBRAIC[279] = ALGEBRAIC[270]+ALGEBRAIC[76];
RATES[11] = (pow(( ALGEBRAIC[64]*ALGEBRAIC[279])/( 500.000*CONSTANTS[105]), CONSTANTS[106]) - STATES[11])/57600.0;
RATES[12] = (pow(( ALGEBRAIC[40]*ALGEBRAIC[279])/( 75.0000*CONSTANTS[107]), CONSTANTS[106]) - STATES[12])/57600.0;
ALGEBRAIC[125] = ALGEBRAIC[124] - STATES[9];
ALGEBRAIC[133] =  pow((ALGEBRAIC[125] - CONSTANTS[83])/CONSTANTS[83], 2.00000)*1.00000;
ALGEBRAIC[130] = pow((CONSTANTS[79]/ALGEBRAIC[124])/5.00000, CONSTANTS[80]);
ALGEBRAIC[134] =  ALGEBRAIC[130]*CONSTANTS[81]+ALGEBRAIC[133];
ALGEBRAIC[127] = STATES[10]/ALGEBRAIC[124];
ALGEBRAIC[129] =  0.280000*ALGEBRAIC[127]+ 0.00190000*pow(ALGEBRAIC[127], 2.00000);
ALGEBRAIC[131] =  ALGEBRAIC[130]*2.00000;
ALGEBRAIC[132] =  ALGEBRAIC[131]*ALGEBRAIC[129]*CONSTANTS[82];
ALGEBRAIC[263] = (ALGEBRAIC[235] - CONSTANTS[74]<0.00000 ? 0.00000 : ALGEBRAIC[235] - CONSTANTS[74]);
ALGEBRAIC[269] = pow( ALGEBRAIC[263]*CONSTANTS[75], CONSTANTS[76]);
ALGEBRAIC[277] =  (((ALGEBRAIC[235] - ALGEBRAIC[30]) - ALGEBRAIC[134])+ALGEBRAIC[132])*CONSTANTS[68]+ALGEBRAIC[269];
ALGEBRAIC[135] = ALGEBRAIC[134] - ALGEBRAIC[131];
ALGEBRAIC[136] = (ALGEBRAIC[135]+CONSTANTS[84]) - ALGEBRAIC[133];
ALGEBRAIC[138] = (ALGEBRAIC[136]>7.00000 ? 7.00000 : ALGEBRAIC[136]);
ALGEBRAIC[139] = (ALGEBRAIC[138]<0.00000 ? 0.00000 :  ALGEBRAIC[138]*0.0200000);
ALGEBRAIC[284] = ((ALGEBRAIC[139] - ALGEBRAIC[277]) - ALGEBRAIC[149])+CONSTANTS[69];
RATES[7] = ALGEBRAIC[284];
ALGEBRAIC[28] = (ALGEBRAIC[27] - CONSTANTS[70]<0.00000 ? 0.00000 : ALGEBRAIC[27] - CONSTANTS[70]);
ALGEBRAIC[29] = CONSTANTS[71] -  pow(ALGEBRAIC[28], CONSTANTS[72])*CONSTANTS[73];
ALGEBRAIC[278] =  ALGEBRAIC[269]*ALGEBRAIC[27]+ (ALGEBRAIC[27] - ALGEBRAIC[127])*0.00104000;
ALGEBRAIC[140] =  ALGEBRAIC[127]*ALGEBRAIC[139];
ALGEBRAIC[285] = (((ALGEBRAIC[29]+ALGEBRAIC[140]) - ALGEBRAIC[278]) - ALGEBRAIC[152])+ CONSTANTS[69]*72.0000;
RATES[8] = ALGEBRAIC[285];
ALGEBRAIC[286] = ALGEBRAIC[278] - ALGEBRAIC[140];
RATES[10] = ALGEBRAIC[286];
ALGEBRAIC[287] = ALGEBRAIC[279] - ALGEBRAIC[212];
RATES[16] = ALGEBRAIC[287];
ALGEBRAIC[288] = ALGEBRAIC[164] - ALGEBRAIC[279];
RATES[17] = ALGEBRAIC[288];
ALGEBRAIC[265] =  STATES[29]*ALGEBRAIC[234];
rootfind_2(VOI, CONSTANTS, RATES, STATES, ALGEBRAIC, pret);
ALGEBRAIC[115] = (ALGEBRAIC[114]>35.0000 ? 35.0000 : ALGEBRAIC[114]);
ALGEBRAIC[118] =  ALGEBRAIC[24]*CONSTANTS[128]*(1.00000 - pow(35.0001 - ALGEBRAIC[115], 3.00000)/42875.0);
ALGEBRAIC[280] = ALGEBRAIC[271] - ALGEBRAIC[118];
ALGEBRAIC[290] = (ALGEBRAIC[113]<6.00000&&ALGEBRAIC[280]<0.00000 ?  ALGEBRAIC[280]*0.100000 : ALGEBRAIC[280]);
RATES[34] = ALGEBRAIC[290];
ALGEBRAIC[276] =  (ALGEBRAIC[268]+STATES[1])*CONSTANTS[18];
ALGEBRAIC[283] = (ALGEBRAIC[276]<1.00000e-05 ? 1.00000e-05 : ALGEBRAIC[276]);
ALGEBRAIC[293] = (CONSTANTS[19]>0.00000 ? CONSTANTS[19] : ALGEBRAIC[283]+CONSTANTS[20]);
RATES[2] = (ALGEBRAIC[293] - STATES[2])/CONSTANTS[21];
ALGEBRAIC[281] = ALGEBRAIC[271]+ALGEBRAIC[229];
ALGEBRAIC[154] = STATES[37]+1.00000;
ALGEBRAIC[291] =  ALGEBRAIC[281]*CONSTANTS[130]*0.0266670*ALGEBRAIC[176]*ALGEBRAIC[154];
ALGEBRAIC[294] = CONSTANTS[131] - (ALGEBRAIC[281]/ALGEBRAIC[291])/0.761000;
ALGEBRAIC[153] = CONSTANTS[132]/(CONSTANTS[133]+ALGEBRAIC[123]);
ALGEBRAIC[296] =  (ALGEBRAIC[294] - ALGEBRAIC[169])*ALGEBRAIC[153];
ALGEBRAIC[298] = (ALGEBRAIC[296] - ALGEBRAIC[281])/( ALGEBRAIC[168]*1.00000);
RATES[29] = ALGEBRAIC[298];
ALGEBRAIC[301] =  ALGEBRAIC[232]*ALGEBRAIC[82]*0.00616190;
ALGEBRAIC[4] =  (ALGEBRAIC[1] - 1.00000)*CONSTANTS[13]+1.00000;
ALGEBRAIC[5] = (ALGEBRAIC[4]<CONSTANTS[15] ? CONSTANTS[15] : ALGEBRAIC[4]>CONSTANTS[14] ? CONSTANTS[14] : ALGEBRAIC[4]);
ALGEBRAIC[13] = CONSTANTS[37] - (CONSTANTS[37] - 1.00000)/( ((CONSTANTS[36] - 1.00000)/(CONSTANTS[36] - CONSTANTS[37]))*(STATES[3] - 1.00000)+1.00000);
ALGEBRAIC[14] = (ALGEBRAIC[13]<CONSTANTS[36] ? CONSTANTS[36] : ALGEBRAIC[13]);
ALGEBRAIC[267] =  ALGEBRAIC[252]*CONSTANTS[185] - CONSTANTS[186];
ALGEBRAIC[274] = (ALGEBRAIC[267]<1.00000 ? 1.00000 : ALGEBRAIC[267]);
ALGEBRAIC[282] =  ( (ALGEBRAIC[250] - 1.20000)*CONSTANTS[188]+1.20000)*CONSTANTS[189];
ALGEBRAIC[292] = ((ALGEBRAIC[252]+CONSTANTS[190]) - ALGEBRAIC[282]) - ALGEBRAIC[274];
ALGEBRAIC[295] =  ALGEBRAIC[292]*CONSTANTS[191];
ALGEBRAIC[297] = ALGEBRAIC[295];
ALGEBRAIC[299] =  (ALGEBRAIC[297] - 1.00000)*CONSTANTS[193]+1.00000;
ALGEBRAIC[300] = (ALGEBRAIC[299]<0.000100000 ? 0.000100000 : ALGEBRAIC[299]);
ALGEBRAIC[302] =  (( ALGEBRAIC[5]*ALGEBRAIC[300]*CONSTANTS[194])/CONSTANTS[195])*( (ALGEBRAIC[14] - 1.00000)*CONSTANTS[196]+1.00000);
ALGEBRAIC[303] = (ALGEBRAIC[302]<CONSTANTS[197] ? CONSTANTS[197] : ALGEBRAIC[302]);
ALGEBRAIC[225] =  ( (ALGEBRAIC[6] - 1.00000)*CONSTANTS[198]+1.00000)*0.100000;
ALGEBRAIC[226] = (ALGEBRAIC[225]<0.00000 ? 0.00000 : ALGEBRAIC[225]);
ALGEBRAIC[309] = (ALGEBRAIC[301] - ALGEBRAIC[303]) - ALGEBRAIC[226];
ALGEBRAIC[310] = (ALGEBRAIC[309]<1.00000e-08 ? 1.00000e-08 : ALGEBRAIC[309]);
ALGEBRAIC[311] =  ALGEBRAIC[310]*CONSTANTS[178];
ALGEBRAIC[313] = ( CONSTANTS[100]*ALGEBRAIC[187] - ALGEBRAIC[311])+ CONSTANTS[101]*142.000;
RATES[23] = ALGEBRAIC[313];
ALGEBRAIC[304] = ( ALGEBRAIC[301]*ALGEBRAIC[84])/ALGEBRAIC[82];
ALGEBRAIC[305] =  (ALGEBRAIC[300] - 1.00000)*CONSTANTS[199];
ALGEBRAIC[306] =  (ALGEBRAIC[232] - 1.00000)*CONSTANTS[200]+1.00000;
ALGEBRAIC[307] = (ALGEBRAIC[306]<0.100000 ? 0.100000 : ALGEBRAIC[306]);
ALGEBRAIC[227] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[201]+1.00000;
ALGEBRAIC[228] = (ALGEBRAIC[227]<CONSTANTS[202] ? CONSTANTS[202] : ALGEBRAIC[227]);
ALGEBRAIC[308] = ( pow(ALGEBRAIC[84]/4.40000, CONSTANTS[203])*ALGEBRAIC[3]*0.0800000*ALGEBRAIC[307])/ALGEBRAIC[228];
ALGEBRAIC[312] = ((ALGEBRAIC[304]+ALGEBRAIC[308]) - STATES[42]) - ALGEBRAIC[305];
ALGEBRAIC[314] = (ALGEBRAIC[312]<0.00000 ? 0.00000 : ALGEBRAIC[312]);
ALGEBRAIC[315] =  ALGEBRAIC[314]*CONSTANTS[178];
ALGEBRAIC[317] = CONSTANTS[103] - ALGEBRAIC[315];
RATES[24] = ALGEBRAIC[317];
ALGEBRAIC[224] = STATES[41]/STATES[22];
ALGEBRAIC[316] =  pow(ALGEBRAIC[251], 2.00000)*ALGEBRAIC[224]*3.84000;
ALGEBRAIC[318] =  ALGEBRAIC[316]*CONSTANTS[178];
RATES[41] = CONSTANTS[187] - ALGEBRAIC[318];
ALGEBRAIC[319] = ALGEBRAIC[316]+ 2.00000*(ALGEBRAIC[310]+ALGEBRAIC[314]);
ALGEBRAIC[320] = (ALGEBRAIC[319]>0.600000 ? 0.600000 : ALGEBRAIC[319]);
ALGEBRAIC[321] = ALGEBRAIC[319] - 0.600000;
ALGEBRAIC[322] = (ALGEBRAIC[321]<0.00000 ? 0.00000 : ALGEBRAIC[321]);
ALGEBRAIC[323] = ALGEBRAIC[320]/( 600.000*ALGEBRAIC[14])+ALGEBRAIC[322]/360.000;
RATES[42] =  ( (ALGEBRAIC[314]/ALGEBRAIC[323])*0.000451800 - STATES[42])*1.00000;
ALGEBRAIC[324] =  ALGEBRAIC[323]*CONSTANTS[178];
RATES[22] = STATES[21] - ALGEBRAIC[324];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[31] = ((((((STATES[7]+STATES[13]) - STATES[16]) - STATES[17]) - STATES[18]) - STATES[19]) - STATES[20])/2.00000;
ALGEBRAIC[45] = STATES[18]+ ALGEBRAIC[31]*0.128000;
ALGEBRAIC[46] = ALGEBRAIC[45] - 0.380000;
ALGEBRAIC[47] = ALGEBRAIC[46]/0.0100000;
ALGEBRAIC[48] = ( (ALGEBRAIC[47] - 1.00000)*1.00000<0.00000 ? 0.00000 :  (ALGEBRAIC[47] - 1.00000)*1.00000);
ALGEBRAIC[32] = STATES[20]+ ALGEBRAIC[31]*0.0574000;
ALGEBRAIC[33] = ALGEBRAIC[32] - 0.100000;
ALGEBRAIC[34] = ALGEBRAIC[33]/0.00500000;
ALGEBRAIC[35] = ( (ALGEBRAIC[34]+1.00000)*2.00000<0.00000 ? 0.00000 :  (ALGEBRAIC[34]+1.00000)*2.00000);
ALGEBRAIC[50] = (ALGEBRAIC[48]+ALGEBRAIC[35])/3.00000;
ALGEBRAIC[52] = (CONSTANTS[38]>0.00000 ? CONSTANTS[38] : ALGEBRAIC[50]+CONSTANTS[39]);
ALGEBRAIC[62] = STATES[17]+ ALGEBRAIC[31]*0.261000;
ALGEBRAIC[63] = ALGEBRAIC[62] - 0.495000;
ALGEBRAIC[64] = ALGEBRAIC[63]/0.00355000;
ALGEBRAIC[65] = (CONSTANTS[42]>1.00000e-07 ? CONSTANTS[42] : ALGEBRAIC[64] - CONSTANTS[43]);
ALGEBRAIC[66] = (ALGEBRAIC[65]<160.000&&ALGEBRAIC[65]>=80.0000 ?  0.0166670*(160.000 - ALGEBRAIC[65]) : ALGEBRAIC[65]<80.0000 ? 1.33360 : 0.00000);
ALGEBRAIC[68] =  (ALGEBRAIC[66] - 1.00000)*CONSTANTS[46]+1.00000;
ALGEBRAIC[70] = ALGEBRAIC[68] - CONSTANTS[49];
ALGEBRAIC[81] = STATES[22] - STATES[25];
ALGEBRAIC[82] = STATES[23]/ALGEBRAIC[81];
ALGEBRAIC[0] = CONSTANTS[8] - (CONSTANTS[8] - 1.00000)/( ((CONSTANTS[9] - 1.00000)/(CONSTANTS[9] - CONSTANTS[8]))*(STATES[0] - 1.00000)*CONSTANTS[10]+1.00000);
ALGEBRAIC[1] =  (ALGEBRAIC[0] - 1.00000)*CONSTANTS[11]+1.00000;
ALGEBRAIC[2] =  (ALGEBRAIC[1] - 1.00000)*CONSTANTS[12]+1.00000;
ALGEBRAIC[3] = (ALGEBRAIC[2]<0.200000 ? 0.200000 : ALGEBRAIC[2]);
ALGEBRAIC[77] =  (ALGEBRAIC[3] - 1.00000)*CONSTANTS[102]+1.00000;
ALGEBRAIC[78] = (STATES[24] - 3000.00)/( ALGEBRAIC[77]*9.33330);
ALGEBRAIC[79] = STATES[24] - ALGEBRAIC[78];
ALGEBRAIC[80] = ALGEBRAIC[79]/STATES[25];
ALGEBRAIC[85] = ALGEBRAIC[80] - ALGEBRAIC[82];
ALGEBRAIC[89] =  ALGEBRAIC[85]*CONSTANTS[104];
ALGEBRAIC[83] = (ALGEBRAIC[82] - CONSTANTS[29])/(142.000 - CONSTANTS[29]);
ALGEBRAIC[87] = (ALGEBRAIC[83]<0.00000 ? 0.00000 : ALGEBRAIC[83]);
ALGEBRAIC[67] = (ALGEBRAIC[65]>CONSTANTS[30] ? CONSTANTS[30] : ALGEBRAIC[65]);
ALGEBRAIC[69] =  pow(CONSTANTS[30] - ALGEBRAIC[67], 2.00000)*CONSTANTS[31];
ALGEBRAIC[91] = ALGEBRAIC[87]+ALGEBRAIC[69]+CONSTANTS[32];
ALGEBRAIC[94] = (ALGEBRAIC[91]<0.00000 ? 0.00000 : ALGEBRAIC[91]);
ALGEBRAIC[6] = CONSTANTS[22] - (CONSTANTS[22] - 1.00000)/( ((CONSTANTS[23] - 1.00000)/(CONSTANTS[23] - CONSTANTS[22]))*(STATES[2] - 1.00000)*CONSTANTS[24]+1.00000);
ALGEBRAIC[7] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[3]+1.00000;
ALGEBRAIC[84] = ALGEBRAIC[78]/ALGEBRAIC[81];
ALGEBRAIC[86] = (ALGEBRAIC[84] - 3.30000)/1.00000;
ALGEBRAIC[90] =  ALGEBRAIC[7]*0.909000*ALGEBRAIC[86];
ALGEBRAIC[93] =  (ALGEBRAIC[90] - 1.00000)*CONSTANTS[4]+1.00000;
ALGEBRAIC[96] = (ALGEBRAIC[93]<0.00000 ? 0.00000 : ALGEBRAIC[93]);
ALGEBRAIC[98] = (CONSTANTS[6]>0.00000 ? CONSTANTS[6] : ALGEBRAIC[96]+CONSTANTS[5]);
ALGEBRAIC[101] = (STATES[30]<0.000100000 ? 0.000100000 : STATES[30]);
ALGEBRAIC[102] =  CONSTANTS[119]*ALGEBRAIC[101];
ALGEBRAIC[104] = ALGEBRAIC[102] - 38.0000;
ALGEBRAIC[106] =  ALGEBRAIC[104]*CONSTANTS[110]+1.00000;
ALGEBRAIC[107] =  CONSTANTS[113]*ALGEBRAIC[104]+1.00000;
ALGEBRAIC[113] = (STATES[34]<0.00000 ? 0.00000 : STATES[34]);
ALGEBRAIC[114] =  ALGEBRAIC[113]*0.486110;
ALGEBRAIC[116] =  (ALGEBRAIC[114] - 10.0000)*CONSTANTS[108];
ALGEBRAIC[117] = ALGEBRAIC[114] - CONSTANTS[120];
ALGEBRAIC[119] =  ALGEBRAIC[117]*CONSTANTS[121]+1.00000;
ALGEBRAIC[120] =  CONSTANTS[123]*ALGEBRAIC[117]+1.00000;
ALGEBRAIC[121] =  CONSTANTS[125]*ALGEBRAIC[117]+1.00000;
ALGEBRAIC[123] = (STATES[35]<0.00100000 ? 0.00100000 : STATES[35]);
ALGEBRAIC[124] = (ALGEBRAIC[81] - STATES[7]) - ALGEBRAIC[123];
ALGEBRAIC[27] = STATES[8]/STATES[7];
ALGEBRAIC[30] =  0.280000*ALGEBRAIC[27]+ 0.00190000*pow(ALGEBRAIC[27], 2.00000);
ALGEBRAIC[38] = STATES[19]+ ALGEBRAIC[31]*0.155000;
ALGEBRAIC[39] = ALGEBRAIC[38] - 0.306250;
ALGEBRAIC[40] = ALGEBRAIC[39]/0.00480000;
ALGEBRAIC[41] =  0.0260000*ALGEBRAIC[40];
ALGEBRAIC[42] = (ALGEBRAIC[41]<1.00000e-05 ? 1.00000e-05 : ALGEBRAIC[41]);
ALGEBRAIC[43] = pow(ALGEBRAIC[42], 0.500000);
ALGEBRAIC[44] = 1.00000/ALGEBRAIC[43];
ALGEBRAIC[51] = ALGEBRAIC[47]+18.0000;
ALGEBRAIC[53] = 1.00000/( ALGEBRAIC[51]*0.0357000);
ALGEBRAIC[122] = ( (ALGEBRAIC[40] - ALGEBRAIC[47])*ALGEBRAIC[53])/(ALGEBRAIC[53]+ALGEBRAIC[44])+ALGEBRAIC[47];
ALGEBRAIC[142] = (STATES[36]<0.0250000 ? 0.0250000 : STATES[36]);
ALGEBRAIC[143] = ALGEBRAIC[142]/ALGEBRAIC[123];
ALGEBRAIC[144] =  ALGEBRAIC[143]*0.400000;
ALGEBRAIC[141] = 2.00000 - 0.150000/ALGEBRAIC[123];
ALGEBRAIC[145] =  (((ALGEBRAIC[122] - ALGEBRAIC[141])+ALGEBRAIC[144]) - ALGEBRAIC[30])*CONSTANTS[129];
ALGEBRAIC[147] =  (ALGEBRAIC[141]+11.0000)*0.000300000;
ALGEBRAIC[148] = ALGEBRAIC[145] - ALGEBRAIC[147];
ALGEBRAIC[149] = ALGEBRAIC[148];
ALGEBRAIC[150] =  ALGEBRAIC[147]*ALGEBRAIC[143];
ALGEBRAIC[146] =  (ALGEBRAIC[27] - ALGEBRAIC[143])*0.000225000;
ALGEBRAIC[151] = ALGEBRAIC[146] - ALGEBRAIC[150];
ALGEBRAIC[152] = ALGEBRAIC[151];
ALGEBRAIC[18] = CONSTANTS[56] - (CONSTANTS[56] - 1.00000)/exp( CONSTANTS[58]*(STATES[6] - 1.00000));
ALGEBRAIC[19] = (ALGEBRAIC[18]<CONSTANTS[57] ? CONSTANTS[57] : ALGEBRAIC[18]);
ALGEBRAIC[21] = ALGEBRAIC[19] - 1.00000;
ALGEBRAIC[22] =  ALGEBRAIC[21]*CONSTANTS[61]+1.00000;
ALGEBRAIC[88] = (STATES[26]>1.00000 ? 1.00000 : STATES[26]);
ALGEBRAIC[57] =  (ALGEBRAIC[47]+4.00000)*( CONSTANTS[85]*(ALGEBRAIC[19] - 1.00000)+1.00000) - 4.00000;
ALGEBRAIC[58] = (ALGEBRAIC[57]<=- 2.00000 ? 0.0100000 : ALGEBRAIC[57]>- 2.00000&&ALGEBRAIC[57]<=1.00000 ? 0.0100000+( (3.60000 - 0.0100000)*(ALGEBRAIC[57] - - 2.00000))/(1.00000 - - 2.00000) : ALGEBRAIC[57]>1.00000&&ALGEBRAIC[57]<=5.00000 ? 3.60000+( (9.40000 - 3.60000)*(ALGEBRAIC[57] - 1.00000))/(5.00000 - 1.00000) : ALGEBRAIC[57]>5.00000&&ALGEBRAIC[57]<=8.00000 ? 9.40000+( (11.6000 - 9.40000)*(ALGEBRAIC[57] - 5.00000))/(8.00000 - 5.00000) : ALGEBRAIC[57]>8.00000&&ALGEBRAIC[57]<=12.0000 ? 11.6000+( (13.5000 - 11.6000)*(ALGEBRAIC[57] - 8.00000))/(12.0000 - 8.00000) : 13.5000);
ALGEBRAIC[155] = STATES[7]+STATES[13];
ALGEBRAIC[156] = STATES[13]/ALGEBRAIC[155];
ALGEBRAIC[157] =  100.000*ALGEBRAIC[156];
ALGEBRAIC[158] = (STATES[29]/ALGEBRAIC[157])/5.25000;
ALGEBRAIC[161] = ALGEBRAIC[64]/( ALGEBRAIC[22]*ALGEBRAIC[158]);
ALGEBRAIC[162] = (ALGEBRAIC[161]<=0.00000 ? 1.04000 : ALGEBRAIC[161]>0.00000&&ALGEBRAIC[161]<=60.0000 ? 1.04000+( (1.02500 - 1.04000)*(ALGEBRAIC[161] - 0.00000))/(60.0000 - 0.00000) : ALGEBRAIC[161]>60.0000&&ALGEBRAIC[161]<=125.000 ? 1.02500+( (0.970000 - 1.02500)*(ALGEBRAIC[161] - 60.0000))/(125.000 - 60.0000) : ALGEBRAIC[161]>125.000&&ALGEBRAIC[161]<=160.000 ? 0.970000+( (0.880000 - 0.970000)*(ALGEBRAIC[161] - 125.000))/(160.000 - 125.000) : ALGEBRAIC[161]>160.000&&ALGEBRAIC[161]<=200.000 ? 0.880000+( (0.590000 - 0.880000)*(ALGEBRAIC[161] - 160.000))/(200.000 - 160.000) : ALGEBRAIC[161]>200.000&&ALGEBRAIC[161]<=240.000 ? 0.590000+( (0.00000 - 0.590000)*(ALGEBRAIC[161] - 200.000))/(240.000 - 200.000) : 0.00000);
ALGEBRAIC[163] =  ALGEBRAIC[162]*ALGEBRAIC[58]*ALGEBRAIC[22]*CONSTANTS[88]*ALGEBRAIC[88]*STATES[11];
ALGEBRAIC[73] = (ALGEBRAIC[47] - ALGEBRAIC[64])/3.00000;
ALGEBRAIC[164] = (ALGEBRAIC[73]>0.00000 ? ALGEBRAIC[163]+ALGEBRAIC[73] : ALGEBRAIC[163]);
ALGEBRAIC[54] = ALGEBRAIC[53]+ALGEBRAIC[44];
ALGEBRAIC[55] = ALGEBRAIC[40] - ALGEBRAIC[47];
ALGEBRAIC[56] = ALGEBRAIC[55]/ALGEBRAIC[54];
ALGEBRAIC[166] = ALGEBRAIC[56] - ALGEBRAIC[164];
ALGEBRAIC[159] = (ALGEBRAIC[40]/ALGEBRAIC[22])/ALGEBRAIC[158];
ALGEBRAIC[160] = (ALGEBRAIC[159]<=0.00000 ? 1.06000 : ALGEBRAIC[159]>0.00000&&ALGEBRAIC[159]<=32.0000 ? 1.06000+( (0.970000 - 1.06000)*(ALGEBRAIC[159] - 0.00000))/(32.0000 - 0.00000) : ALGEBRAIC[159]>32.0000&&ALGEBRAIC[159]<=38.4000 ? 0.970000+( (0.930000 - 0.970000)*(ALGEBRAIC[159] - 32.0000))/(38.4000 - 32.0000) : ALGEBRAIC[159]>38.4000&&ALGEBRAIC[159]<=48.0000 ? 0.930000+( (0.800000 - 0.930000)*(ALGEBRAIC[159] - 38.4000))/(48.0000 - 38.4000) : ALGEBRAIC[159]>48.0000&&ALGEBRAIC[159]<=60.8000 ? 0.800000+( (0.460000 - 0.800000)*(ALGEBRAIC[159] - 48.0000))/(60.8000 - 48.0000) : ALGEBRAIC[159]>60.8000&&ALGEBRAIC[159]<=72.0000 ? 0.460000+( (0.00000 - 0.460000)*(ALGEBRAIC[159] - 60.8000))/(72.0000 - 60.8000) : 0.00000);
ALGEBRAIC[165] =  (1.00000 - CONSTANTS[86])*ALGEBRAIC[22]*ALGEBRAIC[160]*CONSTANTS[87]*ALGEBRAIC[88]*STATES[12]+( CONSTANTS[86]*ALGEBRAIC[164])/ALGEBRAIC[58];
ALGEBRAIC[36] =  (ALGEBRAIC[34]+8.00000)*( CONSTANTS[85]*(ALGEBRAIC[19] - 1.00000)+1.00000) - 8.00000;
ALGEBRAIC[37] = (ALGEBRAIC[36]<=- 8.00000 ? 0.00000 : ALGEBRAIC[36]>- 8.00000&&ALGEBRAIC[36]<=- 6.00000 ? 0.00000+( (0.750000 - 0.00000)*(ALGEBRAIC[36] - - 8.00000))/(- 6.00000 - - 8.00000) : ALGEBRAIC[36]>- 6.00000&&ALGEBRAIC[36]<=- 2.00000 ? 0.750000+( (2.60000 - 0.750000)*(ALGEBRAIC[36] - - 6.00000))/(- 2.00000 - - 6.00000) : ALGEBRAIC[36]>- 2.00000&&ALGEBRAIC[36]<=4.00000 ? 2.60000+( (9.80000 - 2.60000)*(ALGEBRAIC[36] - - 2.00000))/(4.00000 - - 2.00000) : ALGEBRAIC[36]>4.00000&&ALGEBRAIC[36]<=12.0000 ? 9.80000+( (13.5000 - 9.80000)*(ALGEBRAIC[36] - 4.00000))/(12.0000 - 4.00000) : 13.5000);
ALGEBRAIC[168] =  ALGEBRAIC[37]*ALGEBRAIC[165];
ALGEBRAIC[170] = ALGEBRAIC[168] - ALGEBRAIC[56];
ALGEBRAIC[169] = (ALGEBRAIC[158]>1.00000 ? 114.000+ (ALGEBRAIC[158] - 1.00000)*6667.00 : ALGEBRAIC[158]>0.936000&&ALGEBRAIC[158]<=1.00000 ? 74.0000+ (ALGEBRAIC[158] - 0.936000)*625.000 : ALGEBRAIC[158]>0.800000&&ALGEBRAIC[158]<=0.936000 ? 46.0000+ (ALGEBRAIC[158] - 0.800000)*205.882 :  ALGEBRAIC[158]*57.5000);
ALGEBRAIC[172] = (ALGEBRAIC[169] - 67.0000)/30.0000;
ALGEBRAIC[174] = (ALGEBRAIC[172]>1.00000 ? 1.00000 : ALGEBRAIC[172]<0.600000 ? 0.600000 : ALGEBRAIC[172]);
ALGEBRAIC[176] = 1.00000/ALGEBRAIC[174];
ALGEBRAIC[179] =  ( (ALGEBRAIC[176] - 1.00000)*3.00000 - STATES[37])*0.000500000;
ALGEBRAIC[71] = (ALGEBRAIC[65]<80.0000&&ALGEBRAIC[65]>=40.0000 ?  0.00500000*(80.0000 - ALGEBRAIC[65])*CONSTANTS[44] : ALGEBRAIC[65]<40.0000 ?  0.200000*CONSTANTS[44] : 0.00000);
ALGEBRAIC[171] = (ALGEBRAIC[169]<80.0000&&ALGEBRAIC[169]>=40.0000 ?  CONSTANTS[45]*(80.0000 - ALGEBRAIC[169]) : ALGEBRAIC[169]<40.0000 ?  CONSTANTS[45]*40.0000 : 0.00000);
ALGEBRAIC[173] = ALGEBRAIC[71]+ALGEBRAIC[171];
ALGEBRAIC[17] = STATES[5];
ALGEBRAIC[72] = (ALGEBRAIC[65]<40.0000 ?  0.0400000*(40.0000 - ALGEBRAIC[65])*CONSTANTS[50] : 0.00000);
ALGEBRAIC[49] =  (15.0000/(ALGEBRAIC[47]+ALGEBRAIC[34]+ALGEBRAIC[40]) - 1.00000)*CONSTANTS[51]+1.00000;
ALGEBRAIC[175] =  ( CONSTANTS[206]*ALGEBRAIC[49]*(ALGEBRAIC[173]+ALGEBRAIC[17]+ALGEBRAIC[72]) - 1.00000)*CONSTANTS[54]+1.00000;
ALGEBRAIC[177] = (ALGEBRAIC[175]<0.00000 ? 0.00000 : ALGEBRAIC[175]);
ALGEBRAIC[180] = (ALGEBRAIC[177] - STATES[6])/CONSTANTS[55];
ALGEBRAIC[189] =  pow(fabs(ALGEBRAIC[34]), CONSTANTS[153])*CONSTANTS[154];
ALGEBRAIC[192] = (ALGEBRAIC[34]<0.00000 ? - ALGEBRAIC[189] : ALGEBRAIC[189]);
ALGEBRAIC[184] =  (CONSTANTS[208] - 40.0000)*ALGEBRAIC[157];
ALGEBRAIC[186] = (ALGEBRAIC[184]+CONSTANTS[209]<0.00000 ? 0.00000 : ALGEBRAIC[184]+CONSTANTS[209]);
ALGEBRAIC[188] = CONSTANTS[137] - ALGEBRAIC[186];
ALGEBRAIC[190] = ( ALGEBRAIC[188]*CONSTANTS[138]*CONSTANTS[139]+5.00000e-06<0.00000 ? 0.00000 :  ALGEBRAIC[188]*CONSTANTS[138]*CONSTANTS[139]+5.00000e-06);
ALGEBRAIC[178] = ALGEBRAIC[157]/( (CONSTANTS[134] - ALGEBRAIC[157])*CONSTANTS[135]);
ALGEBRAIC[181] = ALGEBRAIC[178]+1.50000;
ALGEBRAIC[182] =  0.333300*ALGEBRAIC[181];
ALGEBRAIC[193] =  STATES[13]*CONSTANTS[140]*ALGEBRAIC[182];
ALGEBRAIC[196] = (ALGEBRAIC[190] - ALGEBRAIC[193])+CONSTANTS[141];
ALGEBRAIC[194] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[150]*0.00100000;
ALGEBRAIC[183] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[146]+1.00000;
ALGEBRAIC[185] =  pow(CONSTANTS[147] - ALGEBRAIC[114], 2.00000)*CONSTANTS[148]*ALGEBRAIC[183];
ALGEBRAIC[187] = (ALGEBRAIC[185]<0.800000 ? 0.800000 : ALGEBRAIC[185]>8.00000 ? 8.00000 : ALGEBRAIC[185]);
ALGEBRAIC[191] =  (STATES[3] - 1.00000)*CONSTANTS[149]+1.00000;
ALGEBRAIC[197] =  ALGEBRAIC[191]*ALGEBRAIC[187]*0.00100000;
ALGEBRAIC[199] = (ALGEBRAIC[197]<0.00000 ? 0.00000 : ALGEBRAIC[197]);
ALGEBRAIC[201] =  (ALGEBRAIC[194]+ALGEBRAIC[199])*1.00000;
ALGEBRAIC[203] = (ALGEBRAIC[201]<0.00000 ? 0.00000 : ALGEBRAIC[201]);
ALGEBRAIC[20] = (CONSTANTS[59] -  ALGEBRAIC[19]*CONSTANTS[60])+CONSTANTS[60];
ALGEBRAIC[195] = ALGEBRAIC[192] - STATES[38];
ALGEBRAIC[205] =  ALGEBRAIC[195]*CONSTANTS[157];
ALGEBRAIC[59] = STATES[16]+ ALGEBRAIC[31]*0.398600;
ALGEBRAIC[8] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[25]+1.00000;
ALGEBRAIC[9] = (ALGEBRAIC[8]<CONSTANTS[26] ? CONSTANTS[26] : ALGEBRAIC[8]);
ALGEBRAIC[60] =  (ALGEBRAIC[9] - 1.00000)*CONSTANTS[89];
ALGEBRAIC[206] = ((((ALGEBRAIC[59] - ALGEBRAIC[20]) - ALGEBRAIC[60]) - STATES[15]) - STATES[14]) - ALGEBRAIC[205];
ALGEBRAIC[207] = (ALGEBRAIC[206]<0.000100000 ? 0.000100000 : ALGEBRAIC[206]);
ALGEBRAIC[208] = 3.70000+(ALGEBRAIC[207] - 0.740000)/CONSTANTS[90];
ALGEBRAIC[209] = (ALGEBRAIC[208]<0.000100000 ? 0.000100000 : ALGEBRAIC[208]);
ALGEBRAIC[210] = 0.740000/pow(ALGEBRAIC[209]/( ALGEBRAIC[182]*3.70000), 0.500000);
ALGEBRAIC[61] = (ALGEBRAIC[34]<CONSTANTS[91] ? CONSTANTS[91] : ALGEBRAIC[34]);
ALGEBRAIC[211] = ALGEBRAIC[209] - ALGEBRAIC[61];
ALGEBRAIC[212] = ALGEBRAIC[211]/ALGEBRAIC[210];
ALGEBRAIC[214] = ALGEBRAIC[212] - ALGEBRAIC[168];
ALGEBRAIC[213] = ALGEBRAIC[64] - ALGEBRAIC[209];
ALGEBRAIC[99] = (STATES[27]<CONSTANTS[112] ? CONSTANTS[112] : STATES[27]);
ALGEBRAIC[100] =  ALGEBRAIC[99]*STATES[28];
ALGEBRAIC[11] = CONSTANTS[34] - (CONSTANTS[34] - 1.00000)/( ((CONSTANTS[35] - 1.00000)/(CONSTANTS[35] - CONSTANTS[34]))*(STATES[3] - 1.00000)+1.00000);
ALGEBRAIC[12] = (ALGEBRAIC[11]<CONSTANTS[35] ? CONSTANTS[35] : ALGEBRAIC[11]);
ALGEBRAIC[25] = pow( ALGEBRAIC[21]*CONSTANTS[64]+1.00000, CONSTANTS[65]);
ALGEBRAIC[198] =  ALGEBRAIC[195]*CONSTANTS[156]+1.00000;
ALGEBRAIC[75] = pow(ALGEBRAIC[64]/100.000, CONSTANTS[95]);
ALGEBRAIC[200] = (( ALGEBRAIC[9]*ALGEBRAIC[12]*ALGEBRAIC[25]*ALGEBRAIC[182]*CONSTANTS[1])/ALGEBRAIC[75])/ALGEBRAIC[198];
ALGEBRAIC[204] =  CONSTANTS[98]*ALGEBRAIC[100]*ALGEBRAIC[200]*CONSTANTS[0]*CONSTANTS[97];
ALGEBRAIC[215] = ALGEBRAIC[213]/ALGEBRAIC[204];
ALGEBRAIC[216] =  STATES[29]*ALGEBRAIC[215];
ALGEBRAIC[24] =  ALGEBRAIC[21]*CONSTANTS[63]+1.00000;
ALGEBRAIC[103] = (ALGEBRAIC[102]>38.0000 ? 38.0000 : ALGEBRAIC[102]);
ALGEBRAIC[105] =  ALGEBRAIC[24]*CONSTANTS[117]*CONSTANTS[115]*(1.00000 - pow(38.0001 - ALGEBRAIC[103], 3.00000)/54872.0);
ALGEBRAIC[261] = ALGEBRAIC[229] - ALGEBRAIC[105];
ALGEBRAIC[10] =  (ALGEBRAIC[9] - 1.00000)*CONSTANTS[27]+1.00000;
ALGEBRAIC[15] = CONSTANTS[41] - CONSTANTS[41]/( 0.555556*(1.00000+STATES[4]));
ALGEBRAIC[16] = (ALGEBRAIC[15]<- 1.00000 ? - 1.00000 : ALGEBRAIC[15]);
ALGEBRAIC[26] =  ALGEBRAIC[21]*CONSTANTS[66]+1.00000;
ALGEBRAIC[108] = (STATES[31]<0.500000 ? 0.500000 : STATES[31]);
ALGEBRAIC[109] = (STATES[32]<0.500000 ? 0.500000 : STATES[32]);
ALGEBRAIC[110] = (STATES[33]<0.300000 ? 0.300000 : STATES[33]);
ALGEBRAIC[111] =  ALGEBRAIC[108]*ALGEBRAIC[109]*ALGEBRAIC[110];
ALGEBRAIC[112] =  (ALGEBRAIC[111] - 1.00000)*CONSTANTS[127]+1.00000;
ALGEBRAIC[202] =  CONSTANTS[96]*ALGEBRAIC[112]*ALGEBRAIC[200]*CONSTANTS[0]*CONSTANTS[97];
ALGEBRAIC[217] = (CONSTANTS[159]>0.00000&&CONSTANTS[160]<=0.00000 ? CONSTANTS[159] : CONSTANTS[160]>0.00000 ? STATES[39] : ALGEBRAIC[64] - CONSTANTS[158]);
ALGEBRAIC[218] =  (ALGEBRAIC[25] - 1.00000)*CONSTANTS[167]+1.00000;
ALGEBRAIC[219] = (ALGEBRAIC[218]<0.800000 ? 0.800000 : ALGEBRAIC[218]);
ALGEBRAIC[220] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[168]+1.00000;
ALGEBRAIC[221] = (ALGEBRAIC[220]<CONSTANTS[169] ? CONSTANTS[169] : ALGEBRAIC[220]);
ALGEBRAIC[222] =  (ALGEBRAIC[219] - 1.00000)*CONSTANTS[173]+1.00000;
ALGEBRAIC[223] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[174]+1.00000;
ALGEBRAIC[262] = ALGEBRAIC[232];
ALGEBRAIC[268] = (ALGEBRAIC[262]>1.00000 ? 1.00000/(1.00000+ (ALGEBRAIC[262] - 1.00000)*72.0000) : 10.0000 - 9.00000/(1.00000+ (1.00000 - ALGEBRAIC[262])*8.00000));
ALGEBRAIC[275] =  (ALGEBRAIC[268] - 1.00000)*CONSTANTS[16];
ALGEBRAIC[74] = ALGEBRAIC[64] - ALGEBRAIC[34];
ALGEBRAIC[76] =  ALGEBRAIC[74]*CONSTANTS[99];
ALGEBRAIC[264] =  CONSTANTS[178]*ALGEBRAIC[250];
ALGEBRAIC[270] = ALGEBRAIC[215]+ALGEBRAIC[234]+ALGEBRAIC[264];
ALGEBRAIC[279] = ALGEBRAIC[270]+ALGEBRAIC[76];
ALGEBRAIC[125] = ALGEBRAIC[124] - STATES[9];
ALGEBRAIC[133] =  pow((ALGEBRAIC[125] - CONSTANTS[83])/CONSTANTS[83], 2.00000)*1.00000;
ALGEBRAIC[130] = pow((CONSTANTS[79]/ALGEBRAIC[124])/5.00000, CONSTANTS[80]);
ALGEBRAIC[134] =  ALGEBRAIC[130]*CONSTANTS[81]+ALGEBRAIC[133];
ALGEBRAIC[127] = STATES[10]/ALGEBRAIC[124];
ALGEBRAIC[129] =  0.280000*ALGEBRAIC[127]+ 0.00190000*pow(ALGEBRAIC[127], 2.00000);
ALGEBRAIC[131] =  ALGEBRAIC[130]*2.00000;
ALGEBRAIC[132] =  ALGEBRAIC[131]*ALGEBRAIC[129]*CONSTANTS[82];
ALGEBRAIC[263] = (ALGEBRAIC[235] - CONSTANTS[74]<0.00000 ? 0.00000 : ALGEBRAIC[235] - CONSTANTS[74]);
ALGEBRAIC[269] = pow( ALGEBRAIC[263]*CONSTANTS[75], CONSTANTS[76]);
ALGEBRAIC[277] =  (((ALGEBRAIC[235] - ALGEBRAIC[30]) - ALGEBRAIC[134])+ALGEBRAIC[132])*CONSTANTS[68]+ALGEBRAIC[269];
ALGEBRAIC[135] = ALGEBRAIC[134] - ALGEBRAIC[131];
ALGEBRAIC[136] = (ALGEBRAIC[135]+CONSTANTS[84]) - ALGEBRAIC[133];
ALGEBRAIC[138] = (ALGEBRAIC[136]>7.00000 ? 7.00000 : ALGEBRAIC[136]);
ALGEBRAIC[139] = (ALGEBRAIC[138]<0.00000 ? 0.00000 :  ALGEBRAIC[138]*0.0200000);
ALGEBRAIC[284] = ((ALGEBRAIC[139] - ALGEBRAIC[277]) - ALGEBRAIC[149])+CONSTANTS[69];
ALGEBRAIC[28] = (ALGEBRAIC[27] - CONSTANTS[70]<0.00000 ? 0.00000 : ALGEBRAIC[27] - CONSTANTS[70]);
ALGEBRAIC[29] = CONSTANTS[71] -  pow(ALGEBRAIC[28], CONSTANTS[72])*CONSTANTS[73];
ALGEBRAIC[278] =  ALGEBRAIC[269]*ALGEBRAIC[27]+ (ALGEBRAIC[27] - ALGEBRAIC[127])*0.00104000;
ALGEBRAIC[140] =  ALGEBRAIC[127]*ALGEBRAIC[139];
ALGEBRAIC[285] = (((ALGEBRAIC[29]+ALGEBRAIC[140]) - ALGEBRAIC[278]) - ALGEBRAIC[152])+ CONSTANTS[69]*72.0000;
ALGEBRAIC[286] = ALGEBRAIC[278] - ALGEBRAIC[140];
ALGEBRAIC[287] = ALGEBRAIC[279] - ALGEBRAIC[212];
ALGEBRAIC[288] = ALGEBRAIC[164] - ALGEBRAIC[279];
ALGEBRAIC[265] =  STATES[29]*ALGEBRAIC[234];
ALGEBRAIC[115] = (ALGEBRAIC[114]>35.0000 ? 35.0000 : ALGEBRAIC[114]);
ALGEBRAIC[118] =  ALGEBRAIC[24]*CONSTANTS[128]*(1.00000 - pow(35.0001 - ALGEBRAIC[115], 3.00000)/42875.0);
ALGEBRAIC[280] = ALGEBRAIC[271] - ALGEBRAIC[118];
ALGEBRAIC[290] = (ALGEBRAIC[113]<6.00000&&ALGEBRAIC[280]<0.00000 ?  ALGEBRAIC[280]*0.100000 : ALGEBRAIC[280]);
ALGEBRAIC[276] =  (ALGEBRAIC[268]+STATES[1])*CONSTANTS[18];
ALGEBRAIC[283] = (ALGEBRAIC[276]<1.00000e-05 ? 1.00000e-05 : ALGEBRAIC[276]);
ALGEBRAIC[293] = (CONSTANTS[19]>0.00000 ? CONSTANTS[19] : ALGEBRAIC[283]+CONSTANTS[20]);
ALGEBRAIC[281] = ALGEBRAIC[271]+ALGEBRAIC[229];
ALGEBRAIC[154] = STATES[37]+1.00000;
ALGEBRAIC[291] =  ALGEBRAIC[281]*CONSTANTS[130]*0.0266670*ALGEBRAIC[176]*ALGEBRAIC[154];
ALGEBRAIC[294] = CONSTANTS[131] - (ALGEBRAIC[281]/ALGEBRAIC[291])/0.761000;
ALGEBRAIC[153] = CONSTANTS[132]/(CONSTANTS[133]+ALGEBRAIC[123]);
ALGEBRAIC[296] =  (ALGEBRAIC[294] - ALGEBRAIC[169])*ALGEBRAIC[153];
ALGEBRAIC[298] = (ALGEBRAIC[296] - ALGEBRAIC[281])/( ALGEBRAIC[168]*1.00000);
ALGEBRAIC[301] =  ALGEBRAIC[232]*ALGEBRAIC[82]*0.00616190;
ALGEBRAIC[4] =  (ALGEBRAIC[1] - 1.00000)*CONSTANTS[13]+1.00000;
ALGEBRAIC[5] = (ALGEBRAIC[4]<CONSTANTS[15] ? CONSTANTS[15] : ALGEBRAIC[4]>CONSTANTS[14] ? CONSTANTS[14] : ALGEBRAIC[4]);
ALGEBRAIC[13] = CONSTANTS[37] - (CONSTANTS[37] - 1.00000)/( ((CONSTANTS[36] - 1.00000)/(CONSTANTS[36] - CONSTANTS[37]))*(STATES[3] - 1.00000)+1.00000);
ALGEBRAIC[14] = (ALGEBRAIC[13]<CONSTANTS[36] ? CONSTANTS[36] : ALGEBRAIC[13]);
ALGEBRAIC[267] =  ALGEBRAIC[252]*CONSTANTS[185] - CONSTANTS[186];
ALGEBRAIC[274] = (ALGEBRAIC[267]<1.00000 ? 1.00000 : ALGEBRAIC[267]);
ALGEBRAIC[282] =  ( (ALGEBRAIC[250] - 1.20000)*CONSTANTS[188]+1.20000)*CONSTANTS[189];
ALGEBRAIC[292] = ((ALGEBRAIC[252]+CONSTANTS[190]) - ALGEBRAIC[282]) - ALGEBRAIC[274];
ALGEBRAIC[295] =  ALGEBRAIC[292]*CONSTANTS[191];
ALGEBRAIC[297] = ALGEBRAIC[295];
ALGEBRAIC[299] =  (ALGEBRAIC[297] - 1.00000)*CONSTANTS[193]+1.00000;
ALGEBRAIC[300] = (ALGEBRAIC[299]<0.000100000 ? 0.000100000 : ALGEBRAIC[299]);
ALGEBRAIC[302] =  (( ALGEBRAIC[5]*ALGEBRAIC[300]*CONSTANTS[194])/CONSTANTS[195])*( (ALGEBRAIC[14] - 1.00000)*CONSTANTS[196]+1.00000);
ALGEBRAIC[303] = (ALGEBRAIC[302]<CONSTANTS[197] ? CONSTANTS[197] : ALGEBRAIC[302]);
ALGEBRAIC[225] =  ( (ALGEBRAIC[6] - 1.00000)*CONSTANTS[198]+1.00000)*0.100000;
ALGEBRAIC[226] = (ALGEBRAIC[225]<0.00000 ? 0.00000 : ALGEBRAIC[225]);
ALGEBRAIC[309] = (ALGEBRAIC[301] - ALGEBRAIC[303]) - ALGEBRAIC[226];
ALGEBRAIC[310] = (ALGEBRAIC[309]<1.00000e-08 ? 1.00000e-08 : ALGEBRAIC[309]);
ALGEBRAIC[311] =  ALGEBRAIC[310]*CONSTANTS[178];
ALGEBRAIC[313] = ( CONSTANTS[100]*ALGEBRAIC[187] - ALGEBRAIC[311])+ CONSTANTS[101]*142.000;
ALGEBRAIC[304] = ( ALGEBRAIC[301]*ALGEBRAIC[84])/ALGEBRAIC[82];
ALGEBRAIC[305] =  (ALGEBRAIC[300] - 1.00000)*CONSTANTS[199];
ALGEBRAIC[306] =  (ALGEBRAIC[232] - 1.00000)*CONSTANTS[200]+1.00000;
ALGEBRAIC[307] = (ALGEBRAIC[306]<0.100000 ? 0.100000 : ALGEBRAIC[306]);
ALGEBRAIC[227] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[201]+1.00000;
ALGEBRAIC[228] = (ALGEBRAIC[227]<CONSTANTS[202] ? CONSTANTS[202] : ALGEBRAIC[227]);
ALGEBRAIC[308] = ( pow(ALGEBRAIC[84]/4.40000, CONSTANTS[203])*ALGEBRAIC[3]*0.0800000*ALGEBRAIC[307])/ALGEBRAIC[228];
ALGEBRAIC[312] = ((ALGEBRAIC[304]+ALGEBRAIC[308]) - STATES[42]) - ALGEBRAIC[305];
ALGEBRAIC[314] = (ALGEBRAIC[312]<0.00000 ? 0.00000 : ALGEBRAIC[312]);
ALGEBRAIC[315] =  ALGEBRAIC[314]*CONSTANTS[178];
ALGEBRAIC[317] = CONSTANTS[103] - ALGEBRAIC[315];
ALGEBRAIC[224] = STATES[41]/STATES[22];
ALGEBRAIC[316] =  pow(ALGEBRAIC[251], 2.00000)*ALGEBRAIC[224]*3.84000;
ALGEBRAIC[318] =  ALGEBRAIC[316]*CONSTANTS[178];
ALGEBRAIC[319] = ALGEBRAIC[316]+ 2.00000*(ALGEBRAIC[310]+ALGEBRAIC[314]);
ALGEBRAIC[320] = (ALGEBRAIC[319]>0.600000 ? 0.600000 : ALGEBRAIC[319]);
ALGEBRAIC[321] = ALGEBRAIC[319] - 0.600000;
ALGEBRAIC[322] = (ALGEBRAIC[321]<0.00000 ? 0.00000 : ALGEBRAIC[321]);
ALGEBRAIC[323] = ALGEBRAIC[320]/( 600.000*ALGEBRAIC[14])+ALGEBRAIC[322]/360.000;
ALGEBRAIC[324] =  ALGEBRAIC[323]*CONSTANTS[178];
ALGEBRAIC[23] =  ALGEBRAIC[21]*CONSTANTS[62]+1.00000;
ALGEBRAIC[92] =  72.0000*ALGEBRAIC[23];
ALGEBRAIC[95] =  (ALGEBRAIC[88] - 1.00000)*0.500000+1.00000;
ALGEBRAIC[97] =  (ALGEBRAIC[92]+CONSTANTS[207])*ALGEBRAIC[95];
ALGEBRAIC[126] = (ALGEBRAIC[124]<=0.00000 ? 0.00000 : ALGEBRAIC[124]>0.00000&&ALGEBRAIC[124]<=12.0000 ? 0.00000+( (11.4000 - 0.00000)*(ALGEBRAIC[124] - 0.00000))/(12.0000 - 0.00000) : ALGEBRAIC[124]>12.0000&&ALGEBRAIC[124]<=15.0000 ? 11.4000+( (14.0000 - 11.4000)*(ALGEBRAIC[124] - 12.0000))/(15.0000 - 12.0000) : ALGEBRAIC[124]>15.0000&&ALGEBRAIC[124]<=18.0000 ? 14.0000+( (16.0000 - 14.0000)*(ALGEBRAIC[124] - 15.0000))/(18.0000 - 15.0000) : ALGEBRAIC[124]>18.0000&&ALGEBRAIC[124]<=21.0000 ? 16.0000+( (17.3000 - 16.0000)*(ALGEBRAIC[124] - 18.0000))/(21.0000 - 18.0000) : ALGEBRAIC[124]>21.0000&&ALGEBRAIC[124]<=24.0000 ? 17.3000+( (18.0000 - 17.3000)*(ALGEBRAIC[124] - 21.0000))/(24.0000 - 21.0000) : 18.0000);
ALGEBRAIC[128] = ALGEBRAIC[124] - ALGEBRAIC[126];
ALGEBRAIC[137] = ALGEBRAIC[133] - ALGEBRAIC[135];
ALGEBRAIC[167] = ALGEBRAIC[164]/ALGEBRAIC[97];
ALGEBRAIC[266] =  ALGEBRAIC[251]*CONSTANTS[178];
ALGEBRAIC[289] = ALGEBRAIC[74]/ALGEBRAIC[279];
}

void objfunc_0(double *p, double *hx, int m, int n, void *adata)
{
  struct rootfind_info* rfi = (struct rootfind_info*)adata;
#define VOI rfi->aVOI
#define CONSTANTS rfi->aCONSTANTS
#define RATES rfi->aRATES
#define STATES rfi->aSTATES
#define ALGEBRAIC rfi->aALGEBRAIC
#define pret rfi->aPRET
  ALGEBRAIC[229] = p[0];
  ALGEBRAIC[230] = p[1];
  ALGEBRAIC[231] = p[2];
  hx[0] = ALGEBRAIC[231] - (ALGEBRAIC[216] - ALGEBRAIC[229])/( ALGEBRAIC[157]*5.25000*ALGEBRAIC[215]);
  hx[1] = ALGEBRAIC[230] -  57.1400*ALGEBRAIC[231]*pow(CONSTANTS[115], CONSTANTS[116]);
  hx[2] = ALGEBRAIC[229] -  (ALGEBRAIC[230] - ALGEBRAIC[102])*CONSTANTS[118]*ALGEBRAIC[215];
#undef VOI
#undef CONSTANTS
#undef RATES
#undef STATES
#undef ALGEBRAIC
#undef pret
}
void rootfind_0(double VOI, double* CONSTANTS, double* RATES,
double* STATES, double* ALGEBRAIC, int* pret)
{
  static double p[3] = {56.8057,39.9793,0.699673};
  double bp[3], work[LM_DIF_WORKSZ(3, 3)];
  struct rootfind_info rfi;
  rfi.aVOI = VOI;
  rfi.aCONSTANTS = CONSTANTS;
  rfi.aRATES = RATES;
  rfi.aSTATES = STATES;
  rfi.aALGEBRAIC = ALGEBRAIC;
  rfi.aPRET = pret;
  do_levmar(objfunc_0, p, bp, work, pret, 3, &rfi);
  ALGEBRAIC[229] = p[0];
  ALGEBRAIC[230] = p[1];
  ALGEBRAIC[231] = p[2];
}

void objfunc_1(double *p, double *hx, int m, int n, void *adata)
{
  struct rootfind_info* rfi = (struct rootfind_info*)adata;
#define VOI rfi->aVOI
#define CONSTANTS rfi->aCONSTANTS
#define RATES rfi->aRATES
#define STATES rfi->aSTATES
#define ALGEBRAIC rfi->aALGEBRAIC
#define pret rfi->aPRET
  ALGEBRAIC[232] = p[0];
  ALGEBRAIC[233] = p[1];
  ALGEBRAIC[234] = p[2];
  ALGEBRAIC[235] = p[3];
  ALGEBRAIC[236] = p[4];
  ALGEBRAIC[237] = p[5];
  ALGEBRAIC[238] = p[6];
  ALGEBRAIC[239] = p[7];
  ALGEBRAIC[240] = p[8];
  ALGEBRAIC[241] = p[9];
  ALGEBRAIC[242] = p[10];
  ALGEBRAIC[243] = p[11];
  ALGEBRAIC[244] = p[12];
  ALGEBRAIC[245] = p[13];
  ALGEBRAIC[246] = p[14];
  ALGEBRAIC[247] = p[15];
  ALGEBRAIC[248] = p[16];
  ALGEBRAIC[249] = p[17];
  ALGEBRAIC[250] = p[18];
  ALGEBRAIC[251] = p[19];
  ALGEBRAIC[252] = p[20];
  ALGEBRAIC[253] = p[21];
  ALGEBRAIC[254] = p[22];
  ALGEBRAIC[255] = p[23];
  ALGEBRAIC[256] = p[24];
  ALGEBRAIC[257] = p[25];
  ALGEBRAIC[258] = p[26];
  ALGEBRAIC[259] = p[27];
  ALGEBRAIC[260] = p[28];
  hx[0] = ALGEBRAIC[235] - ( ALGEBRAIC[233]*1.70000*ALGEBRAIC[234]+ALGEBRAIC[209]);
  hx[1] = ALGEBRAIC[236] -  ( (ALGEBRAIC[235] - 17.0000)*CONSTANTS[92]+17.0000)*CONSTANTS[93];
  hx[2] = ALGEBRAIC[237] - CONSTANTS[94]/ALGEBRAIC[236];
  hx[3] = ALGEBRAIC[233] -  ALGEBRAIC[26]*ALGEBRAIC[237]*ALGEBRAIC[182]*ALGEBRAIC[10];
  hx[4] = ALGEBRAIC[238] -  ALGEBRAIC[233]*1.79000;
  hx[5] = ALGEBRAIC[239] - (ALGEBRAIC[202]+ALGEBRAIC[238]);
  hx[6] = ALGEBRAIC[234] - ALGEBRAIC[213]/ALGEBRAIC[239];
  hx[7] = ALGEBRAIC[242] - ( (ALGEBRAIC[232] - 1.00000)*CONSTANTS[163]+1.00000);
  hx[8] = ALGEBRAIC[241] - (ALGEBRAIC[242]<CONSTANTS[164] ? CONSTANTS[164] : ALGEBRAIC[242]>CONSTANTS[165] ? CONSTANTS[165] : ALGEBRAIC[242]);
  hx[9] = ALGEBRAIC[240] - (ALGEBRAIC[241] - STATES[40]);
  hx[10] = ALGEBRAIC[243] -  CONSTANTS[170]*CONSTANTS[2]*ALGEBRAIC[219]*ALGEBRAIC[240]*ALGEBRAIC[221]*40.0000*CONSTANTS[0];
  hx[11] = ALGEBRAIC[245] - ((ALGEBRAIC[243] -  ALGEBRAIC[16]*CONSTANTS[171])+CONSTANTS[171]);
  hx[12] = ALGEBRAIC[244] - (ALGEBRAIC[245]<CONSTANTS[172] ? CONSTANTS[172] : ALGEBRAIC[245]);
  hx[13] = ALGEBRAIC[246] - ( (ALGEBRAIC[240] - 1.00000)*CONSTANTS[175]+1.00000);
  hx[14] = ALGEBRAIC[248] -  43.3330*CONSTANTS[176]*ALGEBRAIC[223]*ALGEBRAIC[246]*CONSTANTS[0]*ALGEBRAIC[222];
  hx[15] = ALGEBRAIC[247] - (ALGEBRAIC[248]<CONSTANTS[177] ? CONSTANTS[177] : ALGEBRAIC[248]);
  hx[16] = ALGEBRAIC[249] - (ALGEBRAIC[244]+ALGEBRAIC[247]);
  hx[17] = ALGEBRAIC[250] - ALGEBRAIC[217]/ALGEBRAIC[249];
  hx[18] = ALGEBRAIC[254] - ( ALGEBRAIC[250]*(1.00000 - ALGEBRAIC[156]))/( ALGEBRAIC[250]*(1.00000 - ALGEBRAIC[156]) - ALGEBRAIC[251]);
  hx[19] = ALGEBRAIC[253] - (ALGEBRAIC[254]<1.00000 ? 1.00000 : ALGEBRAIC[254]);
  hx[20] = ALGEBRAIC[252] - (CONSTANTS[180]>0.00000 ?  pow(ALGEBRAIC[253], 1.35000)*ALGEBRAIC[30]*0.980000 : ALGEBRAIC[30]+4.00000);
  hx[21] = ALGEBRAIC[256] -  ALGEBRAIC[244]*ALGEBRAIC[250];
  hx[22] = ALGEBRAIC[255] - (ALGEBRAIC[217] - ALGEBRAIC[256]);
  hx[23] = ALGEBRAIC[257] - ((ALGEBRAIC[255] - ALGEBRAIC[252]) - CONSTANTS[181]);
  hx[24] = ALGEBRAIC[258] -  ALGEBRAIC[257]*CONSTANTS[182];
  hx[25] = ALGEBRAIC[251] - (ALGEBRAIC[258]<CONSTANTS[183] ? CONSTANTS[183] : ALGEBRAIC[258]);
  hx[26] = ALGEBRAIC[259] -  ALGEBRAIC[251]*8.00000;
  hx[27] = ALGEBRAIC[260] - ( (ALGEBRAIC[259] - 1.00000)*CONSTANTS[184]+1.00000);
  hx[28] = ALGEBRAIC[232] - (ALGEBRAIC[260]<0.00000 ? 0.00000 : ALGEBRAIC[260]);
#undef VOI
#undef CONSTANTS
#undef RATES
#undef STATES
#undef ALGEBRAIC
#undef pret
}
void rootfind_1(double VOI, double* CONSTANTS, double* RATES,
double* STATES, double* ALGEBRAIC, int* pret)
{
  static double p[29] = {1.00051,2.77751,2.79521,16.9144,0.36019,2.77632,4.97174,35.7072,1.00071,1.00071,1.00071,40,40,40,0.6,42.4737,42.4737,84.8171,1.22057,0.125006,37.8383,1.20569,1.20569,51.842,47.88,6.00368,0.125006,1.00005,1.00051};
  double bp[29], work[LM_DIF_WORKSZ(29, 29)];
  struct rootfind_info rfi;
  rfi.aVOI = VOI;
  rfi.aCONSTANTS = CONSTANTS;
  rfi.aRATES = RATES;
  rfi.aSTATES = STATES;
  rfi.aALGEBRAIC = ALGEBRAIC;
  rfi.aPRET = pret;
  do_levmar(objfunc_1, p, bp, work, pret, 29, &rfi);
  ALGEBRAIC[232] = p[0];
  ALGEBRAIC[233] = p[1];
  ALGEBRAIC[234] = p[2];
  ALGEBRAIC[235] = p[3];
  ALGEBRAIC[236] = p[4];
  ALGEBRAIC[237] = p[5];
  ALGEBRAIC[238] = p[6];
  ALGEBRAIC[239] = p[7];
  ALGEBRAIC[240] = p[8];
  ALGEBRAIC[241] = p[9];
  ALGEBRAIC[242] = p[10];
  ALGEBRAIC[243] = p[11];
  ALGEBRAIC[244] = p[12];
  ALGEBRAIC[245] = p[13];
  ALGEBRAIC[246] = p[14];
  ALGEBRAIC[247] = p[15];
  ALGEBRAIC[248] = p[16];
  ALGEBRAIC[249] = p[17];
  ALGEBRAIC[250] = p[18];
  ALGEBRAIC[251] = p[19];
  ALGEBRAIC[252] = p[20];
  ALGEBRAIC[253] = p[21];
  ALGEBRAIC[254] = p[22];
  ALGEBRAIC[255] = p[23];
  ALGEBRAIC[256] = p[24];
  ALGEBRAIC[257] = p[25];
  ALGEBRAIC[258] = p[26];
  ALGEBRAIC[259] = p[27];
  ALGEBRAIC[260] = p[28];
}

void objfunc_2(double *p, double *hx, int m, int n, void *adata)
{
  struct rootfind_info* rfi = (struct rootfind_info*)adata;
#define VOI rfi->aVOI
#define CONSTANTS rfi->aCONSTANTS
#define RATES rfi->aRATES
#define STATES rfi->aSTATES
#define ALGEBRAIC rfi->aALGEBRAIC
#define pret rfi->aPRET
  ALGEBRAIC[271] = p[0];
  ALGEBRAIC[272] = p[1];
  ALGEBRAIC[273] = p[2];
  hx[0] = ALGEBRAIC[273] - (ALGEBRAIC[265] - ALGEBRAIC[271])/( ALGEBRAIC[157]*5.25000*ALGEBRAIC[234]);
  hx[1] = ALGEBRAIC[272] -  ALGEBRAIC[273]*57.1400;
  hx[2] = ALGEBRAIC[271] -  (ALGEBRAIC[272] - ALGEBRAIC[114])*12.8570*ALGEBRAIC[234];
#undef VOI
#undef CONSTANTS
#undef RATES
#undef STATES
#undef ALGEBRAIC
#undef pret
}
void rootfind_2(double VOI, double* CONSTANTS, double* RATES,
double* STATES, double* ALGEBRAIC, int* pret)
{
  static double p[3] = {163.508,37.6,0.694163};
  double bp[3], work[LM_DIF_WORKSZ(3, 3)];
  struct rootfind_info rfi;
  rfi.aVOI = VOI;
  rfi.aCONSTANTS = CONSTANTS;
  rfi.aRATES = RATES;
  rfi.aSTATES = STATES;
  rfi.aALGEBRAIC = ALGEBRAIC;
  rfi.aPRET = pret;
  do_levmar(objfunc_2, p, bp, work, pret, 3, &rfi);
  ALGEBRAIC[271] = p[0];
  ALGEBRAIC[272] = p[1];
  ALGEBRAIC[273] = p[2];
}