/* There are a total of 132 entries in the algebraic variable array. There are a total of 94 entries in each of the rate and state variable arrays. There are a total of 82 entries in the constant variable array. */ /* * VOI is time in component environment (minute). * STATES[0] is c1 in component c1 (molecules). * ALGEBRAIC[6] is v1 in component v1 (flux). * STATES[1] is c2 in component c2 (molecules). * ALGEBRAIC[10] is v6 in component v6 (flux). * CONSTANTS[80] is v13 in component v13 (flux). * STATES[2] is c3 in component c3 (molecules). * ALGEBRAIC[7] is v2 in component v2 (flux). * STATES[3] is c4 in component c4 (molecules). * ALGEBRAIC[8] is v3 in component v3 (flux). * STATES[4] is c5 in component c5 (molecules). * ALGEBRAIC[11] is v7 in component v7 (flux). * ALGEBRAIC[12] is v8 in component v8 (flux). * STATES[5] is c6 in component c6 (molecules). * ALGEBRAIC[14] is v10 in component v10 (flux). * ALGEBRAIC[63] is v60 in component v60 (flux). * STATES[6] is c7 in component c7 (molecules). * ALGEBRAIC[9] is v4 in component v4 (flux). * ALGEBRAIC[124] is v5 in component v5 (flux). * STATES[7] is c8 in component c8 (molecules). * ALGEBRAIC[16] is v12 in component v12 (flux). * ALGEBRAIC[17] is v14 in component v14 (flux). * ALGEBRAIC[65] is v62 in component v62 (flux). * STATES[8] is c9 in component c9 (molecules). * ALGEBRAIC[18] is v15 in component v15 (flux). * ALGEBRAIC[125] is v107 in component v107 (flux). * ALGEBRAIC[126] is v110 in component v110 (flux). * ALGEBRAIC[127] is v113 in component v113 (flux). * ALGEBRAIC[128] is v116 in component v116 (flux). * ALGEBRAIC[129] is v119 in component v119 (flux). * ALGEBRAIC[130] is v122 in component v122 (flux). * ALGEBRAIC[131] is v125 in component v125 (flux). * STATES[9] is c10 in component c10 (molecules). * ALGEBRAIC[15] is v11 in component v11 (flux). * STATES[10] is c11 in component c11 (molecules). * STATES[11] is c12 in component c12 (molecules). * ALGEBRAIC[109] is v106 in component v106 (flux). * ALGEBRAIC[111] is v109 in component v109 (flux). * ALGEBRAIC[113] is v112 in component v112 (flux). * ALGEBRAIC[115] is v115 in component v115 (flux). * ALGEBRAIC[117] is v118 in component v118 (flux). * ALGEBRAIC[119] is v121 in component v121 (flux). * ALGEBRAIC[121] is v124 in component v124 (flux). * STATES[12] is c13 in component c13 (molecules). * ALGEBRAIC[64] is v61 in component v61 (flux). * STATES[13] is c14 in component c14 (molecules). * STATES[14] is c15 in component c15 (molecules). * ALGEBRAIC[19] is v16 in component v16 (flux). * ALGEBRAIC[25] is v22 in component v22 (flux). * ALGEBRAIC[35] is v32 in component v32 (flux). * ALGEBRAIC[37] is v34 in component v34 (flux). * ALGEBRAIC[40] is v37 in component v37 (flux). * ALGEBRAIC[42] is v39 in component v39 (flux). * ALGEBRAIC[105] is v102 in component v102 (flux). * STATES[15] is c16 in component c16 (molecules). * STATES[16] is c17 in component c17 (molecules). * ALGEBRAIC[66] is v63 in component v63 (flux). * ALGEBRAIC[72] is v69 in component v69 (flux). * ALGEBRAIC[82] is v79 in component v79 (flux). * ALGEBRAIC[83] is v80 in component v80 (flux). * ALGEBRAIC[84] is v81 in component v81 (flux). * ALGEBRAIC[85] is v82 in component v82 (flux). * STATES[17] is c18 in component c18 (molecules). * ALGEBRAIC[13] is v9 in component v9 (flux). * ALGEBRAIC[67] is v64 in component v64 (flux). * STATES[18] is c19 in component c19 (molecules). * ALGEBRAIC[68] is v65 in component v65 (flux). * ALGEBRAIC[69] is v66 in component v66 (flux). * ALGEBRAIC[70] is v67 in component v67 (flux). * ALGEBRAIC[71] is v68 in component v68 (flux). * ALGEBRAIC[108] is v105 in component v105 (flux). * STATES[19] is c20 in component c20 (molecules). * ALGEBRAIC[110] is v108 in component v108 (flux). * STATES[20] is c21 in component c21 (molecules). * ALGEBRAIC[112] is v111 in component v111 (flux). * STATES[21] is c22 in component c22 (molecules). * ALGEBRAIC[27] is v24 in component v24 (flux). * ALGEBRAIC[38] is v35 in component v35 (flux). * ALGEBRAIC[41] is v38 in component v38 (flux). * ALGEBRAIC[74] is v71 in component v71 (flux). * STATES[22] is c23 in component c23 (molecules). * ALGEBRAIC[20] is v17 in component v17 (flux). * STATES[23] is c24 in component c24 (molecules). * ALGEBRAIC[75] is v72 in component v72 (flux). * ALGEBRAIC[28] is v25 in component v25 (flux). * ALGEBRAIC[43] is v40 in component v40 (flux). * STATES[24] is c25 in component c25 (molecules). * ALGEBRAIC[21] is v18 in component v18 (flux). * ALGEBRAIC[22] is v19 in component v19 (flux). * ALGEBRAIC[23] is v20 in component v20 (flux). * ALGEBRAIC[24] is v21 in component v21 (flux). * STATES[25] is c26 in component c26 (molecules). * ALGEBRAIC[29] is v26 in component v26 (flux). * ALGEBRAIC[34] is v31 in component v31 (flux). * ALGEBRAIC[81] is v78 in component v78 (flux). * ALGEBRAIC[76] is v73 in component v73 (flux). * STATES[26] is c27 in component c27 (molecules). * STATES[27] is c28 in component c28 (molecules). * ALGEBRAIC[30] is v27 in component v27 (flux). * ALGEBRAIC[31] is v28 in component v28 (flux). * STATES[28] is c29 in component c29 (molecules). * STATES[29] is c30 in component c30 (molecules). * ALGEBRAIC[36] is v33 in component v33 (flux). * ALGEBRAIC[44] is v41 in component v41 (flux). * ALGEBRAIC[86] is v83 in component v83 (flux). * STATES[30] is c31 in component c31 (molecules). * ALGEBRAIC[39] is v36 in component v36 (flux). * STATES[31] is c32 in component c32 (molecules). * ALGEBRAIC[26] is v23 in component v23 (flux). * ALGEBRAIC[106] is v103 in component v103 (flux). * STATES[32] is c33 in component c33 (molecules). * ALGEBRAIC[107] is v104 in component v104 (flux). * STATES[33] is c34 in component c34 (molecules). * ALGEBRAIC[114] is v114 in component v114 (flux). * STATES[34] is c35 in component c35 (molecules). * ALGEBRAIC[33] is v30 in component v30 (flux). * ALGEBRAIC[116] is v117 in component v117 (flux). * STATES[35] is c36 in component c36 (molecules). * ALGEBRAIC[118] is v120 in component v120 (flux). * STATES[36] is c37 in component c37 (molecules). * ALGEBRAIC[120] is v123 in component v123 (flux). * STATES[37] is c38 in component c38 (molecules). * STATES[38] is c39 in component c39 (molecules). * STATES[39] is c40 in component c40 (molecules). * STATES[40] is c41 in component c41 (molecules). * ALGEBRAIC[46] is v43 in component v43 (flux). * ALGEBRAIC[78] is v75 in component v75 (flux). * ALGEBRAIC[88] is v85 in component v85 (flux). * STATES[41] is c42 in component c42 (molecules). * ALGEBRAIC[32] is v29 in component v29 (flux). * STATES[42] is c43 in component c43 (molecules). * STATES[43] is c44 in component c44 (molecules). * ALGEBRAIC[45] is v42 in component v42 (flux). * ALGEBRAIC[87] is v84 in component v84 (flux). * STATES[44] is c45 in component c45 (molecules). * ALGEBRAIC[47] is v44 in component v44 (flux). * ALGEBRAIC[48] is v45 in component v45 (flux). * ALGEBRAIC[49] is v46 in component v46 (flux). * ALGEBRAIC[50] is v47 in component v47 (flux). * STATES[45] is c46 in component c46 (molecules). * STATES[46] is c47 in component c47 (molecules). * ALGEBRAIC[54] is v51 in component v51 (flux). * ALGEBRAIC[89] is v86 in component v86 (flux). * ALGEBRAIC[96] is v93 in component v93 (flux). * STATES[47] is c48 in component c48 (molecules). * STATES[48] is c49 in component c49 (molecules). * ALGEBRAIC[52] is v49 in component v49 (flux). * ALGEBRAIC[53] is v50 in component v50 (flux). * STATES[49] is c50 in component c50 (molecules). * STATES[50] is c51 in component c51 (molecules). * ALGEBRAIC[51] is v48 in component v48 (flux). * ALGEBRAIC[55] is v52 in component v52 (flux). * ALGEBRAIC[56] is v53 in component v53 (flux). * ALGEBRAIC[57] is v54 in component v54 (flux). * ALGEBRAIC[58] is v55 in component v55 (flux). * STATES[51] is c52 in component c52 (molecules). * STATES[52] is c53 in component c53 (molecules). * ALGEBRAIC[93] is v90 in component v90 (flux). * ALGEBRAIC[94] is v91 in component v91 (flux). * ALGEBRAIC[95] is v92 in component v92 (flux). * STATES[53] is c54 in component c54 (molecules). * STATES[54] is c55 in component c55 (molecules). * ALGEBRAIC[62] is v59 in component v59 (flux). * ALGEBRAIC[97] is v94 in component v94 (flux). * ALGEBRAIC[104] is v101 in component v101 (flux). * STATES[55] is c56 in component c56 (molecules). * STATES[56] is c57 in component c57 (molecules). * ALGEBRAIC[60] is v57 in component v57 (flux). * ALGEBRAIC[61] is v58 in component v58 (flux). * STATES[57] is c58 in component c58 (molecules). * STATES[58] is c59 in component c59 (molecules). * ALGEBRAIC[59] is v56 in component v56 (flux). * STATES[59] is c60 in component c60 (molecules). * ALGEBRAIC[101] is v98 in component v98 (flux). * ALGEBRAIC[102] is v99 in component v99 (flux). * ALGEBRAIC[103] is v100 in component v100 (flux). * STATES[60] is c61 in component c61 (molecules). * STATES[61] is c62 in component c62 (molecules). * STATES[62] is c63 in component c63 (molecules). * ALGEBRAIC[73] is v70 in component v70 (flux). * STATES[63] is c64 in component c64 (molecules). * STATES[64] is c65 in component c65 (molecules). * STATES[65] is c66 in component c66 (molecules). * ALGEBRAIC[77] is v74 in component v74 (flux). * ALGEBRAIC[80] is v77 in component v77 (flux). * STATES[66] is c67 in component c67 (molecules). * STATES[67] is c68 in component c68 (molecules). * STATES[68] is c69 in component c69 (molecules). * STATES[69] is c70 in component c70 (molecules). * ALGEBRAIC[79] is v76 in component v76 (flux). * STATES[70] is c71 in component c71 (molecules). * STATES[71] is c72 in component c72 (molecules). * ALGEBRAIC[90] is v87 in component v87 (flux). * ALGEBRAIC[91] is v88 in component v88 (flux). * ALGEBRAIC[92] is v89 in component v89 (flux). * STATES[72] is c73 in component c73 (molecules). * STATES[73] is c74 in component c74 (molecules). * STATES[74] is c75 in component c75 (molecules). * STATES[75] is c76 in component c76 (molecules). * STATES[76] is c77 in component c77 (molecules). * ALGEBRAIC[98] is v95 in component v95 (flux). * ALGEBRAIC[99] is v96 in component v96 (flux). * ALGEBRAIC[100] is v97 in component v97 (flux). * STATES[77] is c78 in component c78 (molecules). * STATES[78] is c79 in component c79 (molecules). * STATES[79] is c80 in component c80 (molecules). * STATES[80] is c81 in component c81 (molecules). * STATES[81] is c82 in component c82 (molecules). * STATES[82] is c83 in component c83 (molecules). * STATES[83] is c84 in component c84 (molecules). * STATES[84] is c85 in component c85 (molecules). * STATES[85] is c86 in component c86 (molecules). * STATES[86] is c87 in component c87 (molecules). * STATES[87] is c88 in component c88 (molecules). * STATES[88] is c89 in component c89 (molecules). * STATES[89] is c90 in component c90 (molecules). * STATES[90] is c91 in component c91 (molecules). * STATES[91] is c92 in component c92 (molecules). * STATES[92] is c93 in component c93 (molecules). * STATES[93] is c94 in component c94 (molecules). * ALGEBRAIC[0] is Raf_activation in component Raf_activation (molecules). * ALGEBRAIC[1] is MEKPP_total in component MEKPP_total (molecules). * ALGEBRAIC[2] is ERKPP_total in component ERKPP_total (molecules). * ALGEBRAIC[3] is Ras_GTP_total in component Ras_GTP_total (molecules). * ALGEBRAIC[4] is ShcP_total in component ShcP_total (molecules). * ALGEBRAIC[5] is EGF_EGFR_totalx2 in component EGF_EGFR_totalx2 (molecules). * CONSTANTS[0] is k1 in component reaction_rates (second_order_rate_constant). * CONSTANTS[1] is kd1 in component reaction_rates (first_order_rate_constant). * CONSTANTS[2] is k2 in component reaction_rates (second_order_rate_constant). * CONSTANTS[3] is kd2 in component reaction_rates (first_order_rate_constant). * CONSTANTS[4] is k3 in component reaction_rates (first_order_rate_constant). * CONSTANTS[5] is kd3 in component reaction_rates (first_order_rate_constant). * CONSTANTS[6] is k4 in component reaction_rates (second_order_rate_constant). * CONSTANTS[7] is kd4 in component reaction_rates (first_order_rate_constant). * ALGEBRAIC[123] is k5 in component reaction_rates (first_order_rate_constant). * CONSTANTS[8] is k6 in component reaction_rates (first_order_rate_constant). * CONSTANTS[9] is kd6 in component reaction_rates (first_order_rate_constant). * CONSTANTS[10] is k8 in component reaction_rates (second_order_rate_constant). * CONSTANTS[11] is kd8 in component reaction_rates (first_order_rate_constant). * CONSTANTS[12] is k10b in component reaction_rates (second_order_rate_constant). * CONSTANTS[13] is kd10 in component reaction_rates (first_order_rate_constant). * CONSTANTS[14] is k13 in component reaction_rates (first_order_rate_constant). * CONSTANTS[15] is k15 in component reaction_rates (first_order_rate_constant). * CONSTANTS[16] is k16 in component reaction_rates (second_order_rate_constant). * CONSTANTS[17] is kd16 in component reaction_rates (first_order_rate_constant). * CONSTANTS[18] is k17 in component reaction_rates (second_order_rate_constant). * CONSTANTS[19] is kd17 in component reaction_rates (first_order_rate_constant). * CONSTANTS[20] is k18 in component reaction_rates (second_order_rate_constant). * CONSTANTS[21] is kd18 in component reaction_rates (first_order_rate_constant). * CONSTANTS[22] is k19 in component reaction_rates (first_order_rate_constant). * CONSTANTS[23] is kd19 in component reaction_rates (second_order_rate_constant). * CONSTANTS[24] is k20 in component reaction_rates (second_order_rate_constant). * CONSTANTS[25] is kd20 in component reaction_rates (first_order_rate_constant). * CONSTANTS[26] is k21 in component reaction_rates (first_order_rate_constant). * CONSTANTS[27] is kd21 in component reaction_rates (second_order_rate_constant). * CONSTANTS[28] is k22 in component reaction_rates (second_order_rate_constant). * CONSTANTS[29] is kd22 in component reaction_rates (first_order_rate_constant). * CONSTANTS[30] is k23 in component reaction_rates (first_order_rate_constant). * CONSTANTS[31] is kd23 in component reaction_rates (first_order_rate_constant). * CONSTANTS[32] is kd24 in component reaction_rates (first_order_rate_constant). * CONSTANTS[33] is k25 in component reaction_rates (second_order_rate_constant). * CONSTANTS[34] is kd25 in component reaction_rates (first_order_rate_constant). * CONSTANTS[35] is k28 in component reaction_rates (second_order_rate_constant). * CONSTANTS[36] is kd28 in component reaction_rates (first_order_rate_constant). * CONSTANTS[37] is k29 in component reaction_rates (first_order_rate_constant). * CONSTANTS[38] is kd29 in component reaction_rates (second_order_rate_constant). * CONSTANTS[39] is k32 in component reaction_rates (first_order_rate_constant). * CONSTANTS[40] is kd32 in component reaction_rates (second_order_rate_constant). * CONSTANTS[41] is k33 in component reaction_rates (first_order_rate_constant). * CONSTANTS[42] is kd33 in component reaction_rates (second_order_rate_constant). * CONSTANTS[43] is k34 in component reaction_rates (first_order_rate_constant). * CONSTANTS[44] is kd34 in component reaction_rates (second_order_rate_constant). * CONSTANTS[45] is k35 in component reaction_rates (first_order_rate_constant). * CONSTANTS[46] is kd35 in component reaction_rates (second_order_rate_constant). * CONSTANTS[47] is Km36 in component reaction_rates (molecules). * CONSTANTS[48] is Vm36 in component reaction_rates (flux). * CONSTANTS[49] is k37 in component reaction_rates (first_order_rate_constant). * CONSTANTS[50] is kd37 in component reaction_rates (second_order_rate_constant). * CONSTANTS[51] is k40 in component reaction_rates (second_order_rate_constant). * CONSTANTS[52] is kd40 in component reaction_rates (first_order_rate_constant). * CONSTANTS[53] is k41 in component reaction_rates (second_order_rate_constant). * CONSTANTS[54] is kd41 in component reaction_rates (first_order_rate_constant). * CONSTANTS[55] is k42 in component reaction_rates (second_order_rate_constant). * CONSTANTS[56] is kd42 in component reaction_rates (first_order_rate_constant). * CONSTANTS[57] is k43 in component reaction_rates (first_order_rate_constant). * CONSTANTS[58] is k44 in component reaction_rates (second_order_rate_constant). * CONSTANTS[59] is kd44 in component reaction_rates (first_order_rate_constant). * CONSTANTS[60] is k45 in component reaction_rates (first_order_rate_constant). * CONSTANTS[61] is k47 in component reaction_rates (first_order_rate_constant). * CONSTANTS[62] is k48 in component reaction_rates (second_order_rate_constant). * CONSTANTS[63] is kd48 in component reaction_rates (first_order_rate_constant). * CONSTANTS[64] is k49 in component reaction_rates (first_order_rate_constant). * CONSTANTS[65] is k50 in component reaction_rates (second_order_rate_constant). * CONSTANTS[66] is kd50 in component reaction_rates (first_order_rate_constant). * CONSTANTS[67] is k52 in component reaction_rates (second_order_rate_constant). * CONSTANTS[68] is kd52 in component reaction_rates (first_order_rate_constant). * CONSTANTS[69] is k53 in component reaction_rates (first_order_rate_constant). * CONSTANTS[70] is k55 in component reaction_rates (first_order_rate_constant). * CONSTANTS[71] is k56 in component reaction_rates (second_order_rate_constant). * CONSTANTS[72] is kd56 in component reaction_rates (first_order_rate_constant). * CONSTANTS[73] is k57 in component reaction_rates (first_order_rate_constant). * CONSTANTS[74] is k58 in component reaction_rates (second_order_rate_constant). * CONSTANTS[75] is kd58 in component reaction_rates (first_order_rate_constant). * CONSTANTS[76] is k59 in component reaction_rates (first_order_rate_constant). * CONSTANTS[77] is k60 in component reaction_rates (first_order_rate_constant). * CONSTANTS[78] is k61 in component reaction_rates (first_order_rate_constant). * ALGEBRAIC[122] is C in component C (flux). * CONSTANTS[79] is RT in component C (molecules). * RATES[0] is d/dt c1 in component c1 (molecules). * RATES[1] is d/dt c2 in component c2 (molecules). * RATES[2] is d/dt c3 in component c3 (molecules). * RATES[3] is d/dt c4 in component c4 (molecules). * RATES[4] is d/dt c5 in component c5 (molecules). * RATES[5] is d/dt c6 in component c6 (molecules). * RATES[6] is d/dt c7 in component c7 (molecules). * RATES[7] is d/dt c8 in component c8 (molecules). * RATES[8] is d/dt c9 in component c9 (molecules). * RATES[9] is d/dt c10 in component c10 (molecules). * RATES[10] is d/dt c11 in component c11 (molecules). * RATES[11] is d/dt c12 in component c12 (molecules). * RATES[12] is d/dt c13 in component c13 (molecules). * RATES[13] is d/dt c14 in component c14 (molecules). * RATES[14] is d/dt c15 in component c15 (molecules). * RATES[15] is d/dt c16 in component c16 (molecules). * RATES[16] is d/dt c17 in component c17 (molecules). * RATES[17] is d/dt c18 in component c18 (molecules). * RATES[18] is d/dt c19 in component c19 (molecules). * RATES[19] is d/dt c20 in component c20 (molecules). * RATES[20] is d/dt c21 in component c21 (molecules). * RATES[21] is d/dt c22 in component c22 (molecules). * RATES[22] is d/dt c23 in component c23 (molecules). * RATES[23] is d/dt c24 in component c24 (molecules). * RATES[24] is d/dt c25 in component c25 (molecules). * RATES[25] is d/dt c26 in component c26 (molecules). * RATES[26] is d/dt c27 in component c27 (molecules). * RATES[27] is d/dt c28 in component c28 (molecules). * RATES[28] is d/dt c29 in component c29 (molecules). * RATES[29] is d/dt c30 in component c30 (molecules). * RATES[30] is d/dt c31 in component c31 (molecules). * RATES[31] is d/dt c32 in component c32 (molecules). * RATES[32] is d/dt c33 in component c33 (molecules). * RATES[33] is d/dt c34 in component c34 (molecules). * RATES[34] is d/dt c35 in component c35 (molecules). * RATES[35] is d/dt c36 in component c36 (molecules). * RATES[36] is d/dt c37 in component c37 (molecules). * RATES[37] is d/dt c38 in component c38 (molecules). * RATES[38] is d/dt c39 in component c39 (molecules). * RATES[39] is d/dt c40 in component c40 (molecules). * RATES[40] is d/dt c41 in component c41 (molecules). * RATES[41] is d/dt c42 in component c42 (molecules). * RATES[42] is d/dt c43 in component c43 (molecules). * RATES[43] is d/dt c44 in component c44 (molecules). * RATES[44] is d/dt c45 in component c45 (molecules). * RATES[45] is d/dt c46 in component c46 (molecules). * RATES[46] is d/dt c47 in component c47 (molecules). * RATES[47] is d/dt c48 in component c48 (molecules). * RATES[48] is d/dt c49 in component c49 (molecules). * RATES[49] is d/dt c50 in component c50 (molecules). * RATES[50] is d/dt c51 in component c51 (molecules). * RATES[51] is d/dt c52 in component c52 (molecules). * RATES[52] is d/dt c53 in component c53 (molecules). * RATES[53] is d/dt c54 in component c54 (molecules). * RATES[54] is d/dt c55 in component c55 (molecules). * RATES[55] is d/dt c56 in component c56 (molecules). * RATES[56] is d/dt c57 in component c57 (molecules). * RATES[57] is d/dt c58 in component c58 (molecules). * RATES[58] is d/dt c59 in component c59 (molecules). * RATES[59] is d/dt c60 in component c60 (molecules). * RATES[60] is d/dt c61 in component c61 (molecules). * RATES[61] is d/dt c62 in component c62 (molecules). * RATES[62] is d/dt c63 in component c63 (molecules). * RATES[63] is d/dt c64 in component c64 (molecules). * RATES[64] is d/dt c65 in component c65 (molecules). * RATES[65] is d/dt c66 in component c66 (molecules). * RATES[66] is d/dt c67 in component c67 (molecules). * RATES[67] is d/dt c68 in component c68 (molecules). * RATES[68] is d/dt c69 in component c69 (molecules). * RATES[69] is d/dt c70 in component c70 (molecules). * RATES[70] is d/dt c71 in component c71 (molecules). * RATES[71] is d/dt c72 in component c72 (molecules). * RATES[72] is d/dt c73 in component c73 (molecules). * RATES[73] is d/dt c74 in component c74 (molecules). * RATES[74] is d/dt c75 in component c75 (molecules). * RATES[75] is d/dt c76 in component c76 (molecules). * RATES[76] is d/dt c77 in component c77 (molecules). * RATES[77] is d/dt c78 in component c78 (molecules). * RATES[78] is d/dt c79 in component c79 (molecules). * RATES[79] is d/dt c80 in component c80 (molecules). * RATES[80] is d/dt c81 in component c81 (molecules). * RATES[81] is d/dt c82 in component c82 (molecules). * RATES[82] is d/dt c83 in component c83 (molecules). * RATES[83] is d/dt c84 in component c84 (molecules). * RATES[84] is d/dt c85 in component c85 (molecules). * RATES[85] is d/dt c86 in component c86 (molecules). * RATES[86] is d/dt c87 in component c87 (molecules). * RATES[87] is d/dt c88 in component c88 (molecules). * RATES[88] is d/dt c89 in component c89 (molecules). * RATES[89] is d/dt c90 in component c90 (molecules). * RATES[90] is d/dt c91 in component c91 (molecules). * RATES[91] is d/dt c92 in component c92 (molecules). * RATES[92] is d/dt c93 in component c93 (molecules). * RATES[93] is d/dt c94 in component c94 (molecules). * There are a total of 2 condition variables. */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { STATES[0] = 4962; STATES[1] = 5e4; STATES[2] = 0; STATES[3] = 0; STATES[4] = 0; STATES[5] = 0; STATES[6] = 0; STATES[7] = 0; STATES[8] = 0; STATES[9] = 0; STATES[10] = 0; STATES[11] = 8.1e4; STATES[12] = 0; STATES[13] = 1.2e4; STATES[14] = 0; STATES[15] = 0; STATES[16] = 0; STATES[17] = 0; STATES[18] = 0; STATES[19] = 0; STATES[20] = 0; STATES[21] = 1.1e4; STATES[22] = 0; STATES[23] = 2.63e4; STATES[24] = 0; STATES[25] = 7.2e4; STATES[26] = 0; STATES[27] = 0; STATES[28] = 0; STATES[29] = 4e4; STATES[30] = 1.01e5; STATES[31] = 0; STATES[32] = 0; STATES[33] = 0; STATES[34] = 0; STATES[35] = 0; STATES[36] = 0; STATES[37] = 0; STATES[38] = 0; STATES[39] = 0; STATES[40] = 4e4; STATES[41] = 0; STATES[42] = 0; STATES[43] = 4e4; STATES[44] = 0; STATES[45] = 0; STATES[46] = 2.2e7; STATES[47] = 0; STATES[48] = 0; STATES[49] = 0; STATES[50] = 0; STATES[51] = 0; STATES[52] = 4e4; STATES[53] = 0; STATES[54] = 2.1e7; STATES[55] = 0; STATES[56] = 0; STATES[57] = 0; STATES[58] = 0; STATES[59] = 1e7; STATES[60] = 0; STATES[61] = 0; STATES[62] = 0; STATES[63] = 0; STATES[64] = 0; STATES[65] = 0; STATES[66] = 0; STATES[67] = 0; STATES[68] = 0; STATES[69] = 0; STATES[70] = 0; STATES[71] = 0; STATES[72] = 0; STATES[73] = 0; STATES[74] = 0; STATES[75] = 0; STATES[76] = 0; STATES[77] = 0; STATES[78] = 0; STATES[79] = 0; STATES[80] = 0; STATES[81] = 0; STATES[82] = 0; STATES[83] = 0; STATES[84] = 0; STATES[85] = 0; STATES[86] = 0; STATES[87] = 0; STATES[88] = 0; STATES[89] = 0; STATES[90] = 0; STATES[91] = 0; STATES[92] = 0; STATES[93] = 0; CONSTANTS[0] = 3e-3; CONSTANTS[1] = 0.228; CONSTANTS[2] = 1e-3; CONSTANTS[3] = 6; CONSTANTS[4] = 60; CONSTANTS[5] = 0.6; CONSTANTS[6] = 1.038e-5; CONSTANTS[7] = 0.0996; CONSTANTS[8] = 3e-3; CONSTANTS[9] = 0.3; CONSTANTS[10] = 1e-4; CONSTANTS[11] = 12; CONSTANTS[12] = 3.25581; CONSTANTS[13] = 0.66; CONSTANTS[14] = 130.2; CONSTANTS[15] = 600000; CONSTANTS[16] = 0.001; CONSTANTS[17] = 16.5; CONSTANTS[18] = 0.001; CONSTANTS[19] = 3.6; CONSTANTS[20] = 0.0015; CONSTANTS[21] = 78; CONSTANTS[22] = 30; CONSTANTS[23] = 1e-5; CONSTANTS[24] = 2.1e-4; CONSTANTS[25] = 24; CONSTANTS[26] = 1.38; CONSTANTS[27] = 2.2e-5; CONSTANTS[28] = 0.0021; CONSTANTS[29] = 6; CONSTANTS[30] = 360; CONSTANTS[31] = 36; CONSTANTS[32] = 33; CONSTANTS[33] = 0.001; CONSTANTS[34] = 1.284; CONSTANTS[35] = 1e-4; CONSTANTS[36] = 0.318; CONSTANTS[37] = 60; CONSTANTS[38] = 7e-5; CONSTANTS[39] = 6; CONSTANTS[40] = 2.4e-5; CONSTANTS[41] = 12; CONSTANTS[42] = 0.0021; CONSTANTS[43] = 1.8; CONSTANTS[44] = 4.5e-4; CONSTANTS[45] = 0.09; CONSTANTS[46] = 4.5e-4; CONSTANTS[47] = 2e14; CONSTANTS[48] = 61200; CONSTANTS[49] = 18; CONSTANTS[50] = 9e-5; CONSTANTS[51] = 0.003; CONSTANTS[52] = 3.84; CONSTANTS[53] = 0.003; CONSTANTS[54] = 2.574; CONSTANTS[55] = 0.007; CONSTANTS[56] = 12; CONSTANTS[57] = 60; CONSTANTS[58] = 0.00111; CONSTANTS[59] = 1.0998; CONSTANTS[60] = 210; CONSTANTS[61] = 174; CONSTANTS[62] = 0.00143; CONSTANTS[63] = 48; CONSTANTS[64] = 3.48; CONSTANTS[65] = 2.5e-5; CONSTANTS[66] = 30; CONSTANTS[67] = 0.00534; CONSTANTS[68] = 1.98; CONSTANTS[69] = 960; CONSTANTS[70] = 342; CONSTANTS[71] = 0.00145; CONSTANTS[72] = 36; CONSTANTS[73] = 16.2; CONSTANTS[74] = 5e-4; CONSTANTS[75] = 30; CONSTANTS[76] = 18; CONSTANTS[77] = 0.04002; CONSTANTS[78] = 0.01002; CONSTANTS[79] = 50000; CONSTANTS[80] = CONSTANTS[14]*1.00000; CONSTANTS[81] = 0.00000; RATES[1] = 0.1001; RATES[2] = 0.1001; RATES[3] = 0.1001; RATES[4] = 0.1001; RATES[5] = 0.1001; RATES[6] = 0.1001; RATES[7] = 0.1001; RATES[8] = 0.1001; RATES[9] = 0.1001; RATES[10] = 0.1001; RATES[11] = 0.1001; RATES[12] = 0.1001; RATES[13] = 0.1001; RATES[14] = 0.1001; RATES[15] = 0.1001; RATES[16] = 0.1001; RATES[17] = 0.1001; RATES[18] = 0.1001; RATES[19] = 0.1001; RATES[20] = 0.1001; RATES[21] = 0.1001; RATES[22] = 0.1001; RATES[23] = 0.1001; RATES[24] = 0.1001; RATES[25] = 0.1001; RATES[26] = 0.1001; RATES[27] = 0.1001; RATES[28] = 0.1001; RATES[29] = 0.1001; RATES[30] = 0.1001; RATES[31] = 0.1001; RATES[32] = 0.1001; RATES[33] = 0.1001; RATES[34] = 0.1001; RATES[35] = 0.1001; RATES[36] = 0.1001; RATES[37] = 0.1001; RATES[38] = 0.1001; RATES[39] = 0.1001; RATES[40] = 0.1001; RATES[41] = 0.1001; RATES[42] = 0.1001; RATES[43] = 0.1001; RATES[44] = 0.1001; RATES[45] = 0.1001; RATES[46] = 0.1001; RATES[47] = 0.1001; RATES[48] = 0.1001; RATES[49] = 0.1001; RATES[50] = 0.1001; RATES[51] = 0.1001; RATES[52] = 0.1001; RATES[53] = 0.1001; RATES[54] = 0.1001; RATES[55] = 0.1001; RATES[56] = 0.1001; RATES[57] = 0.1001; RATES[58] = 0.1001; RATES[59] = 0.1001; RATES[60] = 0.1001; RATES[61] = 0.1001; RATES[62] = 0.1001; RATES[63] = 0.1001; RATES[64] = 0.1001; RATES[65] = 0.1001; RATES[66] = 0.1001; RATES[67] = 0.1001; RATES[68] = 0.1001; RATES[69] = 0.1001; RATES[70] = 0.1001; RATES[71] = 0.1001; RATES[72] = 0.1001; RATES[73] = 0.1001; RATES[74] = 0.1001; RATES[75] = 0.1001; RATES[76] = 0.1001; RATES[77] = 0.1001; RATES[78] = 0.1001; RATES[79] = 0.1001; RATES[80] = 0.1001; RATES[81] = 0.1001; RATES[82] = 0.1001; RATES[83] = 0.1001; RATES[84] = 0.1001; RATES[85] = 0.1001; RATES[86] = 0.1001; RATES[87] = 0.1001; RATES[88] = 0.1001; RATES[89] = 0.1001; RATES[90] = 0.1001; RATES[91] = 0.1001; RATES[92] = 0.1001; RATES[93] = 0.1001; } void computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { resid[0] = RATES[1] - CONSTANTS[80] - (ALGEBRAIC[6]+ALGEBRAIC[10]); resid[1] = RATES[2] - ALGEBRAIC[6] - 2.00000*ALGEBRAIC[7]; resid[2] = RATES[3] - ALGEBRAIC[7] - ALGEBRAIC[8]; resid[3] = RATES[4] - ALGEBRAIC[8] - (ALGEBRAIC[11]+ALGEBRAIC[12]); resid[4] = RATES[5] - ALGEBRAIC[10] - (ALGEBRAIC[14]+ALGEBRAIC[63]); resid[5] = RATES[6] - ALGEBRAIC[9] - ALGEBRAIC[124]; resid[6] = RATES[7] - (ALGEBRAIC[11]+ALGEBRAIC[16]) - (ALGEBRAIC[17]+ALGEBRAIC[65]); resid[7] = RATES[8] - (ALGEBRAIC[124]+ALGEBRAIC[125]+ALGEBRAIC[126]+ALGEBRAIC[127]+ALGEBRAIC[128]+ALGEBRAIC[129]+ALGEBRAIC[130]+ALGEBRAIC[131]) - ALGEBRAIC[18]; resid[8] = RATES[9] - ALGEBRAIC[14] - 2.00000*ALGEBRAIC[15]; resid[9] = RATES[10] - ALGEBRAIC[15] - ALGEBRAIC[16]; resid[10] = RATES[11] - ALGEBRAIC[18] - (ALGEBRAIC[9]+ALGEBRAIC[109]+ALGEBRAIC[111]+ALGEBRAIC[113]+ALGEBRAIC[115]+ALGEBRAIC[117]+ALGEBRAIC[119]+ALGEBRAIC[121]); resid[11] = RATES[12] - ALGEBRAIC[64]; resid[12] = RATES[13] - - (ALGEBRAIC[12]+ALGEBRAIC[17]); resid[13] = RATES[14] - (ALGEBRAIC[12]+ALGEBRAIC[35]+ALGEBRAIC[37]+ALGEBRAIC[40]+ALGEBRAIC[42]) - (ALGEBRAIC[19]+ALGEBRAIC[25]+ALGEBRAIC[105]); resid[14] = RATES[15] - - (ALGEBRAIC[14]+ALGEBRAIC[64]); resid[15] = RATES[16] - (ALGEBRAIC[17]+ALGEBRAIC[82]+ALGEBRAIC[83]+ALGEBRAIC[84]+ALGEBRAIC[85]+ALGEBRAIC[105]) - (ALGEBRAIC[66]+ALGEBRAIC[72]); resid[16] = RATES[17] - (ALGEBRAIC[124]+ALGEBRAIC[13]+ALGEBRAIC[66]) - ALGEBRAIC[67]; resid[17] = RATES[18] - (ALGEBRAIC[69]+ALGEBRAIC[67]+ALGEBRAIC[71]+ALGEBRAIC[108]+ALGEBRAIC[125]) - (ALGEBRAIC[68]+ALGEBRAIC[70]+ALGEBRAIC[83]); resid[18] = RATES[19] - (ALGEBRAIC[68]+ALGEBRAIC[110]+ALGEBRAIC[126]) - ALGEBRAIC[69]; resid[19] = RATES[20] - (ALGEBRAIC[70]+ALGEBRAIC[112]+ALGEBRAIC[127]) - ALGEBRAIC[71]; resid[20] = RATES[21] - ALGEBRAIC[38] - (ALGEBRAIC[19]+ALGEBRAIC[27]+ALGEBRAIC[41]+ALGEBRAIC[66]+ALGEBRAIC[74]); resid[21] = RATES[22] - ALGEBRAIC[19] - (ALGEBRAIC[9]+ALGEBRAIC[13]+ALGEBRAIC[20]); resid[22] = RATES[23] - ALGEBRAIC[38] - (ALGEBRAIC[20]+ALGEBRAIC[28]+ALGEBRAIC[43]+ALGEBRAIC[67]+ALGEBRAIC[75]); resid[23] = RATES[24] - (ALGEBRAIC[20]+ALGEBRAIC[22]+ALGEBRAIC[24]) - (ALGEBRAIC[21]+ALGEBRAIC[23]+ALGEBRAIC[37]+ALGEBRAIC[108]+ALGEBRAIC[109]); resid[24] = RATES[25] - (ALGEBRAIC[24]+ALGEBRAIC[34]+ALGEBRAIC[71]+ALGEBRAIC[81]) - (ALGEBRAIC[21]+ALGEBRAIC[29]+ALGEBRAIC[68]+ALGEBRAIC[76]); resid[25] = RATES[26] - ALGEBRAIC[21] - (ALGEBRAIC[22]+ALGEBRAIC[110]+ALGEBRAIC[111]); resid[26] = RATES[27] - (ALGEBRAIC[22]+ALGEBRAIC[30]) - ALGEBRAIC[31]; resid[27] = RATES[28] - ALGEBRAIC[23] - (ALGEBRAIC[24]+ALGEBRAIC[112]+ALGEBRAIC[113]); resid[28] = RATES[29] - (ALGEBRAIC[36]+ALGEBRAIC[37]+ALGEBRAIC[83]) - (ALGEBRAIC[38]+ALGEBRAIC[44]+ALGEBRAIC[86]); resid[29] = RATES[30] - ALGEBRAIC[39] - (ALGEBRAIC[25]+ALGEBRAIC[72]); resid[30] = RATES[31] - ALGEBRAIC[25] - (ALGEBRAIC[26]+ALGEBRAIC[106]); resid[31] = RATES[32] - ALGEBRAIC[26] - (ALGEBRAIC[27]+ALGEBRAIC[40]+ALGEBRAIC[44]+ALGEBRAIC[107]); resid[32] = RATES[33] - ALGEBRAIC[27] - (ALGEBRAIC[28]+ALGEBRAIC[42]+ALGEBRAIC[114]+ALGEBRAIC[115]); resid[33] = RATES[34] - (ALGEBRAIC[28]+ALGEBRAIC[30]+ALGEBRAIC[34]+ALGEBRAIC[44]) - (ALGEBRAIC[116]+ALGEBRAIC[117]+ALGEBRAIC[29]+ALGEBRAIC[33]+ALGEBRAIC[35]); resid[34] = RATES[35] - ALGEBRAIC[29] - (ALGEBRAIC[30]+ALGEBRAIC[118]+ALGEBRAIC[119]); resid[35] = RATES[36] - ALGEBRAIC[33] - (ALGEBRAIC[34]+ALGEBRAIC[120]+ALGEBRAIC[121]); resid[36] = RATES[37] - (ALGEBRAIC[35]+ALGEBRAIC[43]+ALGEBRAIC[82]) - ALGEBRAIC[36]; resid[37] = RATES[38] - (ALGEBRAIC[41]+ALGEBRAIC[42]+ALGEBRAIC[85]) - ALGEBRAIC[43]; resid[38] = RATES[39] - (ALGEBRAIC[36]+ALGEBRAIC[40]+ALGEBRAIC[84]) - (ALGEBRAIC[39]+ALGEBRAIC[41]); resid[39] = RATES[40] - (ALGEBRAIC[46]+ALGEBRAIC[88]) - (ALGEBRAIC[31]+ALGEBRAIC[78]); resid[40] = RATES[41] - ALGEBRAIC[31] - ALGEBRAIC[32]; resid[41] = RATES[42] - ALGEBRAIC[32] - (ALGEBRAIC[23]+ALGEBRAIC[33]); resid[42] = RATES[43] - (ALGEBRAIC[46]+ALGEBRAIC[88]) - (ALGEBRAIC[45]+ALGEBRAIC[87]); resid[43] = RATES[44] - (ALGEBRAIC[32]+ALGEBRAIC[48]+ALGEBRAIC[50]) - (ALGEBRAIC[45]+ALGEBRAIC[47]+ALGEBRAIC[49]); resid[44] = RATES[45] - ALGEBRAIC[45] - ALGEBRAIC[46]; resid[45] = RATES[46] - (ALGEBRAIC[54]+ALGEBRAIC[96]) - (ALGEBRAIC[47]+ALGEBRAIC[89]); resid[46] = RATES[47] - ALGEBRAIC[47] - ALGEBRAIC[48]; resid[47] = RATES[48] - (ALGEBRAIC[48]+ALGEBRAIC[52]) - (ALGEBRAIC[49]+ALGEBRAIC[53]); resid[48] = RATES[49] - ALGEBRAIC[49] - ALGEBRAIC[50]; resid[49] = RATES[50] - (ALGEBRAIC[50]+ALGEBRAIC[56]+ALGEBRAIC[58]) - (ALGEBRAIC[51]+ALGEBRAIC[55]+ALGEBRAIC[57]); resid[50] = RATES[51] - ALGEBRAIC[51] - ALGEBRAIC[52]; resid[51] = RATES[52] - (ALGEBRAIC[52]+ALGEBRAIC[54]+ALGEBRAIC[94]+ALGEBRAIC[96]) - (ALGEBRAIC[51]+ALGEBRAIC[53]+ALGEBRAIC[93]+ALGEBRAIC[95]); resid[52] = RATES[53] - ALGEBRAIC[53] - ALGEBRAIC[54]; resid[53] = RATES[54] - (ALGEBRAIC[62]+ALGEBRAIC[104]) - (ALGEBRAIC[55]+ALGEBRAIC[97]); resid[54] = RATES[55] - ALGEBRAIC[55] - ALGEBRAIC[56]; resid[55] = RATES[56] - (ALGEBRAIC[56]+ALGEBRAIC[60]) - (ALGEBRAIC[57]+ALGEBRAIC[61]); resid[56] = RATES[57] - ALGEBRAIC[57] - ALGEBRAIC[58]; resid[57] = RATES[58] - ALGEBRAIC[58] - ALGEBRAIC[59]; resid[58] = RATES[59] - (ALGEBRAIC[60]+ALGEBRAIC[62]+ALGEBRAIC[102]+ALGEBRAIC[104]) - (ALGEBRAIC[59]+ALGEBRAIC[61]+ALGEBRAIC[101]+ALGEBRAIC[103]); resid[59] = RATES[60] - ALGEBRAIC[59] - ALGEBRAIC[60]; resid[60] = RATES[61] - ALGEBRAIC[61] - ALGEBRAIC[62]; resid[61] = RATES[62] - (ALGEBRAIC[72]+ALGEBRAIC[106]) - ALGEBRAIC[73]; resid[62] = RATES[63] - (ALGEBRAIC[73]+ALGEBRAIC[107]) - (ALGEBRAIC[74]+ALGEBRAIC[84]+ALGEBRAIC[86]); resid[63] = RATES[64] - (ALGEBRAIC[74]+ALGEBRAIC[114]+ALGEBRAIC[128]) - (ALGEBRAIC[75]+ALGEBRAIC[85]); resid[64] = RATES[65] - (ALGEBRAIC[75]+ALGEBRAIC[77]+ALGEBRAIC[81]+ALGEBRAIC[86]+ALGEBRAIC[116]+ALGEBRAIC[129]) - (ALGEBRAIC[76]+ALGEBRAIC[80]+ALGEBRAIC[82]); resid[65] = RATES[66] - (ALGEBRAIC[76]+ALGEBRAIC[118]+ALGEBRAIC[130]) - ALGEBRAIC[77]; resid[66] = RATES[67] - (ALGEBRAIC[80]+ALGEBRAIC[120]+ALGEBRAIC[131]) - ALGEBRAIC[81]; resid[67] = RATES[68] - (ALGEBRAIC[69]+ALGEBRAIC[77]) - ALGEBRAIC[78]; resid[68] = RATES[69] - ALGEBRAIC[78] - ALGEBRAIC[79]; resid[69] = RATES[70] - ALGEBRAIC[79] - (ALGEBRAIC[70]+ALGEBRAIC[80]); resid[70] = RATES[71] - (ALGEBRAIC[79]+ALGEBRAIC[90]+ALGEBRAIC[92]) - (ALGEBRAIC[87]+ALGEBRAIC[89]+ALGEBRAIC[91]); resid[71] = RATES[72] - ALGEBRAIC[87] - ALGEBRAIC[88]; resid[72] = RATES[73] - ALGEBRAIC[89] - ALGEBRAIC[90]; resid[73] = RATES[74] - (ALGEBRAIC[90]+ALGEBRAIC[94]) - (ALGEBRAIC[95]+ALGEBRAIC[91]); resid[74] = RATES[75] - ALGEBRAIC[91] - ALGEBRAIC[92]; resid[75] = RATES[76] - (ALGEBRAIC[92]+ALGEBRAIC[98]+ALGEBRAIC[100]) - (ALGEBRAIC[93]+ALGEBRAIC[97]+ALGEBRAIC[99]); resid[76] = RATES[77] - ALGEBRAIC[93] - ALGEBRAIC[94]; resid[77] = RATES[78] - ALGEBRAIC[95] - ALGEBRAIC[96]; resid[78] = RATES[79] - ALGEBRAIC[97] - ALGEBRAIC[98]; resid[79] = RATES[80] - (ALGEBRAIC[98]+ALGEBRAIC[102]) - (ALGEBRAIC[99]+ALGEBRAIC[103]); resid[80] = RATES[81] - ALGEBRAIC[99] - ALGEBRAIC[100]; resid[81] = RATES[82] - ALGEBRAIC[100] - ALGEBRAIC[101]; resid[82] = RATES[83] - ALGEBRAIC[101] - ALGEBRAIC[102]; resid[83] = RATES[84] - ALGEBRAIC[103] - ALGEBRAIC[104]; resid[84] = RATES[85] - ALGEBRAIC[63]; resid[85] = RATES[86] - ALGEBRAIC[65]; resid[86] = RATES[87] - ALGEBRAIC[109] - ALGEBRAIC[125]; resid[87] = RATES[88] - ALGEBRAIC[111] - ALGEBRAIC[126]; resid[88] = RATES[89] - ALGEBRAIC[113] - ALGEBRAIC[127]; resid[89] = RATES[90] - ALGEBRAIC[115] - ALGEBRAIC[128]; resid[90] = RATES[91] - ALGEBRAIC[117] - ALGEBRAIC[129]; resid[91] = RATES[92] - ALGEBRAIC[119] - ALGEBRAIC[130]; resid[92] = RATES[93] - ALGEBRAIC[121] - ALGEBRAIC[131]; resid[93] = RATES[0] - CONSTANTS[81]; } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[0] = STATES[44]+STATES[45]+STATES[47]+STATES[49]+STATES[71]+STATES[72]+STATES[73]+STATES[75]; ALGEBRAIC[1] = STATES[50]+STATES[76]; ALGEBRAIC[2] = STATES[58]+STATES[82]; ALGEBRAIC[3] = STATES[41]+STATES[27]+STATES[69]+STATES[68]; ALGEBRAIC[4] = STATES[32]+STATES[33]+STATES[34]+STATES[35]+STATES[36]+STATES[37]+STATES[38]+STATES[39]+STATES[63]+STATES[64]+STATES[65]+STATES[66]+STATES[67]+STATES[90]+STATES[91]+STATES[92]+STATES[93]; ALGEBRAIC[5] = 2.00000*(STATES[4]+STATES[6]+STATES[7]+STATES[10]+STATES[14]+STATES[16]+STATES[17]+STATES[18]+STATES[19]+STATES[20]+STATES[22]+STATES[24]+STATES[26]+STATES[28]+STATES[31]+STATES[32]+STATES[33]+STATES[34]+STATES[35]+STATES[36]+STATES[62]+STATES[63]+STATES[64]+STATES[65]+STATES[66]+STATES[67]+STATES[87]+STATES[88]+STATES[89]+STATES[90]+STATES[91]+STATES[92]+STATES[93]); } void computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[6] = CONSTANTS[0]*STATES[0]*STATES[1] - CONSTANTS[1]*STATES[2]; ALGEBRAIC[7] = CONSTANTS[2]*STATES[2]*STATES[2] - CONSTANTS[3]*STATES[3]; ALGEBRAIC[8] = CONSTANTS[4]*STATES[3] - CONSTANTS[5]*STATES[4]; ALGEBRAIC[9] = CONSTANTS[6]*STATES[22]*STATES[11] - CONSTANTS[7]*STATES[6]; ALGEBRAIC[10] = CONSTANTS[8]*STATES[1] - CONSTANTS[9]*STATES[5]; ALGEBRAIC[11] = CONSTANTS[8]*STATES[4]; ALGEBRAIC[12] = CONSTANTS[10]*STATES[4]*STATES[13] - CONSTANTS[11]*STATES[14]; ALGEBRAIC[13] = CONSTANTS[8]*STATES[22]; ALGEBRAIC[14] = CONSTANTS[12]*STATES[5]*STATES[15] - CONSTANTS[13]*STATES[9]; ALGEBRAIC[15] = CONSTANTS[2]*STATES[9]*STATES[9] - CONSTANTS[3]*STATES[10]; ALGEBRAIC[16] = CONSTANTS[4]*STATES[10] - CONSTANTS[5]*STATES[7]; ALGEBRAIC[17] = CONSTANTS[10]*STATES[7]*STATES[13] - CONSTANTS[11]*STATES[16]; ALGEBRAIC[18] = CONSTANTS[15]*STATES[8]; ALGEBRAIC[19] = CONSTANTS[16]*STATES[21]*STATES[14] - CONSTANTS[17]*STATES[22]; ALGEBRAIC[20] = CONSTANTS[18]*STATES[23]*STATES[22] - CONSTANTS[19]*STATES[24]; ALGEBRAIC[21] = CONSTANTS[20]*STATES[25]*STATES[24] - CONSTANTS[21]*STATES[26]; ALGEBRAIC[22] = CONSTANTS[22]*STATES[26] - CONSTANTS[23]*STATES[27]*STATES[24]; ALGEBRAIC[23] = CONSTANTS[24]*STATES[24]*STATES[42] - CONSTANTS[25]*STATES[28]; ALGEBRAIC[24] = CONSTANTS[26]*STATES[28] - CONSTANTS[27]*STATES[24]*STATES[25]; ALGEBRAIC[25] = CONSTANTS[28]*STATES[30]*STATES[14] - CONSTANTS[29]*STATES[31]; ALGEBRAIC[26] = CONSTANTS[30]*STATES[31] - CONSTANTS[31]*STATES[32]; ALGEBRAIC[27] = CONSTANTS[16]*STATES[21]*STATES[32] - CONSTANTS[32]*STATES[33]; ALGEBRAIC[28] = CONSTANTS[33]*STATES[23]*STATES[33] - CONSTANTS[34]*STATES[34]; ALGEBRAIC[29] = CONSTANTS[20]*STATES[25]*STATES[34] - CONSTANTS[21]*STATES[35]; ALGEBRAIC[30] = CONSTANTS[22]*STATES[35] - CONSTANTS[23]*STATES[34]*STATES[27]; ALGEBRAIC[31] = CONSTANTS[35]*STATES[27]*STATES[40] - CONSTANTS[36]*STATES[41]; ALGEBRAIC[32] = CONSTANTS[37]*STATES[41] - CONSTANTS[38]*STATES[42]*STATES[44]; ALGEBRAIC[33] = CONSTANTS[24]*STATES[34]*STATES[42] - CONSTANTS[25]*STATES[36]; ALGEBRAIC[34] = CONSTANTS[26]*STATES[36] - CONSTANTS[27]*STATES[34]*STATES[25]; ALGEBRAIC[35] = CONSTANTS[39]*STATES[34] - CONSTANTS[40]*STATES[37]*STATES[14]; ALGEBRAIC[36] = CONSTANTS[41]*STATES[37] - CONSTANTS[42]*STATES[39]*STATES[29]; ALGEBRAIC[37] = CONSTANTS[43]*STATES[24] - CONSTANTS[44]*STATES[14]*STATES[29]; ALGEBRAIC[38] = CONSTANTS[45]*STATES[29] - CONSTANTS[46]*STATES[23]*STATES[21]; ALGEBRAIC[39] = ( CONSTANTS[48]*STATES[39])/(CONSTANTS[47]+STATES[39]); ALGEBRAIC[40] = CONSTANTS[49]*STATES[32] - CONSTANTS[50]*STATES[14]*STATES[39]; ALGEBRAIC[41] = CONSTANTS[16]*STATES[21]*STATES[39] - CONSTANTS[32]*STATES[38]; ALGEBRAIC[42] = CONSTANTS[49]*STATES[33] - CONSTANTS[50]*STATES[14]*STATES[38]; ALGEBRAIC[43] = CONSTANTS[51]*STATES[23]*STATES[38] - CONSTANTS[52]*STATES[37]; ALGEBRAIC[44] = CONSTANTS[53]*STATES[29]*STATES[32] - CONSTANTS[54]*STATES[34]; ALGEBRAIC[45] = CONSTANTS[55]*STATES[43]*STATES[44] - CONSTANTS[56]*STATES[45]; ALGEBRAIC[46] = CONSTANTS[57]*STATES[45]; ALGEBRAIC[47] = CONSTANTS[58]*STATES[46]*STATES[44] - CONSTANTS[59]*STATES[47]; ALGEBRAIC[48] = CONSTANTS[60]*STATES[47]; ALGEBRAIC[49] = CONSTANTS[58]*STATES[48]*STATES[44] - CONSTANTS[59]*STATES[49]; ALGEBRAIC[50] = CONSTANTS[61]*STATES[49]; ALGEBRAIC[51] = CONSTANTS[62]*STATES[50]*STATES[52] - CONSTANTS[63]*STATES[51]; ALGEBRAIC[52] = CONSTANTS[64]*STATES[51]; ALGEBRAIC[53] = CONSTANTS[65]*STATES[52]*STATES[48] - CONSTANTS[66]*STATES[53]; ALGEBRAIC[54] = CONSTANTS[64]*STATES[53]; ALGEBRAIC[55] = CONSTANTS[67]*STATES[54]*STATES[50] - CONSTANTS[68]*STATES[55]; ALGEBRAIC[56] = CONSTANTS[69]*STATES[55]; ALGEBRAIC[57] = CONSTANTS[67]*STATES[50]*STATES[56] - CONSTANTS[68]*STATES[57]; ALGEBRAIC[58] = CONSTANTS[70]*STATES[57]; ALGEBRAIC[59] = CONSTANTS[71]*STATES[58]*STATES[59] - CONSTANTS[72]*STATES[60]; ALGEBRAIC[60] = CONSTANTS[73]*STATES[60]; ALGEBRAIC[61] = CONSTANTS[74]*STATES[59]*STATES[56] - CONSTANTS[75]*STATES[61]; ALGEBRAIC[62] = CONSTANTS[76]*STATES[61]; ALGEBRAIC[63] = CONSTANTS[77]*STATES[5]; ALGEBRAIC[64] = CONSTANTS[78]*STATES[15]; ALGEBRAIC[65] = CONSTANTS[77]*STATES[7]; ALGEBRAIC[66] = CONSTANTS[16]*STATES[16]*STATES[21] - CONSTANTS[17]*STATES[17]; ALGEBRAIC[67] = CONSTANTS[18]*STATES[23]*STATES[17] - CONSTANTS[19]*STATES[18]; ALGEBRAIC[68] = CONSTANTS[20]*STATES[25]*STATES[18] - CONSTANTS[21]*STATES[19]; ALGEBRAIC[69] = CONSTANTS[22]*STATES[19] - CONSTANTS[23]*STATES[68]*STATES[18]; ALGEBRAIC[70] = CONSTANTS[24]*STATES[70]*STATES[18] - CONSTANTS[25]*STATES[20]; ALGEBRAIC[71] = CONSTANTS[26]*STATES[20] - CONSTANTS[27]*STATES[18]*STATES[25]; ALGEBRAIC[72] = CONSTANTS[28]*STATES[30]*STATES[16] - CONSTANTS[29]*STATES[62]; ALGEBRAIC[73] = CONSTANTS[30]*STATES[62] - CONSTANTS[31]*STATES[63]; ALGEBRAIC[74] = CONSTANTS[16]*STATES[21]*STATES[63] - CONSTANTS[32]*STATES[64]; ALGEBRAIC[75] = CONSTANTS[33]*STATES[23]*STATES[64] - CONSTANTS[34]*STATES[65]; ALGEBRAIC[76] = CONSTANTS[20]*STATES[25]*STATES[65] - CONSTANTS[21]*STATES[66]; ALGEBRAIC[77] = CONSTANTS[22]*STATES[66] - CONSTANTS[23]*STATES[65]*STATES[68]; ALGEBRAIC[78] = CONSTANTS[35]*STATES[68]*STATES[40] - CONSTANTS[36]*STATES[69]; ALGEBRAIC[79] = CONSTANTS[37]*STATES[69] - CONSTANTS[38]*STATES[70]*STATES[71]; ALGEBRAIC[80] = CONSTANTS[24]*STATES[70]*STATES[65] - CONSTANTS[25]*STATES[67]; ALGEBRAIC[81] = CONSTANTS[26]*STATES[67] - STATES[25]*CONSTANTS[27]*STATES[65]; ALGEBRAIC[82] = CONSTANTS[39]*STATES[65] - CONSTANTS[40]*STATES[16]*STATES[37]; ALGEBRAIC[83] = CONSTANTS[43]*STATES[18] - CONSTANTS[44]*STATES[16]*STATES[29]; ALGEBRAIC[84] = CONSTANTS[49]*STATES[63] - CONSTANTS[50]*STATES[16]*STATES[39]; ALGEBRAIC[85] = CONSTANTS[49]*STATES[64] - CONSTANTS[50]*STATES[16]*STATES[38]; ALGEBRAIC[86] = CONSTANTS[53]*STATES[29]*STATES[63] - CONSTANTS[54]*STATES[65]; ALGEBRAIC[87] = CONSTANTS[55]*STATES[43]*STATES[71] - CONSTANTS[56]*STATES[72]; ALGEBRAIC[88] = CONSTANTS[57]*STATES[72]; ALGEBRAIC[89] = CONSTANTS[58]*STATES[46]*STATES[71] - CONSTANTS[59]*STATES[73]; ALGEBRAIC[90] = CONSTANTS[60]*STATES[73]; ALGEBRAIC[91] = CONSTANTS[58]*STATES[71]*STATES[74] - CONSTANTS[59]*STATES[75]; ALGEBRAIC[92] = CONSTANTS[61]*STATES[75]; ALGEBRAIC[93] = CONSTANTS[62]*STATES[76]*STATES[52] - CONSTANTS[63]*STATES[77]; ALGEBRAIC[94] = CONSTANTS[64]*STATES[77]; ALGEBRAIC[95] = CONSTANTS[65]*STATES[52]*STATES[74] - CONSTANTS[66]*STATES[78]; ALGEBRAIC[96] = CONSTANTS[64]*STATES[78]; ALGEBRAIC[97] = CONSTANTS[67]*STATES[54]*STATES[76] - CONSTANTS[68]*STATES[79]; ALGEBRAIC[98] = CONSTANTS[69]*STATES[79]; ALGEBRAIC[99] = CONSTANTS[67]*STATES[76]*STATES[80] - CONSTANTS[68]*STATES[81]; ALGEBRAIC[100] = CONSTANTS[70]*STATES[81]; ALGEBRAIC[101] = CONSTANTS[71]*STATES[82]*STATES[59] - CONSTANTS[72]*STATES[83]; ALGEBRAIC[102] = CONSTANTS[73]*STATES[83]; ALGEBRAIC[103] = CONSTANTS[74]*STATES[59]*STATES[80] - CONSTANTS[75]*STATES[84]; ALGEBRAIC[104] = CONSTANTS[76]*STATES[84]; ALGEBRAIC[105] = CONSTANTS[8]*STATES[14] - CONSTANTS[9]*STATES[16]; ALGEBRAIC[106] = CONSTANTS[8]*STATES[31] - CONSTANTS[9]*STATES[62]; ALGEBRAIC[107] = CONSTANTS[8]*STATES[32] - CONSTANTS[9]*STATES[63]; ALGEBRAIC[108] = CONSTANTS[8]*STATES[24] - CONSTANTS[9]*STATES[18]; ALGEBRAIC[109] = CONSTANTS[6]*STATES[24]*STATES[11] - CONSTANTS[7]*STATES[87]; ALGEBRAIC[110] = CONSTANTS[8]*STATES[26] - CONSTANTS[9]*STATES[19]; ALGEBRAIC[111] = CONSTANTS[6]*STATES[26]*STATES[11] - CONSTANTS[7]*STATES[88]; ALGEBRAIC[112] = CONSTANTS[8]*STATES[28] - CONSTANTS[9]*STATES[20]; ALGEBRAIC[113] = CONSTANTS[6]*STATES[28]*STATES[11] - CONSTANTS[7]*STATES[89]; ALGEBRAIC[114] = CONSTANTS[8]*STATES[33] - CONSTANTS[9]*STATES[64]; ALGEBRAIC[115] = CONSTANTS[6]*STATES[33]*STATES[11] - CONSTANTS[7]*STATES[90]; ALGEBRAIC[116] = CONSTANTS[8]*STATES[34] - CONSTANTS[9]*STATES[65]; ALGEBRAIC[117] = CONSTANTS[6]*STATES[34]*STATES[11] - CONSTANTS[7]*STATES[91]; ALGEBRAIC[118] = CONSTANTS[8]*STATES[35] - CONSTANTS[9]*STATES[66]; ALGEBRAIC[119] = CONSTANTS[6]*STATES[35]*STATES[11] - CONSTANTS[7]*STATES[92]; ALGEBRAIC[120] = CONSTANTS[8]*STATES[36] - CONSTANTS[9]*STATES[67]; ALGEBRAIC[121] = CONSTANTS[6]*STATES[36]*STATES[11] - CONSTANTS[7]*STATES[93]; ALGEBRAIC[122] = ( CONSTANTS[79]*1.00000)/(CONSTANTS[1]/( CONSTANTS[0]*STATES[0])+1.00000); ALGEBRAIC[123] = (CONDVAR[0]<0.00000 ? 1.55000 : CONDVAR[1]>0.00000 ? 0.200000 : ( ALGEBRAIC[122]*- 1.35000e-05+1.55000)/1.00000); ALGEBRAIC[124] = ALGEBRAIC[123]*STATES[6]; ALGEBRAIC[125] = ALGEBRAIC[123]*STATES[87]; ALGEBRAIC[126] = ALGEBRAIC[123]*STATES[88]; ALGEBRAIC[127] = ALGEBRAIC[123]*STATES[89]; ALGEBRAIC[128] = ALGEBRAIC[123]*STATES[90]; ALGEBRAIC[129] = ALGEBRAIC[123]*STATES[91]; ALGEBRAIC[130] = ALGEBRAIC[123]*STATES[92]; ALGEBRAIC[131] = ALGEBRAIC[123]*STATES[93]; } void getStateInformation(double* SI) { SI[0] = 1.0; SI[1] = 1.0; SI[2] = 1.0; SI[3] = 1.0; SI[4] = 1.0; SI[5] = 1.0; SI[6] = 1.0; SI[7] = 1.0; SI[8] = 1.0; SI[9] = 1.0; SI[10] = 1.0; SI[11] = 1.0; SI[12] = 1.0; SI[13] = 1.0; SI[14] = 1.0; SI[15] = 1.0; SI[16] = 1.0; SI[17] = 1.0; SI[18] = 1.0; SI[19] = 1.0; SI[20] = 1.0; SI[21] = 1.0; SI[22] = 1.0; SI[23] = 1.0; SI[24] = 1.0; SI[25] = 1.0; SI[26] = 1.0; SI[27] = 1.0; SI[28] = 1.0; SI[29] = 1.0; SI[30] = 1.0; SI[31] = 1.0; SI[32] = 1.0; SI[33] = 1.0; SI[34] = 1.0; SI[35] = 1.0; SI[36] = 1.0; SI[37] = 1.0; SI[38] = 1.0; SI[39] = 1.0; SI[40] = 1.0; SI[41] = 1.0; SI[42] = 1.0; SI[43] = 1.0; SI[44] = 1.0; SI[45] = 1.0; SI[46] = 1.0; SI[47] = 1.0; SI[48] = 1.0; SI[49] = 1.0; SI[50] = 1.0; SI[51] = 1.0; SI[52] = 1.0; SI[53] = 1.0; SI[54] = 1.0; SI[55] = 1.0; SI[56] = 1.0; SI[57] = 1.0; SI[58] = 1.0; SI[59] = 1.0; SI[60] = 1.0; SI[61] = 1.0; SI[62] = 1.0; SI[63] = 1.0; SI[64] = 1.0; SI[65] = 1.0; SI[66] = 1.0; SI[67] = 1.0; SI[68] = 1.0; SI[69] = 1.0; SI[70] = 1.0; SI[71] = 1.0; SI[72] = 1.0; SI[73] = 1.0; SI[74] = 1.0; SI[75] = 1.0; SI[76] = 1.0; SI[77] = 1.0; SI[78] = 1.0; SI[79] = 1.0; SI[80] = 1.0; SI[81] = 1.0; SI[82] = 1.0; SI[83] = 1.0; SI[84] = 1.0; SI[85] = 1.0; SI[86] = 1.0; SI[87] = 1.0; SI[88] = 1.0; SI[89] = 1.0; SI[90] = 1.0; SI[91] = 1.0; SI[92] = 1.0; SI[93] = 1.0; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { CONDVAR[0] = ALGEBRAIC[122] - 3100.00; CONDVAR[1] = ALGEBRAIC[122] - 100000.; }