Generated Code

The following is c_ida code generated by the CellML API from this CellML file. (Back to language selection)

The raw code is available.

/*
   There are a total of 229 entries in the algebraic variable array.
   There are a total of 139 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[166] 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[125] is ANMAL in component angiotensin_control_of_aldosterone_secretion (dimensionless).
 * CONSTANTS[3] is ANMALD in component aldosterone_parameter_values (dimensionless).
 * ALGEBRAIC[194] is OSMAL in component osmotic_control_of_aldosterone_secretion (dimensionless).
 * ALGEBRAIC[198] 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[195] is AMRBSC in component aldosterone_secretion (dimensionless).
 * ALGEBRAIC[196] is AMRT in component aldosterone_secretion (dimensionless).
 * ALGEBRAIC[197] 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).
 * STATES[43] is MDFLW in component proximal_tubular_and_macula_densa_flow (L_per_minute).
 * ALGEBRAIC[9] is ANUVN in component angiotensin_effect_on_venous_constriction (dimensionless).
 * ALGEBRAIC[8] is ANU in component angiotensin_effect_on_circulation (dimensionless).
 * STATES[44] is ANGSCR in component instantaneous_angiotensin_formation (dimensionless).
 * STATES[45] 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).
 * STATES[46] is ANX in component time_delayed_angiotensin_formation (dimensionless).
 * STATES[47] is ANPR in component total_angiotensin_formation (dimensionless).
 * CONSTANTS[18] is REK in component angiotensin_parameter_values (dimensionless).
 * STATES[48] is ANPRT in component total_angiotensin_formation (dimensionless).
 * STATES[49] 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[7] 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[165] is CNA in component extracellular_Na_concentration (monovalent_mEq_per_litre).
 * ALGEBRAIC[130] is PA1 in component pressure_driving_autonomic_receptors (mmHg).
 * ALGEBRAIC[13] is ADHMK in component ADH_effect_on_kidney (dimensionless).
 * ALGEBRAIC[11] is ADHMV in component ADH_effect_on_nonrenal_vascular_resistance (dimensionless).
 * STATES[3] is ADHC in component ADH_in_blood (dimensionless).
 * ALGEBRAIC[200] 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[199] is ADHNA1 in component osmotic_control_of_ADH_secretion (dimensionless).
 * ALGEBRAIC[202] 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[201] is ADHPA in component pressure_control_of_ADH_secretion (mmHg).
 * ALGEBRAIC[204] is ADH in component total_ADH_secretion (dimensionless).
 * CONSTANTS[32] is ADHINF in component antidiuretic_hormone_parameter_values (dimensionless).
 * ALGEBRAIC[203] 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[10] 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[12] is ADHMK1 in component ADH_effect_on_kidney (dimensionless).
 * ALGEBRAIC[45] is PLA in component left_atrial_pressure (mmHg).
 * ALGEBRAIC[33] is PRA in component right_atrial_pressure (mmHg).
 * ALGEBRAIC[15] is ANPX in component ANP_effect_on_renal_afferent_arteriolar_resistance (dimensionless).
 * ALGEBRAIC[128] is ANP in component total_ANP_secreted (dimensionless).
 * ALGEBRAIC[126] is ANPL in component total_ANP_secreted (dimensionless).
 * ALGEBRAIC[127] is ANPR2 in component total_ANP_secreted (dimensionless).
 * ALGEBRAIC[129] 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[14] is ANPX1 in component ANP_effect_on_renal_afferent_arteriolar_resistance (dimensionless).
 * ALGEBRAIC[53] is PA in component arterial_pressure_and_pressure_gradient (mmHg).
 * ALGEBRAIC[189] is PO2ART in component arterial_PO2 (mmHg).
 * ALGEBRAIC[38] is PPA in component pulmonary_vasculature_pressure (mmHg).
 * ALGEBRAIC[19] is VVR in component autonomic_drive_on_target_organs_and_tissues (litre).
 * ALGEBRAIC[21] is AUH in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * ALGEBRAIC[22] is AUR in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * ALGEBRAIC[24] is AUM in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * ALGEBRAIC[23] is AOM in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * ALGEBRAIC[25] is AVE in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * ALGEBRAIC[18] 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[131] is AUC in component chemoreceptors_effect_of_PA (dimensionless).
 * CONSTANTS[44] is AUC1 in component autonomics_parameter_values (dimensionless).
 * ALGEBRAIC[206] is AUC3 in component chemoreceptors_effect_of_art_PO2 (dimensionless).
 * CONSTANTS[45] is O2CHMO in component autonomics_parameter_values (per_mmHg).
 * ALGEBRAIC[205] is AUC2 in component chemoreceptors_effect_of_art_PO2 (dimensionless).
 * ALGEBRAIC[16] 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[132] is AUB in component arterial_baroreceptor_reflex (dimensionless).
 * ALGEBRAIC[133] is A1B in component arterial_baroreceptor_reflex (dimensionless).
 * ALGEBRAIC[134] 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[135] is AUN in component CNS_ischemic_reflex (dimensionless).
 * CONSTANTS[50] is AUN1 in component autonomics_parameter_values (dimensionless).
 * ALGEBRAIC[136] 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[208] is AUTTL in component total_autonomic_stimulation (dimensionless).
 * ALGEBRAIC[207] 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[209] is DAU in component actual_autonomic_stimulation (per_minute).
 * STATES[6] is AU1 in component actual_autonomic_stimulation (dimensionless).
 * ALGEBRAIC[17] 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[20] is AUO in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
 * CONSTANTS[67] is MDMP in component autonomics_parameter_values (dimensionless).
 * ALGEBRAIC[61] is VEC in component extracellular_fluid_volume (litre).
 * ALGEBRAIC[187] is PPD in component concentration_of_protein_in_pulmonary_interstitium (gram_per_minute).
 * STATES[50] is RVS in component venous_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[185] is DFP in component pulmonary_interstitial_free_fluid_volume (L_per_minute).
 * ALGEBRAIC[77] is VPF in component pulmonary_interstitial_free_fluid_volume (litre).
 * STATES[51] is BFN in component blood_flow_through_NM_NR_tissues (L_per_minute).
 * ALGEBRAIC[156] is PVS in component venous_average_pressure (mmHg).
 * STATES[52] is PC in component capillary_pressure (mmHg).
 * ALGEBRAIC[26] is CPP in component plasma_protein_concentration (gram_per_L).
 * ALGEBRAIC[29] is PPC in component plasma_colloid_osmotic_pressure (mmHg).
 * STATES[7] is VP in component plasma_volume (litre).
 * ALGEBRAIC[210] is PGH in component hydrostatic_pressure_of_tissue_gel (mmHg).
 * ALGEBRAIC[143] is PTC in component total_osmotic_pressure_of_tissue_gel (mmHg).
 * STATES[53] is VTCPL in component plasma_leakage (L_per_minute).
 * STATES[54] 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[215] is VTL in component lymph_flow (L_per_minute).
 * CONSTANTS[69] is TRPL in component capillary_dynamics_parameter_values (L_per_minute).
 * STATES[55] is VPD in component plasma_volume (L_per_minute).
 * STATES[8] is PRP in component total_plasma_protein (gram).
 * ALGEBRAIC[28] 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[27] 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).
 * STATES[56] is PRCD in component plasma_leakage (mmHg).
 * ALGEBRAIC[139] is CPI in component interstitial_protein_concentration (gram_per_L).
 * STATES[57] is DPC in component protein_influx_into_interstitium (gram_per_minute).
 * ALGEBRAIC[223] is DPL in component lymph_protein_flow (gram_per_minute).
 * STATES[58] is DPP in component total_plasma_protein (gram_per_minute).
 * ALGEBRAIC[137] is VTS in component total_systemic_fluid_volume (litre).
 * ALGEBRAIC[138] 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).
 * STATES[59] is DPI in component total_interstitial_protein (gram_per_minute).
 * ALGEBRAIC[140] is PTCPR in component interstitial_colloid_osmotic_pressure (mmHg).
 * ALGEBRAIC[144] is PTT in component total_tissue_pressure (mmHg).
 * ALGEBRAIC[141] 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[142] 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[211] is PIF in component interstial_free_fluid_pressure (mmHg).
 * ALGEBRAIC[212] is PTS in component interstitial_solid_tissue_pressure (mmHg).
 * CONSTANTS[84] is PLDF in component capillary_dynamics_parameter_values (mmHg).
 * ALGEBRAIC[213] is PLD1 in component lymph_flow (mmHg).
 * ALGEBRAIC[214] is PLD in component lymph_flow (mmHg).
 * ALGEBRAIC[145] is VG in component interstitial_gel_volume (litre).
 * ALGEBRAIC[146] is VIF in component interstitial_free_fluid_volume (litre).
 * ALGEBRAIC[67] is AMM in component global_M_blood_flow_autoregulation_output (dimensionless).
 * ALGEBRAIC[74] is ARM in component global_NM_blood_flow_autoregulation_output (dimensionless).
 * ALGEBRAIC[111] is ATRRFB in component volume_effect_on_arteries (dimensionless).
 * ALGEBRAIC[112] is ATRVFB in component volume_effect_on_unstressed_venous_volume (litre).
 * ALGEBRAIC[62] 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[188] is OSA in component arterial_PO2 (dimensionless).
 * STATES[60] is RBF in component actual_renal_blood_flow (L_per_minute).
 * ALGEBRAIC[92] 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[164] is BFM in component blood_flow_through_M_tissues (L_per_minute).
 * STATES[61] is QAO in component systemic_blood_flow (L_per_minute).
 * ALGEBRAIC[225] is QRO in component right_ventricular_output (L_per_minute).
 * ALGEBRAIC[221] is QLO in component left_ventricular_output (L_per_minute).
 * ALGEBRAIC[42] is RPA in component pulmonary_arterial_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[148] is RPV in component pulmonary_venous_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[154] 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[30] is VBD in component total_blood_volume_change (litre).
 * ALGEBRAIC[159] is QVO in component rate_of_blood_flow_from_veins_to_right_atrium (L_per_minute).
 * ALGEBRAIC[31] is VRA in component right_atrial_blood_volume (litre).
 * ALGEBRAIC[228] is DRA in component right_atrial_blood_volume (L_per_minute).
 * ALGEBRAIC[32] is VRE in component right_atrial_pressure (litre).
 * ALGEBRAIC[34] is PRA1 in component autonomic_stimulation_effect_on_right_atrial_pressure (mmHg).
 * CONSTANTS[85] is HTAUML in component circulatory_dynamics_parameter_values (dimensionless).
 * ALGEBRAIC[217] is RVM in component pressure_effect_on_right_ventricular_pumping (dimensionless).
 * ALGEBRAIC[216] is PP2 in component pressure_effect_on_right_ventricular_pumping (mmHg).
 * ALGEBRAIC[47] is QLN in component left_ventricular_output (L_per_minute).
 * ALGEBRAIC[224] 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[35] is QRN in component right_ventricular_output (dimensionless).
 * ALGEBRAIC[151] is QPO in component rate_of_blood_flow_from_pulmonary_veins_to_left_atrium (L_per_minute).
 * ALGEBRAIC[36] is VPA in component pulmonary_vasculature_blood_volume (litre).
 * ALGEBRAIC[226] is DPA in component pulmonary_vasculature_blood_volume (L_per_minute).
 * ALGEBRAIC[37] is VPE in component pulmonary_vasculature_pressure (litre).
 * ALGEBRAIC[39] is PP1T in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
 * ALGEBRAIC[40] is PP1 in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
 * ALGEBRAIC[41] is CPA in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
 * ALGEBRAIC[147] is PL1 in component pulmonary_venous_resistance (mmHg).
 * ALGEBRAIC[149] is RPT in component total_pulmonary_vascular_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[150] is PGL in component pressure_gradient_through_the_lungs (mmHg).
 * ALGEBRAIC[43] is VLA in component left_atrial_blood_volume (litre).
 * ALGEBRAIC[227] is DLA in component left_atrial_blood_volume (L_per_minute).
 * ALGEBRAIC[44] is VLE in component left_atrial_pressure (litre).
 * ALGEBRAIC[46] is PLA1 in component autonomic_stimulation_effect_on_left_atrial_pressure (mmHg).
 * ALGEBRAIC[219] is LVM in component pumping_effectiveness_of_left_ventricle (dimensionless).
 * ALGEBRAIC[218] is PA2 in component pumping_effectiveness_of_left_ventricle (mmHg).
 * ALGEBRAIC[220] is QLOT in component left_ventricular_output (L_per_minute).
 * CONSTANTS[88] is HSL in component circulatory_dynamics_parameter_values (dimensionless).
 * ALGEBRAIC[152] is QLO1 in component left_ventricular_output (L_per_minute).
 * ALGEBRAIC[48] is VVS in component venous_blood_volume (litre).
 * STATES[62] is DVS in component venous_blood_volume (L_per_minute).
 * ALGEBRAIC[49] is VVA in component angiotensin_induced_venous_constriction (litre).
 * CONSTANTS[89] is ANY in component circulatory_dynamics_parameter_values (litre).
 * ALGEBRAIC[153] is VVE1 in component venous_excess_volume (litre).
 * CONSTANTS[90] is CV in component circulatory_dynamics_parameter_values (L_per_mmHg).
 * ALGEBRAIC[155] is PVS1 in component venous_average_pressure (mmHg).
 * ALGEBRAIC[50] is PR1 in component venous_outflow_pressure_into_heart (mmHg).
 * CONSTANTS[91] is PR1LL in component circulatory_dynamics_parameter_values (mmHg).
 * ALGEBRAIC[157] is RVG in component resistance_from_veins_to_right_atrium (mmHg_minute_per_L).
 * ALGEBRAIC[158] 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).
 * STATES[63] is CN3 in component venous_resistance (dimensionless).
 * STATES[64] is RV1 in component venous_resistance (mmHg_minute_per_L).
 * STATES[65] is NNRVR in component NM_NR_venous_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[51] is VAS in component arterial_blood_volume (litre).
 * STATES[66] is DAS in component arterial_blood_volume (L_per_minute).
 * ALGEBRAIC[54] is PAG in component arterial_pressure_and_pressure_gradient (mmHg).
 * ALGEBRAIC[52] is VAE in component arterial_pressure_and_pressure_gradient (litre).
 * ALGEBRAIC[55] is PAM in component pressure_effect_on_arterial_distention (dimensionless).
 * CONSTANTS[95] is PAEX in component circulatory_dynamics_parameter_values (dimensionless).
 * ALGEBRAIC[160] is R1 in component non_renal_systemic_arterial_resistance_multiplier (dimensionless).
 * ALGEBRAIC[161] 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[162] is PGS in component pressure_gradient_from_arteries_to_veins (mmHg).
 * ALGEBRAIC[163] 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).
 * STATES[67] is RSN in component total_NM_NR_systemic_resistance (mmHg_minute_per_L).
 * ALGEBRAIC[56] 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).
 * STATES[68] is SYSFLO in component systemic_blood_flow (L_per_minute).
 * STATES[69] is RTP in component total_peripheral_resistance (mmHg_minute_per_L).
 * STATES[21] is TVD in component rate_of_fluid_intake (L_per_minute).
 * STATES[70] is NOD in component actual_Na_excretion_rate (monovalent_mEq_per_minute).
 * ALGEBRAIC[101] is STH in component effect_of_salt_appetite_stimulation_on_thirst (dimensionless).
 * STATES[71] is KOD in component actual_K_excretion_rate (monovalent_mEq_per_minute).
 * STATES[72] 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).
 * STATES[73] is NED in component extracellular_Na_concentration (monovalent_mEq_per_minute).
 * STATES[23] is NAE in component extracellular_Na_concentration (monovalent_mEq).
 * ALGEBRAIC[57] is AMK1 in component aldosterone_effect_on_cellular_K_distribution (dimensionless).
 * CONSTANTS[102] is ALCLK in component electrolytes_parameter_values (dimensionless).
 * ALGEBRAIC[58] 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).
 * STATES[74] is KTOTD in component extracellular_K_concentration (monovalent_mEq_per_minute).
 * STATES[25] is VIC in component intracellular_fluid_volume (litre).
 * ALGEBRAIC[60] is CKI in component intracellular_K_concentration (monovalent_mEq_per_litre).
 * ALGEBRAIC[59] is KI in component intracellular_K_concentration (monovalent_mEq).
 * ALGEBRAIC[168] 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[167] is CCD in component intracellular_fluid_volume (monovalent_mEq_per_litre).
 * ALGEBRAIC[76] 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[169] is DHM in component heart_deterioration (per_minute).
 * ALGEBRAIC[63] 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[64] is HDHR in component effect_of_heart_deterioration_on_HR (dimensionless).
 * ALGEBRAIC[65] is HR in component heart_rate (beats_per_minute).
 * ALGEBRAIC[222] is SVO in component stroke_volume_output (litre).
 * ALGEBRAIC[69] is PMO in component pressure_of_O2_in_M_tissue_cells (mmHg).
 * ALGEBRAIC[170] is PDO in component M_autoregulatory_driving_force (mmHg).
 * ALGEBRAIC[171] is POE in component M_ST_sensitivity_control (mmHg).
 * CONSTANTS[110] is POM in component M_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * ALGEBRAIC[66] 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[172] 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[89] is HM in component hematocrit_fraction (dimensionless).
 * STATES[75] is RMO in component delivery_of_O2_to_M_tissues (mL_per_minute).
 * ALGEBRAIC[173] is O2ARTM in component M_O2_blood_supply (mL_per_minute).
 * STATES[76] 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).
 * STATES[77] is OVS in component M_venous_O2_content (dimensionless).
 * ALGEBRAIC[175] is MMO in component metabolic_O2_consumption_by_M_tissue (mL_per_minute).
 * ALGEBRAIC[174] 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[68] is QOM in component volume_of_O2_in_M_tissue (mL).
 * STATES[78] 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[176] is POD in component NM_autoregulatory_driving_force (mmHg).
 * CONSTANTS[120] is POR in component NM_autoregulatory_local_blood_flow_parameter_values (mmHg).
 * ALGEBRAIC[177] is POB in component NM_ST_sensitivity_control (mmHg).
 * CONSTANTS[121] is POK in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * ALGEBRAIC[70] 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[178] is POA in component NM_I_sensitivity_control (mmHg).
 * CONSTANTS[123] is PON in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * ALGEBRAIC[71] 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[179] is POC in component NM_LT_sensitivity_control (mmHg).
 * CONSTANTS[125] is POZ in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * ALGEBRAIC[72] 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[73] is ARM1 in component total_NM_autoregulation (dimensionless).
 * CONSTANTS[127] is AUTOSN in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
 * STATES[79] is DOB in component delivery_of_O2_to_NM_tissues (mL_per_minute).
 * STATES[80] is O2ARTN in component NM_O2_blood_supply (mL_per_minute).
 * STATES[81] is POV in component NM_venous_O2_content (mmHg).
 * STATES[82] is OSV in component NM_venous_O2_content (dimensionless).
 * ALGEBRAIC[181] 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[180] is P1O in component O2_consumption_by_NM_tissue (mmHg).
 * ALGEBRAIC[75] is QO2 in component volume_of_O2_in_NM_tissue (mL).
 * STATES[83] is DO2N in component volume_of_O2_in_NM_tissue (mL_per_minute).
 * STATES[84] 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[182] is PCP in component pulmonary_capillary_pressure (mmHg).
 * ALGEBRAIC[81] is POS in component colloid_osmotic_pressure_of_pulmonary_interstitium (mmHg).
 * ALGEBRAIC[78] is PPI in component pulmonary_interstitial_fluid_pressure (mmHg).
 * ALGEBRAIC[183] 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[83] is PLF in component lung_lymphatic_protein_flow (L_per_minute).
 * ALGEBRAIC[184] 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[84] is PPO in component lung_lymphatic_protein_flow (gram_per_minute).
 * ALGEBRAIC[82] is PPN in component protein_leakage_into_pulmonary_interstitium (gram_per_minute).
 * ALGEBRAIC[80] is CPN in component concentration_of_protein_in_pulmonary_interstitium (gram_per_L).
 * ALGEBRAIC[186] 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[79] is PPR in component concentration_of_protein_in_pulmonary_interstitium (gram).
 * STATES[85] is O2UTIL in component total_O2_utilization (mL_per_minute).
 * ALGEBRAIC[86] is O2VAD2 in component progressive_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
 * ALGEBRAIC[192] is O2VTS2 in component acute_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
 * STATES[86] is ALVENT in component alveolar_ventilation (L_per_minute).
 * CONSTANTS[130] is VNTSTM in component pulmonary_O2_uptake_parameter_values (dimensionless).
 * STATES[87] is PO2ALV in component alveolar_PO2 (mmHg).
 * CONSTANTS[131] is PO2AMB in component pulmonary_O2_uptake_parameter_values (mmHg).
 * STATES[88] 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[85] is RSPDFC in component respiratory_O2_diffusion_into_capillaries (mL_per_minute_per_mmHg).
 * STATES[89] is DOVA in component O2_volume_of_arterial_blood (mL_per_L_per_minute).
 * ALGEBRAIC[191] is O2VTST in component acute_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
 * ALGEBRAIC[190] is O2VTST1 in component acute_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
 * ALGEBRAIC[193] 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[88] is HM1 in component hematocrit_fraction (dimensionless).
 * ALGEBRAIC[87] is VB in component hematocrit_fraction (litre).
 * ALGEBRAIC[90] 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[91] is VIB in component blood_viscosity (dimensionless).
 * ALGEBRAIC[95] 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[93] is HM3 in component oxygen_stimulation (mmHg).
 * CONSTANTS[209] is HM4 in component oxygen_stimulation (mmHg).
 * ALGEBRAIC[94] is HM5 in component oxygen_stimulation (mmHg).
 * ALGEBRAIC[96] 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[97] 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[98] 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[99] is ANMSML in component effect_of_salt_appetite_stimulation_on_thirst (dimensionless).
 * ALGEBRAIC[100] is STH1 in component effect_of_salt_appetite_stimulation_on_thirst (dimensionless).
 * ALGEBRAIC[102] 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[103] 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[105] is AHTH in component rate_of_fluid_intake (dimensionless).
 * ALGEBRAIC[104] is AHTH1 in component rate_of_fluid_intake (dimensionless).
 * ALGEBRAIC[107] is TVZ in component rate_of_fluid_intake (L_per_minute).
 * ALGEBRAIC[106] is TVZ1 in component rate_of_fluid_intake (L_per_minute).
 * ALGEBRAIC[109] 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[108] 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[110] 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[113] 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).
 * STATES[90] 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).
 * STATES[91] is RNAUG1 in component renal_autoregulatory_feedback_factor (dimensionless).
 * STATES[92] is RNAUG1T in component renal_autoregulatory_feedback_factor (dimensionless).
 * STATES[40] is RNAUG3 in component renal_autoregulatory_feedback_factor (dimensionless).
 * ALGEBRAIC[115] is AUMK in component autonomic_effect_on_AAR (dimensionless).
 * CONSTANTS[167] is ARF in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[114] is AUMKT in component autonomic_effect_on_AAR (dimensionless).
 * ALGEBRAIC[117] 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[116] is ANMAR1 in component angiotensin_effect_on_AAR (dimensionless).
 * STATES[93] is AAR1 in component AAR_calculation (mmHg_minute_per_L).
 * CONSTANTS[170] is AARK in component kidney_parameter_values (mmHg_minute_per_L).
 * STATES[94] 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).
 * STATES[95] is AART in component atrial_natriuretic_peptide_effect_on_AAR (mmHg_minute_per_L).
 * ALGEBRAIC[118] is AUMK2 in component autonomic_effect_on_EAR (dimensionless).
 * CONSTANTS[173] is AUMK1 in component kidney_parameter_values (dimensionless).
 * ALGEBRAIC[119] is ANMER in component angiotensin_effect_on_EAR (dimensionless).
 * CONSTANTS[174] is ANMEM in component kidney_parameter_values (dimensionless).
 * STATES[96] is RNAUG4 in component effect_of_renal_autoregulatory_feedback_on_EAR (dimensionless).
 * CONSTANTS[175] is EFAFR in component kidney_parameter_values (dimensionless).
 * STATES[97] 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).
 * STATES[98] is EAR1 in component EAR_calculation (mmHg_minute_per_L).
 * STATES[99] is RR in component total_renal_resistance (mmHg_minute_per_L).
 * STATES[100] is RFN in component normal_renal_blood_flow (L_per_minute).
 * CONSTANTS[178] is REK in component kidney_parameter_values (dimensionless).
 * STATES[101] is GFN in component glomerular_filtration_rate (L_per_minute).
 * STATES[102] 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).
 * STATES[103] is EFAFPR in component glomerular_colloid_osmotic_pressure (dimensionless).
 * STATES[104] is EFAFPR1 in component glomerular_colloid_osmotic_pressure (dimensionless).
 * STATES[105] is GLP in component glomerular_pressure (mmHg).
 * STATES[106] is APD in component glomerular_pressure (mmHg).
 * STATES[107] 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).
 * STATES[108] is PFL in component glomerular_filtration_rate (mmHg).
 * STATES[109] is GFN1 in component glomerular_filtration_rate (L_per_minute).
 * CONSTANTS[184] is MDFL1 in component kidney_parameter_values (dimensionless).
 * STATES[110] is PTFL in component proximal_tubular_and_macula_densa_flow (L_per_minute).
 * STATES[111] is MDFLWT in component proximal_tubular_and_macula_densa_flow (L_per_minute).
 * STATES[112] 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).
 * STATES[113] is RTSPPC1 in component renal_tissue_osmotic_pressure (mmHg).
 * STATES[114] 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[120] is PLURC in component plasma_urea_concentration (mOsm_per_litre).
 * STATES[115] 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).
 * STATES[116] 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).
 * STATES[117] is RABSPR in component peritubular_capillary_reabsorption_factor (mmHg).
 * STATES[118] is RFAB1 in component peritubular_capillary_reabsorption_factor (dimensionless).
 * STATES[119] is RFAB in component peritubular_capillary_reabsorption_factor (dimensionless).
 * STATES[120] is RFABD1 in component peritubular_capillary_reabsorption_factor (dimensionless).
 * STATES[121] is DTNAI in component distal_tubular_Na_delivery (monovalent_mEq_per_minute).
 * STATES[122] 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).
 * STATES[123] is DTNARA1 in component Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
 * ALGEBRAIC[122] 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[121] is DTNANG1 in component angiotensin_induced_Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
 * STATES[124] is DTKI in component distal_tubular_K_delivery (monovalent_mEq_per_minute).
 * STATES[125] 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).
 * STATES[126] 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).
 * STATES[127] is MDFLK1 in component effect_of_fluid_flow_on_distal_K_reabsorption (monovalent_mEq_per_minute).
 * STATES[128] is KODN in component normal_K_excretion (monovalent_mEq_per_minute).
 * STATES[129] 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).
 * STATES[130] 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[123] is ANMKE1 in component K_secretion_from_distal_tubules (dimensionless).
 * ALGEBRAIC[124] is ANMKE in component K_secretion_from_distal_tubules (dimensionless).
 * STATES[131] is NODN in component normal_Na_excretion (monovalent_mEq_per_minute).
 * STATES[132] is NODN1 in component normal_Na_excretion (monovalent_mEq_per_minute).
 * STATES[133] is KODN1 in component normal_K_excretion (monovalent_mEq_per_minute).
 * STATES[134] is DTURI in component normal_urea_excretion (mOsm_per_minute).
 * STATES[135] is OSMOPN1 in component normal_osmolar_and_water_excretion (mOsm_per_minute).
 * STATES[136] is OSMOPN in component normal_osmolar_and_water_excretion (mOsm_per_minute).
 * STATES[137] is OSMOP1T in component normal_urine_volume (mOsm_per_minute).
 * STATES[138] 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).
 * There are a total of 141 condition variables.
 */
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;
STATES[43] = 1.00051;
STATES[44] = 0.1001;
STATES[45] = 0.1001;
STATES[46] = 0.1001;
STATES[47] = 0.1001;
STATES[48] = 0.1001;
STATES[49] = 0.1001;
STATES[50] = 2.77751;
STATES[51] = 2.79521;
STATES[52] = 16.9144;
STATES[53] = 0.1001;
STATES[54] = 0.1001;
STATES[55] = 0.1001;
STATES[56] = 0.1001;
STATES[57] = 0.1001;
STATES[58] = 0.1001;
STATES[59] = 0.1001;
STATES[60] = 0.1001;
STATES[61] = 0.1001;
STATES[62] = 0.1001;
STATES[63] = 0.36019;
STATES[64] = 2.77632;
STATES[65] = 4.97174;
STATES[66] = 0.1001;
STATES[67] = 35.7072;
STATES[68] = 0.1001;
STATES[69] = 0.1001;
STATES[70] = 0.1001;
STATES[71] = 0.1001;
STATES[72] = 0.1001;
STATES[73] = 0.1001;
STATES[74] = 0.1001;
STATES[75] = 56.8057;
STATES[76] = 39.9793;
STATES[77] = 0.699673;
STATES[78] = 0.1001;
STATES[79] = 163.508;
STATES[80] = 0.1001;
STATES[81] = 37.6;
STATES[82] = 0.694163;
STATES[83] = 0.1001;
STATES[84] = 0.1001;
STATES[85] = 0.1001;
STATES[86] = 0.1001;
STATES[87] = 0.1001;
STATES[88] = 0.1001;
STATES[89] = 0.1001;
STATES[90] = 1.00071;
STATES[91] = 1.00071;
STATES[92] = 1.00071;
STATES[93] = 40;
STATES[94] = 40;
STATES[95] = 40;
STATES[96] = 0.6;
STATES[97] = 42.4737;
STATES[98] = 42.4737;
STATES[99] = 84.8171;
STATES[100] = 1.22057;
STATES[101] = 0.125006;
STATES[102] = 37.8383;
STATES[103] = 1.20569;
STATES[104] = 1.20569;
STATES[105] = 51.842;
STATES[106] = 47.88;
STATES[107] = 0.1001;
STATES[108] = 6.00368;
STATES[109] = 0.125006;
STATES[110] = 1.00005;
STATES[111] = 1.00051;
STATES[112] = 0.1001;
STATES[113] = 0.1001;
STATES[114] = 0.1001;
STATES[115] = 0.1001;
STATES[116] = 0.1001;
STATES[117] = 0.1001;
STATES[118] = 0.1001;
STATES[119] = 0.1001;
STATES[120] = 0.1001;
STATES[121] = 0.1001;
STATES[122] = 0.1001;
STATES[123] = 0.1001;
STATES[124] = 0.1001;
STATES[125] = 0.1001;
STATES[126] = 0.1001;
STATES[127] = 0.1001;
STATES[128] = 0.1001;
STATES[129] = 0.1001;
STATES[130] = 0.1001;
STATES[131] = 0.1001;
STATES[132] = 0.1001;
STATES[133] = 0.1001;
STATES[134] = 0.1001;
STATES[135] = 0.1001;
STATES[136] = 0.1001;
STATES[137] = 0.1001;
STATES[138] = 0.1001;
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[2] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
RATES[6] = 0.1001;
RATES[7] = 0.1001;
RATES[8] = 0.1001;
RATES[9] = 0.1001;
RATES[10] = 0.1001;
RATES[20] = 0.1001;
RATES[19] = 0.1001;
RATES[18] = 0.1001;
RATES[16] = 0.1001;
RATES[17] = 0.1001;
RATES[23] = 0.1001;
RATES[24] = 0.1001;
RATES[25] = 0.1001;
RATES[22] = 0.1001;
RATES[11] = 0.1001;
RATES[12] = 0.1001;
RATES[26] = 0.1001;
RATES[27] = 0.1001;
RATES[28] = 0.1001;
RATES[30] = 0.1001;
RATES[31] = 0.1001;
RATES[32] = 0.1001;
RATES[33] = 0.1001;
RATES[34] = 0.1001;
RATES[35] = 0.1001;
RATES[36] = 0.1001;
RATES[29] = 0.1001;
RATES[37] = 0.1001;
RATES[13] = 0.1001;
RATES[15] = 0.1001;
RATES[14] = 0.1001;
RATES[21] = 0.1001;
RATES[38] = 0.1001;
RATES[39] = 0.1001;
RATES[40] = 0.1001;
RATES[41] = 0.1001;
RATES[42] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - (ALGEBRAIC[198] - STATES[0])/CONSTANTS[7];
resid[1] = STATES[45] - STATES[43];
resid[2] = STATES[44] - (CONDVAR[5]>0.00000 ? 1.00000/(1.00000+ (STATES[45] - 1.00000)*72.0000) : 10.0000 - 9.00000/(1.00000+ (1.00000 - STATES[45])*8.00000));
resid[3] = STATES[46] -  (STATES[44] - 1.00000)*CONSTANTS[16];
resid[4] = RATES[1] - (STATES[46] - STATES[1])/CONSTANTS[17];
resid[5] = STATES[48] -  (STATES[44]+STATES[1])*CONSTANTS[18];
resid[6] = STATES[47] - (CONDVAR[6]<0.00000 ? 1.00000e-05 : STATES[48]);
resid[7] = STATES[49] - (CONDVAR[7]>0.00000 ? CONSTANTS[19] : STATES[47]+CONSTANTS[20]);
resid[8] = RATES[2] - (STATES[49] - STATES[2])/CONSTANTS[21];
resid[9] = RATES[3] - (ALGEBRAIC[204] - STATES[3])/CONSTANTS[33];
resid[10] = RATES[4] - (ALGEBRAIC[129] - STATES[4])/CONSTANTS[40];
resid[11] = RATES[5] - (ALGEBRAIC[134] - STATES[5])/CONSTANTS[48];
resid[12] = RATES[6] - ALGEBRAIC[209];
resid[13] = STATES[52] -  STATES[50]*1.70000*STATES[51]+ALGEBRAIC[156];
resid[14] = STATES[54] -  (((STATES[52] - ALGEBRAIC[29]) - ALGEBRAIC[210])+ALGEBRAIC[143])*CONSTANTS[68]+STATES[53];
resid[15] = STATES[55] - ((ALGEBRAIC[215] - STATES[54]) - ALGEBRAIC[185])+CONSTANTS[69];
resid[16] = RATES[7] - STATES[55];
resid[17] = STATES[56] - (CONDVAR[32]<0.00000 ? 0.00000 : STATES[52] - CONSTANTS[74]);
resid[18] = STATES[53] - pow( STATES[56]*CONSTANTS[75], CONSTANTS[76]);
resid[19] = STATES[57] -  STATES[53]*ALGEBRAIC[26]+ (ALGEBRAIC[26] - ALGEBRAIC[139])*0.00104000;
resid[20] = STATES[58] - (((ALGEBRAIC[28]+ALGEBRAIC[223]) - STATES[57]) - ALGEBRAIC[187])+ CONSTANTS[69]*72.0000;
resid[21] = RATES[8] - STATES[58];
resid[22] = RATES[9] -  ( (ALGEBRAIC[137] - 12.0000)*CONSTANTS[77] - STATES[9])*CONSTANTS[78];
resid[23] = STATES[59] - STATES[57] - ALGEBRAIC[223];
resid[24] = RATES[10] - STATES[59];
resid[25] = RATES[20] - ALGEBRAIC[228];
resid[26] = RATES[19] - ALGEBRAIC[226];
resid[27] = RATES[18] - ALGEBRAIC[227];
resid[28] = STATES[62] - STATES[61] - ALGEBRAIC[159];
resid[29] = RATES[16] - STATES[62];
resid[30] = STATES[63] -  ( (STATES[52] - 17.0000)*CONSTANTS[92]+17.0000)*CONSTANTS[93];
resid[31] = STATES[64] - CONSTANTS[94]/STATES[63];
resid[32] = STATES[50] -  ALGEBRAIC[25]*STATES[64]*ALGEBRAIC[92]*ALGEBRAIC[9];
resid[33] = STATES[65] -  STATES[50]*1.79000;
resid[34] = STATES[66] - ALGEBRAIC[221] - STATES[61];
resid[35] = RATES[17] - STATES[66];
resid[36] = STATES[67] - ALGEBRAIC[161]+STATES[65];
resid[37] = STATES[51] - ALGEBRAIC[162]/STATES[67];
resid[38] = STATES[68] - ALGEBRAIC[164]+STATES[51]+STATES[60];
resid[39] = STATES[61] - STATES[68]+ALGEBRAIC[56];
resid[40] = STATES[69] - ALGEBRAIC[54]/STATES[61];
resid[41] = STATES[73] - ( CONSTANTS[100]*ALGEBRAIC[101] - STATES[70])+ CONSTANTS[101]*142.000;
resid[42] = RATES[23] - STATES[73];
resid[43] = STATES[74] - CONSTANTS[103] - STATES[71];
resid[44] = RATES[24] - STATES[74];
resid[45] = RATES[25] - ALGEBRAIC[168];
resid[46] = RATES[22] - STATES[21] - STATES[72];
resid[47] = RATES[11] - (pow(( ALGEBRAIC[53]*STATES[61])/( 500.000*CONSTANTS[105]), CONSTANTS[106]) - STATES[11])/57600.0;
resid[48] = RATES[12] - (pow(( ALGEBRAIC[38]*STATES[61])/( 75.0000*CONSTANTS[107]), CONSTANTS[106]) - STATES[12])/57600.0;
resid[49] = RATES[26] - ALGEBRAIC[169];
resid[50] = RATES[27] - ( ALGEBRAIC[171]*1.00000 - STATES[27])/CONSTANTS[111];
resid[51] = RATES[28] - ( ALGEBRAIC[172]*1.00000 - STATES[28])/CONSTANTS[114];
resid[52] = STATES[77] - (ALGEBRAIC[173] - STATES[75])/( ALGEBRAIC[89]*5.25000*ALGEBRAIC[164]);
resid[53] = STATES[76] -  57.1400*STATES[77]*pow(CONSTANTS[115], CONSTANTS[116]);
resid[54] = STATES[75] -  (STATES[76] - ALGEBRAIC[69])*CONSTANTS[118]*ALGEBRAIC[164];
resid[55] = STATES[78] - STATES[75] - ALGEBRAIC[175];
resid[56] = RATES[30] - STATES[78];
resid[57] = RATES[31] - ( ALGEBRAIC[177]*1.00000 - STATES[31])/CONSTANTS[122];
resid[58] = RATES[32] - ( ALGEBRAIC[178]*1.00000 - STATES[32])/CONSTANTS[124];
resid[59] = RATES[33] - ( ALGEBRAIC[179]*1.00000 - STATES[33])/CONSTANTS[126];
resid[60] = STATES[80] -  STATES[29]*STATES[51];
resid[61] = STATES[82] - (STATES[80] - STATES[79])/( ALGEBRAIC[89]*5.25000*STATES[51]);
resid[62] = STATES[81] -  STATES[82]*57.1400;
resid[63] = STATES[79] -  (STATES[81] - ALGEBRAIC[76])*12.8570*STATES[51];
resid[64] = STATES[84] - STATES[79] - ALGEBRAIC[181];
resid[65] = STATES[83] - (CONDVAR[99]<0.00000&&CONDVAR[100]<0.00000 ?  STATES[84]*0.100000 : STATES[84]);
resid[66] = RATES[34] - STATES[83];
resid[67] = RATES[35] - ALGEBRAIC[185];
resid[68] = RATES[36] - ALGEBRAIC[187];
resid[69] = STATES[85] - STATES[79]+STATES[75];
resid[70] = STATES[86] -  STATES[85]*CONSTANTS[130]*0.0266670*ALGEBRAIC[192]*ALGEBRAIC[86];
resid[71] = STATES[87] - CONSTANTS[131] - (STATES[85]/STATES[86])/0.761000;
resid[72] = STATES[88] -  (STATES[87] - ALGEBRAIC[189])*ALGEBRAIC[85];
resid[73] = STATES[89] - (STATES[88] - STATES[85])/( ALGEBRAIC[225]*1.00000);
resid[74] = RATES[29] - STATES[89];
resid[75] = RATES[37] - ALGEBRAIC[193];
resid[76] = RATES[13] - ALGEBRAIC[98];
resid[77] = RATES[15] - ( (ALGEBRAIC[154] - 0.740000)*CONSTANTS[142] - STATES[15])/CONSTANTS[143];
resid[78] = RATES[14] - ( (ALGEBRAIC[154] - 0.740000)*CONSTANTS[144] - STATES[14])/CONSTANTS[145];
resid[79] = RATES[21] - ((ALGEBRAIC[107]+CONSTANTS[151]) - STATES[21])/CONSTANTS[152];
resid[80] = RATES[38] - (ALGEBRAIC[109] - STATES[38])/CONSTANTS[155];
resid[81] = RATES[39] - ((100.000+ (ALGEBRAIC[53] - 100.000)*CONSTANTS[161]) - STATES[39])/CONSTANTS[162];
resid[82] = STATES[92] -  (STATES[43] - 1.00000)*CONSTANTS[163]+1.00000;
resid[83] = STATES[91] - (CONDVAR[121]<0.00000 ? CONSTANTS[164] : CONDVAR[122]>0.00000 ? CONSTANTS[165] : STATES[92]);
resid[84] = STATES[90] - STATES[91] - STATES[40];
resid[85] = RATES[40] -  (STATES[90] - 1.00000)*CONSTANTS[166];
resid[86] = STATES[93] -  CONSTANTS[170]*CONSTANTS[2]*ALGEBRAIC[115]*STATES[90]*ALGEBRAIC[117]*40.0000*CONSTANTS[0];
resid[87] = STATES[95] - (STATES[93] -  ALGEBRAIC[15]*CONSTANTS[171])+CONSTANTS[171];
resid[88] = STATES[94] - (CONDVAR[125]<0.00000 ? CONSTANTS[172] : STATES[95]);
resid[89] = STATES[96] -  (STATES[90] - 1.00000)*CONSTANTS[175]+1.00000;
resid[90] = STATES[98] -  43.3330*CONSTANTS[176]*ALGEBRAIC[119]*STATES[96]*CONSTANTS[0]*ALGEBRAIC[118];
resid[91] = STATES[97] - (CONDVAR[126]<0.00000 ? CONSTANTS[177] : STATES[98]);
resid[92] = STATES[99] - STATES[94]+STATES[97];
resid[93] = STATES[100] - ALGEBRAIC[113]/STATES[99];
resid[94] = STATES[60] -  CONSTANTS[178]*STATES[100];
resid[95] = STATES[104] - ( STATES[100]*(1.00000 - ALGEBRAIC[88]))/( STATES[100]*(1.00000 - ALGEBRAIC[88]) - STATES[101]);
resid[96] = STATES[103] - (CONDVAR[127]<0.00000 ? 1.00000 : STATES[104]);
resid[97] = STATES[102] - (CONDVAR[128]>0.00000 ?  pow(STATES[103], 1.35000)*ALGEBRAIC[29]*0.980000 : ALGEBRAIC[29]+4.00000);
resid[98] = STATES[106] -  STATES[94]*STATES[100];
resid[99] = STATES[105] - ALGEBRAIC[113] - STATES[106];
resid[100] = STATES[108] - (STATES[105] - STATES[102]) - CONSTANTS[181];
resid[101] = STATES[109] -  STATES[108]*CONSTANTS[182];
resid[102] = STATES[101] - (CONDVAR[129]<0.00000 ? CONSTANTS[183] : STATES[109]);
resid[103] = STATES[107] -  STATES[101]*CONSTANTS[178];
resid[104] = STATES[110] -  STATES[101]*8.00000;
resid[105] = STATES[111] -  (STATES[110] - 1.00000)*CONSTANTS[184]+1.00000;
resid[106] = STATES[43] - (CONDVAR[130]<0.00000 ? 0.00000 : STATES[111]);
resid[107] = STATES[113] -  STATES[102]*CONSTANTS[185] - CONSTANTS[186];
resid[108] = STATES[112] - (CONDVAR[131]<0.00000 ? 1.00000 : STATES[113]);
resid[109] = RATES[41] - CONSTANTS[187] - STATES[114];
resid[110] = STATES[115] -  ( (STATES[100] - 1.20000)*CONSTANTS[188]+1.20000)*CONSTANTS[189];
resid[111] = STATES[117] - ((STATES[102]+CONSTANTS[190]) - STATES[115]) - STATES[112];
resid[112] = STATES[118] -  STATES[117]*CONSTANTS[191];
resid[113] = STATES[119] - STATES[118];
resid[114] = STATES[120] -  (STATES[119] - 1.00000)*CONSTANTS[193]+1.00000;
resid[115] = STATES[116] - (CONDVAR[132]<0.00000 ? 0.000100000 : STATES[120]);
resid[116] = STATES[121] -  STATES[43]*ALGEBRAIC[165]*0.00616190;
resid[117] = STATES[123] -  (( ALGEBRAIC[5]*STATES[116]*CONSTANTS[194])/CONSTANTS[195])*( (ALGEBRAIC[13] - 1.00000)*CONSTANTS[196]+1.00000);
resid[118] = STATES[122] - (CONDVAR[133]<0.00000 ? CONSTANTS[197] : STATES[123]);
resid[119] = STATES[124] - ( STATES[121]*ALGEBRAIC[166])/ALGEBRAIC[165];
resid[120] = STATES[125] -  (STATES[116] - 1.00000)*CONSTANTS[199];
resid[121] = STATES[127] -  (STATES[43] - 1.00000)*CONSTANTS[200]+1.00000;
resid[122] = STATES[126] - (CONDVAR[135]<0.00000 ? 0.100000 : STATES[127]);
resid[123] = RATES[42] -  ( (STATES[128]/STATES[129])*0.000451800 - STATES[42])*1.00000;
resid[124] = STATES[130] - ( pow(ALGEBRAIC[166]/4.40000, CONSTANTS[203])*ALGEBRAIC[3]*0.0800000*STATES[126])/ALGEBRAIC[124];
resid[125] = STATES[132] - (STATES[121] - STATES[122]) - ALGEBRAIC[122];
resid[126] = STATES[131] - (CONDVAR[137]<0.00000 ? 1.00000e-08 : STATES[132]);
resid[127] = STATES[133] - ((STATES[124]+STATES[130]) - STATES[42]) - STATES[125];
resid[128] = STATES[128] - (CONDVAR[138]<0.00000 ? 0.00000 : STATES[133]);
resid[129] = STATES[134] -  pow(STATES[101], 2.00000)*ALGEBRAIC[120]*3.84000;
resid[130] = STATES[135] - STATES[134]+ 2.00000*(STATES[131]+STATES[128]);
resid[131] = STATES[136] - (CONDVAR[139]>0.00000 ? 0.600000 : STATES[135]);
resid[132] = STATES[137] - STATES[135] - 0.600000;
resid[133] = STATES[138] - (CONDVAR[140]<0.00000 ? 0.00000 : STATES[137]);
resid[134] = STATES[129] - STATES[136]/( 600.000*ALGEBRAIC[13])+STATES[138]/360.000;
resid[135] = STATES[70] -  STATES[131]*CONSTANTS[178];
resid[136] = STATES[71] -  STATES[128]*CONSTANTS[178];
resid[137] = STATES[114] -  STATES[134]*CONSTANTS[178];
resid[138] = STATES[72] -  STATES[129]*CONSTANTS[178];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[22] =  ALGEBRAIC[20]*CONSTANTS[62]+1.00000;
ALGEBRAIC[63] =  72.0000*ALGEBRAIC[22];
ALGEBRAIC[64] =  (ALGEBRAIC[62] - 1.00000)*0.500000+1.00000;
ALGEBRAIC[65] =  (ALGEBRAIC[63]+CONSTANTS[207])*ALGEBRAIC[64];
ALGEBRAIC[145] = (CONDVAR[35]<=0.00000 ? 0.00000 : CONDVAR[36]>0.00000&&CONDVAR[37]<=0.00000 ? 0.00000+( (11.4000 - 0.00000)*(ALGEBRAIC[137] - 0.00000))/(12.0000 - 0.00000) : CONDVAR[38]>0.00000&&CONDVAR[39]<=0.00000 ? 11.4000+( (14.0000 - 11.4000)*(ALGEBRAIC[137] - 12.0000))/(15.0000 - 12.0000) : CONDVAR[40]>0.00000&&CONDVAR[41]<=0.00000 ? 14.0000+( (16.0000 - 14.0000)*(ALGEBRAIC[137] - 15.0000))/(18.0000 - 15.0000) : CONDVAR[42]>0.00000&&CONDVAR[43]<=0.00000 ? 16.0000+( (17.3000 - 16.0000)*(ALGEBRAIC[137] - 18.0000))/(21.0000 - 18.0000) : CONDVAR[44]>0.00000&&CONDVAR[45]<=0.00000 ? 17.3000+( (18.0000 - 17.3000)*(ALGEBRAIC[137] - 21.0000))/(24.0000 - 21.0000) : 18.0000);
ALGEBRAIC[146] = ALGEBRAIC[137] - ALGEBRAIC[145];
ALGEBRAIC[212] = ALGEBRAIC[144] - ALGEBRAIC[211];
ALGEBRAIC[222] = ALGEBRAIC[221]/ALGEBRAIC[65];
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
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] = (CONDVAR[2]<0.00000 ? 0.200000 : ALGEBRAIC[2]);
ALGEBRAIC[4] =  (ALGEBRAIC[1] - 1.00000)*CONSTANTS[13]+1.00000;
ALGEBRAIC[5] = (CONDVAR[3]<0.00000 ? CONSTANTS[15] : CONDVAR[4]>0.00000 ? CONSTANTS[14] : ALGEBRAIC[4]);
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[25]+1.00000;
ALGEBRAIC[8] = (CONDVAR[8]<0.00000 ? CONSTANTS[26] : ALGEBRAIC[7]);
ALGEBRAIC[9] =  (ALGEBRAIC[8] - 1.00000)*CONSTANTS[27]+1.00000;
ALGEBRAIC[12] = CONSTANTS[37] - (CONSTANTS[37] - 1.00000)/( ((CONSTANTS[36] - 1.00000)/(CONSTANTS[36] - CONSTANTS[37]))*(STATES[3] - 1.00000)+1.00000);
ALGEBRAIC[13] = (CONDVAR[13]<0.00000 ? CONSTANTS[36] : ALGEBRAIC[12]);
ALGEBRAIC[14] = CONSTANTS[41] - CONSTANTS[41]/( 0.555556*(1.00000+STATES[4]));
ALGEBRAIC[15] = (CONDVAR[17]<0.00000 ? - 1.00000 : ALGEBRAIC[14]);
ALGEBRAIC[17] = CONSTANTS[56] - (CONSTANTS[56] - 1.00000)/exp( CONSTANTS[58]*(STATES[6] - 1.00000));
ALGEBRAIC[18] = (CONDVAR[30]<0.00000 ? CONSTANTS[57] : ALGEBRAIC[17]);
ALGEBRAIC[20] = ALGEBRAIC[18] - 1.00000;
ALGEBRAIC[25] =  ALGEBRAIC[20]*CONSTANTS[66]+1.00000;
ALGEBRAIC[26] = STATES[8]/STATES[7];
ALGEBRAIC[27] = (CONDVAR[31]<0.00000 ? 0.00000 : ALGEBRAIC[26] - CONSTANTS[70]);
ALGEBRAIC[28] = CONSTANTS[71] -  pow(ALGEBRAIC[27], CONSTANTS[72])*CONSTANTS[73];
ALGEBRAIC[29] =  0.280000*ALGEBRAIC[26]+ 0.00190000*pow(ALGEBRAIC[26], 2.00000);
ALGEBRAIC[30] = ((((((STATES[7]+STATES[13]) - STATES[16]) - STATES[17]) - STATES[18]) - STATES[19]) - STATES[20])/2.00000;
ALGEBRAIC[36] = STATES[19]+ ALGEBRAIC[30]*0.155000;
ALGEBRAIC[37] = ALGEBRAIC[36] - 0.306250;
ALGEBRAIC[38] = ALGEBRAIC[37]/0.00480000;
ALGEBRAIC[51] = STATES[17]+ ALGEBRAIC[30]*0.261000;
ALGEBRAIC[52] = ALGEBRAIC[51] - 0.495000;
ALGEBRAIC[53] = ALGEBRAIC[52]/0.00355000;
ALGEBRAIC[31] = STATES[20]+ ALGEBRAIC[30]*0.0574000;
ALGEBRAIC[32] = ALGEBRAIC[31] - 0.100000;
ALGEBRAIC[33] = ALGEBRAIC[32]/0.00500000;
ALGEBRAIC[54] = ALGEBRAIC[53] - ALGEBRAIC[33];
ALGEBRAIC[56] =  ALGEBRAIC[54]*CONSTANTS[99];
ALGEBRAIC[68] = (CONDVAR[94]<0.00000 ? 0.000100000 : STATES[30]);
ALGEBRAIC[69] =  CONSTANTS[119]*ALGEBRAIC[68];
ALGEBRAIC[75] = (CONDVAR[101]<0.00000 ? 0.00000 : STATES[34]);
ALGEBRAIC[76] =  ALGEBRAIC[75]*0.486110;
ALGEBRAIC[77] = (CONDVAR[102]<0.00000 ? 0.00100000 : STATES[35]);
ALGEBRAIC[85] = CONSTANTS[132]/(CONSTANTS[133]+ALGEBRAIC[77]);
ALGEBRAIC[86] = STATES[37]+1.00000;
ALGEBRAIC[87] = STATES[7]+STATES[13];
ALGEBRAIC[88] = STATES[13]/ALGEBRAIC[87];
ALGEBRAIC[89] =  100.000*ALGEBRAIC[88];
ALGEBRAIC[90] = ALGEBRAIC[89]/( (CONSTANTS[134] - ALGEBRAIC[89])*CONSTANTS[135]);
ALGEBRAIC[91] = ALGEBRAIC[90]+1.50000;
ALGEBRAIC[92] =  0.333300*ALGEBRAIC[91];
ALGEBRAIC[93] =  (CONSTANTS[208] - 40.0000)*ALGEBRAIC[89];
ALGEBRAIC[94] = (CONDVAR[111]<0.00000 ? 0.00000 : ALGEBRAIC[93]+CONSTANTS[209]);
ALGEBRAIC[95] = CONSTANTS[137] - ALGEBRAIC[94];
ALGEBRAIC[96] = (CONDVAR[112]<0.00000 ? 0.00000 :  ALGEBRAIC[95]*CONSTANTS[138]*CONSTANTS[139]+5.00000e-06);
ALGEBRAIC[97] =  STATES[13]*CONSTANTS[140]*ALGEBRAIC[92];
ALGEBRAIC[98] = (ALGEBRAIC[96] - ALGEBRAIC[97])+CONSTANTS[141];
ALGEBRAIC[99] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[146]+1.00000;
ALGEBRAIC[100] =  pow(CONSTANTS[147] - ALGEBRAIC[76], 2.00000)*CONSTANTS[148]*ALGEBRAIC[99];
ALGEBRAIC[101] = (CONDVAR[113]<0.00000 ? 0.800000 : CONDVAR[114]>0.00000 ? 8.00000 : ALGEBRAIC[100]);
ALGEBRAIC[103] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[150]*0.00100000;
ALGEBRAIC[102] =  (STATES[3] - 1.00000)*CONSTANTS[149]+1.00000;
ALGEBRAIC[104] =  ALGEBRAIC[102]*ALGEBRAIC[101]*0.00100000;
ALGEBRAIC[105] = (CONDVAR[115]<0.00000 ? 0.00000 : ALGEBRAIC[104]);
ALGEBRAIC[106] =  (ALGEBRAIC[103]+ALGEBRAIC[105])*1.00000;
ALGEBRAIC[107] = (CONDVAR[116]<0.00000 ? 0.00000 : ALGEBRAIC[106]);
ALGEBRAIC[108] =  pow(fabs(ALGEBRAIC[33]), CONSTANTS[153])*CONSTANTS[154];
ALGEBRAIC[109] = (CONDVAR[117]<0.00000 ? - ALGEBRAIC[108] : ALGEBRAIC[108]);
ALGEBRAIC[113] = (CONDVAR[118]>0.00000&&CONDVAR[119]<=0.00000 ? CONSTANTS[159] : CONDVAR[120]>0.00000 ? STATES[39] : ALGEBRAIC[53] - CONSTANTS[158]);
ALGEBRAIC[24] = pow( ALGEBRAIC[20]*CONSTANTS[64]+1.00000, CONSTANTS[65]);
ALGEBRAIC[114] =  (ALGEBRAIC[24] - 1.00000)*CONSTANTS[167]+1.00000;
ALGEBRAIC[115] = (CONDVAR[123]<0.00000 ? 0.800000 : ALGEBRAIC[114]);
ALGEBRAIC[116] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[168]+1.00000;
ALGEBRAIC[117] = (CONDVAR[124]<0.00000 ? CONSTANTS[169] : ALGEBRAIC[116]);
ALGEBRAIC[118] =  (ALGEBRAIC[115] - 1.00000)*CONSTANTS[173]+1.00000;
ALGEBRAIC[119] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[174]+1.00000;
ALGEBRAIC[120] = STATES[41]/STATES[22];
ALGEBRAIC[121] =  ( (ALGEBRAIC[6] - 1.00000)*CONSTANTS[198]+1.00000)*0.100000;
ALGEBRAIC[122] = (CONDVAR[134]<0.00000 ? 0.00000 : ALGEBRAIC[121]);
ALGEBRAIC[123] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[201]+1.00000;
ALGEBRAIC[124] = (CONDVAR[136]<0.00000 ? CONSTANTS[202] : ALGEBRAIC[123]);
ALGEBRAIC[43] = STATES[18]+ ALGEBRAIC[30]*0.128000;
ALGEBRAIC[44] = ALGEBRAIC[43] - 0.380000;
ALGEBRAIC[45] = ALGEBRAIC[44]/0.0100000;
ALGEBRAIC[126] = (CONDVAR[14]<0.00000 ? 0.00000 :  (ALGEBRAIC[45] - 1.00000)*1.00000);
ALGEBRAIC[127] = (CONDVAR[15]<0.00000 ? 0.00000 :  (ALGEBRAIC[33]+1.00000)*2.00000);
ALGEBRAIC[128] = (ALGEBRAIC[126]+ALGEBRAIC[127])/3.00000;
ALGEBRAIC[129] = (CONDVAR[16]>0.00000 ? CONSTANTS[38] : ALGEBRAIC[128]+CONSTANTS[39]);
ALGEBRAIC[130] = (CONDVAR[18]>0.00000 ? CONSTANTS[42] : ALGEBRAIC[53] - CONSTANTS[43]);
ALGEBRAIC[132] = (CONDVAR[25]<0.00000&&CONDVAR[26]>=0.00000 ?  0.0166670*(160.000 - ALGEBRAIC[130]) : CONDVAR[27]<0.00000 ? 1.33360 : 0.00000);
ALGEBRAIC[133] =  (ALGEBRAIC[132] - 1.00000)*CONSTANTS[46]+1.00000;
ALGEBRAIC[134] = ALGEBRAIC[133] - CONSTANTS[49];
ALGEBRAIC[61] = STATES[22] - STATES[25];
ALGEBRAIC[137] = (ALGEBRAIC[61] - STATES[7]) - ALGEBRAIC[77];
ALGEBRAIC[139] = STATES[10]/ALGEBRAIC[137];
ALGEBRAIC[140] =  0.280000*ALGEBRAIC[139]+ 0.00190000*pow(ALGEBRAIC[139], 2.00000);
ALGEBRAIC[141] = pow((CONSTANTS[79]/ALGEBRAIC[137])/5.00000, CONSTANTS[80]);
ALGEBRAIC[142] =  ALGEBRAIC[141]*2.00000;
ALGEBRAIC[143] =  ALGEBRAIC[142]*ALGEBRAIC[140]*CONSTANTS[82];
ALGEBRAIC[19] = (CONSTANTS[59] -  ALGEBRAIC[18]*CONSTANTS[60])+CONSTANTS[60];
ALGEBRAIC[110] = ALGEBRAIC[109] - STATES[38];
ALGEBRAIC[112] =  ALGEBRAIC[110]*CONSTANTS[157];
ALGEBRAIC[48] = STATES[16]+ ALGEBRAIC[30]*0.398600;
ALGEBRAIC[49] =  (ALGEBRAIC[8] - 1.00000)*CONSTANTS[89];
ALGEBRAIC[153] = ((((ALGEBRAIC[48] - ALGEBRAIC[19]) - ALGEBRAIC[49]) - STATES[15]) - STATES[14]) - ALGEBRAIC[112];
ALGEBRAIC[154] = (CONDVAR[88]<0.00000 ? 0.000100000 : ALGEBRAIC[153]);
ALGEBRAIC[155] = 3.70000+(ALGEBRAIC[154] - 0.740000)/CONSTANTS[90];
ALGEBRAIC[156] = (CONDVAR[89]<0.00000 ? 0.000100000 : ALGEBRAIC[155]);
ALGEBRAIC[157] = 0.740000/pow(ALGEBRAIC[156]/( ALGEBRAIC[92]*3.70000), 0.500000);
ALGEBRAIC[50] = (CONDVAR[90]<0.00000 ? CONSTANTS[91] : ALGEBRAIC[33]);
ALGEBRAIC[158] = ALGEBRAIC[156] - ALGEBRAIC[50];
ALGEBRAIC[159] = ALGEBRAIC[158]/ALGEBRAIC[157];
ALGEBRAIC[70] = (CONDVAR[95]<0.00000 ? 0.500000 : STATES[31]);
ALGEBRAIC[71] = (CONDVAR[96]<0.00000 ? 0.500000 : STATES[32]);
ALGEBRAIC[72] = (CONDVAR[97]<0.00000 ? 0.300000 : STATES[33]);
ALGEBRAIC[73] =  ALGEBRAIC[70]*ALGEBRAIC[71]*ALGEBRAIC[72];
ALGEBRAIC[74] =  (ALGEBRAIC[73] - 1.00000)*CONSTANTS[127]+1.00000;
ALGEBRAIC[10] = CONSTANTS[34] - (CONSTANTS[34] - 1.00000)/( ((CONSTANTS[35] - 1.00000)/(CONSTANTS[35] - CONSTANTS[34]))*(STATES[3] - 1.00000)+1.00000);
ALGEBRAIC[11] = (CONDVAR[12]<0.00000 ? CONSTANTS[35] : ALGEBRAIC[10]);
ALGEBRAIC[111] =  ALGEBRAIC[110]*CONSTANTS[156]+1.00000;
ALGEBRAIC[55] = pow(ALGEBRAIC[53]/100.000, CONSTANTS[95]);
ALGEBRAIC[160] = (( ALGEBRAIC[8]*ALGEBRAIC[11]*ALGEBRAIC[24]*ALGEBRAIC[92]*CONSTANTS[1])/ALGEBRAIC[55])/ALGEBRAIC[111];
ALGEBRAIC[161] =  CONSTANTS[96]*ALGEBRAIC[74]*ALGEBRAIC[160]*CONSTANTS[0]*CONSTANTS[97];
ALGEBRAIC[162] = ALGEBRAIC[53] - ALGEBRAIC[156];
ALGEBRAIC[66] = (CONDVAR[92]<0.00000 ? CONSTANTS[112] : STATES[27]);
ALGEBRAIC[67] =  ALGEBRAIC[66]*STATES[28];
ALGEBRAIC[163] =  CONSTANTS[98]*ALGEBRAIC[67]*ALGEBRAIC[160]*CONSTANTS[0]*CONSTANTS[97];
ALGEBRAIC[164] = ALGEBRAIC[162]/ALGEBRAIC[163];
ALGEBRAIC[165] = STATES[23]/ALGEBRAIC[61];
ALGEBRAIC[57] =  (ALGEBRAIC[3] - 1.00000)*CONSTANTS[102]+1.00000;
ALGEBRAIC[58] = (STATES[24] - 3000.00)/( ALGEBRAIC[57]*9.33330);
ALGEBRAIC[166] = ALGEBRAIC[58]/ALGEBRAIC[61];
ALGEBRAIC[59] = STATES[24] - ALGEBRAIC[58];
ALGEBRAIC[60] = ALGEBRAIC[59]/STATES[25];
ALGEBRAIC[167] = ALGEBRAIC[60] - ALGEBRAIC[165];
ALGEBRAIC[168] =  ALGEBRAIC[167]*CONSTANTS[104];
ALGEBRAIC[169] =  (ALGEBRAIC[76] - 10.0000)*CONSTANTS[108];
ALGEBRAIC[170] = ALGEBRAIC[69] - 38.0000;
ALGEBRAIC[171] =  ALGEBRAIC[170]*CONSTANTS[110]+1.00000;
ALGEBRAIC[172] =  CONSTANTS[113]*ALGEBRAIC[170]+1.00000;
ALGEBRAIC[173] =  STATES[29]*ALGEBRAIC[164];
ALGEBRAIC[23] =  ALGEBRAIC[20]*CONSTANTS[63]+1.00000;
ALGEBRAIC[174] = (CONDVAR[93]>0.00000 ? 38.0000 : ALGEBRAIC[69]);
ALGEBRAIC[175] =  ALGEBRAIC[23]*CONSTANTS[117]*CONSTANTS[115]*(1.00000 - pow(38.0001 - ALGEBRAIC[174], 3.00000)/54872.0);
ALGEBRAIC[176] = ALGEBRAIC[76] - CONSTANTS[120];
ALGEBRAIC[177] =  ALGEBRAIC[176]*CONSTANTS[121]+1.00000;
ALGEBRAIC[178] =  CONSTANTS[123]*ALGEBRAIC[176]+1.00000;
ALGEBRAIC[179] =  CONSTANTS[125]*ALGEBRAIC[176]+1.00000;
ALGEBRAIC[180] = (CONDVAR[98]>0.00000 ? 35.0000 : ALGEBRAIC[76]);
ALGEBRAIC[181] =  ALGEBRAIC[23]*CONSTANTS[128]*(1.00000 - pow(35.0001 - ALGEBRAIC[180], 3.00000)/42875.0);
ALGEBRAIC[39] =  0.0260000*ALGEBRAIC[38];
ALGEBRAIC[40] = (CONDVAR[66]<0.00000 ? 1.00000e-05 : ALGEBRAIC[39]);
ALGEBRAIC[41] = pow(ALGEBRAIC[40], 0.500000);
ALGEBRAIC[42] = 1.00000/ALGEBRAIC[41];
ALGEBRAIC[147] = ALGEBRAIC[45]+18.0000;
ALGEBRAIC[148] = 1.00000/( ALGEBRAIC[147]*0.0357000);
ALGEBRAIC[182] = ( (ALGEBRAIC[38] - ALGEBRAIC[45])*ALGEBRAIC[148])/(ALGEBRAIC[148]+ALGEBRAIC[42])+ALGEBRAIC[45];
ALGEBRAIC[79] = (CONDVAR[103]<0.00000 ? 0.0250000 : STATES[36]);
ALGEBRAIC[80] = ALGEBRAIC[79]/ALGEBRAIC[77];
ALGEBRAIC[81] =  ALGEBRAIC[80]*0.400000;
ALGEBRAIC[78] = 2.00000 - 0.150000/ALGEBRAIC[77];
ALGEBRAIC[183] =  (((ALGEBRAIC[182] - ALGEBRAIC[78])+ALGEBRAIC[81]) - ALGEBRAIC[29])*CONSTANTS[129];
ALGEBRAIC[83] =  (ALGEBRAIC[78]+11.0000)*0.000300000;
ALGEBRAIC[184] = ALGEBRAIC[183] - ALGEBRAIC[83];
ALGEBRAIC[185] = ALGEBRAIC[184];
ALGEBRAIC[84] =  ALGEBRAIC[83]*ALGEBRAIC[80];
ALGEBRAIC[82] =  (ALGEBRAIC[26] - ALGEBRAIC[80])*0.000225000;
ALGEBRAIC[186] = ALGEBRAIC[82] - ALGEBRAIC[84];
ALGEBRAIC[187] = ALGEBRAIC[186];
ALGEBRAIC[188] = (STATES[29]/ALGEBRAIC[89])/5.25000;
ALGEBRAIC[189] = (CONDVAR[104]>0.00000 ? 114.000+ (ALGEBRAIC[188] - 1.00000)*6667.00 : CONDVAR[105]>0.00000&&CONDVAR[106]<=0.00000 ? 74.0000+ (ALGEBRAIC[188] - 0.936000)*625.000 : CONDVAR[107]>0.00000&&CONDVAR[108]<=0.00000 ? 46.0000+ (ALGEBRAIC[188] - 0.800000)*205.882 :  ALGEBRAIC[188]*57.5000);
ALGEBRAIC[190] = (ALGEBRAIC[189] - 67.0000)/30.0000;
ALGEBRAIC[191] = (CONDVAR[109]>0.00000 ? 1.00000 : CONDVAR[110]<0.00000 ? 0.600000 : ALGEBRAIC[190]);
ALGEBRAIC[192] = 1.00000/ALGEBRAIC[191];
ALGEBRAIC[193] =  ( (ALGEBRAIC[192] - 1.00000)*3.00000 - STATES[37])*0.000500000;
ALGEBRAIC[125] =  (ALGEBRAIC[6] - 1.00000)*CONSTANTS[3]+1.00000;
ALGEBRAIC[194] = (ALGEBRAIC[166] - 3.30000)/1.00000;
ALGEBRAIC[195] =  ALGEBRAIC[125]*0.909000*ALGEBRAIC[194];
ALGEBRAIC[196] =  (ALGEBRAIC[195] - 1.00000)*CONSTANTS[4]+1.00000;
ALGEBRAIC[197] = (CONDVAR[0]<0.00000 ? 0.00000 : ALGEBRAIC[196]);
ALGEBRAIC[198] = (CONDVAR[1]>0.00000 ? CONSTANTS[6] : ALGEBRAIC[197]+CONSTANTS[5]);
ALGEBRAIC[199] = (ALGEBRAIC[165] - CONSTANTS[29])/(142.000 - CONSTANTS[29]);
ALGEBRAIC[200] = (CONDVAR[9]<0.00000 ? 0.00000 : ALGEBRAIC[199]);
ALGEBRAIC[201] = (CONDVAR[10]>0.00000 ? CONSTANTS[30] : ALGEBRAIC[130]);
ALGEBRAIC[202] =  pow(CONSTANTS[30] - ALGEBRAIC[201], 2.00000)*CONSTANTS[31];
ALGEBRAIC[203] = ALGEBRAIC[200]+ALGEBRAIC[202]+CONSTANTS[32];
ALGEBRAIC[204] = (CONDVAR[11]<0.00000 ? 0.00000 : ALGEBRAIC[203]);
ALGEBRAIC[131] = (CONDVAR[19]<0.00000&&CONDVAR[20]>=0.00000 ?  0.00500000*(80.0000 - ALGEBRAIC[130])*CONSTANTS[44] : CONDVAR[21]<0.00000 ?  0.200000*CONSTANTS[44] : 0.00000);
ALGEBRAIC[205] = (CONDVAR[22]<0.00000&&CONDVAR[23]>=0.00000 ?  CONSTANTS[45]*(80.0000 - ALGEBRAIC[189]) : CONDVAR[24]<0.00000 ?  CONSTANTS[45]*40.0000 : 0.00000);
ALGEBRAIC[206] = ALGEBRAIC[131]+ALGEBRAIC[205];
ALGEBRAIC[16] = STATES[5];
ALGEBRAIC[135] = (CONDVAR[28]<0.00000 ?  0.0400000*(40.0000 - ALGEBRAIC[130])*CONSTANTS[50] : 0.00000);
ALGEBRAIC[136] =  (15.0000/(ALGEBRAIC[45]+ALGEBRAIC[33]+ALGEBRAIC[38]) - 1.00000)*CONSTANTS[51]+1.00000;
ALGEBRAIC[207] =  ( CONSTANTS[206]*ALGEBRAIC[136]*(ALGEBRAIC[206]+ALGEBRAIC[16]+ALGEBRAIC[135]) - 1.00000)*CONSTANTS[54]+1.00000;
ALGEBRAIC[208] = (CONDVAR[29]<0.00000 ? 0.00000 : ALGEBRAIC[207]);
ALGEBRAIC[209] = (ALGEBRAIC[208] - STATES[6])/CONSTANTS[55];
ALGEBRAIC[138] = ALGEBRAIC[137] - STATES[9];
ALGEBRAIC[144] =  pow((ALGEBRAIC[138] - CONSTANTS[83])/CONSTANTS[83], 2.00000)*1.00000;
ALGEBRAIC[210] =  ALGEBRAIC[141]*CONSTANTS[81]+ALGEBRAIC[144];
ALGEBRAIC[211] = ALGEBRAIC[210] - ALGEBRAIC[142];
ALGEBRAIC[213] = (ALGEBRAIC[211]+CONSTANTS[84]) - ALGEBRAIC[144];
ALGEBRAIC[214] = (CONDVAR[33]>0.00000 ? 7.00000 : ALGEBRAIC[213]);
ALGEBRAIC[215] = (CONDVAR[34]<0.00000 ? 0.00000 :  ALGEBRAIC[214]*0.0200000);
ALGEBRAIC[21] =  ALGEBRAIC[20]*CONSTANTS[61]+1.00000;
ALGEBRAIC[62] = (CONDVAR[91]>0.00000 ? 1.00000 : STATES[26]);
ALGEBRAIC[46] =  (ALGEBRAIC[45]+4.00000)*( CONSTANTS[85]*(ALGEBRAIC[18] - 1.00000)+1.00000) - 4.00000;
ALGEBRAIC[47] = (CONDVAR[78]<=0.00000 ? 0.0100000 : CONDVAR[79]>0.00000&&CONDVAR[80]<=0.00000 ? 0.0100000+( (3.60000 - 0.0100000)*(ALGEBRAIC[46] - - 2.00000))/(1.00000 - - 2.00000) : CONDVAR[81]>0.00000&&CONDVAR[82]<=0.00000 ? 3.60000+( (9.40000 - 3.60000)*(ALGEBRAIC[46] - 1.00000))/(5.00000 - 1.00000) : CONDVAR[83]>0.00000&&CONDVAR[84]<=0.00000 ? 9.40000+( (11.6000 - 9.40000)*(ALGEBRAIC[46] - 5.00000))/(8.00000 - 5.00000) : CONDVAR[85]>0.00000&&CONDVAR[86]<=0.00000 ? 11.6000+( (13.5000 - 11.6000)*(ALGEBRAIC[46] - 8.00000))/(12.0000 - 8.00000) : 13.5000);
ALGEBRAIC[218] = ALGEBRAIC[53]/( ALGEBRAIC[21]*ALGEBRAIC[188]);
ALGEBRAIC[219] = (CONDVAR[67]<=0.00000 ? 1.04000 : CONDVAR[68]>0.00000&&CONDVAR[69]<=0.00000 ? 1.04000+( (1.02500 - 1.04000)*(ALGEBRAIC[218] - 0.00000))/(60.0000 - 0.00000) : CONDVAR[70]>0.00000&&CONDVAR[71]<=0.00000 ? 1.02500+( (0.970000 - 1.02500)*(ALGEBRAIC[218] - 60.0000))/(125.000 - 60.0000) : CONDVAR[72]>0.00000&&CONDVAR[73]<=0.00000 ? 0.970000+( (0.880000 - 0.970000)*(ALGEBRAIC[218] - 125.000))/(160.000 - 125.000) : CONDVAR[74]>0.00000&&CONDVAR[75]<=0.00000 ? 0.880000+( (0.590000 - 0.880000)*(ALGEBRAIC[218] - 160.000))/(200.000 - 160.000) : CONDVAR[76]>0.00000&&CONDVAR[77]<=0.00000 ? 0.590000+( (0.00000 - 0.590000)*(ALGEBRAIC[218] - 200.000))/(240.000 - 200.000) : 0.00000);
ALGEBRAIC[220] =  ALGEBRAIC[219]*ALGEBRAIC[47]*ALGEBRAIC[21]*CONSTANTS[88]*ALGEBRAIC[62]*STATES[11];
ALGEBRAIC[152] = (ALGEBRAIC[45] - ALGEBRAIC[53])/3.00000;
ALGEBRAIC[221] = (CONDVAR[87]>0.00000 ? ALGEBRAIC[220]+ALGEBRAIC[152] : ALGEBRAIC[220]);
ALGEBRAIC[223] =  ALGEBRAIC[139]*ALGEBRAIC[215];
ALGEBRAIC[216] = (ALGEBRAIC[38]/ALGEBRAIC[21])/ALGEBRAIC[188];
ALGEBRAIC[217] = (CONDVAR[46]<=0.00000 ? 1.06000 : CONDVAR[47]>0.00000&&CONDVAR[48]<=0.00000 ? 1.06000+( (0.970000 - 1.06000)*(ALGEBRAIC[216] - 0.00000))/(32.0000 - 0.00000) : CONDVAR[49]>0.00000&&CONDVAR[50]<=0.00000 ? 0.970000+( (0.930000 - 0.970000)*(ALGEBRAIC[216] - 32.0000))/(38.4000 - 32.0000) : CONDVAR[51]>0.00000&&CONDVAR[52]<=0.00000 ? 0.930000+( (0.800000 - 0.930000)*(ALGEBRAIC[216] - 38.4000))/(48.0000 - 38.4000) : CONDVAR[53]>0.00000&&CONDVAR[54]<=0.00000 ? 0.800000+( (0.460000 - 0.800000)*(ALGEBRAIC[216] - 48.0000))/(60.8000 - 48.0000) : CONDVAR[55]>0.00000&&CONDVAR[56]<=0.00000 ? 0.460000+( (0.00000 - 0.460000)*(ALGEBRAIC[216] - 60.8000))/(72.0000 - 60.8000) : 0.00000);
ALGEBRAIC[224] =  (1.00000 - CONSTANTS[86])*ALGEBRAIC[21]*ALGEBRAIC[217]*CONSTANTS[87]*ALGEBRAIC[62]*STATES[12]+( CONSTANTS[86]*ALGEBRAIC[221])/ALGEBRAIC[47];
ALGEBRAIC[34] =  (ALGEBRAIC[33]+8.00000)*( CONSTANTS[85]*(ALGEBRAIC[18] - 1.00000)+1.00000) - 8.00000;
ALGEBRAIC[35] = (CONDVAR[57]<=0.00000 ? 0.00000 : CONDVAR[58]>0.00000&&CONDVAR[59]<=0.00000 ? 0.00000+( (0.750000 - 0.00000)*(ALGEBRAIC[34] - - 8.00000))/(- 6.00000 - - 8.00000) : CONDVAR[60]>0.00000&&CONDVAR[61]<=0.00000 ? 0.750000+( (2.60000 - 0.750000)*(ALGEBRAIC[34] - - 6.00000))/(- 2.00000 - - 6.00000) : CONDVAR[62]>0.00000&&CONDVAR[63]<=0.00000 ? 2.60000+( (9.80000 - 2.60000)*(ALGEBRAIC[34] - - 2.00000))/(4.00000 - - 2.00000) : CONDVAR[64]>0.00000&&CONDVAR[65]<=0.00000 ? 9.80000+( (13.5000 - 9.80000)*(ALGEBRAIC[34] - 4.00000))/(12.0000 - 4.00000) : 13.5000);
ALGEBRAIC[225] =  ALGEBRAIC[35]*ALGEBRAIC[224];
ALGEBRAIC[149] = ALGEBRAIC[148]+ALGEBRAIC[42];
ALGEBRAIC[150] = ALGEBRAIC[38] - ALGEBRAIC[45];
ALGEBRAIC[151] = ALGEBRAIC[150]/ALGEBRAIC[149];
ALGEBRAIC[226] = ALGEBRAIC[225] - ALGEBRAIC[151];
ALGEBRAIC[227] = ALGEBRAIC[151] - ALGEBRAIC[221];
ALGEBRAIC[228] = ALGEBRAIC[159] - ALGEBRAIC[225];
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[43] = 0.0;
SI[44] = 0.0;
SI[45] = 0.0;
SI[1] = 1.0;
SI[46] = 0.0;
SI[47] = 0.0;
SI[48] = 0.0;
SI[49] = 0.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
SI[5] = 1.0;
SI[6] = 1.0;
SI[50] = 0.0;
SI[51] = 0.0;
SI[52] = 0.0;
SI[7] = 1.0;
SI[53] = 0.0;
SI[54] = 0.0;
SI[55] = 0.0;
SI[8] = 1.0;
SI[56] = 0.0;
SI[57] = 0.0;
SI[58] = 0.0;
SI[9] = 1.0;
SI[10] = 1.0;
SI[59] = 0.0;
SI[11] = 1.0;
SI[12] = 1.0;
SI[60] = 0.0;
SI[13] = 1.0;
SI[14] = 1.0;
SI[15] = 1.0;
SI[61] = 0.0;
SI[16] = 1.0;
SI[17] = 1.0;
SI[18] = 1.0;
SI[19] = 1.0;
SI[20] = 1.0;
SI[62] = 0.0;
SI[63] = 0.0;
SI[64] = 0.0;
SI[65] = 0.0;
SI[66] = 0.0;
SI[67] = 0.0;
SI[68] = 0.0;
SI[69] = 0.0;
SI[21] = 1.0;
SI[70] = 0.0;
SI[71] = 0.0;
SI[72] = 0.0;
SI[22] = 1.0;
SI[73] = 0.0;
SI[23] = 1.0;
SI[24] = 1.0;
SI[74] = 0.0;
SI[25] = 1.0;
SI[26] = 1.0;
SI[27] = 1.0;
SI[28] = 1.0;
SI[29] = 1.0;
SI[75] = 0.0;
SI[76] = 0.0;
SI[77] = 0.0;
SI[78] = 0.0;
SI[30] = 1.0;
SI[31] = 1.0;
SI[32] = 1.0;
SI[33] = 1.0;
SI[79] = 0.0;
SI[80] = 0.0;
SI[81] = 0.0;
SI[82] = 0.0;
SI[83] = 0.0;
SI[84] = 0.0;
SI[34] = 1.0;
SI[35] = 1.0;
SI[36] = 1.0;
SI[85] = 0.0;
SI[86] = 0.0;
SI[87] = 0.0;
SI[88] = 0.0;
SI[89] = 0.0;
SI[37] = 1.0;
SI[38] = 1.0;
SI[39] = 1.0;
SI[90] = 0.0;
SI[91] = 0.0;
SI[92] = 0.0;
SI[40] = 1.0;
SI[93] = 0.0;
SI[94] = 0.0;
SI[95] = 0.0;
SI[96] = 0.0;
SI[97] = 0.0;
SI[98] = 0.0;
SI[99] = 0.0;
SI[100] = 0.0;
SI[101] = 0.0;
SI[102] = 0.0;
SI[103] = 0.0;
SI[104] = 0.0;
SI[105] = 0.0;
SI[106] = 0.0;
SI[107] = 0.0;
SI[108] = 0.0;
SI[109] = 0.0;
SI[110] = 0.0;
SI[111] = 0.0;
SI[112] = 0.0;
SI[113] = 0.0;
SI[114] = 0.0;
SI[41] = 1.0;
SI[115] = 0.0;
SI[116] = 0.0;
SI[117] = 0.0;
SI[118] = 0.0;
SI[119] = 0.0;
SI[120] = 0.0;
SI[121] = 0.0;
SI[122] = 0.0;
SI[123] = 0.0;
SI[124] = 0.0;
SI[125] = 0.0;
SI[126] = 0.0;
SI[127] = 0.0;
SI[128] = 0.0;
SI[129] = 0.0;
SI[42] = 1.0;
SI[130] = 0.0;
SI[131] = 0.0;
SI[132] = 0.0;
SI[133] = 0.0;
SI[134] = 0.0;
SI[135] = 0.0;
SI[136] = 0.0;
SI[137] = 0.0;
SI[138] = 0.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = ALGEBRAIC[196] - 0.00000;
CONDVAR[1] = CONSTANTS[6] - 0.00000;
CONDVAR[2] = ALGEBRAIC[2] - 0.200000;
CONDVAR[3] = ALGEBRAIC[4] - CONSTANTS[15];
CONDVAR[4] = ALGEBRAIC[4] - CONSTANTS[14];
CONDVAR[5] = STATES[45] - 1.00000;
CONDVAR[6] = STATES[48] - 1.00000e-05;
CONDVAR[7] = CONSTANTS[19] - 0.00000;
CONDVAR[8] = ALGEBRAIC[7] - CONSTANTS[26];
CONDVAR[9] = ALGEBRAIC[199] - 0.00000;
CONDVAR[10] = ALGEBRAIC[130] - CONSTANTS[30];
CONDVAR[11] = ALGEBRAIC[203] - 0.00000;
CONDVAR[12] = ALGEBRAIC[10] - CONSTANTS[35];
CONDVAR[13] = ALGEBRAIC[12] - CONSTANTS[36];
CONDVAR[14] =  (ALGEBRAIC[45] - 1.00000)*1.00000 - 0.00000;
CONDVAR[15] =  (ALGEBRAIC[33]+1.00000)*2.00000 - 0.00000;
CONDVAR[16] = CONSTANTS[38] - 0.00000;
CONDVAR[17] = ALGEBRAIC[14] - - 1.00000;
CONDVAR[18] = CONSTANTS[42] - 1.00000e-07;
CONDVAR[19] = ALGEBRAIC[130] - 80.0000;
CONDVAR[20] = ALGEBRAIC[130] - 40.0000;
CONDVAR[21] = ALGEBRAIC[130] - 40.0000;
CONDVAR[22] = ALGEBRAIC[189] - 80.0000;
CONDVAR[23] = ALGEBRAIC[189] - 40.0000;
CONDVAR[24] = ALGEBRAIC[189] - 40.0000;
CONDVAR[25] = ALGEBRAIC[130] - 160.000;
CONDVAR[26] = ALGEBRAIC[130] - 80.0000;
CONDVAR[27] = ALGEBRAIC[130] - 80.0000;
CONDVAR[28] = ALGEBRAIC[130] - 40.0000;
CONDVAR[29] = ALGEBRAIC[207] - 0.00000;
CONDVAR[30] = ALGEBRAIC[17] - CONSTANTS[57];
CONDVAR[31] = (ALGEBRAIC[26] - CONSTANTS[70]) - 0.00000;
CONDVAR[32] = (STATES[52] - CONSTANTS[74]) - 0.00000;
CONDVAR[33] = ALGEBRAIC[213] - 7.00000;
CONDVAR[34] = ALGEBRAIC[214] - 0.00000;
CONDVAR[35] = ALGEBRAIC[137] - 0.00000;
CONDVAR[36] = ALGEBRAIC[137] - 0.00000;
CONDVAR[37] = ALGEBRAIC[137] - 12.0000;
CONDVAR[38] = ALGEBRAIC[137] - 12.0000;
CONDVAR[39] = ALGEBRAIC[137] - 15.0000;
CONDVAR[40] = ALGEBRAIC[137] - 15.0000;
CONDVAR[41] = ALGEBRAIC[137] - 18.0000;
CONDVAR[42] = ALGEBRAIC[137] - 18.0000;
CONDVAR[43] = ALGEBRAIC[137] - 21.0000;
CONDVAR[44] = ALGEBRAIC[137] - 21.0000;
CONDVAR[45] = ALGEBRAIC[137] - 24.0000;
CONDVAR[46] = ALGEBRAIC[216] - 0.00000;
CONDVAR[47] = ALGEBRAIC[216] - 0.00000;
CONDVAR[48] = ALGEBRAIC[216] - 32.0000;
CONDVAR[49] = ALGEBRAIC[216] - 32.0000;
CONDVAR[50] = ALGEBRAIC[216] - 38.4000;
CONDVAR[51] = ALGEBRAIC[216] - 38.4000;
CONDVAR[52] = ALGEBRAIC[216] - 48.0000;
CONDVAR[53] = ALGEBRAIC[216] - 48.0000;
CONDVAR[54] = ALGEBRAIC[216] - 60.8000;
CONDVAR[55] = ALGEBRAIC[216] - 60.8000;
CONDVAR[56] = ALGEBRAIC[216] - 72.0000;
CONDVAR[57] = ALGEBRAIC[34] - - 8.00000;
CONDVAR[58] = ALGEBRAIC[34] - - 8.00000;
CONDVAR[59] = ALGEBRAIC[34] - - 6.00000;
CONDVAR[60] = ALGEBRAIC[34] - - 6.00000;
CONDVAR[61] = ALGEBRAIC[34] - - 2.00000;
CONDVAR[62] = ALGEBRAIC[34] - - 2.00000;
CONDVAR[63] = ALGEBRAIC[34] - 4.00000;
CONDVAR[64] = ALGEBRAIC[34] - 4.00000;
CONDVAR[65] = ALGEBRAIC[34] - 12.0000;
CONDVAR[66] = ALGEBRAIC[39] - 1.00000e-05;
CONDVAR[67] = ALGEBRAIC[218] - 0.00000;
CONDVAR[68] = ALGEBRAIC[218] - 0.00000;
CONDVAR[69] = ALGEBRAIC[218] - 60.0000;
CONDVAR[70] = ALGEBRAIC[218] - 60.0000;
CONDVAR[71] = ALGEBRAIC[218] - 125.000;
CONDVAR[72] = ALGEBRAIC[218] - 125.000;
CONDVAR[73] = ALGEBRAIC[218] - 160.000;
CONDVAR[74] = ALGEBRAIC[218] - 160.000;
CONDVAR[75] = ALGEBRAIC[218] - 200.000;
CONDVAR[76] = ALGEBRAIC[218] - 200.000;
CONDVAR[77] = ALGEBRAIC[218] - 240.000;
CONDVAR[78] = ALGEBRAIC[46] - - 2.00000;
CONDVAR[79] = ALGEBRAIC[46] - - 2.00000;
CONDVAR[80] = ALGEBRAIC[46] - 1.00000;
CONDVAR[81] = ALGEBRAIC[46] - 1.00000;
CONDVAR[82] = ALGEBRAIC[46] - 5.00000;
CONDVAR[83] = ALGEBRAIC[46] - 5.00000;
CONDVAR[84] = ALGEBRAIC[46] - 8.00000;
CONDVAR[85] = ALGEBRAIC[46] - 8.00000;
CONDVAR[86] = ALGEBRAIC[46] - 12.0000;
CONDVAR[87] = ALGEBRAIC[152] - 0.00000;
CONDVAR[88] = ALGEBRAIC[153] - 0.000100000;
CONDVAR[89] = ALGEBRAIC[155] - 0.000100000;
CONDVAR[90] = ALGEBRAIC[33] - CONSTANTS[91];
CONDVAR[91] = STATES[26] - 1.00000;
CONDVAR[92] = STATES[27] - CONSTANTS[112];
CONDVAR[93] = ALGEBRAIC[69] - 38.0000;
CONDVAR[94] = STATES[30] - 0.000100000;
CONDVAR[95] = STATES[31] - 0.500000;
CONDVAR[96] = STATES[32] - 0.500000;
CONDVAR[97] = STATES[33] - 0.300000;
CONDVAR[98] = ALGEBRAIC[76] - 35.0000;
CONDVAR[99] = ALGEBRAIC[75] - 6.00000;
CONDVAR[100] = STATES[84] - 0.00000;
CONDVAR[101] = STATES[34] - 0.00000;
CONDVAR[102] = STATES[35] - 0.00100000;
CONDVAR[103] = STATES[36] - 0.0250000;
CONDVAR[104] = ALGEBRAIC[188] - 1.00000;
CONDVAR[105] = ALGEBRAIC[188] - 0.936000;
CONDVAR[106] = ALGEBRAIC[188] - 1.00000;
CONDVAR[107] = ALGEBRAIC[188] - 0.800000;
CONDVAR[108] = ALGEBRAIC[188] - 0.936000;
CONDVAR[109] = ALGEBRAIC[190] - 1.00000;
CONDVAR[110] = ALGEBRAIC[190] - 0.600000;
CONDVAR[111] = (ALGEBRAIC[93]+CONSTANTS[209]) - 0.00000;
CONDVAR[112] = ( ALGEBRAIC[95]*CONSTANTS[138]*CONSTANTS[139]+5.00000e-06) - 0.00000;
CONDVAR[113] = ALGEBRAIC[100] - 0.800000;
CONDVAR[114] = ALGEBRAIC[100] - 8.00000;
CONDVAR[115] = ALGEBRAIC[104] - 0.00000;
CONDVAR[116] = ALGEBRAIC[106] - 0.00000;
CONDVAR[117] = ALGEBRAIC[33] - 0.00000;
CONDVAR[118] = CONSTANTS[159] - 0.00000;
CONDVAR[119] = CONSTANTS[160] - 0.00000;
CONDVAR[120] = CONSTANTS[160] - 0.00000;
CONDVAR[121] = STATES[92] - CONSTANTS[164];
CONDVAR[122] = STATES[92] - CONSTANTS[165];
CONDVAR[123] = ALGEBRAIC[114] - 0.800000;
CONDVAR[124] = ALGEBRAIC[116] - CONSTANTS[169];
CONDVAR[125] = STATES[95] - CONSTANTS[172];
CONDVAR[126] = STATES[98] - CONSTANTS[177];
CONDVAR[127] = STATES[104] - 1.00000;
CONDVAR[128] = CONSTANTS[180] - 0.00000;
CONDVAR[129] = STATES[109] - CONSTANTS[183];
CONDVAR[130] = STATES[111] - 0.00000;
CONDVAR[131] = STATES[113] - 1.00000;
CONDVAR[132] = STATES[120] - 0.000100000;
CONDVAR[133] = STATES[123] - CONSTANTS[197];
CONDVAR[134] = ALGEBRAIC[121] - 0.00000;
CONDVAR[135] = STATES[127] - 0.100000;
CONDVAR[136] = ALGEBRAIC[123] - CONSTANTS[202];
CONDVAR[137] = STATES[132] - 1.00000e-08;
CONDVAR[138] = STATES[133] - 0.00000;
CONDVAR[139] = STATES[135] - 0.600000;
CONDVAR[140] = STATES[137] - 0.00000;
}