Generated Code
The following is c code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/*
There are a total of 86 entries in the algebraic variable array.
There are a total of 44 entries in each of the rate and state variable arrays.
There are a total of 96 entries in the constant variable array.
*/
/*
* VOI is time in component Circ_Environment (ms).
* CONSTANTS[0] is HR in component Circ_Environment (ratepm).
* CONSTANTS[91] is hrf in component Circ_Environment (Hz).
* ALGEBRAIC[0] is beattime in component LATiming (ms).
* CONSTANTS[1] is Esys in component LAElastanceFunction (elastance).
* CONSTANTS[2] is Edia in component LAElastanceFunction (elastance).
* CONSTANTS[3] is TsK in component LAElastanceFunction (ms).
* CONSTANTS[92] is Ts in component LAElastanceFunction (ms).
* ALGEBRAIC[1] is E_LA in component LAElastanceFunction (elastance).
* ALGEBRAIC[82] is P_PV in component PulmonaryVein (mmHg).
* ALGEBRAIC[9] is F_LV in component LeftVentricle (flow).
* CONSTANTS[4] is P_LA_ext in component LeftAtrium (mmHg).
* STATES[0] is V_LA in component LeftAtrium (mL).
* CONSTANTS[5] is V_LA_0 in component LeftAtrium (mL).
* ALGEBRAIC[3] is P_LA in component LeftAtrium (mmHg).
* ALGEBRAIC[83] is F_LA in component LeftAtrium (flow).
* CONSTANTS[6] is R_LA in component LeftAtrium (resistance).
* ALGEBRAIC[2] is conc_X_LA in component LeftAtrium (concentration).
* ALGEBRAIC[81] is conc_X_PV in component PulmonaryVein (concentration).
* STATES[1] is sub_X_LA in component LeftAtrium (mL).
* CONSTANTS[7] is PRint in component LVTiming (ms).
* ALGEBRAIC[4] is beattime in component LVTiming (ms).
* CONSTANTS[8] is Esys in component LVElastanceFunction (elastance).
* CONSTANTS[9] is Edia in component LVElastanceFunction (elastance).
* CONSTANTS[10] is TsK in component LVElastanceFunction (second).
* CONSTANTS[93] is Ts in component LVElastanceFunction (second).
* ALGEBRAIC[5] is E_LV in component LVElastanceFunction (elastance).
* ALGEBRAIC[12] is F_AO in component Aorta (flow).
* ALGEBRAIC[15] is F_CO in component CoronaryCirc (flow).
* CONSTANTS[11] is pLVPext in component LeftVentricle (mmHg).
* STATES[2] is V_LV in component LeftVentricle (mL).
* CONSTANTS[12] is V_LV_0 in component LeftVentricle (mL).
* ALGEBRAIC[7] is P_LV in component LeftVentricle (mmHg).
* CONSTANTS[13] is R_LV in component LeftVentricle (resistance).
* ALGEBRAIC[6] is conc_X_LV in component LeftVentricle (concentration).
* STATES[3] is sub_X_LV in component LeftVentricle (mL).
* ALGEBRAIC[8] is fiber_stress in component LeftVentricle (mmHg).
* ALGEBRAIC[19] is F_CR in component CarotidCirc (flow).
* ALGEBRAIC[22] is F_AR in component Arteries (flow).
* CONSTANTS[14] is C_AO in component Aorta (capacitance).
* CONSTANTS[15] is L_AO in component Aorta (inductance).
* CONSTANTS[16] is R_AO in component Aorta (resistance).
* CONSTANTS[17] is V_AO_0 in component Aorta (mL).
* STATES[4] is V_AO in component Aorta (mL).
* STATES[5] is F_AOalways in component Aorta (flow).
* ALGEBRAIC[11] is P_AO in component Aorta (mmHg).
* ALGEBRAIC[10] is conc_X_AO in component Aorta (concentration).
* STATES[6] is sub_X_AO in component Aorta (mL).
* ALGEBRAIC[62] is P_RA in component RightAtrium (mmHg).
* CONSTANTS[18] is R_CO in component CoronaryCirc (resistance).
* CONSTANTS[19] is R_COV in component CoronaryCirc (resistance).
* CONSTANTS[20] is V_CO_0 in component CoronaryCirc (mL).
* CONSTANTS[21] is C_CO in component CoronaryCirc (capacitance).
* ALGEBRAIC[63] is F_COV in component CoronaryCirc (flow).
* STATES[7] is V_CO in component CoronaryCirc (mL).
* ALGEBRAIC[14] is P_CO in component CoronaryCirc (mmHg).
* ALGEBRAIC[13] is conc_X_CO in component CoronaryCirc (concentration).
* STATES[8] is sub_X_CO in component CoronaryCirc (mL).
* ALGEBRAIC[55] is P_VC in component VenaCava (mmHg).
* CONSTANTS[22] is R_CR in component CarotidCirc (resistance).
* CONSTANTS[23] is R_CRV in component CarotidCirc (resistance).
* CONSTANTS[24] is V_CR_0 in component CarotidCirc (mL).
* CONSTANTS[25] is C_CR in component CarotidCirc (capacitance).
* ALGEBRAIC[56] is F_CRV in component CarotidCirc (flow).
* STATES[9] is V_CR in component CarotidCirc (mL).
* ALGEBRAIC[18] is P_CR in component CarotidCirc (mmHg).
* ALGEBRAIC[17] is conc_X_CR in component CarotidCirc (concentration).
* STATES[10] is sub_X_CR in component CarotidCirc (mL).
* ALGEBRAIC[16] is conc_X_CRtissue in component CarotidCirc (concentration).
* CONSTANTS[26] is D_CRtissue in component CarotidCirc (flow).
* CONSTANTS[27] is V_CRtissue in component CarotidCirc (mL).
* STATES[11] is sub_X_CRtissue in component CarotidCirc (mL).
* ALGEBRAIC[25] is F_AD in component Adipose (flow).
* ALGEBRAIC[28] is F_MU in component Muscle (flow).
* ALGEBRAIC[31] is F_GI in component GutIntestine (flow).
* ALGEBRAIC[35] is F_LI in component Liver (flow).
* ALGEBRAIC[38] is F_KI in component Kidney (flow).
* ALGEBRAIC[44] is F_OT in component OtherTissue (flow).
* ALGEBRAIC[41] is F_SK in component Skin (flow).
* CONSTANTS[28] is R_AR in component Arteries (resistance).
* CONSTANTS[29] is C_AR in component Arteries (capacitance).
* CONSTANTS[30] is V_AR_0 in component Arteries (mL).
* STATES[12] is V_AR in component Arteries (mL).
* ALGEBRAIC[21] is P_AR in component Arteries (mmHg).
* ALGEBRAIC[20] is conc_X_AR in component Arteries (concentration).
* STATES[13] is sub_X_AR in component Arteries (mL).
* ALGEBRAIC[46] is P_VE in component Veins (mmHg).
* CONSTANTS[31] is C_AD in component Adipose (capacitance).
* CONSTANTS[32] is R_AD in component Adipose (resistance).
* CONSTANTS[33] is R_ADV in component Adipose (resistance).
* CONSTANTS[34] is V_AD_0 in component Adipose (mL).
* STATES[14] is V_AD in component Adipose (mL).
* ALGEBRAIC[47] is F_ADV in component Adipose (flow).
* ALGEBRAIC[24] is P_AD in component Adipose (mmHg).
* ALGEBRAIC[23] is conc_X_AD in component Adipose (concentration).
* STATES[15] is sub_X_AD in component Adipose (mL).
* CONSTANTS[35] is C_MU in component Muscle (capacitance).
* CONSTANTS[36] is R_MU in component Muscle (resistance).
* CONSTANTS[37] is R_MUV in component Muscle (resistance).
* CONSTANTS[38] is V_MU_0 in component Muscle (mL).
* STATES[16] is V_MU in component Muscle (mL).
* ALGEBRAIC[48] is F_MUV in component Muscle (flow).
* ALGEBRAIC[27] is P_MU in component Muscle (mmHg).
* ALGEBRAIC[26] is conc_X_MU in component Muscle (concentration).
* STATES[17] is sub_X_MU in component Muscle (mL).
* ALGEBRAIC[33] is P_LI in component Liver (mmHg).
* CONSTANTS[39] is C_GI in component GutIntestine (capacitance).
* CONSTANTS[40] is R_GI in component GutIntestine (resistance).
* CONSTANTS[41] is R_GIV in component GutIntestine (resistance).
* CONSTANTS[42] is V_GI_0 in component GutIntestine (mL).
* STATES[18] is V_GI in component GutIntestine (mL).
* ALGEBRAIC[34] is F_GIV in component GutIntestine (flow).
* ALGEBRAIC[30] is P_GI in component GutIntestine (mmHg).
* ALGEBRAIC[29] is conc_X_GI in component GutIntestine (concentration).
* STATES[19] is sub_X_GI in component GutIntestine (mL).
* CONSTANTS[43] is C_LI in component Liver (capacitance).
* CONSTANTS[44] is R_LI in component Liver (resistance).
* CONSTANTS[45] is R_LIV in component Liver (resistance).
* CONSTANTS[46] is V_LI_0 in component Liver (mL).
* STATES[20] is V_LI in component Liver (mL).
* ALGEBRAIC[49] is F_LIV in component Liver (flow).
* ALGEBRAIC[32] is conc_X_LI in component Liver (concentration).
* STATES[21] is sub_X_LI in component Liver (mL).
* CONSTANTS[47] is metabolism in component Liver (flow).
* CONSTANTS[48] is C_KI in component Kidney (capacitance).
* CONSTANTS[49] is R_KI in component Kidney (resistance).
* CONSTANTS[50] is V_KI_0 in component Kidney (mL).
* CONSTANTS[51] is R_KIV in component Kidney (resistance).
* STATES[22] is V_KI in component Kidney (mL).
* ALGEBRAIC[50] is F_KIV in component Kidney (flow).
* ALGEBRAIC[37] is P_KI in component Kidney (mmHg).
* ALGEBRAIC[36] is conc_X_KI in component Kidney (concentration).
* STATES[23] is sub_X_KI in component Kidney (mL).
* CONSTANTS[52] is C_SK in component Skin (capacitance).
* CONSTANTS[53] is R_SK in component Skin (resistance).
* CONSTANTS[54] is R_SKV in component Skin (resistance).
* CONSTANTS[55] is V_SK_0 in component Skin (mL).
* STATES[24] is V_SK in component Skin (mL).
* ALGEBRAIC[51] is F_SKV in component Skin (flow).
* ALGEBRAIC[40] is P_SK in component Skin (mmHg).
* ALGEBRAIC[39] is conc_X_SK in component Skin (concentration).
* STATES[25] is sub_X_SK in component Skin (mL).
* CONSTANTS[56] is C_OT in component OtherTissue (capacitance).
* CONSTANTS[57] is R_OT in component OtherTissue (resistance).
* CONSTANTS[58] is R_OTV in component OtherTissue (resistance).
* CONSTANTS[59] is V_OT_0 in component OtherTissue (mL).
* STATES[26] is V_OT in component OtherTissue (mL).
* ALGEBRAIC[52] is F_OTV in component OtherTissue (flow).
* ALGEBRAIC[43] is P_OT in component OtherTissue (mmHg).
* ALGEBRAIC[42] is conc_X_OT in component OtherTissue (concentration).
* STATES[27] is sub_X_OT in component OtherTissue (mL).
* CONSTANTS[60] is R_VEV in component Veins (resistance).
* STATES[28] is V_VE in component Veins (mL).
* CONSTANTS[61] is C_VE in component Veins (capacitance).
* CONSTANTS[62] is V_VE_0 in component Veins (mL).
* ALGEBRAIC[53] is F_VE in component Veins (flow).
* ALGEBRAIC[57] is F_VEV in component Veins (flow).
* ALGEBRAIC[45] is conc_X_VE in component Veins (concentration).
* STATES[29] is sub_X_VE in component Veins (mL).
* CONSTANTS[63] is R_VCV in component VenaCava (resistance).
* STATES[30] is V_VC in component VenaCava (mL).
* CONSTANTS[64] is V_VC_0 in component VenaCava (mL).
* CONSTANTS[65] is C_VC in component VenaCava (capacitance).
* ALGEBRAIC[58] is F_VC in component VenaCava (flow).
* ALGEBRAIC[64] is F_VCV in component VenaCava (flow).
* ALGEBRAIC[54] is conc_X_VC in component VenaCava (concentration).
* STATES[31] is sub_X_VC in component VenaCava (mL).
* ALGEBRAIC[59] is beattime in component RATiming (ms).
* CONSTANTS[66] is Esys in component RAElastanceFunction (elastance).
* CONSTANTS[67] is Edia in component RAElastanceFunction (elastance).
* CONSTANTS[68] is TsK in component RAElastanceFunction (second).
* CONSTANTS[94] is Ts in component RAElastanceFunction (second).
* ALGEBRAIC[60] is E_RA in component RAElastanceFunction (elastance).
* ALGEBRAIC[70] is F_RV in component RightVentricle (flow).
* CONSTANTS[69] is P_RA_ext in component RightAtrium (mmHg).
* STATES[32] is V_RA in component RightAtrium (mL).
* CONSTANTS[70] is V_RA_0 in component RightAtrium (mL).
* ALGEBRAIC[65] is F_RA in component RightAtrium (flow).
* ALGEBRAIC[61] is conc_X_RA in component RightAtrium (concentration).
* STATES[33] is sub_X_RA in component RightAtrium (mL).
* CONSTANTS[71] is PRinterval in component RVTiming (second).
* ALGEBRAIC[66] is beattime in component RVTiming (ms).
* CONSTANTS[72] is Esys in component RVElastanceFunction (elastance).
* CONSTANTS[73] is Edia in component RVElastanceFunction (elastance).
* CONSTANTS[74] is TsK in component RVElastanceFunction (second).
* CONSTANTS[95] is Ts in component RVElastanceFunction (second).
* ALGEBRAIC[67] is E_RV in component RVElastanceFunction (elastance).
* ALGEBRAIC[73] is F_PA in component PulmonaryArtery (flow).
* CONSTANTS[75] is P_RV_ext in component RightVentricle (mmHg).
* STATES[34] is V_RV in component RightVentricle (mL).
* CONSTANTS[76] is V_RV_0 in component RightVentricle (mL).
* ALGEBRAIC[69] is P_RV in component RightVentricle (mmHg).
* CONSTANTS[77] is R_RV in component RightVentricle (resistance).
* ALGEBRAIC[68] is conc_X_RV in component RightVentricle (concentration).
* STATES[35] is sub_X_RV in component RightVentricle (mL).
* ALGEBRAIC[76] is F_Pa in component PulmonaryArteries (flow).
* CONSTANTS[78] is R_PA in component PulmonaryArtery (resistance).
* CONSTANTS[79] is C_PA in component PulmonaryArtery (capacitance).
* CONSTANTS[80] is V_PA_0 in component PulmonaryArtery (mL).
* STATES[36] is V_PA in component PulmonaryArtery (mL).
* ALGEBRAIC[72] is P_PA in component PulmonaryArtery (mmHg).
* ALGEBRAIC[71] is conc_X_PA in component PulmonaryArtery (concentration).
* STATES[37] is sub_X_PA in component PulmonaryArtery (mL).
* ALGEBRAIC[80] is F_PC in component PulmonaryCapillaries (flow).
* ALGEBRAIC[79] is F_SH in component Shunt (flow).
* STATES[38] is V_Pa in component PulmonaryArteries (mL).
* CONSTANTS[81] is C_Pa in component PulmonaryArteries (capacitance).
* CONSTANTS[82] is R_Pa in component PulmonaryArteries (resistance).
* CONSTANTS[83] is V_Pa_0 in component PulmonaryArteries (mL).
* ALGEBRAIC[75] is P_Pa in component PulmonaryArteries (mmHg).
* ALGEBRAIC[74] is conc_X_Pa in component PulmonaryArteries (concentration).
* STATES[39] is sub_X_Pa in component PulmonaryArteries (mL).
* CONSTANTS[84] is R_SH in component Shunt (resistance).
* ALGEBRAIC[78] is P_PC in component PulmonaryCapillaries (mmHg).
* CONSTANTS[85] is C_PC in component PulmonaryCapillaries (capacitance).
* CONSTANTS[86] is R_PC in component PulmonaryCapillaries (resistance).
* CONSTANTS[87] is R_PCV in component PulmonaryCapillaries (resistance).
* CONSTANTS[88] is V_PC_0 in component PulmonaryCapillaries (mL).
* STATES[40] is V_PC in component PulmonaryCapillaries (mL).
* ALGEBRAIC[84] is F_PCV in component PulmonaryCapillaries (flow).
* ALGEBRAIC[77] is conc_X_PC in component PulmonaryCapillaries (concentration).
* STATES[41] is sub_X_PC in component PulmonaryCapillaries (mL).
* CONSTANTS[89] is C_PV in component PulmonaryVein (capacitance).
* CONSTANTS[90] is V_PV_0 in component PulmonaryVein (mL).
* STATES[42] is V_PV in component PulmonaryVein (mL).
* ALGEBRAIC[85] is F_PV in component PulmonaryVein (flow).
* STATES[43] is sub_X_PV in component PulmonaryVein (mL).
* RATES[1] is d/dt sub_X_LA in component LeftAtrium (mL).
* RATES[0] is d/dt V_LA in component LeftAtrium (mL).
* RATES[3] is d/dt sub_X_LV in component LeftVentricle (mL).
* RATES[2] is d/dt V_LV in component LeftVentricle (mL).
* RATES[6] is d/dt sub_X_AO in component Aorta (mL).
* RATES[5] is d/dt F_AOalways in component Aorta (flow).
* RATES[4] is d/dt V_AO in component Aorta (mL).
* RATES[8] is d/dt sub_X_CO in component CoronaryCirc (mL).
* RATES[7] is d/dt V_CO in component CoronaryCirc (mL).
* RATES[11] is d/dt sub_X_CRtissue in component CarotidCirc (mL).
* RATES[10] is d/dt sub_X_CR in component CarotidCirc (mL).
* RATES[9] is d/dt V_CR in component CarotidCirc (mL).
* RATES[13] is d/dt sub_X_AR in component Arteries (mL).
* RATES[12] is d/dt V_AR in component Arteries (mL).
* RATES[15] is d/dt sub_X_AD in component Adipose (mL).
* RATES[14] is d/dt V_AD in component Adipose (mL).
* RATES[17] is d/dt sub_X_MU in component Muscle (mL).
* RATES[16] is d/dt V_MU in component Muscle (mL).
* RATES[19] is d/dt sub_X_GI in component GutIntestine (mL).
* RATES[18] is d/dt V_GI in component GutIntestine (mL).
* RATES[21] is d/dt sub_X_LI in component Liver (mL).
* RATES[20] is d/dt V_LI in component Liver (mL).
* RATES[23] is d/dt sub_X_KI in component Kidney (mL).
* RATES[22] is d/dt V_KI in component Kidney (mL).
* RATES[25] is d/dt sub_X_SK in component Skin (mL).
* RATES[24] is d/dt V_SK in component Skin (mL).
* RATES[27] is d/dt sub_X_OT in component OtherTissue (mL).
* RATES[26] is d/dt V_OT in component OtherTissue (mL).
* RATES[29] is d/dt sub_X_VE in component Veins (mL).
* RATES[28] is d/dt V_VE in component Veins (mL).
* RATES[31] is d/dt sub_X_VC in component VenaCava (mL).
* RATES[30] is d/dt V_VC in component VenaCava (mL).
* RATES[33] is d/dt sub_X_RA in component RightAtrium (mL).
* RATES[32] is d/dt V_RA in component RightAtrium (mL).
* RATES[35] is d/dt sub_X_RV in component RightVentricle (mL).
* RATES[34] is d/dt V_RV in component RightVentricle (mL).
* RATES[37] is d/dt sub_X_PA in component PulmonaryArtery (mL).
* RATES[36] is d/dt V_PA in component PulmonaryArtery (mL).
* RATES[39] is d/dt sub_X_Pa in component PulmonaryArteries (mL).
* RATES[38] is d/dt V_Pa in component PulmonaryArteries (mL).
* RATES[41] is d/dt sub_X_PC in component PulmonaryCapillaries (mL).
* RATES[40] is d/dt V_PC in component PulmonaryCapillaries (mL).
* RATES[43] is d/dt sub_X_PV in component PulmonaryVein (mL).
* RATES[42] is d/dt V_PV in component PulmonaryVein (mL).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 77;
CONSTANTS[1] = 0.145;
CONSTANTS[2] = 0.131;
CONSTANTS[3] = 0.2;
CONSTANTS[4] = 0;
STATES[0] = 144.038898121062;
CONSTANTS[5] = 15;
CONSTANTS[6] = 0.01;
STATES[1] = 0;
CONSTANTS[7] = 0.12;
CONSTANTS[8] = 5.6;
CONSTANTS[9] = 0.19;
CONSTANTS[10] = 0.35;
CONSTANTS[11] = 0;
STATES[2] = 118.75315693072;
CONSTANTS[12] = 30;
CONSTANTS[13] = 0.001;
STATES[3] = 0;
CONSTANTS[14] = 1.2;
CONSTANTS[15] = 5e-5;
CONSTANTS[16] = 0.01;
CONSTANTS[17] = 200;
STATES[4] = 302.714548324312;
STATES[5] = -6898.97886842467;
STATES[6] = 0;
CONSTANTS[18] = 1;
CONSTANTS[19] = 0.2;
CONSTANTS[20] = 10;
CONSTANTS[21] = 0.1;
STATES[7] = 11.4012164759732;
STATES[8] = 0;
CONSTANTS[22] = 5;
CONSTANTS[23] = 0.02;
CONSTANTS[24] = 10;
CONSTANTS[25] = 1;
STATES[9] = 27.2147903053333;
STATES[10] = 0;
CONSTANTS[26] = 1;
CONSTANTS[27] = 1;
STATES[11] = 0;
CONSTANTS[28] = 1.5;
CONSTANTS[29] = 0.5;
CONSTANTS[30] = 370;
STATES[12] = 383.557125288364;
STATES[13] = 0;
CONSTANTS[31] = 1.81;
CONSTANTS[32] = 1;
CONSTANTS[33] = 0.33;
CONSTANTS[34] = 75;
STATES[14] = 111.639088723426;
STATES[15] = 0;
CONSTANTS[35] = 1.81;
CONSTANTS[36] = 1;
CONSTANTS[37] = 0.33;
CONSTANTS[38] = 75;
STATES[16] = 111.639088723426;
STATES[17] = 0;
CONSTANTS[39] = 1.81;
CONSTANTS[40] = 1;
CONSTANTS[41] = 0.33;
CONSTANTS[42] = 75;
STATES[18] = 111.639088723426;
STATES[19] = 0;
CONSTANTS[43] = 1.81;
CONSTANTS[44] = 1;
CONSTANTS[45] = 0.33;
CONSTANTS[46] = 75;
STATES[20] = 111.639088723426;
STATES[21] = 0;
CONSTANTS[47] = 0;
CONSTANTS[48] = 1.81;
CONSTANTS[49] = 1;
CONSTANTS[50] = 75;
CONSTANTS[51] = 0.33;
STATES[22] = 111.639088723426;
STATES[23] = 0;
CONSTANTS[52] = 1.81;
CONSTANTS[53] = 1;
CONSTANTS[54] = 0.33;
CONSTANTS[55] = 75;
STATES[24] = 111.639088723426;
STATES[25] = 0;
CONSTANTS[56] = 1.81;
CONSTANTS[57] = 1;
CONSTANTS[58] = 0.33;
CONSTANTS[59] = 75;
STATES[26] = 111.639088723426;
STATES[27] = 0;
CONSTANTS[60] = 0.0223;
STATES[28] = 833.291684587657;
CONSTANTS[61] = 13.24;
CONSTANTS[62] = 596;
STATES[29] = 5;
CONSTANTS[63] = 0.0267;
STATES[30] = 1849.93788284599;
CONSTANTS[64] = 600;
CONSTANTS[65] = 73.88;
STATES[31] = 0;
CONSTANTS[66] = 0.11;
CONSTANTS[67] = 0.099;
CONSTANTS[68] = 0.2;
CONSTANTS[69] = 0;
STATES[32] = 152.150765432928;
CONSTANTS[70] = 15;
STATES[33] = 0;
CONSTANTS[71] = 0.12;
CONSTANTS[72] = 0.67;
CONSTANTS[73] = 0.104;
CONSTANTS[74] = 0.35;
CONSTANTS[75] = 0;
STATES[34] = 154.337916417774;
CONSTANTS[76] = 25;
CONSTANTS[77] = 0.002;
STATES[35] = 0;
CONSTANTS[78] = 0.0227;
CONSTANTS[79] = 2.222;
CONSTANTS[80] = 50;
STATES[36] = 112.059036699839;
STATES[37] = 0;
STATES[38] = 66.9910537203609;
CONSTANTS[81] = 1.481;
CONSTANTS[82] = 0.053;
CONSTANTS[83] = 30;
STATES[39] = 0;
CONSTANTS[84] = 2;
CONSTANTS[85] = 1.778;
CONSTANTS[86] = 0.0379;
CONSTANTS[87] = 0.0379;
CONSTANTS[88] = 53;
STATES[40] = 91.7784213531534;
STATES[41] = 0;
CONSTANTS[89] = 5;
CONSTANTS[90] = 150;
STATES[42] = 239.25897115602;
STATES[43] = 0;
CONSTANTS[91] = CONSTANTS[0]/0.0600000;
CONSTANTS[92] = CONSTANTS[3]* pow((CONSTANTS[91]/1000.00), 1.0 / 2);
CONSTANTS[93] = CONSTANTS[10]* pow((CONSTANTS[91]/1000.00), 1.0 / 2);
CONSTANTS[94] = CONSTANTS[68]* pow((CONSTANTS[91]/1000.00), 1.0 / 2);
CONSTANTS[95] = CONSTANTS[74]* pow((CONSTANTS[91]/1000.00), 1.0 / 2);
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[4] = (VOI - floor(VOI/CONSTANTS[91])*CONSTANTS[91])/1000.00 - CONSTANTS[7];
ALGEBRAIC[5] = (ALGEBRAIC[4]*0.00100000>=0.00000&&ALGEBRAIC[4]*0.00100000<=CONSTANTS[93] ? CONSTANTS[9]+( (CONSTANTS[8] - CONSTANTS[9])*(1.00000 - cos(( 3.14159265358979*(ALGEBRAIC[4]*0.00100000))/CONSTANTS[93])))/2.00000 : ALGEBRAIC[4]*0.00100000< 1.50000*CONSTANTS[93]&&ALGEBRAIC[4]*0.00100000>=CONSTANTS[93] ? CONSTANTS[9]+( (CONSTANTS[8] - CONSTANTS[9])*(1.00000+cos(( 2.00000* 3.14159265358979*(ALGEBRAIC[4]*0.00100000 - CONSTANTS[93]))/CONSTANTS[93])))/2.00000 : CONSTANTS[9]);
ALGEBRAIC[7] = ALGEBRAIC[5]*(STATES[2] - CONSTANTS[12])+CONSTANTS[11];
ALGEBRAIC[11] = (1.00000/CONSTANTS[14])*(STATES[4] - CONSTANTS[17]);
RATES[5] = ( (1.00000/CONSTANTS[15])*(ALGEBRAIC[7] - ALGEBRAIC[11]) - (CONSTANTS[16]/CONSTANTS[15])*STATES[5])/1000.00;
ALGEBRAIC[0] = (VOI - floor(VOI/CONSTANTS[91])*CONSTANTS[91])/1000.00;
ALGEBRAIC[1] = (ALGEBRAIC[0]>=0.00000&&ALGEBRAIC[0]<=CONSTANTS[92] ? CONSTANTS[2]+( (CONSTANTS[1] - CONSTANTS[2])*(1.00000 - cos(( 3.14159265358979*ALGEBRAIC[0])/CONSTANTS[92])))/2.00000 : ALGEBRAIC[0]< 1.50000*CONSTANTS[92]&&ALGEBRAIC[0]>=CONSTANTS[92] ? CONSTANTS[2]+( (CONSTANTS[1] - CONSTANTS[2])*(1.00000+cos(( 2.00000* 3.14159265358979*(ALGEBRAIC[0] - CONSTANTS[92]))/CONSTANTS[92])))/2.00000 : CONSTANTS[2]);
ALGEBRAIC[3] = ALGEBRAIC[1]*(STATES[0] - CONSTANTS[5])+CONSTANTS[4];
ALGEBRAIC[9] = (ALGEBRAIC[3]>ALGEBRAIC[7] ? (ALGEBRAIC[3] - ALGEBRAIC[7])/CONSTANTS[13] : 0.00000);
ALGEBRAIC[2] = STATES[1]/STATES[0];
ALGEBRAIC[12] = (ALGEBRAIC[7]>ALGEBRAIC[11] ? STATES[5] : 0.00000);
ALGEBRAIC[14] = (1.00000/CONSTANTS[21])*(STATES[7] - CONSTANTS[20]);
ALGEBRAIC[15] = (ALGEBRAIC[7]>ALGEBRAIC[14] ? (ALGEBRAIC[7] - ALGEBRAIC[14])/CONSTANTS[18] : 0.00000);
ALGEBRAIC[6] = STATES[3]/STATES[2];
RATES[3] = ( ALGEBRAIC[2]*ALGEBRAIC[9] - ALGEBRAIC[6]*(ALGEBRAIC[12]+ALGEBRAIC[15]))/1000.00;
RATES[2] = ((ALGEBRAIC[9] - ALGEBRAIC[12]) - ALGEBRAIC[15])/1000.00;
ALGEBRAIC[17] = STATES[10]/STATES[9];
ALGEBRAIC[16] = STATES[11]/CONSTANTS[27];
RATES[11] = ( (ALGEBRAIC[17] - ALGEBRAIC[16])*CONSTANTS[26])/1000.00;
ALGEBRAIC[18] = (1.00000/CONSTANTS[25])*(STATES[9] - CONSTANTS[24]);
ALGEBRAIC[19] = (ALGEBRAIC[11] - ALGEBRAIC[18])/CONSTANTS[22];
ALGEBRAIC[21] = (1.00000/CONSTANTS[29])*(STATES[12] - CONSTANTS[30]);
ALGEBRAIC[22] = (ALGEBRAIC[11] - ALGEBRAIC[21])/CONSTANTS[28];
ALGEBRAIC[10] = STATES[6]/STATES[4];
RATES[6] = ( ALGEBRAIC[6]*ALGEBRAIC[12] - ALGEBRAIC[10]*(ALGEBRAIC[19]+ALGEBRAIC[22]))/1000.00;
RATES[4] = ((ALGEBRAIC[12] - ALGEBRAIC[22]) - ALGEBRAIC[19])/1000.00;
ALGEBRAIC[30] = (1.00000/CONSTANTS[39])*(STATES[18] - CONSTANTS[42]);
ALGEBRAIC[31] = (ALGEBRAIC[21] - ALGEBRAIC[30])/CONSTANTS[40];
ALGEBRAIC[20] = STATES[13]/STATES[12];
ALGEBRAIC[33] = (1.00000/CONSTANTS[43])*(STATES[20] - CONSTANTS[46]);
ALGEBRAIC[34] = (ALGEBRAIC[30] - ALGEBRAIC[33])/CONSTANTS[41];
ALGEBRAIC[29] = STATES[19]/STATES[18];
RATES[19] = (( ALGEBRAIC[20]*ALGEBRAIC[31] - ALGEBRAIC[29]*ALGEBRAIC[34])/1000.00)*0.00100000;
RATES[18] = ((ALGEBRAIC[31] - ALGEBRAIC[34])/1000.00)*0.00100000;
ALGEBRAIC[24] = (1.00000/CONSTANTS[31])*(STATES[14] - CONSTANTS[34]);
ALGEBRAIC[25] = (ALGEBRAIC[21] - ALGEBRAIC[24])/CONSTANTS[32];
ALGEBRAIC[27] = (1.00000/CONSTANTS[35])*(STATES[16] - CONSTANTS[38]);
ALGEBRAIC[28] = (ALGEBRAIC[21] - ALGEBRAIC[27])/CONSTANTS[36];
ALGEBRAIC[35] = (ALGEBRAIC[21] - ALGEBRAIC[33])/CONSTANTS[44];
ALGEBRAIC[37] = (1.00000/CONSTANTS[48])*(STATES[22] - CONSTANTS[50]);
ALGEBRAIC[38] = (ALGEBRAIC[21] - ALGEBRAIC[37])/CONSTANTS[49];
ALGEBRAIC[43] = (1.00000/CONSTANTS[56])*(STATES[26] - CONSTANTS[59]);
ALGEBRAIC[44] = (ALGEBRAIC[21] - ALGEBRAIC[43])/CONSTANTS[57];
ALGEBRAIC[40] = (1.00000/CONSTANTS[52])*(STATES[24] - CONSTANTS[55]);
ALGEBRAIC[41] = (ALGEBRAIC[21] - ALGEBRAIC[40])/CONSTANTS[53];
RATES[13] = ( ALGEBRAIC[10]*ALGEBRAIC[22] - ALGEBRAIC[20]*(ALGEBRAIC[41]+ALGEBRAIC[25]+ALGEBRAIC[28]+ALGEBRAIC[31]+ALGEBRAIC[35]+ALGEBRAIC[38]+ALGEBRAIC[44]))/1000.00;
RATES[12] = (((((((ALGEBRAIC[22] - ALGEBRAIC[25]) - ALGEBRAIC[28]) - ALGEBRAIC[31]) - ALGEBRAIC[35]) - ALGEBRAIC[38]) - ALGEBRAIC[44]) - ALGEBRAIC[41])/1000.00;
ALGEBRAIC[46] = (1.00000/CONSTANTS[61])*(STATES[28] - CONSTANTS[62]);
ALGEBRAIC[47] = (ALGEBRAIC[24] - ALGEBRAIC[46])/CONSTANTS[33];
ALGEBRAIC[23] = STATES[15]/STATES[14];
RATES[15] = (( ALGEBRAIC[20]*ALGEBRAIC[25] - ALGEBRAIC[23]*ALGEBRAIC[47])/1000.00)*0.00100000;
RATES[14] = ((ALGEBRAIC[25] - ALGEBRAIC[47])/1000.00)*0.00100000;
ALGEBRAIC[48] = (ALGEBRAIC[27] - ALGEBRAIC[46])/CONSTANTS[37];
ALGEBRAIC[26] = STATES[17]/STATES[16];
RATES[17] = (( ALGEBRAIC[20]*ALGEBRAIC[28] - ALGEBRAIC[26]*ALGEBRAIC[48])/1000.00)*0.00100000;
RATES[16] = ((ALGEBRAIC[28] - ALGEBRAIC[48])/1000.00)*0.00100000;
ALGEBRAIC[49] = (ALGEBRAIC[33] - ALGEBRAIC[46])/CONSTANTS[45];
ALGEBRAIC[32] = STATES[21]/STATES[20];
RATES[21] = (((( ALGEBRAIC[29]*ALGEBRAIC[34]+ ALGEBRAIC[20]*ALGEBRAIC[35]) - ALGEBRAIC[32]*ALGEBRAIC[49])+ ALGEBRAIC[32]*CONSTANTS[47])/1000.00)*0.00100000;
RATES[20] = (((ALGEBRAIC[34]+ALGEBRAIC[35]) - ALGEBRAIC[49])/1000.00)*0.00100000;
ALGEBRAIC[50] = (ALGEBRAIC[37] - ALGEBRAIC[46])/CONSTANTS[51];
ALGEBRAIC[36] = STATES[23]/STATES[22];
RATES[23] = (( ALGEBRAIC[20]*ALGEBRAIC[38] - ALGEBRAIC[36]*ALGEBRAIC[50])/1000.00)*0.00100000;
RATES[22] = ((ALGEBRAIC[38] - ALGEBRAIC[50])/1000.00)*0.00100000;
ALGEBRAIC[51] = (ALGEBRAIC[40] - ALGEBRAIC[46])/CONSTANTS[54];
ALGEBRAIC[39] = STATES[25]/STATES[24];
RATES[25] = (( ALGEBRAIC[20]*ALGEBRAIC[41] - ALGEBRAIC[39]*ALGEBRAIC[51])/1000.00)*0.00100000;
RATES[24] = ((ALGEBRAIC[41] - ALGEBRAIC[51])/1000.00)*0.00100000;
ALGEBRAIC[52] = (ALGEBRAIC[43] - ALGEBRAIC[46])/CONSTANTS[58];
ALGEBRAIC[42] = STATES[27]/STATES[26];
RATES[27] = (( ALGEBRAIC[20]*ALGEBRAIC[44] - ALGEBRAIC[42]*ALGEBRAIC[52])/1000.00)*0.00100000;
RATES[26] = ((ALGEBRAIC[44] - ALGEBRAIC[52])/1000.00)*0.00100000;
ALGEBRAIC[55] = (1.00000/CONSTANTS[65])*(STATES[30] - CONSTANTS[64]);
ALGEBRAIC[56] = (ALGEBRAIC[18] - ALGEBRAIC[55])/CONSTANTS[23];
RATES[10] = (( ALGEBRAIC[10]*ALGEBRAIC[19] - ALGEBRAIC[17]*ALGEBRAIC[56]) - (ALGEBRAIC[17] - ALGEBRAIC[16])*CONSTANTS[26])/1000.00;
RATES[9] = (ALGEBRAIC[19] - ALGEBRAIC[56])/1000.00;
ALGEBRAIC[57] = (ALGEBRAIC[46] - ALGEBRAIC[55])/CONSTANTS[60];
ALGEBRAIC[45] = STATES[29]/STATES[28];
RATES[29] = (( ALGEBRAIC[39]*ALGEBRAIC[51]+ ALGEBRAIC[23]*ALGEBRAIC[47]+ ALGEBRAIC[26]*ALGEBRAIC[48]+ ALGEBRAIC[32]*ALGEBRAIC[49]+ ALGEBRAIC[36]*ALGEBRAIC[50]+ ALGEBRAIC[42]*ALGEBRAIC[52]) - ALGEBRAIC[45]*ALGEBRAIC[57])/1000.00;
ALGEBRAIC[53] = ALGEBRAIC[51]+ALGEBRAIC[47]+ALGEBRAIC[48]+ALGEBRAIC[49]+ALGEBRAIC[50]+ALGEBRAIC[52];
RATES[28] = (ALGEBRAIC[53] - ALGEBRAIC[57])/1000.00;
ALGEBRAIC[59] = (VOI - floor(VOI/CONSTANTS[91])*CONSTANTS[91])/1000.00;
ALGEBRAIC[60] = (ALGEBRAIC[59]*0.00100000>=0.00000&&ALGEBRAIC[59]*0.00100000<=CONSTANTS[94] ? CONSTANTS[67]+( (CONSTANTS[66] - CONSTANTS[67])*(1.00000 - cos(( 3.14159265358979*(ALGEBRAIC[59]*0.00100000))/CONSTANTS[94])))/2.00000 : ALGEBRAIC[59]*0.00100000< 1.50000*CONSTANTS[94]&&ALGEBRAIC[59]*0.00100000>=CONSTANTS[94] ? CONSTANTS[67]+( (CONSTANTS[66] - CONSTANTS[67])*(1.00000+cos(( 2.00000* 3.14159265358979*(ALGEBRAIC[59]*0.00100000 - CONSTANTS[94]))/CONSTANTS[94])))/2.00000 : CONSTANTS[67]);
ALGEBRAIC[62] = ALGEBRAIC[60]*(STATES[32] - CONSTANTS[70])+CONSTANTS[69];
ALGEBRAIC[63] = (ALGEBRAIC[14] - ALGEBRAIC[62])/CONSTANTS[19];
ALGEBRAIC[13] = STATES[8]/STATES[7];
RATES[8] = ( ALGEBRAIC[6]*ALGEBRAIC[15] - ALGEBRAIC[13]*ALGEBRAIC[63])/1000.00;
RATES[7] = (ALGEBRAIC[15] - ALGEBRAIC[63])/1000.00;
ALGEBRAIC[64] = (ALGEBRAIC[55] - ALGEBRAIC[62])/CONSTANTS[63];
ALGEBRAIC[54] = STATES[31]/STATES[30];
RATES[31] = (( ALGEBRAIC[45]*ALGEBRAIC[57]+ ALGEBRAIC[17]*ALGEBRAIC[56]) - ALGEBRAIC[54]*ALGEBRAIC[64])/1000.00;
ALGEBRAIC[58] = ALGEBRAIC[56]+ALGEBRAIC[57];
RATES[30] = (ALGEBRAIC[58] - ALGEBRAIC[64])/1000.00;
ALGEBRAIC[66] = (VOI - floor(VOI/CONSTANTS[91])*CONSTANTS[91])/1000.00 - CONSTANTS[71];
ALGEBRAIC[67] = (ALGEBRAIC[66]*0.00100000>=0.00000&&ALGEBRAIC[66]*0.00100000<=CONSTANTS[95] ? CONSTANTS[73]+( (CONSTANTS[72] - CONSTANTS[73])*(1.00000 - cos(( 3.14159265358979*(ALGEBRAIC[66]*0.00100000))/CONSTANTS[95])))/2.00000 : ALGEBRAIC[66]*0.00100000< 1.50000*CONSTANTS[95]&&ALGEBRAIC[66]*0.00100000>=CONSTANTS[95] ? CONSTANTS[73]+( (CONSTANTS[72] - CONSTANTS[73])*(1.00000+cos(( 2.00000* 3.14159265358979*(ALGEBRAIC[66]*0.00100000 - CONSTANTS[95]))/CONSTANTS[95])))/2.00000 : CONSTANTS[73]);
ALGEBRAIC[69] = ALGEBRAIC[67]*(STATES[34] - CONSTANTS[76])+CONSTANTS[75];
ALGEBRAIC[70] = (ALGEBRAIC[62]>ALGEBRAIC[69] ? (ALGEBRAIC[62] - ALGEBRAIC[69])/CONSTANTS[77] : 0.00000);
ALGEBRAIC[61] = STATES[33]/STATES[32];
RATES[33] = (( ALGEBRAIC[54]*ALGEBRAIC[64]+ ALGEBRAIC[13]*ALGEBRAIC[63]) - ALGEBRAIC[61]*ALGEBRAIC[70])/1000.00;
ALGEBRAIC[65] = ALGEBRAIC[64]+ALGEBRAIC[63];
RATES[32] = (ALGEBRAIC[65] - ALGEBRAIC[70])/1000.00;
ALGEBRAIC[72] = (1.00000/CONSTANTS[79])*(STATES[36] - CONSTANTS[80]);
ALGEBRAIC[73] = (ALGEBRAIC[69]>ALGEBRAIC[72] ? (ALGEBRAIC[69] - ALGEBRAIC[72])/CONSTANTS[78] : 0.00000);
ALGEBRAIC[68] = STATES[35]/STATES[34];
RATES[35] = ( ALGEBRAIC[61]*ALGEBRAIC[70] - ALGEBRAIC[68]*ALGEBRAIC[73])/1000.00;
RATES[34] = (ALGEBRAIC[70] - ALGEBRAIC[73])/1000.00;
ALGEBRAIC[75] = (1.00000/CONSTANTS[81])*(STATES[38] - CONSTANTS[83]);
ALGEBRAIC[76] = (ALGEBRAIC[72] - ALGEBRAIC[75])/CONSTANTS[82];
ALGEBRAIC[71] = STATES[37]/STATES[36];
RATES[37] = ( ALGEBRAIC[68]*ALGEBRAIC[73] - ALGEBRAIC[71]*ALGEBRAIC[76])/1000.00;
RATES[36] = (ALGEBRAIC[73] - ALGEBRAIC[76])/1000.00;
ALGEBRAIC[78] = (1.00000/CONSTANTS[85])*(STATES[40] - CONSTANTS[88]);
ALGEBRAIC[80] = (ALGEBRAIC[75] - ALGEBRAIC[78])/CONSTANTS[86];
ALGEBRAIC[79] = (ALGEBRAIC[75] - ALGEBRAIC[78])/CONSTANTS[84];
ALGEBRAIC[74] = STATES[39]/STATES[38];
RATES[39] = ( ALGEBRAIC[71]*ALGEBRAIC[76] - ALGEBRAIC[74]*(ALGEBRAIC[80]+ALGEBRAIC[79]))/1000.00;
RATES[38] = ((ALGEBRAIC[76] - ALGEBRAIC[80]) - ALGEBRAIC[79])/1000.00;
ALGEBRAIC[82] = (1.00000/CONSTANTS[89])*(STATES[42] - CONSTANTS[90]);
ALGEBRAIC[83] = (ALGEBRAIC[82] - ALGEBRAIC[3])/CONSTANTS[6];
ALGEBRAIC[81] = STATES[43]/STATES[42];
RATES[1] = ( ALGEBRAIC[81]*ALGEBRAIC[83] - ALGEBRAIC[2]*ALGEBRAIC[9])/1000.00;
RATES[0] = (ALGEBRAIC[83] - ALGEBRAIC[9])/1000.00;
ALGEBRAIC[84] = (ALGEBRAIC[78] - ALGEBRAIC[82])/CONSTANTS[87];
ALGEBRAIC[77] = STATES[41]/STATES[40];
RATES[41] = ( ALGEBRAIC[74]*ALGEBRAIC[80] - ALGEBRAIC[77]*ALGEBRAIC[84])/1000.00;
RATES[40] = (ALGEBRAIC[80] - ALGEBRAIC[84])/1000.00;
RATES[43] = (( ALGEBRAIC[77]*ALGEBRAIC[84]+ ALGEBRAIC[74]*ALGEBRAIC[79]) - ALGEBRAIC[81]*ALGEBRAIC[83])/1000.00;
ALGEBRAIC[85] = ALGEBRAIC[84]+ALGEBRAIC[79];
RATES[42] = (ALGEBRAIC[85] - ALGEBRAIC[83])/1000.00;
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[4] = (VOI - floor(VOI/CONSTANTS[91])*CONSTANTS[91])/1000.00 - CONSTANTS[7];
ALGEBRAIC[5] = (ALGEBRAIC[4]*0.00100000>=0.00000&&ALGEBRAIC[4]*0.00100000<=CONSTANTS[93] ? CONSTANTS[9]+( (CONSTANTS[8] - CONSTANTS[9])*(1.00000 - cos(( 3.14159265358979*(ALGEBRAIC[4]*0.00100000))/CONSTANTS[93])))/2.00000 : ALGEBRAIC[4]*0.00100000< 1.50000*CONSTANTS[93]&&ALGEBRAIC[4]*0.00100000>=CONSTANTS[93] ? CONSTANTS[9]+( (CONSTANTS[8] - CONSTANTS[9])*(1.00000+cos(( 2.00000* 3.14159265358979*(ALGEBRAIC[4]*0.00100000 - CONSTANTS[93]))/CONSTANTS[93])))/2.00000 : CONSTANTS[9]);
ALGEBRAIC[7] = ALGEBRAIC[5]*(STATES[2] - CONSTANTS[12])+CONSTANTS[11];
ALGEBRAIC[11] = (1.00000/CONSTANTS[14])*(STATES[4] - CONSTANTS[17]);
ALGEBRAIC[0] = (VOI - floor(VOI/CONSTANTS[91])*CONSTANTS[91])/1000.00;
ALGEBRAIC[1] = (ALGEBRAIC[0]>=0.00000&&ALGEBRAIC[0]<=CONSTANTS[92] ? CONSTANTS[2]+( (CONSTANTS[1] - CONSTANTS[2])*(1.00000 - cos(( 3.14159265358979*ALGEBRAIC[0])/CONSTANTS[92])))/2.00000 : ALGEBRAIC[0]< 1.50000*CONSTANTS[92]&&ALGEBRAIC[0]>=CONSTANTS[92] ? CONSTANTS[2]+( (CONSTANTS[1] - CONSTANTS[2])*(1.00000+cos(( 2.00000* 3.14159265358979*(ALGEBRAIC[0] - CONSTANTS[92]))/CONSTANTS[92])))/2.00000 : CONSTANTS[2]);
ALGEBRAIC[3] = ALGEBRAIC[1]*(STATES[0] - CONSTANTS[5])+CONSTANTS[4];
ALGEBRAIC[9] = (ALGEBRAIC[3]>ALGEBRAIC[7] ? (ALGEBRAIC[3] - ALGEBRAIC[7])/CONSTANTS[13] : 0.00000);
ALGEBRAIC[2] = STATES[1]/STATES[0];
ALGEBRAIC[12] = (ALGEBRAIC[7]>ALGEBRAIC[11] ? STATES[5] : 0.00000);
ALGEBRAIC[14] = (1.00000/CONSTANTS[21])*(STATES[7] - CONSTANTS[20]);
ALGEBRAIC[15] = (ALGEBRAIC[7]>ALGEBRAIC[14] ? (ALGEBRAIC[7] - ALGEBRAIC[14])/CONSTANTS[18] : 0.00000);
ALGEBRAIC[6] = STATES[3]/STATES[2];
ALGEBRAIC[17] = STATES[10]/STATES[9];
ALGEBRAIC[16] = STATES[11]/CONSTANTS[27];
ALGEBRAIC[18] = (1.00000/CONSTANTS[25])*(STATES[9] - CONSTANTS[24]);
ALGEBRAIC[19] = (ALGEBRAIC[11] - ALGEBRAIC[18])/CONSTANTS[22];
ALGEBRAIC[21] = (1.00000/CONSTANTS[29])*(STATES[12] - CONSTANTS[30]);
ALGEBRAIC[22] = (ALGEBRAIC[11] - ALGEBRAIC[21])/CONSTANTS[28];
ALGEBRAIC[10] = STATES[6]/STATES[4];
ALGEBRAIC[30] = (1.00000/CONSTANTS[39])*(STATES[18] - CONSTANTS[42]);
ALGEBRAIC[31] = (ALGEBRAIC[21] - ALGEBRAIC[30])/CONSTANTS[40];
ALGEBRAIC[20] = STATES[13]/STATES[12];
ALGEBRAIC[33] = (1.00000/CONSTANTS[43])*(STATES[20] - CONSTANTS[46]);
ALGEBRAIC[34] = (ALGEBRAIC[30] - ALGEBRAIC[33])/CONSTANTS[41];
ALGEBRAIC[29] = STATES[19]/STATES[18];
ALGEBRAIC[24] = (1.00000/CONSTANTS[31])*(STATES[14] - CONSTANTS[34]);
ALGEBRAIC[25] = (ALGEBRAIC[21] - ALGEBRAIC[24])/CONSTANTS[32];
ALGEBRAIC[27] = (1.00000/CONSTANTS[35])*(STATES[16] - CONSTANTS[38]);
ALGEBRAIC[28] = (ALGEBRAIC[21] - ALGEBRAIC[27])/CONSTANTS[36];
ALGEBRAIC[35] = (ALGEBRAIC[21] - ALGEBRAIC[33])/CONSTANTS[44];
ALGEBRAIC[37] = (1.00000/CONSTANTS[48])*(STATES[22] - CONSTANTS[50]);
ALGEBRAIC[38] = (ALGEBRAIC[21] - ALGEBRAIC[37])/CONSTANTS[49];
ALGEBRAIC[43] = (1.00000/CONSTANTS[56])*(STATES[26] - CONSTANTS[59]);
ALGEBRAIC[44] = (ALGEBRAIC[21] - ALGEBRAIC[43])/CONSTANTS[57];
ALGEBRAIC[40] = (1.00000/CONSTANTS[52])*(STATES[24] - CONSTANTS[55]);
ALGEBRAIC[41] = (ALGEBRAIC[21] - ALGEBRAIC[40])/CONSTANTS[53];
ALGEBRAIC[46] = (1.00000/CONSTANTS[61])*(STATES[28] - CONSTANTS[62]);
ALGEBRAIC[47] = (ALGEBRAIC[24] - ALGEBRAIC[46])/CONSTANTS[33];
ALGEBRAIC[23] = STATES[15]/STATES[14];
ALGEBRAIC[48] = (ALGEBRAIC[27] - ALGEBRAIC[46])/CONSTANTS[37];
ALGEBRAIC[26] = STATES[17]/STATES[16];
ALGEBRAIC[49] = (ALGEBRAIC[33] - ALGEBRAIC[46])/CONSTANTS[45];
ALGEBRAIC[32] = STATES[21]/STATES[20];
ALGEBRAIC[50] = (ALGEBRAIC[37] - ALGEBRAIC[46])/CONSTANTS[51];
ALGEBRAIC[36] = STATES[23]/STATES[22];
ALGEBRAIC[51] = (ALGEBRAIC[40] - ALGEBRAIC[46])/CONSTANTS[54];
ALGEBRAIC[39] = STATES[25]/STATES[24];
ALGEBRAIC[52] = (ALGEBRAIC[43] - ALGEBRAIC[46])/CONSTANTS[58];
ALGEBRAIC[42] = STATES[27]/STATES[26];
ALGEBRAIC[55] = (1.00000/CONSTANTS[65])*(STATES[30] - CONSTANTS[64]);
ALGEBRAIC[56] = (ALGEBRAIC[18] - ALGEBRAIC[55])/CONSTANTS[23];
ALGEBRAIC[57] = (ALGEBRAIC[46] - ALGEBRAIC[55])/CONSTANTS[60];
ALGEBRAIC[45] = STATES[29]/STATES[28];
ALGEBRAIC[53] = ALGEBRAIC[51]+ALGEBRAIC[47]+ALGEBRAIC[48]+ALGEBRAIC[49]+ALGEBRAIC[50]+ALGEBRAIC[52];
ALGEBRAIC[59] = (VOI - floor(VOI/CONSTANTS[91])*CONSTANTS[91])/1000.00;
ALGEBRAIC[60] = (ALGEBRAIC[59]*0.00100000>=0.00000&&ALGEBRAIC[59]*0.00100000<=CONSTANTS[94] ? CONSTANTS[67]+( (CONSTANTS[66] - CONSTANTS[67])*(1.00000 - cos(( 3.14159265358979*(ALGEBRAIC[59]*0.00100000))/CONSTANTS[94])))/2.00000 : ALGEBRAIC[59]*0.00100000< 1.50000*CONSTANTS[94]&&ALGEBRAIC[59]*0.00100000>=CONSTANTS[94] ? CONSTANTS[67]+( (CONSTANTS[66] - CONSTANTS[67])*(1.00000+cos(( 2.00000* 3.14159265358979*(ALGEBRAIC[59]*0.00100000 - CONSTANTS[94]))/CONSTANTS[94])))/2.00000 : CONSTANTS[67]);
ALGEBRAIC[62] = ALGEBRAIC[60]*(STATES[32] - CONSTANTS[70])+CONSTANTS[69];
ALGEBRAIC[63] = (ALGEBRAIC[14] - ALGEBRAIC[62])/CONSTANTS[19];
ALGEBRAIC[13] = STATES[8]/STATES[7];
ALGEBRAIC[64] = (ALGEBRAIC[55] - ALGEBRAIC[62])/CONSTANTS[63];
ALGEBRAIC[54] = STATES[31]/STATES[30];
ALGEBRAIC[58] = ALGEBRAIC[56]+ALGEBRAIC[57];
ALGEBRAIC[66] = (VOI - floor(VOI/CONSTANTS[91])*CONSTANTS[91])/1000.00 - CONSTANTS[71];
ALGEBRAIC[67] = (ALGEBRAIC[66]*0.00100000>=0.00000&&ALGEBRAIC[66]*0.00100000<=CONSTANTS[95] ? CONSTANTS[73]+( (CONSTANTS[72] - CONSTANTS[73])*(1.00000 - cos(( 3.14159265358979*(ALGEBRAIC[66]*0.00100000))/CONSTANTS[95])))/2.00000 : ALGEBRAIC[66]*0.00100000< 1.50000*CONSTANTS[95]&&ALGEBRAIC[66]*0.00100000>=CONSTANTS[95] ? CONSTANTS[73]+( (CONSTANTS[72] - CONSTANTS[73])*(1.00000+cos(( 2.00000* 3.14159265358979*(ALGEBRAIC[66]*0.00100000 - CONSTANTS[95]))/CONSTANTS[95])))/2.00000 : CONSTANTS[73]);
ALGEBRAIC[69] = ALGEBRAIC[67]*(STATES[34] - CONSTANTS[76])+CONSTANTS[75];
ALGEBRAIC[70] = (ALGEBRAIC[62]>ALGEBRAIC[69] ? (ALGEBRAIC[62] - ALGEBRAIC[69])/CONSTANTS[77] : 0.00000);
ALGEBRAIC[61] = STATES[33]/STATES[32];
ALGEBRAIC[65] = ALGEBRAIC[64]+ALGEBRAIC[63];
ALGEBRAIC[72] = (1.00000/CONSTANTS[79])*(STATES[36] - CONSTANTS[80]);
ALGEBRAIC[73] = (ALGEBRAIC[69]>ALGEBRAIC[72] ? (ALGEBRAIC[69] - ALGEBRAIC[72])/CONSTANTS[78] : 0.00000);
ALGEBRAIC[68] = STATES[35]/STATES[34];
ALGEBRAIC[75] = (1.00000/CONSTANTS[81])*(STATES[38] - CONSTANTS[83]);
ALGEBRAIC[76] = (ALGEBRAIC[72] - ALGEBRAIC[75])/CONSTANTS[82];
ALGEBRAIC[71] = STATES[37]/STATES[36];
ALGEBRAIC[78] = (1.00000/CONSTANTS[85])*(STATES[40] - CONSTANTS[88]);
ALGEBRAIC[80] = (ALGEBRAIC[75] - ALGEBRAIC[78])/CONSTANTS[86];
ALGEBRAIC[79] = (ALGEBRAIC[75] - ALGEBRAIC[78])/CONSTANTS[84];
ALGEBRAIC[74] = STATES[39]/STATES[38];
ALGEBRAIC[82] = (1.00000/CONSTANTS[89])*(STATES[42] - CONSTANTS[90]);
ALGEBRAIC[83] = (ALGEBRAIC[82] - ALGEBRAIC[3])/CONSTANTS[6];
ALGEBRAIC[81] = STATES[43]/STATES[42];
ALGEBRAIC[84] = (ALGEBRAIC[78] - ALGEBRAIC[82])/CONSTANTS[87];
ALGEBRAIC[77] = STATES[41]/STATES[40];
ALGEBRAIC[85] = ALGEBRAIC[84]+ALGEBRAIC[79];
ALGEBRAIC[8] = ALGEBRAIC[7]*(1.00000+( 3.00000*STATES[2])/200.000);
}