Generated Code
The following is c_ida code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/*
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.;
}