/* There are a total of 126 entries in the algebraic variable array. There are a total of 46 entries in each of the rate and state variable arrays. There are a total of 175 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[174] 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). * STATES[0] is q_LV in component main (litre). * CONSTANTS[6] is q_LV_us in component main (litre). * ALGEBRAIC[6] is u_LV in component main (kPa). * ALGEBRAIC[4] is E_LV in component main (kPa_per_L). * CONSTANTS[7] is E_LVmax in component main (kPa_per_L). * CONSTANTS[8] is E_LVmin in component main (kPa_per_L). * STATES[1] is q_RV in component main (litre). * CONSTANTS[9] is q_RV_us in component main (litre). * ALGEBRAIC[10] is u_RV in component main (kPa). * ALGEBRAIC[8] is E_RV in component main (kPa_per_L). * CONSTANTS[10] is E_RVmax in component main (kPa_per_L). * CONSTANTS[11] is E_RVmin in component main (kPa_per_L). * STATES[2] is q_LA in component main (litre). * CONSTANTS[12] is q_LA_us in component main (litre). * ALGEBRAIC[12] is u_LA in component main (kPa). * CONSTANTS[13] is E_LA in component main (kPa_per_L). * STATES[3] is q_RA in component main (litre). * CONSTANTS[14] is q_RA_us in component main (litre). * ALGEBRAIC[14] is u_RA in component main (kPa). * CONSTANTS[15] is E_RA in component main (kPa_per_L). * ALGEBRAIC[24] is v_AV in component main (L_per_s). * CONSTANTS[16] is R_AV in component main (kPa_s_per_L). * ALGEBRAIC[23] is v_PV in component main (L_per_s). * CONSTANTS[17] is R_PV in component main (kPa_s_per_L). * ALGEBRAIC[16] is v_TV in component main (L_per_s). * CONSTANTS[18] is R_TV in component main (kPa_s_per_L). * ALGEBRAIC[18] is v_MV in component main (L_per_s). * CONSTANTS[19] is R_MV in component main (kPa_s_per_L). * STATES[4] is q_pulmA in component main (litre). * CONSTANTS[20] is q_pulmA_us in component main (litre). * ALGEBRAIC[20] is u_pulmA in component main (kPa). * CONSTANTS[21] is E_pulmA in component main (kPa_per_L). * ALGEBRAIC[27] is v_lungA in component main (L_per_s). * CONSTANTS[22] is R_lungA in component main (kPa_s_per_L). * STATES[5] is q_lungC in component main (litre). * CONSTANTS[23] is q_lungC_us in component main (litre). * ALGEBRAIC[25] is u_lungC in component main (kPa). * CONSTANTS[24] is E_lungC in component main (kPa_per_L). * ALGEBRAIC[29] is v_pulmV in component main (L_per_s). * CONSTANTS[25] is R_pulmV in component main (kPa_s_per_L). * STATES[6] is q_AA in component main (litre). * CONSTANTS[26] is q_AA_us in component main (litre). * ALGEBRAIC[21] is u_AA in component main (kPa). * CONSTANTS[27] is E_AA in component main (kPa_per_L). * ALGEBRAIC[28] is v_brachiocephalicA in component main (L_per_s). * CONSTANTS[28] is R_brachiocephalicA in component main (kPa_s_per_L). * STATES[7] is q_brachiocephalicA in component main (litre). * CONSTANTS[29] is q_brachiocephalicA_us in component main (litre). * ALGEBRAIC[26] is u_brachiocephalicA in component main (kPa). * CONSTANTS[30] is E_brachiocephalicA in component main (kPa_per_L). * ALGEBRAIC[31] is v_rSubclavianA in component main (L_per_s). * CONSTANTS[31] is R_rSubclavianA in component main (kPa_s_per_L). * ALGEBRAIC[35] is v_rCCA in component main (L_per_s). * CONSTANTS[32] is R_rCCA in component main (kPa_s_per_L). * STATES[8] is q_rSubclavianA in component main (litre). * CONSTANTS[33] is q_rSubclavianA_us in component main (litre). * ALGEBRAIC[30] is u_rSubclavianA in component main (kPa). * CONSTANTS[34] is E_rSubclavianA in component main (kPa_per_L). * ALGEBRAIC[44] is v_rBrachialA in component main (L_per_s). * CONSTANTS[35] is R_rBrachialA in component main (kPa_s_per_L). * ALGEBRAIC[36] is v_rVertebralA in component main (L_per_s). * CONSTANTS[36] is R_rVertebralA in component main (kPa_s_per_L). * ALGEBRAIC[37] is v_lCCA in component main (L_per_s). * CONSTANTS[37] is R_lCCA in component main (kPa_s_per_L). * ALGEBRAIC[33] is v_lSubclavianA in component main (L_per_s). * CONSTANTS[38] is R_lSubclavianA in component main (kPa_s_per_L). * STATES[9] is q_lSubclavianA in component main (litre). * CONSTANTS[39] is q_lSubclavianA_us in component main (litre). * ALGEBRAIC[32] is u_lSubclavianA in component main (kPa). * CONSTANTS[40] is E_lSubclavianA in component main (kPa_per_L). * ALGEBRAIC[41] is v_lBrachialA in component main (L_per_s). * CONSTANTS[41] is R_lBrachialA in component main (kPa_s_per_L). * ALGEBRAIC[38] is v_lVertebralA in component main (L_per_s). * CONSTANTS[42] is R_lVertebralA in component main (kPa_s_per_L). * STATES[10] is q_headC in component main (litre). * CONSTANTS[43] is q_headC_us in component main (litre). * ALGEBRAIC[34] is u_headC in component main (kPa). * CONSTANTS[44] is E_headC in component main (kPa_per_L). * ALGEBRAIC[39] is v_headV in component main (L_per_s). * CONSTANTS[45] is R_headV in component main (kPa_s_per_L). * STATES[11] is q_lArmC in component main (litre). * CONSTANTS[46] is q_lArmC_us in component main (litre). * ALGEBRAIC[40] is u_lArmC in component main (kPa). * CONSTANTS[47] is E_lArmC in component main (kPa_per_L). * ALGEBRAIC[42] is v_lArmV in component main (L_per_s). * CONSTANTS[48] is R_lArmV in component main (kPa_s_per_L). * STATES[12] is q_rArmC in component main (litre). * CONSTANTS[49] is q_rArmC_us in component main (litre). * ALGEBRAIC[43] is u_rArmC in component main (kPa). * CONSTANTS[50] is E_rArmC in component main (kPa_per_L). * ALGEBRAIC[45] is v_rArmV in component main (L_per_s). * CONSTANTS[51] is R_rArmV in component main (kPa_s_per_L). * STATES[13] is q_A1 in component main (litre). * CONSTANTS[52] is q_A1_us in component main (litre). * STATES[14] is q_A2 in component main (litre). * CONSTANTS[53] is q_A2_us in component main (litre). * STATES[15] is q_A3 in component main (litre). * CONSTANTS[54] is q_A3_us in component main (litre). * STATES[16] is q_A4 in component main (litre). * CONSTANTS[55] is q_A4_us in component main (litre). * STATES[17] is q_A5 in component main (litre). * CONSTANTS[56] is q_A5_us in component main (litre). * STATES[18] is q_A6 in component main (litre). * CONSTANTS[57] is q_A6_us in component main (litre). * STATES[19] is q_A7 in component main (litre). * CONSTANTS[58] is q_A7_us in component main (litre). * STATES[20] is q_A8 in component main (litre). * CONSTANTS[59] is q_A8_us in component main (litre). * ALGEBRAIC[54] is v_A1 in component main (L_per_s). * ALGEBRAIC[55] is v_A2 in component main (L_per_s). * ALGEBRAIC[56] is v_A3 in component main (L_per_s). * ALGEBRAIC[57] is v_A4 in component main (L_per_s). * ALGEBRAIC[58] is v_A5 in component main (L_per_s). * ALGEBRAIC[59] is v_A6 in component main (L_per_s). * ALGEBRAIC[60] is v_A7 in component main (L_per_s). * ALGEBRAIC[61] is v_A8 in component main (L_per_s). * ALGEBRAIC[46] is u_A1 in component main (kPa). * ALGEBRAIC[47] is u_A2 in component main (kPa). * ALGEBRAIC[48] is u_A3 in component main (kPa). * ALGEBRAIC[49] is u_A4 in component main (kPa). * ALGEBRAIC[50] is u_A5 in component main (kPa). * ALGEBRAIC[51] is u_A6 in component main (kPa). * ALGEBRAIC[52] is u_A7 in component main (kPa). * ALGEBRAIC[53] is u_A8 in component main (kPa). * CONSTANTS[60] is E_A1 in component main (kPa_per_L). * CONSTANTS[61] is E_A2 in component main (kPa_per_L). * CONSTANTS[62] is E_A3 in component main (kPa_per_L). * CONSTANTS[63] is E_A4 in component main (kPa_per_L). * CONSTANTS[64] is E_A5 in component main (kPa_per_L). * CONSTANTS[65] is E_A6 in component main (kPa_per_L). * CONSTANTS[66] is E_A7 in component main (kPa_per_L). * CONSTANTS[67] is E_A8 in component main (kPa_per_L). * CONSTANTS[68] is R_A1 in component main (kPa_s_per_L). * CONSTANTS[69] is R_A2 in component main (kPa_s_per_L). * CONSTANTS[70] is R_A3 in component main (kPa_s_per_L). * CONSTANTS[71] is R_A4 in component main (kPa_s_per_L). * CONSTANTS[72] is R_A5 in component main (kPa_s_per_L). * CONSTANTS[73] is R_A6 in component main (kPa_s_per_L). * CONSTANTS[74] is R_A7 in component main (kPa_s_per_L). * CONSTANTS[75] is R_A8 in component main (kPa_s_per_L). * ALGEBRAIC[63] is v_lInfPhrenicA in component main (L_per_s). * CONSTANTS[76] is R_lInfPhrenicA in component main (kPa_s_per_L). * STATES[21] is q_lAdrenalGlandC in component main (litre). * CONSTANTS[77] is q_lAdrenalGlandC_us in component main (litre). * ALGEBRAIC[62] is u_lAdrenalGlandC in component main (kPa). * CONSTANTS[78] is E_lAdrenalGlandC in component main (kPa_per_L). * ALGEBRAIC[65] is v_lAdrenalV in component main (L_per_s). * CONSTANTS[79] is R_lAdrenalV in component main (kPa_s_per_L). * STATES[22] is q_lInfPhrenicV in component main (litre). * CONSTANTS[80] is q_lInfPhrenicV_us in component main (litre). * ALGEBRAIC[64] is u_lInfPhrenicV in component main (kPa). * CONSTANTS[81] is E_lInfPhrenicV in component main (kPa_per_L). * ALGEBRAIC[104] is v_lInfPhrenicV in component main (L_per_s). * CONSTANTS[82] is R_lInfPhrenicV in component main (kPa_s_per_L). * ALGEBRAIC[67] is v_rInfPhrenicA in component main (L_per_s). * CONSTANTS[83] is R_rInfPhrenicA in component main (kPa_s_per_L). * STATES[23] is q_rAdrenalGlandC in component main (litre). * CONSTANTS[84] is q_rAdrenalGlandC_us in component main (litre). * ALGEBRAIC[66] is u_rAdrenalGlandC in component main (kPa). * CONSTANTS[85] is E_rAdrenalGlandC in component main (kPa_per_L). * ALGEBRAIC[69] is v_rAdrenalV in component main (L_per_s). * CONSTANTS[86] is R_rAdrenalV in component main (kPa_s_per_L). * STATES[24] is q_rInfPhrenicV in component main (litre). * CONSTANTS[87] is q_rInfPhrenicV_us in component main (litre). * ALGEBRAIC[68] is u_rInfPhrenicV in component main (kPa). * CONSTANTS[88] is E_rInfPhrenicV in component main (kPa_per_L). * ALGEBRAIC[106] is v_rInfPhrenicV in component main (L_per_s). * CONSTANTS[89] is R_rInfPhrenicV in component main (kPa_s_per_L). * ALGEBRAIC[71] is v_celiacA in component main (L_per_s). * CONSTANTS[90] is R_celiacA in component main (kPa_s_per_L). * STATES[25] is q_celiacA in component main (litre). * CONSTANTS[91] is q_celiacA_us in component main (litre). * ALGEBRAIC[70] is u_celiacA in component main (kPa). * CONSTANTS[92] is E_celiacA in component main (kPa_per_L). * CONSTANTS[93] is z_celiacA in component main (meter). * ALGEBRAIC[73] is v_supMesA in component main (L_per_s). * CONSTANTS[94] is R_supMesA in component main (kPa_s_per_L). * STATES[26] is q_supMesA in component main (litre). * CONSTANTS[95] is q_supMesA_us in component main (litre). * ALGEBRAIC[72] is u_supMesA in component main (kPa). * CONSTANTS[96] is E_supMesA in component main (kPa_per_L). * CONSTANTS[97] is z_supMesA in component main (meter). * ALGEBRAIC[75] is v_lRenalA in component main (L_per_s). * CONSTANTS[98] is R_lRenalA in component main (kPa_s_per_L). * STATES[27] is q_lGlomerulusC in component main (litre). * CONSTANTS[99] is q_lGlomerulusC_us in component main (litre). * ALGEBRAIC[74] is u_lGlomerulusC in component main (kPa). * CONSTANTS[100] is E_lGlomerulusC in component main (kPa_per_L). * ALGEBRAIC[77] is v_lEffArteriol in component main (L_per_s). * CONSTANTS[101] is R_lEffArteriol in component main (kPa_s_per_L). * STATES[28] is q_lPeritubularC in component main (litre). * CONSTANTS[102] is q_lPeritubularC_us in component main (litre). * ALGEBRAIC[76] is u_lPeritubularC in component main (kPa). * CONSTANTS[103] is E_lPeritubularC in component main (kPa_per_L). * ALGEBRAIC[110] is v_lRenalV in component main (L_per_s). * CONSTANTS[104] is R_lRenalV in component main (kPa_s_per_L). * ALGEBRAIC[79] is v_rRenalA in component main (L_per_s). * CONSTANTS[105] is R_rRenalA in component main (kPa_s_per_L). * STATES[29] is q_rGlomerulusC in component main (litre). * CONSTANTS[106] is q_rGlomerulusC_us in component main (litre). * ALGEBRAIC[78] is u_rGlomerulusC in component main (kPa). * CONSTANTS[107] is E_rGlomerulusC in component main (kPa_per_L). * ALGEBRAIC[81] is v_rEffArteriol in component main (L_per_s). * CONSTANTS[108] is R_rEffArteriol in component main (kPa_s_per_L). * STATES[30] is q_rPeritubularC in component main (litre). * CONSTANTS[109] is q_rPeritubularC_us in component main (litre). * ALGEBRAIC[80] is u_rPeritubularC in component main (kPa). * CONSTANTS[110] is E_rPeritubularC in component main (kPa_per_L). * ALGEBRAIC[112] is v_rRenalV in component main (L_per_s). * CONSTANTS[111] is R_rRenalV in component main (kPa_s_per_L). * ALGEBRAIC[92] is v_infMesA in component main (L_per_s). * CONSTANTS[112] is R_infMesA in component main (kPa_s_per_L). * ALGEBRAIC[83] is v_stomachA in component main (L_per_s). * CONSTANTS[113] is R_stomachA in component main (kPa_s_per_L). * STATES[31] is q_stomachC in component main (litre). * CONSTANTS[114] is q_stomachC_us in component main (litre). * ALGEBRAIC[82] is u_stomachC in component main (kPa). * CONSTANTS[115] is E_stomachC in component main (kPa_per_L). * CONSTANTS[116] is z_stomachC in component main (meter). * ALGEBRAIC[95] is v_stomachV in component main (L_per_s). * CONSTANTS[117] is R_stomachV in component main (kPa_s_per_L). * ALGEBRAIC[85] is v_spleenA in component main (L_per_s). * CONSTANTS[118] is R_spleenA in component main (kPa_s_per_L). * STATES[32] is q_spleenC in component main (litre). * CONSTANTS[119] is q_spleenC_us in component main (litre). * ALGEBRAIC[84] is u_spleenC in component main (kPa). * CONSTANTS[120] is E_spleenC in component main (kPa_per_L). * CONSTANTS[121] is z_spleenC in component main (meter). * ALGEBRAIC[96] is v_spleenV in component main (L_per_s). * CONSTANTS[122] is R_spleenV in component main (kPa_s_per_L). * ALGEBRAIC[87] is v_pancreas1A in component main (L_per_s). * ALGEBRAIC[88] is v_pancreas2A in component main (L_per_s). * CONSTANTS[123] is R_pancreas1A in component main (kPa_s_per_L). * CONSTANTS[124] is R_pancreas2A in component main (kPa_s_per_L). * STATES[33] is q_pancreasC in component main (litre). * CONSTANTS[125] is q_pancreasC_us in component main (litre). * ALGEBRAIC[86] is u_pancreasC in component main (kPa). * CONSTANTS[126] is E_pancreasC in component main (kPa_per_L). * CONSTANTS[127] is z_pancreasC in component main (meter). * ALGEBRAIC[97] is v_pancreasV in component main (L_per_s). * CONSTANTS[128] is R_pancreasV in component main (kPa_s_per_L). * ALGEBRAIC[90] is v_intestineA in component main (L_per_s). * CONSTANTS[129] is R_intestineA in component main (kPa_s_per_L). * STATES[34] is q_intestineC in component main (litre). * CONSTANTS[130] is q_intestineC_us in component main (litre). * ALGEBRAIC[89] is u_intestineC in component main (kPa). * CONSTANTS[131] is E_intestineC in component main (kPa_per_L). * CONSTANTS[132] is z_ in component main (meter). * ALGEBRAIC[98] is v_intestineV in component main (L_per_s). * CONSTANTS[133] is R_intestineV in component main (kPa_s_per_L). * ALGEBRAIC[93] is v_colonA in component main (L_per_s). * CONSTANTS[134] is R_colonA in component main (kPa_s_per_L). * STATES[35] is q_colonC in component main (litre). * CONSTANTS[135] is q_colonC_us in component main (litre). * ALGEBRAIC[91] is u_colonC in component main (kPa). * CONSTANTS[136] is E_colonC in component main (kPa_per_L). * CONSTANTS[137] is z_colonC in component main (meter). * ALGEBRAIC[99] is v_colonV in component main (L_per_s). * CONSTANTS[138] is R_colonV in component main (kPa_s_per_L). * STATES[36] is q_portalV in component main (litre). * CONSTANTS[139] is q_portalV_us in component main (litre). * ALGEBRAIC[94] is u_portalV in component main (kPa). * CONSTANTS[140] is E_portalV in component main (kPa_per_L). * CONSTANTS[141] is z_portalV in component main (meter). * ALGEBRAIC[101] is v_portalV in component main (L_per_s). * CONSTANTS[142] is R_portalV in component main (kPa_s_per_L). * ALGEBRAIC[102] is v_liverA in component main (L_per_s). * CONSTANTS[143] is R_liverA in component main (kPa_s_per_L). * STATES[37] is q_liverC in component main (litre). * CONSTANTS[144] is q_liverC_us in component main (litre). * ALGEBRAIC[100] is u_liverC in component main (kPa). * CONSTANTS[145] is E_liverC in component main (kPa_per_L). * CONSTANTS[146] is z_liverC in component main (meter). * ALGEBRAIC[108] is v_liverV in component main (L_per_s). * CONSTANTS[147] is R_liverV in component main (kPa_s_per_L). * STATES[38] is q_V1 in component main (litre). * CONSTANTS[148] is q_V1_us in component main (litre). * STATES[39] is q_V2 in component main (litre). * CONSTANTS[149] is q_V2_us in component main (litre). * STATES[40] is q_V3 in component main (litre). * CONSTANTS[150] is q_V3_us in component main (litre). * STATES[41] is q_V4 in component main (litre). * CONSTANTS[151] is q_V4_us in component main (litre). * STATES[42] is q_V5 in component main (litre). * CONSTANTS[152] is q_V5_us in component main (litre). * STATES[43] is q_V6 in component main (litre). * CONSTANTS[153] is q_V6_us in component main (litre). * ALGEBRAIC[114] is v_V1 in component main (L_per_s). * ALGEBRAIC[115] is v_V2 in component main (L_per_s). * ALGEBRAIC[116] is v_V3 in component main (L_per_s). * ALGEBRAIC[117] is v_V4 in component main (L_per_s). * ALGEBRAIC[118] is v_V5 in component main (L_per_s). * ALGEBRAIC[119] is v_V6 in component main (L_per_s). * ALGEBRAIC[103] is u_V1 in component main (kPa). * ALGEBRAIC[105] is u_V2 in component main (kPa). * ALGEBRAIC[107] is u_V3 in component main (kPa). * ALGEBRAIC[109] is u_V4 in component main (kPa). * ALGEBRAIC[111] is u_V5 in component main (kPa). * ALGEBRAIC[113] is u_V6 in component main (kPa). * CONSTANTS[154] is E_V1 in component main (kPa_per_L). * CONSTANTS[155] is E_V2 in component main (kPa_per_L). * CONSTANTS[156] is E_V3 in component main (kPa_per_L). * CONSTANTS[157] is E_V4 in component main (kPa_per_L). * CONSTANTS[158] is E_V5 in component main (kPa_per_L). * CONSTANTS[159] is E_V6 in component main (kPa_per_L). * CONSTANTS[160] is R_V1 in component main (kPa_s_per_L). * CONSTANTS[161] is R_V2 in component main (kPa_s_per_L). * CONSTANTS[162] is R_V3 in component main (kPa_s_per_L). * CONSTANTS[163] is R_V4 in component main (kPa_s_per_L). * CONSTANTS[164] is R_V5 in component main (kPa_s_per_L). * CONSTANTS[165] is R_V6 in component main (kPa_s_per_L). * ALGEBRAIC[121] is v_lCommonIliacA in component main (L_per_s). * CONSTANTS[166] is R_lCommonIliacA in component main (kPa_s_per_L). * STATES[44] is q_lLegC in component main (litre). * CONSTANTS[167] is q_lLegC_us in component main (litre). * ALGEBRAIC[120] is u_lLegC in component main (kPa). * CONSTANTS[168] is E_lLegC in component main (kPa_per_L). * ALGEBRAIC[122] is v_lCommonIliacV in component main (L_per_s). * CONSTANTS[169] is R_lCommonIliacV in component main (kPa_s_per_L). * ALGEBRAIC[124] is v_rCommonIliacA in component main (L_per_s). * CONSTANTS[170] is R_rCommonIliacA in component main (kPa_s_per_L). * STATES[45] is q_rLegC in component main (litre). * CONSTANTS[171] is q_rLegC_us in component main (litre). * ALGEBRAIC[123] is u_rLegC in component main (kPa). * CONSTANTS[172] is E_rLegC in component main (kPa_per_L). * ALGEBRAIC[125] is v_rCommonIliacV in component main (L_per_s). * CONSTANTS[173] is R_rCommonIliacV in component main (kPa_s_per_L). * ALGEBRAIC[3] is q_TotHeart in component main (litre). * ALGEBRAIC[5] is q_TotLung in component main (litre). * ALGEBRAIC[7] is q_TotHead in component main (litre). * ALGEBRAIC[9] is q_TotThoracic in component main (litre). * ALGEBRAIC[11] is q_TotAorta in component main (litre). * ALGEBRAIC[13] is q_TotRenal in component main (litre). * ALGEBRAIC[15] is q_TotGut in component main (litre). * ALGEBRAIC[17] is q_TotVenaCava in component main (litre). * ALGEBRAIC[19] is q_TotLimbs in component main (litre). * ALGEBRAIC[22] is q_TotBlood in component main (litre). * RATES[0] is d/dt q_LV in component main (litre). * RATES[1] is d/dt q_RV in component main (litre). * RATES[2] is d/dt q_LA in component main (litre). * RATES[3] is d/dt q_RA in component main (litre). * RATES[4] is d/dt q_pulmA in component main (litre). * RATES[5] is d/dt q_lungC in component main (litre). * RATES[6] is d/dt q_AA in component main (litre). * RATES[7] is d/dt q_brachiocephalicA in component main (litre). * RATES[8] is d/dt q_rSubclavianA in component main (litre). * RATES[9] is d/dt q_lSubclavianA in component main (litre). * RATES[10] is d/dt q_headC in component main (litre). * RATES[11] is d/dt q_lArmC in component main (litre). * RATES[12] is d/dt q_rArmC in component main (litre). * RATES[13] is d/dt q_A1 in component main (litre). * RATES[14] is d/dt q_A2 in component main (litre). * RATES[15] is d/dt q_A3 in component main (litre). * RATES[16] is d/dt q_A4 in component main (litre). * RATES[17] is d/dt q_A5 in component main (litre). * RATES[18] is d/dt q_A6 in component main (litre). * RATES[19] is d/dt q_A7 in component main (litre). * RATES[20] is d/dt q_A8 in component main (litre). * RATES[21] is d/dt q_lAdrenalGlandC in component main (litre). * RATES[22] is d/dt q_lInfPhrenicV in component main (litre). * RATES[23] is d/dt q_rAdrenalGlandC in component main (litre). * RATES[24] is d/dt q_rInfPhrenicV in component main (litre). * RATES[25] is d/dt q_celiacA in component main (litre). * RATES[26] is d/dt q_supMesA in component main (litre). * RATES[27] is d/dt q_lGlomerulusC in component main (litre). * RATES[28] is d/dt q_lPeritubularC in component main (litre). * RATES[29] is d/dt q_rGlomerulusC in component main (litre). * RATES[30] is d/dt q_rPeritubularC in component main (litre). * RATES[31] is d/dt q_stomachC in component main (litre). * RATES[32] is d/dt q_spleenC in component main (litre). * RATES[33] is d/dt q_pancreasC in component main (litre). * RATES[34] is d/dt q_intestineC in component main (litre). * RATES[35] is d/dt q_colonC in component main (litre). * RATES[36] is d/dt q_portalV in component main (litre). * RATES[37] is d/dt q_liverC in component main (litre). * RATES[38] is d/dt q_V1 in component main (litre). * RATES[39] is d/dt q_V2 in component main (litre). * RATES[40] is d/dt q_V3 in component main (litre). * RATES[41] is d/dt q_V4 in component main (litre). * RATES[42] is d/dt q_V5 in component main (litre). * RATES[43] is d/dt q_V6 in component main (litre). * RATES[44] is d/dt q_lLegC in component main (litre). * RATES[45] is d/dt q_rLegC in component main (litre). */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { CONSTANTS[0] = 1000; CONSTANTS[1] = 9.81; CONSTANTS[2] = 90; CONSTANTS[3] = 1; CONSTANTS[4] = 0.1; CONSTANTS[5] = 0.1; STATES[0] = 2.24894867772612; CONSTANTS[6] = 0.005; CONSTANTS[7] = 320; CONSTANTS[8] = 6.667; STATES[1] = 2.24212874273495; CONSTANTS[9] = 0.01; CONSTANTS[10] = 70; CONSTANTS[11] = 2; STATES[2] = 1.50038322871398; CONSTANTS[12] = 0.004; CONSTANTS[13] = 10; STATES[3] = 0.582873822068968; CONSTANTS[14] = 0.004; CONSTANTS[15] = 10; CONSTANTS[16] = 1; CONSTANTS[17] = 1; CONSTANTS[18] = 1; CONSTANTS[19] = 1; STATES[4] = 0.274822042097049; CONSTANTS[20] = 0.2; CONSTANTS[21] = 200; CONSTANTS[22] = 1; STATES[5] = 0.274821056256821; CONSTANTS[23] = 0.2; CONSTANTS[24] = 200; CONSTANTS[25] = 1; STATES[6] = 1.43716360225001; CONSTANTS[26] = 0.12; CONSTANTS[27] = 80; CONSTANTS[28] = 100; STATES[7] = 0.658321340279823; CONSTANTS[29] = 0.2; CONSTANTS[30] = 200; CONSTANTS[31] = 100; CONSTANTS[32] = 100; STATES[8] = 0.578637981993128; CONSTANTS[33] = 0.2; CONSTANTS[34] = 200; CONSTANTS[35] = 100; CONSTANTS[36] = 100; CONSTANTS[37] = 100; CONSTANTS[38] = 100; STATES[9] = 0.609560299960373; CONSTANTS[39] = 0.2; CONSTANTS[40] = 200; CONSTANTS[41] = 100; CONSTANTS[42] = 100; STATES[10] = 0.575700923676085; CONSTANTS[43] = 0.2; CONSTANTS[44] = 200; CONSTANTS[45] = 100; STATES[11] = 0.380415025779693; CONSTANTS[46] = 0.06; CONSTANTS[47] = 150; CONSTANTS[48] = 100; STATES[12] = 0.380415025779692; CONSTANTS[49] = 0.06; CONSTANTS[50] = 150; CONSTANTS[51] = 100; STATES[13] = 0.0726453866278764; CONSTANTS[52] = 0.06; STATES[14] = 0.0722595904606296; CONSTANTS[53] = 0.06; STATES[15] = 0.0717727262354559; CONSTANTS[54] = 0.06; STATES[16] = 0.0715328891139091; CONSTANTS[55] = 0.06; STATES[17] = 0.0714293557557069; CONSTANTS[56] = 0.06; STATES[18] = 0.0713497500336812; CONSTANTS[57] = 0.06; STATES[19] = 0.0712949341825069; CONSTANTS[58] = 0.06; STATES[20] = 0.0712201889071758; CONSTANTS[59] = 0.06; CONSTANTS[60] = 400; CONSTANTS[61] = 400; CONSTANTS[62] = 400; CONSTANTS[63] = 400; CONSTANTS[64] = 400; CONSTANTS[65] = 400; CONSTANTS[66] = 400; CONSTANTS[67] = 400; CONSTANTS[68] = 10000000; CONSTANTS[69] = 100; CONSTANTS[70] = 100; CONSTANTS[71] = 100; CONSTANTS[72] = 100; CONSTANTS[73] = 100; CONSTANTS[74] = 100; CONSTANTS[75] = 100; CONSTANTS[76] = 100; STATES[21] = 0.226013520134949; CONSTANTS[77] = 0.2; CONSTANTS[78] = 200; CONSTANTS[79] = 100; STATES[22] = 0.226425525942081; CONSTANTS[80] = 0.2; CONSTANTS[81] = 200; CONSTANTS[82] = 100; CONSTANTS[83] = 100; STATES[23] = 0.224729694488936; CONSTANTS[84] = 0.2; CONSTANTS[85] = 200; CONSTANTS[86] = 100; STATES[24] = 0.224899115493557; CONSTANTS[87] = 0.2; CONSTANTS[88] = 200; CONSTANTS[89] = 100; CONSTANTS[90] = 100; STATES[25] = 0.0715428919998289; CONSTANTS[91] = 0.06; CONSTANTS[92] = 400; CONSTANTS[93] = 0.2; CONSTANTS[94] = 100; STATES[26] = 0.071414457768405; CONSTANTS[95] = 0.06; CONSTANTS[96] = 400; CONSTANTS[97] = 0.2; CONSTANTS[98] = 100; STATES[27] = 0.222843940557744; CONSTANTS[99] = 0.2; CONSTANTS[100] = 200; CONSTANTS[101] = 100; STATES[28] = 0.222892809988734; CONSTANTS[102] = 0.2; CONSTANTS[103] = 200; CONSTANTS[104] = 100; CONSTANTS[105] = 100; STATES[29] = 0.222680450943523; CONSTANTS[106] = 0.2; CONSTANTS[107] = 200; CONSTANTS[108] = 100; STATES[30] = 0.222721287850087; CONSTANTS[109] = 0.2; CONSTANTS[110] = 200; CONSTANTS[111] = 100; CONSTANTS[112] = 100; CONSTANTS[113] = 100; STATES[31] = 0.105910755847991; CONSTANTS[114] = 0.06; CONSTANTS[115] = 100; CONSTANTS[116] = 0.2; CONSTANTS[117] = 2400; CONSTANTS[118] = 100; STATES[32] = 0.105910034847806; CONSTANTS[119] = 0.06; CONSTANTS[120] = 100; CONSTANTS[121] = 0.2; CONSTANTS[122] = 3000; CONSTANTS[123] = 100; CONSTANTS[124] = 100; STATES[33] = 0.105787155209054; CONSTANTS[125] = 0.06; CONSTANTS[126] = 100; CONSTANTS[127] = 0.2; CONSTANTS[128] = 4000; CONSTANTS[129] = 100; STATES[34] = 0.105459895332545; CONSTANTS[130] = 0.06; CONSTANTS[131] = 100; CONSTANTS[132] = 0.2; CONSTANTS[133] = 1200; CONSTANTS[134] = 100; STATES[35] = 0.105317375726849; CONSTANTS[135] = 0.06; CONSTANTS[136] = 100; CONSTANTS[137] = 0.2; CONSTANTS[138] = 2000; STATES[36] = 0.152013945600315; CONSTANTS[139] = 0.06; CONSTANTS[140] = 50; CONSTANTS[141] = 0.2; CONSTANTS[142] = 100; CONSTANTS[143] = 800; STATES[37] = 0.526188169704937; CONSTANTS[144] = 0.06; CONSTANTS[145] = 10; CONSTANTS[146] = -0.5; CONSTANTS[147] = 1; STATES[38] = 0.0732527312969717; CONSTANTS[148] = 0.06; STATES[39] = 0.0724355014489091; CONSTANTS[149] = 0.06; STATES[40] = 0.0716607456049039; CONSTANTS[150] = 0.06; STATES[41] = 0.071501504166401; CONSTANTS[151] = 0.06; STATES[42] = 0.0714111536200381; CONSTANTS[152] = 0.06; STATES[43] = 0.0713861207574313; CONSTANTS[153] = 0.06; CONSTANTS[154] = 400; CONSTANTS[155] = 400; CONSTANTS[156] = 400; CONSTANTS[157] = 400; CONSTANTS[158] = 400; CONSTANTS[159] = 400; CONSTANTS[160] = 100; CONSTANTS[161] = 100; CONSTANTS[162] = 100; CONSTANTS[163] = 100; CONSTANTS[164] = 100; CONSTANTS[165] = 100; CONSTANTS[166] = 100; STATES[44] = 0.0898390162613648; CONSTANTS[167] = 0.06; CONSTANTS[168] = 150; CONSTANTS[169] = 4000; CONSTANTS[170] = 100; STATES[45] = 0.0898390162613648; CONSTANTS[171] = 0.06; CONSTANTS[172] = 150; CONSTANTS[173] = 4000; CONSTANTS[174] = ( 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[4] = (ALGEBRAIC[2]<0.400000 ? (( (CONSTANTS[7] - CONSTANTS[8])*ALGEBRAIC[2])/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(ALGEBRAIC[2]/CONSTANTS[5], 2.00000)))+CONSTANTS[8] : ALGEBRAIC[2]<0.600000 ? (( (CONSTANTS[7] - CONSTANTS[8])*0.400000)/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(0.400000/CONSTANTS[5], 2.00000)))+CONSTANTS[8] : CONSTANTS[8]); ALGEBRAIC[6] = ALGEBRAIC[4]*(STATES[0] - CONSTANTS[6]); ALGEBRAIC[21] = CONSTANTS[27]*(STATES[6] - CONSTANTS[26]); ALGEBRAIC[24] = (ALGEBRAIC[6]>ALGEBRAIC[21] ? (ALGEBRAIC[6] - ALGEBRAIC[21])/CONSTANTS[16] : 0.00000); ALGEBRAIC[12] = CONSTANTS[13]*(STATES[2] - CONSTANTS[12]); ALGEBRAIC[18] = (ALGEBRAIC[12]>ALGEBRAIC[6] ? (ALGEBRAIC[12] - ALGEBRAIC[6])/CONSTANTS[19] : 0.00000); RATES[0] = ALGEBRAIC[18] - ALGEBRAIC[24]; ALGEBRAIC[8] = (ALGEBRAIC[2]<0.400000 ? (( (CONSTANTS[10] - CONSTANTS[11])*ALGEBRAIC[2])/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(ALGEBRAIC[2]/CONSTANTS[5], 2.00000)))+CONSTANTS[11] : ALGEBRAIC[2]<0.600000 ? (( (CONSTANTS[10] - CONSTANTS[11])*0.400000)/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(0.400000/CONSTANTS[5], 2.00000)))+CONSTANTS[11] : CONSTANTS[11]); ALGEBRAIC[10] = ALGEBRAIC[8]*(STATES[1] - CONSTANTS[9]); ALGEBRAIC[20] = CONSTANTS[21]*(STATES[4] - CONSTANTS[20]); ALGEBRAIC[23] = (ALGEBRAIC[10]>ALGEBRAIC[20] ? (ALGEBRAIC[10] - ALGEBRAIC[20])/CONSTANTS[17] : 0.00000); ALGEBRAIC[14] = CONSTANTS[15]*(STATES[3] - CONSTANTS[14]); ALGEBRAIC[16] = (ALGEBRAIC[14]>ALGEBRAIC[10] ? (ALGEBRAIC[14] - ALGEBRAIC[10])/CONSTANTS[18] : 0.00000); RATES[1] = ALGEBRAIC[16] - ALGEBRAIC[23]; ALGEBRAIC[25] = CONSTANTS[24]*(STATES[5] - CONSTANTS[23]); ALGEBRAIC[27] = (ALGEBRAIC[20] - ALGEBRAIC[25])/CONSTANTS[22]; RATES[4] = ALGEBRAIC[23] - ALGEBRAIC[27]; ALGEBRAIC[29] = (ALGEBRAIC[25] - ALGEBRAIC[12])/CONSTANTS[25]; RATES[2] = ALGEBRAIC[29] - ALGEBRAIC[18]; RATES[5] = ALGEBRAIC[27] - ALGEBRAIC[29]; ALGEBRAIC[26] = CONSTANTS[30]*(STATES[7] - CONSTANTS[29]); ALGEBRAIC[28] = (ALGEBRAIC[21] - ALGEBRAIC[26])/CONSTANTS[28]; ALGEBRAIC[30] = CONSTANTS[34]*(STATES[8] - CONSTANTS[33]); ALGEBRAIC[31] = (ALGEBRAIC[26] - ALGEBRAIC[30])/CONSTANTS[31]; ALGEBRAIC[34] = CONSTANTS[44]*(STATES[10] - CONSTANTS[43]); ALGEBRAIC[35] = (ALGEBRAIC[26] - ALGEBRAIC[34])/CONSTANTS[32]; RATES[7] = (ALGEBRAIC[28] - ALGEBRAIC[31]) - ALGEBRAIC[35]; ALGEBRAIC[36] = (ALGEBRAIC[30] - ALGEBRAIC[34])/CONSTANTS[36]; ALGEBRAIC[37] = (ALGEBRAIC[21] - ALGEBRAIC[34])/CONSTANTS[37]; ALGEBRAIC[32] = CONSTANTS[40]*(STATES[9] - CONSTANTS[39]); ALGEBRAIC[38] = (ALGEBRAIC[32] - ALGEBRAIC[34])/CONSTANTS[42]; ALGEBRAIC[39] = (ALGEBRAIC[34] - ALGEBRAIC[14])/CONSTANTS[45]; RATES[10] = (ALGEBRAIC[36]+ALGEBRAIC[35]+ALGEBRAIC[37]+ALGEBRAIC[38]) - ALGEBRAIC[39]; ALGEBRAIC[33] = (ALGEBRAIC[21] - ALGEBRAIC[32])/CONSTANTS[38]; ALGEBRAIC[40] = CONSTANTS[47]*(STATES[11] - CONSTANTS[46]); ALGEBRAIC[41] = (ALGEBRAIC[32] - ALGEBRAIC[40])/CONSTANTS[41]; RATES[9] = (ALGEBRAIC[33] - ALGEBRAIC[41]) - ALGEBRAIC[38]; ALGEBRAIC[42] = (ALGEBRAIC[40] - ALGEBRAIC[14])/CONSTANTS[48]; RATES[11] = ALGEBRAIC[41] - ALGEBRAIC[42]; ALGEBRAIC[43] = CONSTANTS[50]*(STATES[12] - CONSTANTS[49]); ALGEBRAIC[44] = (ALGEBRAIC[30] - ALGEBRAIC[43])/CONSTANTS[35]; RATES[8] = (ALGEBRAIC[31] - ALGEBRAIC[44]) - ALGEBRAIC[36]; ALGEBRAIC[45] = (ALGEBRAIC[43] - ALGEBRAIC[14])/CONSTANTS[51]; RATES[12] = ALGEBRAIC[44] - ALGEBRAIC[45]; ALGEBRAIC[46] = CONSTANTS[60]*(STATES[13] - CONSTANTS[52]); ALGEBRAIC[54] = (ALGEBRAIC[21] - ALGEBRAIC[46])/CONSTANTS[68]; RATES[6] = (((ALGEBRAIC[24] - ALGEBRAIC[54]) - ALGEBRAIC[28]) - ALGEBRAIC[37]) - ALGEBRAIC[33]; ALGEBRAIC[47] = CONSTANTS[61]*(STATES[14] - CONSTANTS[53]); ALGEBRAIC[55] = (ALGEBRAIC[46] - ALGEBRAIC[47])/CONSTANTS[69]; ALGEBRAIC[62] = CONSTANTS[78]*(STATES[21] - CONSTANTS[77]); ALGEBRAIC[63] = (ALGEBRAIC[46] - ALGEBRAIC[62])/CONSTANTS[76]; RATES[13] = (ALGEBRAIC[54] - ALGEBRAIC[55]) - ALGEBRAIC[63]; ALGEBRAIC[64] = CONSTANTS[81]*(STATES[22] - CONSTANTS[80]); ALGEBRAIC[65] = (ALGEBRAIC[62] - ALGEBRAIC[64])/CONSTANTS[79]; RATES[21] = ALGEBRAIC[63] - ALGEBRAIC[65]; ALGEBRAIC[48] = CONSTANTS[62]*(STATES[15] - CONSTANTS[54]); ALGEBRAIC[56] = (ALGEBRAIC[47] - ALGEBRAIC[48])/CONSTANTS[70]; ALGEBRAIC[66] = CONSTANTS[85]*(STATES[23] - CONSTANTS[84]); ALGEBRAIC[67] = (ALGEBRAIC[47] - ALGEBRAIC[66])/CONSTANTS[83]; RATES[14] = (ALGEBRAIC[55] - ALGEBRAIC[56]) - ALGEBRAIC[67]; ALGEBRAIC[68] = CONSTANTS[88]*(STATES[24] - CONSTANTS[87]); ALGEBRAIC[69] = (ALGEBRAIC[66] - ALGEBRAIC[68])/CONSTANTS[86]; RATES[23] = ALGEBRAIC[67] - ALGEBRAIC[69]; ALGEBRAIC[49] = CONSTANTS[63]*(STATES[16] - CONSTANTS[55]); ALGEBRAIC[57] = (ALGEBRAIC[48] - ALGEBRAIC[49])/CONSTANTS[71]; ALGEBRAIC[70] = CONSTANTS[92]*(STATES[25] - CONSTANTS[91]); ALGEBRAIC[71] = (ALGEBRAIC[48] - ALGEBRAIC[70])/CONSTANTS[90]; RATES[15] = (ALGEBRAIC[56] - ALGEBRAIC[57]) - ALGEBRAIC[71]; ALGEBRAIC[50] = CONSTANTS[64]*(STATES[17] - CONSTANTS[56]); ALGEBRAIC[58] = (ALGEBRAIC[49] - ALGEBRAIC[50])/CONSTANTS[72]; ALGEBRAIC[72] = CONSTANTS[96]*(STATES[26] - CONSTANTS[95]); ALGEBRAIC[73] = (ALGEBRAIC[49] - ALGEBRAIC[72])/CONSTANTS[94]; RATES[16] = (ALGEBRAIC[57] - ALGEBRAIC[58]) - ALGEBRAIC[73]; ALGEBRAIC[51] = CONSTANTS[65]*(STATES[18] - CONSTANTS[57]); ALGEBRAIC[59] = (ALGEBRAIC[50] - ALGEBRAIC[51])/CONSTANTS[73]; ALGEBRAIC[74] = CONSTANTS[100]*(STATES[27] - CONSTANTS[99]); ALGEBRAIC[75] = (ALGEBRAIC[50] - ALGEBRAIC[74])/CONSTANTS[98]; RATES[17] = (ALGEBRAIC[58] - ALGEBRAIC[59]) - ALGEBRAIC[75]; ALGEBRAIC[76] = CONSTANTS[103]*(STATES[28] - CONSTANTS[102]); ALGEBRAIC[77] = (ALGEBRAIC[74] - ALGEBRAIC[76])/CONSTANTS[101]; RATES[27] = ALGEBRAIC[75] - ALGEBRAIC[77]; ALGEBRAIC[52] = CONSTANTS[66]*(STATES[19] - CONSTANTS[58]); ALGEBRAIC[60] = (ALGEBRAIC[51] - ALGEBRAIC[52])/CONSTANTS[74]; ALGEBRAIC[78] = CONSTANTS[107]*(STATES[29] - CONSTANTS[106]); ALGEBRAIC[79] = (ALGEBRAIC[51] - ALGEBRAIC[78])/CONSTANTS[105]; RATES[18] = (ALGEBRAIC[59] - ALGEBRAIC[60]) - ALGEBRAIC[79]; ALGEBRAIC[80] = CONSTANTS[110]*(STATES[30] - CONSTANTS[109]); ALGEBRAIC[81] = (ALGEBRAIC[78] - ALGEBRAIC[80])/CONSTANTS[108]; RATES[29] = ALGEBRAIC[79] - ALGEBRAIC[81]; ALGEBRAIC[53] = CONSTANTS[67]*(STATES[20] - CONSTANTS[59]); ALGEBRAIC[61] = (ALGEBRAIC[52] - ALGEBRAIC[53])/CONSTANTS[75]; ALGEBRAIC[91] = CONSTANTS[136]*(STATES[35] - CONSTANTS[135]); ALGEBRAIC[92] = (ALGEBRAIC[52] - ALGEBRAIC[91])/CONSTANTS[112]; RATES[19] = (ALGEBRAIC[60] - ALGEBRAIC[61]) - ALGEBRAIC[92]; ALGEBRAIC[86] = CONSTANTS[126]*(STATES[33] - CONSTANTS[125]); ALGEBRAIC[88] = (ALGEBRAIC[72] - ALGEBRAIC[86])/CONSTANTS[124]; ALGEBRAIC[89] = CONSTANTS[131]*(STATES[34] - CONSTANTS[130]); ALGEBRAIC[90] = (ALGEBRAIC[72] - ALGEBRAIC[89])/CONSTANTS[129]; ALGEBRAIC[93] = (ALGEBRAIC[72] - ALGEBRAIC[91])/CONSTANTS[134]; RATES[26] = ((ALGEBRAIC[73] - ALGEBRAIC[88]) - ALGEBRAIC[90]) - ALGEBRAIC[93]; ALGEBRAIC[82] = CONSTANTS[115]*(STATES[31] - CONSTANTS[114]); ALGEBRAIC[83] = (ALGEBRAIC[70] - ALGEBRAIC[82])/CONSTANTS[113]; ALGEBRAIC[94] = CONSTANTS[140]*(STATES[36] - CONSTANTS[139]); ALGEBRAIC[95] = (ALGEBRAIC[82] - ALGEBRAIC[94])/CONSTANTS[117]; RATES[31] = ALGEBRAIC[83] - ALGEBRAIC[95]; ALGEBRAIC[84] = CONSTANTS[120]*(STATES[32] - CONSTANTS[119]); ALGEBRAIC[85] = (ALGEBRAIC[70] - ALGEBRAIC[84])/CONSTANTS[118]; ALGEBRAIC[96] = (ALGEBRAIC[84] - ALGEBRAIC[94])/CONSTANTS[122]; RATES[32] = ALGEBRAIC[85] - ALGEBRAIC[96]; ALGEBRAIC[87] = (ALGEBRAIC[70] - ALGEBRAIC[86])/CONSTANTS[123]; ALGEBRAIC[97] = (ALGEBRAIC[86] - ALGEBRAIC[94])/CONSTANTS[128]; RATES[33] = (ALGEBRAIC[87]+ALGEBRAIC[88]) - ALGEBRAIC[97]; ALGEBRAIC[98] = (ALGEBRAIC[89] - ALGEBRAIC[94])/CONSTANTS[133]; RATES[34] = ALGEBRAIC[90] - ALGEBRAIC[98]; ALGEBRAIC[99] = (ALGEBRAIC[91] - ALGEBRAIC[94])/CONSTANTS[138]; RATES[35] = (ALGEBRAIC[93]+ALGEBRAIC[92]) - ALGEBRAIC[99]; ALGEBRAIC[100] = CONSTANTS[145]*(STATES[37] - CONSTANTS[144]); ALGEBRAIC[101] = (ALGEBRAIC[94] - ALGEBRAIC[100])/CONSTANTS[142]; RATES[36] = (ALGEBRAIC[95]+ALGEBRAIC[96]+ALGEBRAIC[97]+ALGEBRAIC[98]+ALGEBRAIC[99]) - ALGEBRAIC[101]; ALGEBRAIC[102] = (ALGEBRAIC[70] - ALGEBRAIC[100])/CONSTANTS[143]; RATES[25] = (((ALGEBRAIC[71] - ALGEBRAIC[102]) - ALGEBRAIC[83]) - ALGEBRAIC[85]) - ALGEBRAIC[87]; ALGEBRAIC[103] = CONSTANTS[154]*(STATES[38] - CONSTANTS[148]); ALGEBRAIC[104] = (ALGEBRAIC[64] - ALGEBRAIC[103])/CONSTANTS[82]; RATES[22] = ALGEBRAIC[65] - ALGEBRAIC[104]; ALGEBRAIC[105] = CONSTANTS[155]*(STATES[39] - CONSTANTS[149]); ALGEBRAIC[106] = (ALGEBRAIC[68] - ALGEBRAIC[105])/CONSTANTS[89]; RATES[24] = ALGEBRAIC[69] - ALGEBRAIC[106]; ALGEBRAIC[107] = CONSTANTS[156]*(STATES[40] - CONSTANTS[150]); ALGEBRAIC[108] = (ALGEBRAIC[100] - ALGEBRAIC[107])/CONSTANTS[147]; RATES[37] = (ALGEBRAIC[102]+ALGEBRAIC[101]) - ALGEBRAIC[108]; ALGEBRAIC[109] = CONSTANTS[157]*(STATES[41] - CONSTANTS[151]); ALGEBRAIC[110] = (ALGEBRAIC[76] - ALGEBRAIC[109])/CONSTANTS[104]; RATES[28] = ALGEBRAIC[77] - ALGEBRAIC[110]; ALGEBRAIC[111] = CONSTANTS[158]*(STATES[42] - CONSTANTS[152]); ALGEBRAIC[112] = (ALGEBRAIC[80] - ALGEBRAIC[111])/CONSTANTS[111]; RATES[30] = ALGEBRAIC[81] - ALGEBRAIC[112]; ALGEBRAIC[114] = (ALGEBRAIC[103] - ALGEBRAIC[14])/CONSTANTS[160]; RATES[3] = (ALGEBRAIC[39]+ALGEBRAIC[42]+ALGEBRAIC[45]+ALGEBRAIC[114]) - ALGEBRAIC[16]; ALGEBRAIC[115] = (ALGEBRAIC[105] - ALGEBRAIC[103])/CONSTANTS[161]; RATES[38] = (ALGEBRAIC[115]+ALGEBRAIC[104]) - ALGEBRAIC[114]; ALGEBRAIC[116] = (ALGEBRAIC[107] - ALGEBRAIC[105])/CONSTANTS[162]; RATES[39] = (ALGEBRAIC[116]+ALGEBRAIC[106]) - ALGEBRAIC[115]; ALGEBRAIC[117] = (ALGEBRAIC[109] - ALGEBRAIC[107])/CONSTANTS[163]; RATES[40] = (ALGEBRAIC[117]+ALGEBRAIC[108]) - ALGEBRAIC[116]; ALGEBRAIC[118] = (ALGEBRAIC[111] - ALGEBRAIC[109])/CONSTANTS[164]; RATES[41] = (ALGEBRAIC[118]+ALGEBRAIC[110]) - ALGEBRAIC[117]; ALGEBRAIC[113] = CONSTANTS[159]*(STATES[43] - CONSTANTS[153]); ALGEBRAIC[119] = (ALGEBRAIC[113] - ALGEBRAIC[111])/CONSTANTS[165]; RATES[42] = (ALGEBRAIC[119]+ALGEBRAIC[112]) - ALGEBRAIC[118]; ALGEBRAIC[120] = CONSTANTS[168]*(STATES[44] - CONSTANTS[167]); ALGEBRAIC[121] = (ALGEBRAIC[53] - ALGEBRAIC[120])/CONSTANTS[166]; ALGEBRAIC[122] = (ALGEBRAIC[120] - ALGEBRAIC[113])/CONSTANTS[169]; RATES[44] = ALGEBRAIC[121] - ALGEBRAIC[122]; ALGEBRAIC[123] = CONSTANTS[172]*(STATES[45] - CONSTANTS[171]); ALGEBRAIC[124] = (ALGEBRAIC[53] - ALGEBRAIC[123])/CONSTANTS[170]; RATES[20] = (ALGEBRAIC[61] - ALGEBRAIC[121]) - ALGEBRAIC[124]; ALGEBRAIC[125] = (ALGEBRAIC[123] - ALGEBRAIC[113])/CONSTANTS[173]; RATES[43] = (ALGEBRAIC[122]+ALGEBRAIC[125]) - ALGEBRAIC[119]; RATES[45] = ALGEBRAIC[124] - ALGEBRAIC[125]; } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[2] = VOI/CONSTANTS[3] - floor(VOI/CONSTANTS[3]); ALGEBRAIC[4] = (ALGEBRAIC[2]<0.400000 ? (( (CONSTANTS[7] - CONSTANTS[8])*ALGEBRAIC[2])/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(ALGEBRAIC[2]/CONSTANTS[5], 2.00000)))+CONSTANTS[8] : ALGEBRAIC[2]<0.600000 ? (( (CONSTANTS[7] - CONSTANTS[8])*0.400000)/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(0.400000/CONSTANTS[5], 2.00000)))+CONSTANTS[8] : CONSTANTS[8]); ALGEBRAIC[6] = ALGEBRAIC[4]*(STATES[0] - CONSTANTS[6]); ALGEBRAIC[21] = CONSTANTS[27]*(STATES[6] - CONSTANTS[26]); ALGEBRAIC[24] = (ALGEBRAIC[6]>ALGEBRAIC[21] ? (ALGEBRAIC[6] - ALGEBRAIC[21])/CONSTANTS[16] : 0.00000); ALGEBRAIC[12] = CONSTANTS[13]*(STATES[2] - CONSTANTS[12]); ALGEBRAIC[18] = (ALGEBRAIC[12]>ALGEBRAIC[6] ? (ALGEBRAIC[12] - ALGEBRAIC[6])/CONSTANTS[19] : 0.00000); ALGEBRAIC[8] = (ALGEBRAIC[2]<0.400000 ? (( (CONSTANTS[10] - CONSTANTS[11])*ALGEBRAIC[2])/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(ALGEBRAIC[2]/CONSTANTS[5], 2.00000)))+CONSTANTS[11] : ALGEBRAIC[2]<0.600000 ? (( (CONSTANTS[10] - CONSTANTS[11])*0.400000)/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(0.400000/CONSTANTS[5], 2.00000)))+CONSTANTS[11] : CONSTANTS[11]); ALGEBRAIC[10] = ALGEBRAIC[8]*(STATES[1] - CONSTANTS[9]); ALGEBRAIC[20] = CONSTANTS[21]*(STATES[4] - CONSTANTS[20]); ALGEBRAIC[23] = (ALGEBRAIC[10]>ALGEBRAIC[20] ? (ALGEBRAIC[10] - ALGEBRAIC[20])/CONSTANTS[17] : 0.00000); ALGEBRAIC[14] = CONSTANTS[15]*(STATES[3] - CONSTANTS[14]); ALGEBRAIC[16] = (ALGEBRAIC[14]>ALGEBRAIC[10] ? (ALGEBRAIC[14] - ALGEBRAIC[10])/CONSTANTS[18] : 0.00000); ALGEBRAIC[25] = CONSTANTS[24]*(STATES[5] - CONSTANTS[23]); ALGEBRAIC[27] = (ALGEBRAIC[20] - ALGEBRAIC[25])/CONSTANTS[22]; ALGEBRAIC[29] = (ALGEBRAIC[25] - ALGEBRAIC[12])/CONSTANTS[25]; ALGEBRAIC[26] = CONSTANTS[30]*(STATES[7] - CONSTANTS[29]); ALGEBRAIC[28] = (ALGEBRAIC[21] - ALGEBRAIC[26])/CONSTANTS[28]; ALGEBRAIC[30] = CONSTANTS[34]*(STATES[8] - CONSTANTS[33]); ALGEBRAIC[31] = (ALGEBRAIC[26] - ALGEBRAIC[30])/CONSTANTS[31]; ALGEBRAIC[34] = CONSTANTS[44]*(STATES[10] - CONSTANTS[43]); ALGEBRAIC[35] = (ALGEBRAIC[26] - ALGEBRAIC[34])/CONSTANTS[32]; ALGEBRAIC[36] = (ALGEBRAIC[30] - ALGEBRAIC[34])/CONSTANTS[36]; ALGEBRAIC[37] = (ALGEBRAIC[21] - ALGEBRAIC[34])/CONSTANTS[37]; ALGEBRAIC[32] = CONSTANTS[40]*(STATES[9] - CONSTANTS[39]); ALGEBRAIC[38] = (ALGEBRAIC[32] - ALGEBRAIC[34])/CONSTANTS[42]; ALGEBRAIC[39] = (ALGEBRAIC[34] - ALGEBRAIC[14])/CONSTANTS[45]; ALGEBRAIC[33] = (ALGEBRAIC[21] - ALGEBRAIC[32])/CONSTANTS[38]; ALGEBRAIC[40] = CONSTANTS[47]*(STATES[11] - CONSTANTS[46]); ALGEBRAIC[41] = (ALGEBRAIC[32] - ALGEBRAIC[40])/CONSTANTS[41]; ALGEBRAIC[42] = (ALGEBRAIC[40] - ALGEBRAIC[14])/CONSTANTS[48]; ALGEBRAIC[43] = CONSTANTS[50]*(STATES[12] - CONSTANTS[49]); ALGEBRAIC[44] = (ALGEBRAIC[30] - ALGEBRAIC[43])/CONSTANTS[35]; ALGEBRAIC[45] = (ALGEBRAIC[43] - ALGEBRAIC[14])/CONSTANTS[51]; ALGEBRAIC[46] = CONSTANTS[60]*(STATES[13] - CONSTANTS[52]); ALGEBRAIC[54] = (ALGEBRAIC[21] - ALGEBRAIC[46])/CONSTANTS[68]; ALGEBRAIC[47] = CONSTANTS[61]*(STATES[14] - CONSTANTS[53]); ALGEBRAIC[55] = (ALGEBRAIC[46] - ALGEBRAIC[47])/CONSTANTS[69]; ALGEBRAIC[62] = CONSTANTS[78]*(STATES[21] - CONSTANTS[77]); ALGEBRAIC[63] = (ALGEBRAIC[46] - ALGEBRAIC[62])/CONSTANTS[76]; ALGEBRAIC[64] = CONSTANTS[81]*(STATES[22] - CONSTANTS[80]); ALGEBRAIC[65] = (ALGEBRAIC[62] - ALGEBRAIC[64])/CONSTANTS[79]; ALGEBRAIC[48] = CONSTANTS[62]*(STATES[15] - CONSTANTS[54]); ALGEBRAIC[56] = (ALGEBRAIC[47] - ALGEBRAIC[48])/CONSTANTS[70]; ALGEBRAIC[66] = CONSTANTS[85]*(STATES[23] - CONSTANTS[84]); ALGEBRAIC[67] = (ALGEBRAIC[47] - ALGEBRAIC[66])/CONSTANTS[83]; ALGEBRAIC[68] = CONSTANTS[88]*(STATES[24] - CONSTANTS[87]); ALGEBRAIC[69] = (ALGEBRAIC[66] - ALGEBRAIC[68])/CONSTANTS[86]; ALGEBRAIC[49] = CONSTANTS[63]*(STATES[16] - CONSTANTS[55]); ALGEBRAIC[57] = (ALGEBRAIC[48] - ALGEBRAIC[49])/CONSTANTS[71]; ALGEBRAIC[70] = CONSTANTS[92]*(STATES[25] - CONSTANTS[91]); ALGEBRAIC[71] = (ALGEBRAIC[48] - ALGEBRAIC[70])/CONSTANTS[90]; ALGEBRAIC[50] = CONSTANTS[64]*(STATES[17] - CONSTANTS[56]); ALGEBRAIC[58] = (ALGEBRAIC[49] - ALGEBRAIC[50])/CONSTANTS[72]; ALGEBRAIC[72] = CONSTANTS[96]*(STATES[26] - CONSTANTS[95]); ALGEBRAIC[73] = (ALGEBRAIC[49] - ALGEBRAIC[72])/CONSTANTS[94]; ALGEBRAIC[51] = CONSTANTS[65]*(STATES[18] - CONSTANTS[57]); ALGEBRAIC[59] = (ALGEBRAIC[50] - ALGEBRAIC[51])/CONSTANTS[73]; ALGEBRAIC[74] = CONSTANTS[100]*(STATES[27] - CONSTANTS[99]); ALGEBRAIC[75] = (ALGEBRAIC[50] - ALGEBRAIC[74])/CONSTANTS[98]; ALGEBRAIC[76] = CONSTANTS[103]*(STATES[28] - CONSTANTS[102]); ALGEBRAIC[77] = (ALGEBRAIC[74] - ALGEBRAIC[76])/CONSTANTS[101]; ALGEBRAIC[52] = CONSTANTS[66]*(STATES[19] - CONSTANTS[58]); ALGEBRAIC[60] = (ALGEBRAIC[51] - ALGEBRAIC[52])/CONSTANTS[74]; ALGEBRAIC[78] = CONSTANTS[107]*(STATES[29] - CONSTANTS[106]); ALGEBRAIC[79] = (ALGEBRAIC[51] - ALGEBRAIC[78])/CONSTANTS[105]; ALGEBRAIC[80] = CONSTANTS[110]*(STATES[30] - CONSTANTS[109]); ALGEBRAIC[81] = (ALGEBRAIC[78] - ALGEBRAIC[80])/CONSTANTS[108]; ALGEBRAIC[53] = CONSTANTS[67]*(STATES[20] - CONSTANTS[59]); ALGEBRAIC[61] = (ALGEBRAIC[52] - ALGEBRAIC[53])/CONSTANTS[75]; ALGEBRAIC[91] = CONSTANTS[136]*(STATES[35] - CONSTANTS[135]); ALGEBRAIC[92] = (ALGEBRAIC[52] - ALGEBRAIC[91])/CONSTANTS[112]; ALGEBRAIC[86] = CONSTANTS[126]*(STATES[33] - CONSTANTS[125]); ALGEBRAIC[88] = (ALGEBRAIC[72] - ALGEBRAIC[86])/CONSTANTS[124]; ALGEBRAIC[89] = CONSTANTS[131]*(STATES[34] - CONSTANTS[130]); ALGEBRAIC[90] = (ALGEBRAIC[72] - ALGEBRAIC[89])/CONSTANTS[129]; ALGEBRAIC[93] = (ALGEBRAIC[72] - ALGEBRAIC[91])/CONSTANTS[134]; ALGEBRAIC[82] = CONSTANTS[115]*(STATES[31] - CONSTANTS[114]); ALGEBRAIC[83] = (ALGEBRAIC[70] - ALGEBRAIC[82])/CONSTANTS[113]; ALGEBRAIC[94] = CONSTANTS[140]*(STATES[36] - CONSTANTS[139]); ALGEBRAIC[95] = (ALGEBRAIC[82] - ALGEBRAIC[94])/CONSTANTS[117]; ALGEBRAIC[84] = CONSTANTS[120]*(STATES[32] - CONSTANTS[119]); ALGEBRAIC[85] = (ALGEBRAIC[70] - ALGEBRAIC[84])/CONSTANTS[118]; ALGEBRAIC[96] = (ALGEBRAIC[84] - ALGEBRAIC[94])/CONSTANTS[122]; ALGEBRAIC[87] = (ALGEBRAIC[70] - ALGEBRAIC[86])/CONSTANTS[123]; ALGEBRAIC[97] = (ALGEBRAIC[86] - ALGEBRAIC[94])/CONSTANTS[128]; ALGEBRAIC[98] = (ALGEBRAIC[89] - ALGEBRAIC[94])/CONSTANTS[133]; ALGEBRAIC[99] = (ALGEBRAIC[91] - ALGEBRAIC[94])/CONSTANTS[138]; ALGEBRAIC[100] = CONSTANTS[145]*(STATES[37] - CONSTANTS[144]); ALGEBRAIC[101] = (ALGEBRAIC[94] - ALGEBRAIC[100])/CONSTANTS[142]; ALGEBRAIC[102] = (ALGEBRAIC[70] - ALGEBRAIC[100])/CONSTANTS[143]; ALGEBRAIC[103] = CONSTANTS[154]*(STATES[38] - CONSTANTS[148]); ALGEBRAIC[104] = (ALGEBRAIC[64] - ALGEBRAIC[103])/CONSTANTS[82]; ALGEBRAIC[105] = CONSTANTS[155]*(STATES[39] - CONSTANTS[149]); ALGEBRAIC[106] = (ALGEBRAIC[68] - ALGEBRAIC[105])/CONSTANTS[89]; ALGEBRAIC[107] = CONSTANTS[156]*(STATES[40] - CONSTANTS[150]); ALGEBRAIC[108] = (ALGEBRAIC[100] - ALGEBRAIC[107])/CONSTANTS[147]; ALGEBRAIC[109] = CONSTANTS[157]*(STATES[41] - CONSTANTS[151]); ALGEBRAIC[110] = (ALGEBRAIC[76] - ALGEBRAIC[109])/CONSTANTS[104]; ALGEBRAIC[111] = CONSTANTS[158]*(STATES[42] - CONSTANTS[152]); ALGEBRAIC[112] = (ALGEBRAIC[80] - ALGEBRAIC[111])/CONSTANTS[111]; ALGEBRAIC[114] = (ALGEBRAIC[103] - ALGEBRAIC[14])/CONSTANTS[160]; ALGEBRAIC[115] = (ALGEBRAIC[105] - ALGEBRAIC[103])/CONSTANTS[161]; ALGEBRAIC[116] = (ALGEBRAIC[107] - ALGEBRAIC[105])/CONSTANTS[162]; ALGEBRAIC[117] = (ALGEBRAIC[109] - ALGEBRAIC[107])/CONSTANTS[163]; ALGEBRAIC[118] = (ALGEBRAIC[111] - ALGEBRAIC[109])/CONSTANTS[164]; ALGEBRAIC[113] = CONSTANTS[159]*(STATES[43] - CONSTANTS[153]); ALGEBRAIC[119] = (ALGEBRAIC[113] - ALGEBRAIC[111])/CONSTANTS[165]; ALGEBRAIC[120] = CONSTANTS[168]*(STATES[44] - CONSTANTS[167]); ALGEBRAIC[121] = (ALGEBRAIC[53] - ALGEBRAIC[120])/CONSTANTS[166]; ALGEBRAIC[122] = (ALGEBRAIC[120] - ALGEBRAIC[113])/CONSTANTS[169]; ALGEBRAIC[123] = CONSTANTS[172]*(STATES[45] - CONSTANTS[171]); ALGEBRAIC[124] = (ALGEBRAIC[53] - ALGEBRAIC[123])/CONSTANTS[170]; ALGEBRAIC[125] = (ALGEBRAIC[123] - ALGEBRAIC[113])/CONSTANTS[173]; 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[3]+STATES[1]+STATES[2]+STATES[0]; ALGEBRAIC[5] = STATES[5]+STATES[4]; ALGEBRAIC[7] = STATES[7]+STATES[8]+STATES[9]+STATES[10]; ALGEBRAIC[9] = STATES[21]+STATES[22]+STATES[23]+STATES[24]; ALGEBRAIC[11] = STATES[6]+STATES[13]+STATES[14]+STATES[15]+STATES[16]+STATES[17]+STATES[18]+STATES[19]+STATES[20]; ALGEBRAIC[13] = STATES[27]+STATES[28]+STATES[29]+STATES[30]; ALGEBRAIC[15] = STATES[25]+STATES[26]+STATES[31]+STATES[32]+STATES[33]+STATES[34]+STATES[35]+STATES[36]+STATES[37]; ALGEBRAIC[17] = STATES[38]+STATES[39]+STATES[40]+STATES[41]+STATES[42]+STATES[43]; ALGEBRAIC[19] = STATES[11]+STATES[12]+STATES[44]+STATES[45]; ALGEBRAIC[22] = ALGEBRAIC[3]+ALGEBRAIC[5]+(ALGEBRAIC[7])+ALGEBRAIC[9]+ALGEBRAIC[11]+ALGEBRAIC[13]+ALGEBRAIC[15]+ALGEBRAIC[17]+ALGEBRAIC[19]; }