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 93 entries in the algebraic variable array.
There are a total of 33 entries in each of the rate and state variable arrays.
There are a total of 139 entries in the constant variable array.
*/
/*
* VOI is t in component main (second).
* CONSTANTS[0] is rho in component main (Js2_per_m5).
* CONSTANTS[1] is g in component main (m_per_s2).
* CONSTANTS[2] is theta_deg in component main (dim).
* CONSTANTS[138] is theta_rad in component main (radian).
* ALGEBRAIC[0] is z in component main (dim).
* CONSTANTS[3] is T in component main (second).
* ALGEBRAIC[2] is mt in component main (dim).
* ALGEBRAIC[1] is mta in component main (dim).
* CONSTANTS[4] is delay in component main (dim).
* CONSTANTS[5] is To_pk in component main (dim).
* ALGEBRAIC[4] is E_B_LV in component main (kPa_per_L).
* CONSTANTS[6] is E_B_LVmax in component main (kPa_per_L).
* CONSTANTS[7] is E_B_LVmin in component main (kPa_per_L).
* STATES[0] is q_B_LV in component main (litre).
* CONSTANTS[8] is q_B_LV_us in component main (litre).
* ALGEBRAIC[6] is u_B_LV in component main (kPa).
* ALGEBRAIC[8] is E_B_RV in component main (kPa_per_L).
* CONSTANTS[9] is E_B_RVmax in component main (kPa_per_L).
* CONSTANTS[10] is E_B_RVmin in component main (kPa_per_L).
* STATES[1] is q_B_RV in component main (litre).
* CONSTANTS[11] is q_B_RV_us in component main (litre).
* ALGEBRAIC[10] is u_B_RV in component main (kPa).
* CONSTANTS[12] is E_B_LA in component main (kPa_per_L).
* STATES[2] is q_B_LA in component main (litre).
* CONSTANTS[13] is q_B_LA_us in component main (litre).
* ALGEBRAIC[12] is u_B_LA in component main (kPa).
* CONSTANTS[14] is E_B_RA in component main (kPa_per_L).
* STATES[3] is q_B_RA in component main (litre).
* CONSTANTS[15] is q_B_RA_us in component main (litre).
* ALGEBRAIC[14] is u_B_RA in component main (kPa).
* CONSTANTS[16] is R_B_AV in component main (kPa_s_per_L).
* ALGEBRAIC[29] is v_B_AV in component main (L_per_s).
* CONSTANTS[17] is R_B_PV in component main (kPa_s_per_L).
* ALGEBRAIC[20] is v_B_PV in component main (L_per_s).
* CONSTANTS[18] is R_B_TV in component main (kPa_s_per_L).
* ALGEBRAIC[16] is v_B_TV in component main (L_per_s).
* CONSTANTS[19] is R_B_MV in component main (kPa_s_per_L).
* ALGEBRAIC[17] is v_B_MV in component main (L_per_s).
* STATES[4] is q_B_PA in component main (litre).
* CONSTANTS[20] is q_B_PA_us in component main (litre).
* ALGEBRAIC[18] is u_B_PA in component main (kPa).
* CONSTANTS[21] is E_B_PA in component main (kPa_per_L).
* STATES[5] is q_B_lung in component main (litre).
* CONSTANTS[22] is q_B_lung_us in component main (litre).
* ALGEBRAIC[22] is u_B_lung in component main (kPa).
* CONSTANTS[23] is E_B_lung in component main (kPa_per_L).
* CONSTANTS[24] is z_lung in component main (meter).
* ALGEBRAIC[25] is v_B_lung1 in component main (L_per_s).
* ALGEBRAIC[28] is v_B_lung2 in component main (L_per_s).
* CONSTANTS[25] is R_B_lung1 in component main (kPa_s_per_L).
* CONSTANTS[26] is R_B_lung2 in component main (kPa_s_per_L).
* CONSTANTS[27] is E_B_pulmVein in component main (kPa_per_L).
* STATES[6] is q_B_pulmVein in component main (litre).
* CONSTANTS[28] is q_B_pulmVein_us in component main (litre).
* ALGEBRAIC[27] is u_B_pulmVein in component main (kPa).
* ALGEBRAIC[30] is v_B_pulmVein in component main (L_per_s).
* CONSTANTS[29] is R_B_pulmVein in component main (kPa_s_per_L).
* STATES[7] is q_B_brain in component main (litre).
* CONSTANTS[30] is q_B_brain_us in component main (litre).
* STATES[8] is q_B_brainVein in component main (litre).
* CONSTANTS[31] is q_B_brainVein_us in component main (litre).
* ALGEBRAIC[19] is u_B_brain in component main (kPa).
* ALGEBRAIC[21] is u_B_brainVein in component main (kPa).
* CONSTANTS[32] is E_B_brain in component main (kPa_per_L).
* CONSTANTS[33] is E_B_brainVein in component main (kPa_per_L).
* ALGEBRAIC[31] is v_B_brain1 in component main (L_per_s).
* ALGEBRAIC[23] is v_B_brain2 in component main (L_per_s).
* ALGEBRAIC[24] is v_B_brain3 in component main (L_per_s).
* CONSTANTS[34] is R_B_brain1 in component main (kPa_s_per_L).
* CONSTANTS[35] is R_B_brain2 in component main (kPa_s_per_L).
* CONSTANTS[36] is R_B_brain3 in component main (kPa_s_per_L).
* STATES[9] is q_B_AA in component main (litre).
* CONSTANTS[37] is q_B_AA_us in component main (litre).
* ALGEBRAIC[26] is u_B_AA in component main (kPa).
* CONSTANTS[38] is E_B_AA in component main (kPa_per_L).
* STATES[10] is q_B_celiac in component main (litre).
* CONSTANTS[39] is q_B_celiac_us in component main (litre).
* ALGEBRAIC[32] is u_B_celiac in component main (kPa).
* CONSTANTS[40] is E_B_celiac in component main (kPa_per_L).
* CONSTANTS[41] is z_celiac in component main (meter).
* ALGEBRAIC[33] is v_B_celiac in component main (L_per_s).
* CONSTANTS[42] is R_B_celiac in component main (kPa_s_per_L).
* STATES[11] is q_B_supMes in component main (litre).
* CONSTANTS[43] is q_B_supMes_us in component main (litre).
* ALGEBRAIC[34] is u_B_supMes in component main (kPa).
* CONSTANTS[44] is E_B_supMes in component main (kPa_per_L).
* CONSTANTS[45] is z_supMes in component main (meter).
* ALGEBRAIC[35] is v_B_supMes in component main (L_per_s).
* CONSTANTS[46] is R_B_supMes in component main (kPa_s_per_L).
* ALGEBRAIC[46] is v_B_infMes in component main (L_per_s).
* CONSTANTS[47] is R_B_infMes in component main (kPa_s_per_L).
* STATES[12] is q_B_stomach in component main (litre).
* CONSTANTS[48] is q_B_stomach_us in component main (litre).
* ALGEBRAIC[36] is u_B_stomach in component main (kPa).
* CONSTANTS[49] is E_B_stomach in component main (kPa_per_L).
* CONSTANTS[50] is z_stomach in component main (meter).
* ALGEBRAIC[37] is v_B_stomach1 in component main (L_per_s).
* ALGEBRAIC[49] is v_B_stomach2 in component main (L_per_s).
* CONSTANTS[51] is R_B_stomach1 in component main (kPa_s_per_L).
* CONSTANTS[52] is R_B_stomach2 in component main (kPa_s_per_L).
* STATES[13] is q_B_spleen in component main (litre).
* CONSTANTS[53] is q_B_spleen_us in component main (litre).
* ALGEBRAIC[38] is u_B_spleen in component main (kPa).
* CONSTANTS[54] is E_B_spleen in component main (kPa_per_L).
* CONSTANTS[55] is z_spleen in component main (meter).
* ALGEBRAIC[39] is v_B_spleen1 in component main (L_per_s).
* ALGEBRAIC[50] is v_B_spleen2 in component main (L_per_s).
* CONSTANTS[56] is R_B_spleen1 in component main (kPa_s_per_L).
* CONSTANTS[57] is R_B_spleen2 in component main (kPa_s_per_L).
* STATES[14] is q_B_pancreas in component main (litre).
* CONSTANTS[58] is q_B_pancreas_us in component main (litre).
* ALGEBRAIC[40] is u_B_pancreas in component main (kPa).
* CONSTANTS[59] is E_B_pancreas in component main (kPa_per_L).
* CONSTANTS[60] is z_pancreas in component main (meter).
* ALGEBRAIC[41] is v_B_pancreas1 in component main (L_per_s).
* ALGEBRAIC[42] is v_B_pancreas2 in component main (L_per_s).
* ALGEBRAIC[51] is v_B_pancreas3 in component main (L_per_s).
* CONSTANTS[61] is R_B_pancreas1 in component main (kPa_s_per_L).
* CONSTANTS[62] is R_B_pancreas2 in component main (kPa_s_per_L).
* CONSTANTS[63] is R_B_pancreas3 in component main (kPa_s_per_L).
* STATES[15] is q_B_intestine in component main (litre).
* CONSTANTS[64] is q_B_intestine_us in component main (litre).
* ALGEBRAIC[43] is u_B_intestine in component main (kPa).
* CONSTANTS[65] is E_B_intestine in component main (kPa_per_L).
* CONSTANTS[66] is z_ in component main (meter).
* ALGEBRAIC[44] is v_B_intestine1 in component main (L_per_s).
* ALGEBRAIC[52] is v_B_intestine2 in component main (L_per_s).
* CONSTANTS[67] is R_B_intestine1 in component main (kPa_s_per_L).
* CONSTANTS[68] is R_B_intestine2 in component main (kPa_s_per_L).
* STATES[16] is q_B_colon in component main (litre).
* CONSTANTS[69] is q_B_colon_us in component main (litre).
* ALGEBRAIC[45] is u_B_colon in component main (kPa).
* CONSTANTS[70] is E_B_colon in component main (kPa_per_L).
* CONSTANTS[71] is z_colon in component main (meter).
* ALGEBRAIC[47] is v_B_colon1 in component main (L_per_s).
* ALGEBRAIC[53] is v_B_colon2 in component main (L_per_s).
* CONSTANTS[72] is R_B_colon1 in component main (kPa_s_per_L).
* CONSTANTS[73] is R_B_colon2 in component main (kPa_s_per_L).
* STATES[17] is q_B_portalVein in component main (litre).
* CONSTANTS[74] is q_B_portalVein_us in component main (litre).
* ALGEBRAIC[48] is u_B_portalVein in component main (kPa).
* CONSTANTS[75] is E_B_portalVein in component main (kPa_per_L).
* CONSTANTS[76] is z_portalVein in component main (meter).
* ALGEBRAIC[55] is v_B_portalVein in component main (L_per_s).
* CONSTANTS[77] is R_B_portalVein in component main (kPa_s_per_L).
* STATES[18] is q_B_liver in component main (litre).
* CONSTANTS[78] is q_B_liver_us in component main (litre).
* ALGEBRAIC[54] is u_B_liver in component main (kPa).
* CONSTANTS[79] is E_B_liver in component main (kPa_per_L).
* CONSTANTS[80] is z_liver in component main (meter).
* ALGEBRAIC[56] is v_B_liver1 in component main (L_per_s).
* ALGEBRAIC[57] is v_B_liver2 in component main (L_per_s).
* CONSTANTS[81] is R_B_liver1 in component main (kPa_s_per_L).
* CONSTANTS[82] is R_B_liver2 in component main (kPa_s_per_L).
* STATES[19] is q_B_lGlomerulus in component main (litre).
* STATES[20] is q_B_rGlomerulus in component main (litre).
* STATES[21] is q_B_lPeritubular in component main (litre).
* STATES[22] is q_B_rPeritubular in component main (litre).
* CONSTANTS[83] is q_B_lGlomerulus_us in component main (litre).
* CONSTANTS[84] is q_B_rGlomerulus_us in component main (litre).
* CONSTANTS[85] is q_B_lPeritubular_us in component main (litre).
* CONSTANTS[86] is q_B_rPeritubular_us in component main (litre).
* ALGEBRAIC[58] is u_B_lGlomerulus in component main (kPa).
* ALGEBRAIC[59] is u_B_rGlomerulus in component main (kPa).
* ALGEBRAIC[60] is u_B_lPeritubular in component main (kPa).
* ALGEBRAIC[61] is u_B_rPeritubular in component main (kPa).
* CONSTANTS[87] is E_B_lGlomerulus in component main (kPa_per_L).
* CONSTANTS[88] is E_B_rGlomerulus in component main (kPa_per_L).
* CONSTANTS[89] is E_B_lPeritubular in component main (kPa_per_L).
* CONSTANTS[90] is E_B_rPeritubular in component main (kPa_per_L).
* CONSTANTS[91] is z_kidney in component main (meter).
* ALGEBRAIC[63] is v_B_lRenal1 in component main (L_per_s).
* ALGEBRAIC[64] is v_B_rRenal1 in component main (L_per_s).
* ALGEBRAIC[65] is v_B_lRenal2 in component main (L_per_s).
* ALGEBRAIC[62] is v_B_rRenal2 in component main (L_per_s).
* ALGEBRAIC[66] is v_B_lRenal3 in component main (L_per_s).
* ALGEBRAIC[67] is v_B_rRenal3 in component main (L_per_s).
* CONSTANTS[92] is R_B_lRenal1 in component main (kPa_s_per_L).
* CONSTANTS[93] is R_B_rRenal1 in component main (kPa_s_per_L).
* CONSTANTS[94] is R_B_lRenal2 in component main (kPa_s_per_L).
* CONSTANTS[95] is R_B_rRenal2 in component main (kPa_s_per_L).
* CONSTANTS[96] is R_B_lRenal3 in component main (kPa_s_per_L).
* CONSTANTS[97] is R_B_rRenal3 in component main (kPa_s_per_L).
* STATES[23] is q_B_musTorso in component main (litre).
* STATES[24] is q_B_musTorsoVein in component main (litre).
* CONSTANTS[98] is q_B_musTorso_us in component main (litre).
* CONSTANTS[99] is q_B_musTorsoVein_us in component main (litre).
* ALGEBRAIC[68] is u_B_musTorso in component main (kPa).
* ALGEBRAIC[69] is u_B_musTorsoVein in component main (kPa).
* CONSTANTS[100] is E_B_musTorso in component main (kPa_per_L).
* CONSTANTS[101] is E_B_musTorsoVein in component main (kPa_per_L).
* CONSTANTS[102] is z_musTorso in component main (meter).
* ALGEBRAIC[70] is v_B_musTorso1 in component main (L_per_s).
* ALGEBRAIC[71] is v_B_musTorso2 in component main (L_per_s).
* ALGEBRAIC[72] is v_B_musTorso3 in component main (L_per_s).
* CONSTANTS[103] is R_B_musTorso1 in component main (kPa_s_per_L).
* CONSTANTS[104] is R_B_musTorso2 in component main (kPa_s_per_L).
* CONSTANTS[105] is R_B_musTorso3 in component main (kPa_s_per_L).
* STATES[25] is q_B_lArm in component main (litre).
* STATES[26] is q_B_lArmVein in component main (litre).
* CONSTANTS[106] is q_B_lArm_us in component main (litre).
* CONSTANTS[107] is q_B_lArmVein_us in component main (litre).
* ALGEBRAIC[73] is u_B_lArm in component main (kPa).
* ALGEBRAIC[74] is u_B_lArmVein in component main (kPa).
* CONSTANTS[108] is E_B_lArm in component main (kPa_per_L).
* CONSTANTS[109] is E_B_lArmVein in component main (kPa_per_L).
* CONSTANTS[110] is z_lArm in component main (meter).
* ALGEBRAIC[75] is v_B_lArm1 in component main (L_per_s).
* ALGEBRAIC[76] is v_B_lArm2 in component main (L_per_s).
* ALGEBRAIC[77] is v_B_lArm3 in component main (L_per_s).
* CONSTANTS[111] is R_B_lArm1 in component main (kPa_s_per_L).
* CONSTANTS[112] is R_B_lArm2 in component main (kPa_s_per_L).
* CONSTANTS[113] is R_B_lArm3 in component main (kPa_s_per_L).
* STATES[27] is q_B_rArm in component main (litre).
* STATES[28] is q_B_rArmVein in component main (litre).
* CONSTANTS[114] is q_B_rArm_us in component main (litre).
* CONSTANTS[115] is q_B_rArmVein_us in component main (litre).
* ALGEBRAIC[78] is u_B_rArm in component main (kPa).
* ALGEBRAIC[79] is u_B_rArmVein in component main (kPa).
* CONSTANTS[116] is E_B_rArm in component main (kPa_per_L).
* CONSTANTS[117] is E_B_rArmVein in component main (kPa_per_L).
* CONSTANTS[118] is z_rArm in component main (meter).
* ALGEBRAIC[80] is v_B_rArm1 in component main (L_per_s).
* ALGEBRAIC[81] is v_B_rArm2 in component main (L_per_s).
* ALGEBRAIC[82] is v_B_rArm3 in component main (L_per_s).
* CONSTANTS[119] is R_B_rArm1 in component main (kPa_s_per_L).
* CONSTANTS[120] is R_B_rArm2 in component main (kPa_s_per_L).
* CONSTANTS[121] is R_B_rArm3 in component main (kPa_s_per_L).
* STATES[29] is q_B_lLeg in component main (litre).
* STATES[30] is q_B_lLegVein in component main (litre).
* CONSTANTS[122] is q_B_lLeg_us in component main (litre).
* CONSTANTS[123] is q_B_lLegVein_us in component main (litre).
* ALGEBRAIC[83] is u_B_lLeg in component main (kPa).
* ALGEBRAIC[84] is u_B_lLegVein in component main (kPa).
* CONSTANTS[124] is E_B_lLeg in component main (kPa_per_L).
* CONSTANTS[125] is E_B_lLegVein in component main (kPa_per_L).
* CONSTANTS[126] is z_lLeg in component main (meter).
* ALGEBRAIC[85] is v_B_lLeg1 in component main (L_per_s).
* ALGEBRAIC[86] is v_B_lLeg2 in component main (L_per_s).
* ALGEBRAIC[87] is v_B_lLeg3 in component main (L_per_s).
* CONSTANTS[127] is R_B_lLeg1 in component main (kPa_s_per_L).
* CONSTANTS[128] is R_B_lLeg2 in component main (kPa_s_per_L).
* CONSTANTS[129] is R_B_lLeg3 in component main (kPa_s_per_L).
* STATES[31] is q_B_rLeg in component main (litre).
* STATES[32] is q_B_rLegVein in component main (litre).
* CONSTANTS[130] is q_B_rLeg_us in component main (litre).
* CONSTANTS[131] is q_B_rLegVein_us in component main (litre).
* ALGEBRAIC[88] is u_B_rLeg in component main (kPa).
* ALGEBRAIC[89] is u_B_rLegVein in component main (kPa).
* CONSTANTS[132] is E_B_rLeg in component main (kPa_per_L).
* CONSTANTS[133] is E_B_rLegVein in component main (kPa_per_L).
* CONSTANTS[134] is z_rLeg in component main (meter).
* ALGEBRAIC[90] is v_B_rLeg1 in component main (L_per_s).
* ALGEBRAIC[91] is v_B_rLeg2 in component main (L_per_s).
* ALGEBRAIC[92] is v_B_rLeg3 in component main (L_per_s).
* CONSTANTS[135] is R_B_rLeg1 in component main (kPa_s_per_L).
* CONSTANTS[136] is R_B_rLeg2 in component main (kPa_s_per_L).
* CONSTANTS[137] is R_B_rLeg3 in component main (kPa_s_per_L).
* ALGEBRAIC[3] is q_B_heartTot in component main (litre).
* ALGEBRAIC[5] is q_B_lungTot in component main (litre).
* ALGEBRAIC[7] is q_B_renalTot in component main (litre).
* ALGEBRAIC[9] is q_B_gutTot in component main (litre).
* ALGEBRAIC[11] is q_B_brainTot in component main (litre).
* ALGEBRAIC[13] is q_B_limbTot in component main (litre).
* ALGEBRAIC[15] is q_Blood in component main (litre).
* RATES[0] is d/dt q_B_LV in component main (litre).
* RATES[1] is d/dt q_B_RV in component main (litre).
* RATES[2] is d/dt q_B_LA in component main (litre).
* RATES[3] is d/dt q_B_RA in component main (litre).
* RATES[4] is d/dt q_B_PA in component main (litre).
* RATES[5] is d/dt q_B_lung in component main (litre).
* RATES[6] is d/dt q_B_pulmVein in component main (litre).
* RATES[7] is d/dt q_B_brain in component main (litre).
* RATES[8] is d/dt q_B_brainVein in component main (litre).
* RATES[9] is d/dt q_B_AA in component main (litre).
* RATES[10] is d/dt q_B_celiac in component main (litre).
* RATES[11] is d/dt q_B_supMes in component main (litre).
* RATES[12] is d/dt q_B_stomach in component main (litre).
* RATES[13] is d/dt q_B_spleen in component main (litre).
* RATES[14] is d/dt q_B_pancreas in component main (litre).
* RATES[15] is d/dt q_B_intestine in component main (litre).
* RATES[16] is d/dt q_B_colon in component main (litre).
* RATES[17] is d/dt q_B_portalVein in component main (litre).
* RATES[18] is d/dt q_B_liver in component main (litre).
* RATES[19] is d/dt q_B_lGlomerulus in component main (litre).
* RATES[20] is d/dt q_B_rGlomerulus in component main (litre).
* RATES[21] is d/dt q_B_lPeritubular in component main (litre).
* RATES[22] is d/dt q_B_rPeritubular in component main (litre).
* RATES[23] is d/dt q_B_musTorso in component main (litre).
* RATES[24] is d/dt q_B_musTorsoVein in component main (litre).
* RATES[25] is d/dt q_B_lArm in component main (litre).
* RATES[26] is d/dt q_B_lArmVein in component main (litre).
* RATES[27] is d/dt q_B_rArm in component main (litre).
* RATES[28] is d/dt q_B_rArmVein in component main (litre).
* RATES[29] is d/dt q_B_lLeg in component main (litre).
* RATES[30] is d/dt q_B_lLegVein in component main (litre).
* RATES[31] is d/dt q_B_rLeg in component main (litre).
* RATES[32] is d/dt q_B_rLegVein in component main (litre).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1e3;
CONSTANTS[1] = 9.81;
CONSTANTS[2] = 90;
CONSTANTS[3] = 1;
CONSTANTS[4] = 0.1;
CONSTANTS[5] = 0.1;
CONSTANTS[6] = 320;
CONSTANTS[7] = 6.667;
STATES[0] = 1;
CONSTANTS[8] = 0.005;
CONSTANTS[9] = 70;
CONSTANTS[10] = 2;
STATES[1] = 0.2;
CONSTANTS[11] = 0.01;
CONSTANTS[12] = 10;
STATES[2] = 0.02;
CONSTANTS[13] = 0.004;
CONSTANTS[14] = 10;
STATES[3] = 0.02;
CONSTANTS[15] = 0.004;
CONSTANTS[16] = 1;
CONSTANTS[17] = 1;
CONSTANTS[18] = 1;
CONSTANTS[19] = 1;
STATES[4] = 0.3;
CONSTANTS[20] = 0.2;
CONSTANTS[21] = 200;
STATES[5] = 0.2;
CONSTANTS[22] = 0.1;
CONSTANTS[23] = 100;
CONSTANTS[24] = 0.20;
CONSTANTS[25] = 3;
CONSTANTS[26] = 10;
CONSTANTS[27] = 20;
STATES[6] = 0.02;
CONSTANTS[28] = 0.012;
CONSTANTS[29] = 1;
STATES[7] = 0.1;
CONSTANTS[30] = 0.06;
STATES[8] = 0.1;
CONSTANTS[31] = 0.06;
CONSTANTS[32] = 20;
CONSTANTS[33] = 20;
CONSTANTS[34] = 100;
CONSTANTS[35] = 600;
CONSTANTS[36] = 1;
STATES[9] = 0.2;
CONSTANTS[37] = 0.12;
CONSTANTS[38] = 80;
STATES[10] = 0.1;
CONSTANTS[39] = 0.06;
CONSTANTS[40] = 400;
CONSTANTS[41] = 0.20;
CONSTANTS[42] = 100;
STATES[11] = 0.1;
CONSTANTS[43] = 0.06;
CONSTANTS[44] = 400;
CONSTANTS[45] = 0.20;
CONSTANTS[46] = 100;
CONSTANTS[47] = 100;
STATES[12] = 0.1;
CONSTANTS[48] = 0.06;
CONSTANTS[49] = 100;
CONSTANTS[50] = 0.20;
CONSTANTS[51] = 100;
CONSTANTS[52] = 2400;
STATES[13] = 0.1;
CONSTANTS[53] = 0.06;
CONSTANTS[54] = 100;
CONSTANTS[55] = 0.20;
CONSTANTS[56] = 100;
CONSTANTS[57] = 3000;
STATES[14] = 0.1;
CONSTANTS[58] = 0.06;
CONSTANTS[59] = 100;
CONSTANTS[60] = 0.20;
CONSTANTS[61] = 100;
CONSTANTS[62] = 100;
CONSTANTS[63] = 4000;
STATES[15] = 0.1;
CONSTANTS[64] = 0.06;
CONSTANTS[65] = 100;
CONSTANTS[66] = 0.20;
CONSTANTS[67] = 100;
CONSTANTS[68] = 1200;
STATES[16] = 0.1;
CONSTANTS[69] = 0.06;
CONSTANTS[70] = 100;
CONSTANTS[71] = 0.20;
CONSTANTS[72] = 100;
CONSTANTS[73] = 2000;
STATES[17] = 0.1;
CONSTANTS[74] = 0.06;
CONSTANTS[75] = 50;
CONSTANTS[76] = 0.20;
CONSTANTS[77] = 100;
STATES[18] = 0.2;
CONSTANTS[78] = 0.06;
CONSTANTS[79] = 10;
CONSTANTS[80] = -0.50;
CONSTANTS[81] = 800;
CONSTANTS[82] = 1;
STATES[19] = 0.1;
STATES[20] = 0.1;
STATES[21] = 0.1;
STATES[22] = 0.1;
CONSTANTS[83] = 0.06;
CONSTANTS[84] = 0.06;
CONSTANTS[85] = 0.06;
CONSTANTS[86] = 0.06;
CONSTANTS[87] = 150;
CONSTANTS[88] = 150;
CONSTANTS[89] = 50;
CONSTANTS[90] = 50;
CONSTANTS[91] = -0.25;
CONSTANTS[92] = 100;
CONSTANTS[93] = 100;
CONSTANTS[94] = 2000;
CONSTANTS[95] = 2000;
CONSTANTS[96] = 1;
CONSTANTS[97] = 1;
STATES[23] = 0.2;
STATES[24] = 0.2;
CONSTANTS[98] = 0.12;
CONSTANTS[99] = 0.12;
CONSTANTS[100] = 60;
CONSTANTS[101] = 10;
CONSTANTS[102] = -0.50;
CONSTANTS[103] = 100;
CONSTANTS[104] = 1000;
CONSTANTS[105] = 1;
STATES[25] = 0.1;
STATES[26] = 0.1;
CONSTANTS[106] = 0.06;
CONSTANTS[107] = 0.06;
CONSTANTS[108] = 150;
CONSTANTS[109] = 60;
CONSTANTS[110] = -0.50;
CONSTANTS[111] = 100;
CONSTANTS[112] = 6000;
CONSTANTS[113] = 1;
STATES[27] = 0.1;
STATES[28] = 0.1;
CONSTANTS[114] = 0.06;
CONSTANTS[115] = 0.06;
CONSTANTS[116] = 150;
CONSTANTS[117] = 60;
CONSTANTS[118] = -0.50;
CONSTANTS[119] = 100;
CONSTANTS[120] = 6000;
CONSTANTS[121] = 1;
STATES[29] = 0.1;
STATES[30] = 0.1;
CONSTANTS[122] = 0.06;
CONSTANTS[123] = 0.06;
CONSTANTS[124] = 150;
CONSTANTS[125] = 60;
CONSTANTS[126] = -0.50;
CONSTANTS[127] = 100;
CONSTANTS[128] = 4000;
CONSTANTS[129] = 1;
STATES[31] = 0.1;
STATES[32] = 0.1;
CONSTANTS[130] = 0.06;
CONSTANTS[131] = 0.06;
CONSTANTS[132] = 150;
CONSTANTS[133] = 60;
CONSTANTS[134] = -0.50;
CONSTANTS[135] = 100;
CONSTANTS[136] = 4000;
CONSTANTS[137] = 1;
CONSTANTS[138] = ( CONSTANTS[2]* 3.14159265358979)/180.000;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[2] = VOI/CONSTANTS[3] - floor(VOI/CONSTANTS[3]);
ALGEBRAIC[8] = (ALGEBRAIC[2]<0.400000 ? (( (CONSTANTS[9] - CONSTANTS[10])*ALGEBRAIC[2])/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(ALGEBRAIC[2]/CONSTANTS[5], 2.00000)))+CONSTANTS[10] : ALGEBRAIC[2]<0.600000 ? (( (CONSTANTS[9] - CONSTANTS[10])*0.400000)/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(0.400000/CONSTANTS[5], 2.00000)))+CONSTANTS[10] : CONSTANTS[10]);
ALGEBRAIC[10] = ALGEBRAIC[8]*(STATES[1] - CONSTANTS[11]);
ALGEBRAIC[18] = CONSTANTS[21]*(STATES[4] - CONSTANTS[20]);
ALGEBRAIC[20] = (ALGEBRAIC[10]>ALGEBRAIC[18] ? (ALGEBRAIC[10] - ALGEBRAIC[18])/CONSTANTS[17] : 0.00000);
ALGEBRAIC[14] = CONSTANTS[14]*(STATES[3] - CONSTANTS[15]);
ALGEBRAIC[16] = (ALGEBRAIC[14]>ALGEBRAIC[10] ? (ALGEBRAIC[14] - ALGEBRAIC[10])/CONSTANTS[18] : 0.00000);
RATES[1] = ALGEBRAIC[16] - ALGEBRAIC[20];
ALGEBRAIC[22] = CONSTANTS[23]*(STATES[5] - CONSTANTS[22]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[24]*cos(CONSTANTS[138]);
ALGEBRAIC[25] = (ALGEBRAIC[18] - ALGEBRAIC[22])/CONSTANTS[25];
RATES[4] = ALGEBRAIC[20] - ALGEBRAIC[25];
ALGEBRAIC[19] = CONSTANTS[32]*(STATES[7] - CONSTANTS[30]);
ALGEBRAIC[21] = CONSTANTS[33]*(STATES[8] - CONSTANTS[31]);
ALGEBRAIC[23] = (ALGEBRAIC[19] - ALGEBRAIC[21])/CONSTANTS[35];
ALGEBRAIC[24] = (ALGEBRAIC[21] - ALGEBRAIC[14])/CONSTANTS[36];
RATES[8] = ALGEBRAIC[23] - ALGEBRAIC[24];
ALGEBRAIC[4] = (ALGEBRAIC[2]<0.400000 ? (( (CONSTANTS[6] - CONSTANTS[7])*ALGEBRAIC[2])/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(ALGEBRAIC[2]/CONSTANTS[5], 2.00000)))+CONSTANTS[7] : ALGEBRAIC[2]<0.600000 ? (( (CONSTANTS[6] - CONSTANTS[7])*0.400000)/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(0.400000/CONSTANTS[5], 2.00000)))+CONSTANTS[7] : CONSTANTS[7]);
ALGEBRAIC[6] = ALGEBRAIC[4]*(STATES[0] - CONSTANTS[8]);
ALGEBRAIC[26] = CONSTANTS[38]*(STATES[9] - CONSTANTS[37]);
ALGEBRAIC[29] = (ALGEBRAIC[6]>ALGEBRAIC[26] ? (ALGEBRAIC[6] - ALGEBRAIC[26])/CONSTANTS[16] : 0.00000);
ALGEBRAIC[12] = CONSTANTS[12]*(STATES[2] - CONSTANTS[13]);
ALGEBRAIC[17] = (ALGEBRAIC[12]>ALGEBRAIC[6] ? (ALGEBRAIC[12] - ALGEBRAIC[6])/CONSTANTS[19] : 0.00000);
RATES[0] = ALGEBRAIC[17] - ALGEBRAIC[29];
ALGEBRAIC[27] = CONSTANTS[27]*(STATES[6] - CONSTANTS[28]);
ALGEBRAIC[28] = (ALGEBRAIC[22] - ALGEBRAIC[27])/CONSTANTS[26];
RATES[5] = ALGEBRAIC[25] - ALGEBRAIC[28];
ALGEBRAIC[30] = (ALGEBRAIC[27] - ALGEBRAIC[12])/CONSTANTS[29];
RATES[2] = ALGEBRAIC[30] - ALGEBRAIC[17];
RATES[6] = ALGEBRAIC[28] - ALGEBRAIC[30];
ALGEBRAIC[31] = (ALGEBRAIC[26] - ALGEBRAIC[19])/CONSTANTS[34];
RATES[7] = ALGEBRAIC[31] - ALGEBRAIC[23];
ALGEBRAIC[34] = CONSTANTS[44]*(STATES[11] - CONSTANTS[43]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[45]*cos(CONSTANTS[138]);
ALGEBRAIC[35] = (ALGEBRAIC[26] - ALGEBRAIC[34])/CONSTANTS[46];
ALGEBRAIC[40] = CONSTANTS[59]*(STATES[14] - CONSTANTS[58]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[60]*cos(CONSTANTS[138]);
ALGEBRAIC[42] = (ALGEBRAIC[34] - ALGEBRAIC[40])/CONSTANTS[62];
ALGEBRAIC[43] = CONSTANTS[65]*(STATES[15] - CONSTANTS[64]);
ALGEBRAIC[44] = (ALGEBRAIC[34] - ALGEBRAIC[43])/CONSTANTS[67];
ALGEBRAIC[45] = CONSTANTS[70]*(STATES[16] - CONSTANTS[69]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[71]*cos(CONSTANTS[138]);
ALGEBRAIC[47] = (ALGEBRAIC[34] - ALGEBRAIC[45])/CONSTANTS[72];
RATES[11] = ((ALGEBRAIC[35] - ALGEBRAIC[42]) - ALGEBRAIC[44]) - ALGEBRAIC[47];
ALGEBRAIC[32] = CONSTANTS[40]*(STATES[10] - CONSTANTS[39]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[41]*cos(CONSTANTS[138]);
ALGEBRAIC[36] = CONSTANTS[49]*(STATES[12] - CONSTANTS[48]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[50]*cos(CONSTANTS[138]);
ALGEBRAIC[37] = (ALGEBRAIC[32] - ALGEBRAIC[36])/CONSTANTS[51];
ALGEBRAIC[48] = CONSTANTS[75]*(STATES[17] - CONSTANTS[74]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[76]*cos(CONSTANTS[138]);
ALGEBRAIC[49] = (ALGEBRAIC[36] - ALGEBRAIC[48])/CONSTANTS[52];
RATES[12] = ALGEBRAIC[37] - ALGEBRAIC[49];
ALGEBRAIC[38] = CONSTANTS[70]*(STATES[13] - CONSTANTS[53]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[55]*cos(CONSTANTS[138]);
ALGEBRAIC[39] = (ALGEBRAIC[32] - ALGEBRAIC[38])/CONSTANTS[56];
ALGEBRAIC[50] = (ALGEBRAIC[38] - ALGEBRAIC[48])/CONSTANTS[57];
RATES[13] = ALGEBRAIC[39] - ALGEBRAIC[50];
ALGEBRAIC[41] = (ALGEBRAIC[32] - ALGEBRAIC[40])/CONSTANTS[61];
ALGEBRAIC[51] = (ALGEBRAIC[40] - ALGEBRAIC[48])/CONSTANTS[63];
RATES[14] = (ALGEBRAIC[41]+ALGEBRAIC[42]) - ALGEBRAIC[51];
ALGEBRAIC[52] = (ALGEBRAIC[43] - ALGEBRAIC[48])/CONSTANTS[68];
RATES[15] = ALGEBRAIC[44] - ALGEBRAIC[52];
ALGEBRAIC[46] = (ALGEBRAIC[26] - ALGEBRAIC[45])/CONSTANTS[47];
ALGEBRAIC[53] = (ALGEBRAIC[45] - ALGEBRAIC[48])/CONSTANTS[73];
RATES[16] = (ALGEBRAIC[47]+ALGEBRAIC[46]) - ALGEBRAIC[53];
ALGEBRAIC[54] = CONSTANTS[79]*(STATES[18] - CONSTANTS[78]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[80]*cos(CONSTANTS[138]);
ALGEBRAIC[55] = (ALGEBRAIC[48] - ALGEBRAIC[54])/CONSTANTS[77];
RATES[17] = (ALGEBRAIC[49]+ALGEBRAIC[50]+ALGEBRAIC[51]+ALGEBRAIC[52]+ALGEBRAIC[53]) - ALGEBRAIC[55];
ALGEBRAIC[33] = (ALGEBRAIC[26] - ALGEBRAIC[32])/CONSTANTS[42];
ALGEBRAIC[56] = (ALGEBRAIC[32] - ALGEBRAIC[54])/CONSTANTS[81];
RATES[10] = (((ALGEBRAIC[33] - ALGEBRAIC[56]) - ALGEBRAIC[37]) - ALGEBRAIC[39]) - ALGEBRAIC[41];
ALGEBRAIC[57] = (ALGEBRAIC[54] - ALGEBRAIC[14])/CONSTANTS[82];
RATES[18] = (ALGEBRAIC[56]+ALGEBRAIC[55]) - ALGEBRAIC[57];
ALGEBRAIC[58] = CONSTANTS[87]*(STATES[19] - CONSTANTS[83]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[63] = (ALGEBRAIC[26] - ALGEBRAIC[58])/CONSTANTS[92];
ALGEBRAIC[60] = CONSTANTS[89]*(STATES[21] - CONSTANTS[85]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[65] = (ALGEBRAIC[58] - ALGEBRAIC[60])/CONSTANTS[94];
RATES[19] = ALGEBRAIC[63] - ALGEBRAIC[65];
RATES[20] = ALGEBRAIC[63] - ALGEBRAIC[65];
ALGEBRAIC[66] = (ALGEBRAIC[60] - ALGEBRAIC[14])/CONSTANTS[96];
RATES[21] = ALGEBRAIC[65] - ALGEBRAIC[66];
RATES[22] = ALGEBRAIC[65] - ALGEBRAIC[66];
ALGEBRAIC[68] = CONSTANTS[100]*(STATES[23] - CONSTANTS[98]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[102]*cos(CONSTANTS[138]);
ALGEBRAIC[70] = (ALGEBRAIC[26] - ALGEBRAIC[68])/CONSTANTS[103];
ALGEBRAIC[69] = CONSTANTS[101]*(STATES[24] - CONSTANTS[99]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[102]*cos(CONSTANTS[138]);
ALGEBRAIC[71] = (ALGEBRAIC[68] - ALGEBRAIC[69])/CONSTANTS[104];
RATES[23] = ALGEBRAIC[70] - ALGEBRAIC[71];
ALGEBRAIC[72] = (ALGEBRAIC[69] - ALGEBRAIC[14])/CONSTANTS[105];
RATES[24] = ALGEBRAIC[71] - ALGEBRAIC[72];
ALGEBRAIC[73] = CONSTANTS[108]*(STATES[25] - CONSTANTS[106]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[110]*cos(CONSTANTS[138]);
ALGEBRAIC[75] = (ALGEBRAIC[26] - ALGEBRAIC[73])/CONSTANTS[111];
ALGEBRAIC[74] = CONSTANTS[109]*(STATES[26] - CONSTANTS[107]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[110]*cos(CONSTANTS[138]);
ALGEBRAIC[76] = (ALGEBRAIC[73] - ALGEBRAIC[74])/CONSTANTS[112];
RATES[25] = ALGEBRAIC[75] - ALGEBRAIC[76];
ALGEBRAIC[77] = (ALGEBRAIC[74] - ALGEBRAIC[14])/CONSTANTS[113];
RATES[26] = ALGEBRAIC[76] - ALGEBRAIC[77];
ALGEBRAIC[78] = CONSTANTS[116]*(STATES[27] - CONSTANTS[114]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[118]*cos(CONSTANTS[138]);
ALGEBRAIC[80] = (ALGEBRAIC[26] - ALGEBRAIC[78])/CONSTANTS[119];
ALGEBRAIC[79] = CONSTANTS[117]*(STATES[28] - CONSTANTS[114]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[118]*cos(CONSTANTS[138]);
ALGEBRAIC[81] = (ALGEBRAIC[78] - ALGEBRAIC[79])/CONSTANTS[120];
RATES[27] = ALGEBRAIC[80] - ALGEBRAIC[81];
ALGEBRAIC[82] = (ALGEBRAIC[79] - ALGEBRAIC[14])/CONSTANTS[121];
RATES[28] = ALGEBRAIC[81] - ALGEBRAIC[82];
ALGEBRAIC[83] = CONSTANTS[124]*(STATES[29] - CONSTANTS[122]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[126]*cos(CONSTANTS[138]);
ALGEBRAIC[85] = (ALGEBRAIC[26] - ALGEBRAIC[83])/CONSTANTS[127];
ALGEBRAIC[84] = CONSTANTS[125]*(STATES[30] - CONSTANTS[122]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[126]*cos(CONSTANTS[138]);
ALGEBRAIC[86] = (ALGEBRAIC[83] - ALGEBRAIC[84])/CONSTANTS[128];
RATES[29] = ALGEBRAIC[85] - ALGEBRAIC[86];
ALGEBRAIC[87] = (ALGEBRAIC[84] - ALGEBRAIC[14])/CONSTANTS[129];
RATES[30] = ALGEBRAIC[86] - ALGEBRAIC[87];
ALGEBRAIC[59] = CONSTANTS[88]*(STATES[20] - CONSTANTS[84]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[64] = (ALGEBRAIC[26] - ALGEBRAIC[59])/CONSTANTS[93];
ALGEBRAIC[88] = CONSTANTS[132]*(STATES[31] - CONSTANTS[130]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[134]*cos(CONSTANTS[138]);
ALGEBRAIC[90] = (ALGEBRAIC[26] - ALGEBRAIC[88])/CONSTANTS[135];
RATES[9] = ((((((((((ALGEBRAIC[29] - ALGEBRAIC[31]) - ALGEBRAIC[33]) - ALGEBRAIC[35]) - ALGEBRAIC[46]) - ALGEBRAIC[63]) - ALGEBRAIC[64]) - ALGEBRAIC[70]) - ALGEBRAIC[75]) - ALGEBRAIC[80]) - ALGEBRAIC[85]) - ALGEBRAIC[90];
ALGEBRAIC[89] = CONSTANTS[133]*(STATES[32] - CONSTANTS[130]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[134]*cos(CONSTANTS[138]);
ALGEBRAIC[91] = (ALGEBRAIC[88] - ALGEBRAIC[89])/CONSTANTS[136];
RATES[31] = ALGEBRAIC[90] - ALGEBRAIC[91];
ALGEBRAIC[61] = CONSTANTS[90]*(STATES[22] - CONSTANTS[86]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[67] = (ALGEBRAIC[61] - ALGEBRAIC[14])/CONSTANTS[97];
ALGEBRAIC[92] = (ALGEBRAIC[89] - ALGEBRAIC[14])/CONSTANTS[137];
RATES[3] = (ALGEBRAIC[24]+ALGEBRAIC[57]+ALGEBRAIC[66]+ALGEBRAIC[67]+ALGEBRAIC[72]+ALGEBRAIC[77]+ALGEBRAIC[82]+ALGEBRAIC[87]+ALGEBRAIC[92]) - ALGEBRAIC[16];
RATES[32] = ALGEBRAIC[91] - ALGEBRAIC[92];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[2] = VOI/CONSTANTS[3] - floor(VOI/CONSTANTS[3]);
ALGEBRAIC[8] = (ALGEBRAIC[2]<0.400000 ? (( (CONSTANTS[9] - CONSTANTS[10])*ALGEBRAIC[2])/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(ALGEBRAIC[2]/CONSTANTS[5], 2.00000)))+CONSTANTS[10] : ALGEBRAIC[2]<0.600000 ? (( (CONSTANTS[9] - CONSTANTS[10])*0.400000)/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(0.400000/CONSTANTS[5], 2.00000)))+CONSTANTS[10] : CONSTANTS[10]);
ALGEBRAIC[10] = ALGEBRAIC[8]*(STATES[1] - CONSTANTS[11]);
ALGEBRAIC[18] = CONSTANTS[21]*(STATES[4] - CONSTANTS[20]);
ALGEBRAIC[20] = (ALGEBRAIC[10]>ALGEBRAIC[18] ? (ALGEBRAIC[10] - ALGEBRAIC[18])/CONSTANTS[17] : 0.00000);
ALGEBRAIC[14] = CONSTANTS[14]*(STATES[3] - CONSTANTS[15]);
ALGEBRAIC[16] = (ALGEBRAIC[14]>ALGEBRAIC[10] ? (ALGEBRAIC[14] - ALGEBRAIC[10])/CONSTANTS[18] : 0.00000);
ALGEBRAIC[22] = CONSTANTS[23]*(STATES[5] - CONSTANTS[22]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[24]*cos(CONSTANTS[138]);
ALGEBRAIC[25] = (ALGEBRAIC[18] - ALGEBRAIC[22])/CONSTANTS[25];
ALGEBRAIC[19] = CONSTANTS[32]*(STATES[7] - CONSTANTS[30]);
ALGEBRAIC[21] = CONSTANTS[33]*(STATES[8] - CONSTANTS[31]);
ALGEBRAIC[23] = (ALGEBRAIC[19] - ALGEBRAIC[21])/CONSTANTS[35];
ALGEBRAIC[24] = (ALGEBRAIC[21] - ALGEBRAIC[14])/CONSTANTS[36];
ALGEBRAIC[4] = (ALGEBRAIC[2]<0.400000 ? (( (CONSTANTS[6] - CONSTANTS[7])*ALGEBRAIC[2])/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(ALGEBRAIC[2]/CONSTANTS[5], 2.00000)))+CONSTANTS[7] : ALGEBRAIC[2]<0.600000 ? (( (CONSTANTS[6] - CONSTANTS[7])*0.400000)/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(0.400000/CONSTANTS[5], 2.00000)))+CONSTANTS[7] : CONSTANTS[7]);
ALGEBRAIC[6] = ALGEBRAIC[4]*(STATES[0] - CONSTANTS[8]);
ALGEBRAIC[26] = CONSTANTS[38]*(STATES[9] - CONSTANTS[37]);
ALGEBRAIC[29] = (ALGEBRAIC[6]>ALGEBRAIC[26] ? (ALGEBRAIC[6] - ALGEBRAIC[26])/CONSTANTS[16] : 0.00000);
ALGEBRAIC[12] = CONSTANTS[12]*(STATES[2] - CONSTANTS[13]);
ALGEBRAIC[17] = (ALGEBRAIC[12]>ALGEBRAIC[6] ? (ALGEBRAIC[12] - ALGEBRAIC[6])/CONSTANTS[19] : 0.00000);
ALGEBRAIC[27] = CONSTANTS[27]*(STATES[6] - CONSTANTS[28]);
ALGEBRAIC[28] = (ALGEBRAIC[22] - ALGEBRAIC[27])/CONSTANTS[26];
ALGEBRAIC[30] = (ALGEBRAIC[27] - ALGEBRAIC[12])/CONSTANTS[29];
ALGEBRAIC[31] = (ALGEBRAIC[26] - ALGEBRAIC[19])/CONSTANTS[34];
ALGEBRAIC[34] = CONSTANTS[44]*(STATES[11] - CONSTANTS[43]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[45]*cos(CONSTANTS[138]);
ALGEBRAIC[35] = (ALGEBRAIC[26] - ALGEBRAIC[34])/CONSTANTS[46];
ALGEBRAIC[40] = CONSTANTS[59]*(STATES[14] - CONSTANTS[58]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[60]*cos(CONSTANTS[138]);
ALGEBRAIC[42] = (ALGEBRAIC[34] - ALGEBRAIC[40])/CONSTANTS[62];
ALGEBRAIC[43] = CONSTANTS[65]*(STATES[15] - CONSTANTS[64]);
ALGEBRAIC[44] = (ALGEBRAIC[34] - ALGEBRAIC[43])/CONSTANTS[67];
ALGEBRAIC[45] = CONSTANTS[70]*(STATES[16] - CONSTANTS[69]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[71]*cos(CONSTANTS[138]);
ALGEBRAIC[47] = (ALGEBRAIC[34] - ALGEBRAIC[45])/CONSTANTS[72];
ALGEBRAIC[32] = CONSTANTS[40]*(STATES[10] - CONSTANTS[39]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[41]*cos(CONSTANTS[138]);
ALGEBRAIC[36] = CONSTANTS[49]*(STATES[12] - CONSTANTS[48]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[50]*cos(CONSTANTS[138]);
ALGEBRAIC[37] = (ALGEBRAIC[32] - ALGEBRAIC[36])/CONSTANTS[51];
ALGEBRAIC[48] = CONSTANTS[75]*(STATES[17] - CONSTANTS[74]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[76]*cos(CONSTANTS[138]);
ALGEBRAIC[49] = (ALGEBRAIC[36] - ALGEBRAIC[48])/CONSTANTS[52];
ALGEBRAIC[38] = CONSTANTS[70]*(STATES[13] - CONSTANTS[53]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[55]*cos(CONSTANTS[138]);
ALGEBRAIC[39] = (ALGEBRAIC[32] - ALGEBRAIC[38])/CONSTANTS[56];
ALGEBRAIC[50] = (ALGEBRAIC[38] - ALGEBRAIC[48])/CONSTANTS[57];
ALGEBRAIC[41] = (ALGEBRAIC[32] - ALGEBRAIC[40])/CONSTANTS[61];
ALGEBRAIC[51] = (ALGEBRAIC[40] - ALGEBRAIC[48])/CONSTANTS[63];
ALGEBRAIC[52] = (ALGEBRAIC[43] - ALGEBRAIC[48])/CONSTANTS[68];
ALGEBRAIC[46] = (ALGEBRAIC[26] - ALGEBRAIC[45])/CONSTANTS[47];
ALGEBRAIC[53] = (ALGEBRAIC[45] - ALGEBRAIC[48])/CONSTANTS[73];
ALGEBRAIC[54] = CONSTANTS[79]*(STATES[18] - CONSTANTS[78]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[80]*cos(CONSTANTS[138]);
ALGEBRAIC[55] = (ALGEBRAIC[48] - ALGEBRAIC[54])/CONSTANTS[77];
ALGEBRAIC[33] = (ALGEBRAIC[26] - ALGEBRAIC[32])/CONSTANTS[42];
ALGEBRAIC[56] = (ALGEBRAIC[32] - ALGEBRAIC[54])/CONSTANTS[81];
ALGEBRAIC[57] = (ALGEBRAIC[54] - ALGEBRAIC[14])/CONSTANTS[82];
ALGEBRAIC[58] = CONSTANTS[87]*(STATES[19] - CONSTANTS[83]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[63] = (ALGEBRAIC[26] - ALGEBRAIC[58])/CONSTANTS[92];
ALGEBRAIC[60] = CONSTANTS[89]*(STATES[21] - CONSTANTS[85]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[65] = (ALGEBRAIC[58] - ALGEBRAIC[60])/CONSTANTS[94];
ALGEBRAIC[66] = (ALGEBRAIC[60] - ALGEBRAIC[14])/CONSTANTS[96];
ALGEBRAIC[68] = CONSTANTS[100]*(STATES[23] - CONSTANTS[98]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[102]*cos(CONSTANTS[138]);
ALGEBRAIC[70] = (ALGEBRAIC[26] - ALGEBRAIC[68])/CONSTANTS[103];
ALGEBRAIC[69] = CONSTANTS[101]*(STATES[24] - CONSTANTS[99]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[102]*cos(CONSTANTS[138]);
ALGEBRAIC[71] = (ALGEBRAIC[68] - ALGEBRAIC[69])/CONSTANTS[104];
ALGEBRAIC[72] = (ALGEBRAIC[69] - ALGEBRAIC[14])/CONSTANTS[105];
ALGEBRAIC[73] = CONSTANTS[108]*(STATES[25] - CONSTANTS[106]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[110]*cos(CONSTANTS[138]);
ALGEBRAIC[75] = (ALGEBRAIC[26] - ALGEBRAIC[73])/CONSTANTS[111];
ALGEBRAIC[74] = CONSTANTS[109]*(STATES[26] - CONSTANTS[107]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[110]*cos(CONSTANTS[138]);
ALGEBRAIC[76] = (ALGEBRAIC[73] - ALGEBRAIC[74])/CONSTANTS[112];
ALGEBRAIC[77] = (ALGEBRAIC[74] - ALGEBRAIC[14])/CONSTANTS[113];
ALGEBRAIC[78] = CONSTANTS[116]*(STATES[27] - CONSTANTS[114]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[118]*cos(CONSTANTS[138]);
ALGEBRAIC[80] = (ALGEBRAIC[26] - ALGEBRAIC[78])/CONSTANTS[119];
ALGEBRAIC[79] = CONSTANTS[117]*(STATES[28] - CONSTANTS[114]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[118]*cos(CONSTANTS[138]);
ALGEBRAIC[81] = (ALGEBRAIC[78] - ALGEBRAIC[79])/CONSTANTS[120];
ALGEBRAIC[82] = (ALGEBRAIC[79] - ALGEBRAIC[14])/CONSTANTS[121];
ALGEBRAIC[83] = CONSTANTS[124]*(STATES[29] - CONSTANTS[122]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[126]*cos(CONSTANTS[138]);
ALGEBRAIC[85] = (ALGEBRAIC[26] - ALGEBRAIC[83])/CONSTANTS[127];
ALGEBRAIC[84] = CONSTANTS[125]*(STATES[30] - CONSTANTS[122]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[126]*cos(CONSTANTS[138]);
ALGEBRAIC[86] = (ALGEBRAIC[83] - ALGEBRAIC[84])/CONSTANTS[128];
ALGEBRAIC[87] = (ALGEBRAIC[84] - ALGEBRAIC[14])/CONSTANTS[129];
ALGEBRAIC[59] = CONSTANTS[88]*(STATES[20] - CONSTANTS[84]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[64] = (ALGEBRAIC[26] - ALGEBRAIC[59])/CONSTANTS[93];
ALGEBRAIC[88] = CONSTANTS[132]*(STATES[31] - CONSTANTS[130]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[134]*cos(CONSTANTS[138]);
ALGEBRAIC[90] = (ALGEBRAIC[26] - ALGEBRAIC[88])/CONSTANTS[135];
ALGEBRAIC[89] = CONSTANTS[133]*(STATES[32] - CONSTANTS[130]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[134]*cos(CONSTANTS[138]);
ALGEBRAIC[91] = (ALGEBRAIC[88] - ALGEBRAIC[89])/CONSTANTS[136];
ALGEBRAIC[61] = CONSTANTS[90]*(STATES[22] - CONSTANTS[86]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[67] = (ALGEBRAIC[61] - ALGEBRAIC[14])/CONSTANTS[97];
ALGEBRAIC[92] = (ALGEBRAIC[89] - ALGEBRAIC[14])/CONSTANTS[137];
ALGEBRAIC[0] = (VOI>100.000&&VOI<101.000 ? 1.00000 : 0.00000);
ALGEBRAIC[1] = (VOI/CONSTANTS[3] - CONSTANTS[4]) - floor((VOI - CONSTANTS[4]*CONSTANTS[3])/CONSTANTS[3]);
ALGEBRAIC[3] = STATES[0]+STATES[1]+STATES[2]+STATES[3]+STATES[9]+STATES[4];
ALGEBRAIC[5] = STATES[5]+STATES[6];
ALGEBRAIC[7] = STATES[19]+STATES[20]+STATES[21]+STATES[22];
ALGEBRAIC[9] = STATES[10]+STATES[11]+STATES[12]+STATES[13]+STATES[14]+STATES[15]+STATES[16]+STATES[17]+STATES[18]+STATES[23]+STATES[24];
ALGEBRAIC[11] = STATES[7]+STATES[8];
ALGEBRAIC[13] = STATES[25]+STATES[26]+STATES[27]+STATES[28]+STATES[29]+STATES[30]+STATES[31]+STATES[32];
ALGEBRAIC[15] = ALGEBRAIC[3]+ALGEBRAIC[5]+ALGEBRAIC[7]+ALGEBRAIC[9]+ALGEBRAIC[11]+ALGEBRAIC[13];
ALGEBRAIC[62] = (ALGEBRAIC[59] - ALGEBRAIC[61])/CONSTANTS[95];
}