/* 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]; }