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