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 148 entries in the algebraic variable array.
There are a total of 103 entries in each of the rate and state variable arrays.
There are a total of 96 entries in the constant variable array.
*/
/*
* VOI is time in component environment (second).
* STATES[0] is c1 in component c1 (molar).
* ALGEBRAIC[0] is v1 in component v1 (flux).
* STATES[1] is c2 in component c2 (molar).
* ALGEBRAIC[5] is v6 in component v6 (flux).
* ALGEBRAIC[12] is v13 in component v13 (flux).
* STATES[2] is c3 in component c3 (molar).
* ALGEBRAIC[1] is v2 in component v2 (flux).
* STATES[3] is c4 in component c4 (molar).
* ALGEBRAIC[2] is v3 in component v3 (flux).
* STATES[4] is c5 in component c5 (molar).
* ALGEBRAIC[6] is v7 in component v7 (flux).
* ALGEBRAIC[7] is v8 in component v8 (flux).
* STATES[5] is c6 in component c6 (molar).
* ALGEBRAIC[9] is v10 in component v10 (flux).
* ALGEBRAIC[59] is v60 in component v60 (flux).
* STATES[6] is c7 in component c7 (molar).
* ALGEBRAIC[3] is v4 in component v4 (flux).
* ALGEBRAIC[4] is v5 in component v5 (flux).
* STATES[7] is c8 in component c8 (molar).
* ALGEBRAIC[11] is v12 in component v12 (flux).
* ALGEBRAIC[13] is v14 in component v14 (flux).
* ALGEBRAIC[61] is v62 in component v62 (flux).
* STATES[8] is c9 in component c9 (molar).
* ALGEBRAIC[14] is v15 in component v15 (flux).
* ALGEBRAIC[106] is v107 in component v107 (flux).
* ALGEBRAIC[109] is v110 in component v110 (flux).
* ALGEBRAIC[112] is v113 in component v113 (flux).
* ALGEBRAIC[115] is v116 in component v116 (flux).
* ALGEBRAIC[118] is v119 in component v119 (flux).
* ALGEBRAIC[121] is v122 in component v122 (flux).
* ALGEBRAIC[124] is v125 in component v125 (flux).
* STATES[9] is c10 in component c10 (molar).
* ALGEBRAIC[10] is v11 in component v11 (flux).
* STATES[10] is c11 in component c11 (molar).
* STATES[11] is c12 in component c12 (molar).
* ALGEBRAIC[105] is v106 in component v106 (flux).
* ALGEBRAIC[108] is v109 in component v109 (flux).
* ALGEBRAIC[111] is v112 in component v112 (flux).
* ALGEBRAIC[114] is v115 in component v115 (flux).
* ALGEBRAIC[117] is v118 in component v118 (flux).
* ALGEBRAIC[120] is v121 in component v121 (flux).
* ALGEBRAIC[123] is v124 in component v124 (flux).
* STATES[12] is c13 in component c13 (molar).
* ALGEBRAIC[60] is v61 in component v61 (flux).
* STATES[13] is c14 in component c14 (molar).
* STATES[14] is c15 in component c15 (molar).
* ALGEBRAIC[15] is v16 in component v16 (flux).
* ALGEBRAIC[21] is v22 in component v22 (flux).
* ALGEBRAIC[31] is v32 in component v32 (flux).
* ALGEBRAIC[33] is v34 in component v34 (flux).
* ALGEBRAIC[36] is v37 in component v37 (flux).
* ALGEBRAIC[38] is v39 in component v39 (flux).
* ALGEBRAIC[101] is v102 in component v102 (flux).
* STATES[15] is c16 in component c16 (molar).
* STATES[16] is c17 in component c17 (molar).
* ALGEBRAIC[62] is v63 in component v63 (flux).
* ALGEBRAIC[68] is v69 in component v69 (flux).
* ALGEBRAIC[78] is v79 in component v79 (flux).
* ALGEBRAIC[79] is v80 in component v80 (flux).
* ALGEBRAIC[80] is v81 in component v81 (flux).
* ALGEBRAIC[81] is v82 in component v82 (flux).
* ALGEBRAIC[131] is v132 in component v132 (flux).
* STATES[17] is c18 in component c18 (molar).
* ALGEBRAIC[8] is v9 in component v9 (flux).
* ALGEBRAIC[63] is v64 in component v64 (flux).
* ALGEBRAIC[132] is v133 in component v133 (flux).
* STATES[18] is c19 in component c19 (molar).
* ALGEBRAIC[64] is v65 in component v65 (flux).
* ALGEBRAIC[65] is v66 in component v66 (flux).
* ALGEBRAIC[66] is v67 in component v67 (flux).
* ALGEBRAIC[67] is v68 in component v68 (flux).
* ALGEBRAIC[104] is v105 in component v105 (flux).
* ALGEBRAIC[126] is v127 in component v127 (flux).
* ALGEBRAIC[133] is v134 in component v134 (flux).
* STATES[19] is c20 in component c20 (molar).
* ALGEBRAIC[107] is v108 in component v108 (flux).
* ALGEBRAIC[134] is v135 in component v135 (flux).
* STATES[20] is c21 in component c21 (molar).
* ALGEBRAIC[110] is v111 in component v111 (flux).
* ALGEBRAIC[135] is v136 in component v136 (flux).
* STATES[21] is c22 in component c22 (molar).
* ALGEBRAIC[23] is v24 in component v24 (flux).
* ALGEBRAIC[34] is v35 in component v35 (flux).
* ALGEBRAIC[37] is v38 in component v38 (flux).
* ALGEBRAIC[70] is v71 in component v71 (flux).
* STATES[22] is c23 in component c23 (molar).
* ALGEBRAIC[16] is v17 in component v17 (flux).
* STATES[23] is c24 in component c24 (molar).
* ALGEBRAIC[71] is v72 in component v72 (flux).
* ALGEBRAIC[129] is v130 in component v130 (flux).
* ALGEBRAIC[130] is v131 in component v131 (flux).
* ALGEBRAIC[24] is v25 in component v25 (flux).
* ALGEBRAIC[39] is v40 in component v40 (flux).
* STATES[24] is c25 in component c25 (molar).
* ALGEBRAIC[17] is v18 in component v18 (flux).
* ALGEBRAIC[18] is v19 in component v19 (flux).
* ALGEBRAIC[19] is v20 in component v20 (flux).
* ALGEBRAIC[20] is v21 in component v21 (flux).
* ALGEBRAIC[125] is v126 in component v126 (flux).
* STATES[25] is c26 in component c26 (molar).
* ALGEBRAIC[25] is v26 in component v26 (flux).
* ALGEBRAIC[30] is v31 in component v31 (flux).
* ALGEBRAIC[77] is v78 in component v78 (flux).
* ALGEBRAIC[72] is v73 in component v73 (flux).
* STATES[26] is c27 in component c27 (molar).
* STATES[27] is c28 in component c28 (molar).
* ALGEBRAIC[26] is v27 in component v27 (flux).
* ALGEBRAIC[27] is v28 in component v28 (flux).
* STATES[28] is c29 in component c29 (molar).
* STATES[29] is c30 in component c30 (molar).
* ALGEBRAIC[32] is v33 in component v33 (flux).
* ALGEBRAIC[40] is v41 in component v41 (flux).
* ALGEBRAIC[82] is v83 in component v83 (flux).
* STATES[30] is c31 in component c31 (molar).
* ALGEBRAIC[35] is v36 in component v36 (flux).
* STATES[31] is c32 in component c32 (molar).
* ALGEBRAIC[22] is v23 in component v23 (flux).
* ALGEBRAIC[102] is v103 in component v103 (flux).
* STATES[32] is c33 in component c33 (molar).
* ALGEBRAIC[103] is v104 in component v104 (flux).
* STATES[33] is c34 in component c34 (molar).
* ALGEBRAIC[113] is v114 in component v114 (flux).
* STATES[34] is c35 in component c35 (molar).
* ALGEBRAIC[29] is v30 in component v30 (flux).
* ALGEBRAIC[116] is v117 in component v117 (flux).
* ALGEBRAIC[127] is v128 in component v128 (flux).
* STATES[35] is c36 in component c36 (molar).
* ALGEBRAIC[119] is v120 in component v120 (flux).
* STATES[36] is c37 in component c37 (molar).
* ALGEBRAIC[122] is v123 in component v123 (flux).
* STATES[37] is c38 in component c38 (molar).
* STATES[38] is c39 in component c39 (molar).
* STATES[39] is c40 in component c40 (molar).
* STATES[40] is c41 in component c41 (molar).
* ALGEBRAIC[42] is v43 in component v43 (flux).
* ALGEBRAIC[74] is v75 in component v75 (flux).
* ALGEBRAIC[84] is v85 in component v85 (flux).
* STATES[41] is c42 in component c42 (molar).
* ALGEBRAIC[28] is v29 in component v29 (flux).
* STATES[42] is c43 in component c43 (molar).
* STATES[43] is c44 in component c44 (molar).
* ALGEBRAIC[41] is v42 in component v42 (flux).
* ALGEBRAIC[83] is v84 in component v84 (flux).
* STATES[44] is c45 in component c45 (molar).
* ALGEBRAIC[43] is v44 in component v44 (flux).
* ALGEBRAIC[44] is v45 in component v45 (flux).
* ALGEBRAIC[45] is v46 in component v46 (flux).
* ALGEBRAIC[46] is v47 in component v47 (flux).
* STATES[45] is c46 in component c46 (molar).
* STATES[46] is c47 in component c47 (molar).
* ALGEBRAIC[50] is v51 in component v51 (flux).
* ALGEBRAIC[85] is v86 in component v86 (flux).
* ALGEBRAIC[92] is v93 in component v93 (flux).
* STATES[47] is c48 in component c48 (molar).
* STATES[48] is c49 in component c49 (molar).
* ALGEBRAIC[48] is v49 in component v49 (flux).
* ALGEBRAIC[49] is v50 in component v50 (flux).
* STATES[49] is c50 in component c50 (molar).
* STATES[50] is c51 in component c51 (molar).
* ALGEBRAIC[47] is v48 in component v48 (flux).
* ALGEBRAIC[51] is v52 in component v52 (flux).
* ALGEBRAIC[52] is v53 in component v53 (flux).
* ALGEBRAIC[53] is v54 in component v54 (flux).
* ALGEBRAIC[54] is v55 in component v55 (flux).
* STATES[51] is c52 in component c52 (molar).
* STATES[52] is c53 in component c53 (molar).
* ALGEBRAIC[89] is v90 in component v90 (flux).
* ALGEBRAIC[90] is v91 in component v91 (flux).
* ALGEBRAIC[91] is v92 in component v92 (flux).
* STATES[53] is c54 in component c54 (molar).
* STATES[54] is c55 in component c55 (molar).
* ALGEBRAIC[58] is v59 in component v59 (flux).
* ALGEBRAIC[93] is v94 in component v94 (flux).
* ALGEBRAIC[100] is v101 in component v101 (flux).
* STATES[55] is c56 in component c56 (molar).
* STATES[56] is c57 in component c57 (molar).
* ALGEBRAIC[56] is v57 in component v57 (flux).
* ALGEBRAIC[57] is v58 in component v58 (flux).
* STATES[57] is c58 in component c58 (molar).
* STATES[58] is c59 in component c59 (molar).
* ALGEBRAIC[55] is v56 in component v56 (flux).
* ALGEBRAIC[142] is v143 in component v143 (flux).
* ALGEBRAIC[143] is v144 in component v144 (flux).
* ALGEBRAIC[144] is v145 in component v145 (flux).
* STATES[59] is c60 in component c60 (molar).
* ALGEBRAIC[97] is v98 in component v98 (flux).
* ALGEBRAIC[98] is v99 in component v99 (flux).
* ALGEBRAIC[99] is v100 in component v100 (flux).
* STATES[60] is c61 in component c61 (molar).
* STATES[61] is c62 in component c62 (molar).
* STATES[62] is c63 in component c63 (molar).
* ALGEBRAIC[69] is v70 in component v70 (flux).
* ALGEBRAIC[136] is v137 in component v137 (flux).
* STATES[63] is c64 in component c64 (molar).
* ALGEBRAIC[137] is v138 in component v138 (flux).
* STATES[64] is c65 in component c65 (molar).
* ALGEBRAIC[138] is v139 in component v139 (flux).
* STATES[65] is c66 in component c66 (molar).
* ALGEBRAIC[73] is v74 in component v74 (flux).
* ALGEBRAIC[76] is v77 in component v77 (flux).
* ALGEBRAIC[128] is v129 in component v129 (flux).
* ALGEBRAIC[139] is v140 in component v140 (flux).
* STATES[66] is c67 in component c67 (molar).
* ALGEBRAIC[140] is v141 in component v141 (flux).
* STATES[67] is c68 in component c68 (molar).
* ALGEBRAIC[141] is v142 in component v142 (flux).
* STATES[68] is c69 in component c69 (molar).
* STATES[69] is c70 in component c70 (molar).
* ALGEBRAIC[75] is v76 in component v76 (flux).
* STATES[70] is c71 in component c71 (molar).
* STATES[71] is c72 in component c72 (molar).
* ALGEBRAIC[86] is v87 in component v87 (flux).
* ALGEBRAIC[87] is v88 in component v88 (flux).
* ALGEBRAIC[88] is v89 in component v89 (flux).
* STATES[72] is c73 in component c73 (molar).
* STATES[73] is c74 in component c74 (molar).
* STATES[74] is c75 in component c75 (molar).
* STATES[75] is c76 in component c76 (molar).
* STATES[76] is c77 in component c77 (molar).
* ALGEBRAIC[94] is v95 in component v95 (flux).
* ALGEBRAIC[95] is v96 in component v96 (flux).
* ALGEBRAIC[96] is v97 in component v97 (flux).
* STATES[77] is c78 in component c78 (molar).
* STATES[78] is c79 in component c79 (molar).
* STATES[79] is c80 in component c80 (molar).
* STATES[80] is c81 in component c81 (molar).
* STATES[81] is c82 in component c82 (molar).
* STATES[82] is c83 in component c83 (molar).
* ALGEBRAIC[145] is v146 in component v146 (flux).
* ALGEBRAIC[146] is v147 in component v147 (flux).
* ALGEBRAIC[147] is v148 in component v148 (flux).
* STATES[83] is c84 in component c84 (molar).
* STATES[84] is c85 in component c85 (molar).
* STATES[85] is c86 in component c86 (molar).
* STATES[86] is c87 in component c87 (molar).
* STATES[87] is c88 in component c88 (molar).
* STATES[88] is c89 in component c89 (molar).
* STATES[89] is c90 in component c90 (molar).
* STATES[90] is c91 in component c91 (molar).
* STATES[91] is c92 in component c92 (molar).
* STATES[92] is c93 in component c93 (molar).
* STATES[93] is c94 in component c94 (molar).
* STATES[94] is c95 in component c95 (molar).
* STATES[95] is c96 in component c96 (molar).
* STATES[96] is c97 in component c97 (molar).
* STATES[97] is c98 in component c98 (molar).
* STATES[98] is c99 in component c99 (molar).
* STATES[99] is c100 in component c100 (molar).
* STATES[100] is c101 in component c101 (molar).
* STATES[101] is c102 in component c102 (molar).
* STATES[102] is c103 in component c103 (molar).
* 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).
* CONSTANTS[8] is k5 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[9] is kd5 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[10] is k6 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[11] is kd6 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[12] is k8 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[13] is kd8 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[14] is k10b in component reaction_rates (second_order_rate_constant).
* CONSTANTS[15] is kd10 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[16] is k13 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[17] is kd13 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[18] is k15 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[19] is kd15 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[20] is k16 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[21] is kd63 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[22] is k17 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[23] is kd17 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[24] is k18 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[25] is kd18 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[26] is k19 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[27] is kd19 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[28] is k20 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[29] is kd20 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[30] is k21 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[31] is kd21 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[32] is k22 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[33] is kd22 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[34] is k23 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[35] is kd23 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[36] is kd24 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[37] is k25 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[38] is kd25 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[39] is k28 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[40] is kd28 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[41] is k29 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[42] is kd29 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[43] is k32 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[44] is kd32 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[45] is k33 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[46] is kd33 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[47] is k34 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[48] is kd34 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[49] is k35 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[50] is kd35 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[51] is k36 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[52] is kd36 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[53] is k37 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[54] is kd37 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[55] is k40 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[56] is kd40 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[57] is k41 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[58] is kd41 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[59] is k42 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[60] is kd42 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[61] is k43 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[62] is kd43 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[63] is k44 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[64] is kd52 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[65] is k45 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[66] is kd45 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[67] is k47 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[68] is kd47 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[69] is k48 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[70] is kd48 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[71] is k49 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[72] is kd49 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[73] is k50 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[74] is kd50 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[75] is k52 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[76] is kd44 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[77] is k53 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[78] is kd53 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[79] is k55 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[80] is kd55 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[81] is k56 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[82] is kd56 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[83] is k57 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[84] is kd57 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[85] is k58 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[86] is kd58 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[87] is k60 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[88] is kd60 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[89] is k61 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[90] is kd61 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[91] is k126 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[92] is kd126 in component reaction_rates (first_order_rate_constant).
* CONSTANTS[93] is k127 in component reaction_rates (second_order_rate_constant).
* CONSTANTS[94] is kd127 in component reaction_rates (first_order_rate_constant).
* RATES[0] is d/dt c1 in component c1 (molar).
* RATES[1] is d/dt c2 in component c2 (molar).
* RATES[2] is d/dt c3 in component c3 (molar).
* RATES[3] is d/dt c4 in component c4 (molar).
* RATES[4] is d/dt c5 in component c5 (molar).
* RATES[5] is d/dt c6 in component c6 (molar).
* RATES[6] is d/dt c7 in component c7 (molar).
* RATES[7] is d/dt c8 in component c8 (molar).
* RATES[8] is d/dt c9 in component c9 (molar).
* RATES[9] is d/dt c10 in component c10 (molar).
* RATES[10] is d/dt c11 in component c11 (molar).
* RATES[11] is d/dt c12 in component c12 (molar).
* RATES[12] is d/dt c13 in component c13 (molar).
* RATES[13] is d/dt c14 in component c14 (molar).
* RATES[14] is d/dt c15 in component c15 (molar).
* RATES[15] is d/dt c16 in component c16 (molar).
* RATES[16] is d/dt c17 in component c17 (molar).
* RATES[17] is d/dt c18 in component c18 (molar).
* RATES[18] is d/dt c19 in component c19 (molar).
* RATES[19] is d/dt c20 in component c20 (molar).
* RATES[20] is d/dt c21 in component c21 (molar).
* RATES[21] is d/dt c22 in component c22 (molar).
* RATES[22] is d/dt c23 in component c23 (molar).
* RATES[23] is d/dt c24 in component c24 (molar).
* RATES[24] is d/dt c25 in component c25 (molar).
* RATES[25] is d/dt c26 in component c26 (molar).
* RATES[26] is d/dt c27 in component c27 (molar).
* RATES[27] is d/dt c28 in component c28 (molar).
* RATES[28] is d/dt c29 in component c29 (molar).
* RATES[29] is d/dt c30 in component c30 (molar).
* RATES[30] is d/dt c31 in component c31 (molar).
* RATES[31] is d/dt c32 in component c32 (molar).
* RATES[32] is d/dt c33 in component c33 (molar).
* RATES[33] is d/dt c34 in component c34 (molar).
* RATES[34] is d/dt c35 in component c35 (molar).
* RATES[35] is d/dt c36 in component c36 (molar).
* RATES[36] is d/dt c37 in component c37 (molar).
* RATES[37] is d/dt c38 in component c38 (molar).
* RATES[38] is d/dt c39 in component c39 (molar).
* RATES[39] is d/dt c40 in component c40 (molar).
* RATES[40] is d/dt c41 in component c41 (molar).
* RATES[41] is d/dt c42 in component c42 (molar).
* RATES[42] is d/dt c43 in component c43 (molar).
* RATES[43] is d/dt c44 in component c44 (molar).
* RATES[44] is d/dt c45 in component c45 (molar).
* RATES[45] is d/dt c46 in component c46 (molar).
* RATES[46] is d/dt c47 in component c47 (molar).
* RATES[47] is d/dt c48 in component c48 (molar).
* RATES[48] is d/dt c49 in component c49 (molar).
* RATES[49] is d/dt c50 in component c50 (molar).
* RATES[50] is d/dt c51 in component c51 (molar).
* RATES[51] is d/dt c52 in component c52 (molar).
* RATES[52] is d/dt c53 in component c53 (molar).
* RATES[53] is d/dt c54 in component c54 (molar).
* RATES[54] is d/dt c55 in component c55 (molar).
* RATES[55] is d/dt c56 in component c56 (molar).
* RATES[56] is d/dt c57 in component c57 (molar).
* RATES[57] is d/dt c58 in component c58 (molar).
* RATES[58] is d/dt c59 in component c59 (molar).
* RATES[59] is d/dt c60 in component c60 (molar).
* RATES[60] is d/dt c61 in component c61 (molar).
* RATES[61] is d/dt c62 in component c62 (molar).
* RATES[62] is d/dt c63 in component c63 (molar).
* RATES[63] is d/dt c64 in component c64 (molar).
* RATES[64] is d/dt c65 in component c65 (molar).
* RATES[65] is d/dt c66 in component c66 (molar).
* RATES[66] is d/dt c67 in component c67 (molar).
* RATES[67] is d/dt c68 in component c68 (molar).
* RATES[68] is d/dt c69 in component c69 (molar).
* RATES[69] is d/dt c70 in component c70 (molar).
* RATES[70] is d/dt c71 in component c71 (molar).
* RATES[71] is d/dt c72 in component c72 (molar).
* RATES[72] is d/dt c73 in component c73 (molar).
* RATES[73] is d/dt c74 in component c74 (molar).
* RATES[74] is d/dt c75 in component c75 (molar).
* RATES[75] is d/dt c76 in component c76 (molar).
* RATES[76] is d/dt c77 in component c77 (molar).
* RATES[77] is d/dt c78 in component c78 (molar).
* RATES[78] is d/dt c79 in component c79 (molar).
* RATES[79] is d/dt c80 in component c80 (molar).
* RATES[80] is d/dt c81 in component c81 (molar).
* RATES[81] is d/dt c82 in component c82 (molar).
* RATES[82] is d/dt c83 in component c83 (molar).
* RATES[83] is d/dt c84 in component c84 (molar).
* RATES[84] is d/dt c85 in component c85 (molar).
* RATES[85] is d/dt c86 in component c86 (molar).
* RATES[86] is d/dt c87 in component c87 (molar).
* RATES[87] is d/dt c88 in component c88 (molar).
* RATES[88] is d/dt c89 in component c89 (molar).
* RATES[89] is d/dt c90 in component c90 (molar).
* RATES[90] is d/dt c91 in component c91 (molar).
* RATES[91] is d/dt c92 in component c92 (molar).
* RATES[92] is d/dt c93 in component c93 (molar).
* RATES[93] is d/dt c94 in component c94 (molar).
* RATES[94] is d/dt c95 in component c95 (molar).
* RATES[95] is d/dt c96 in component c96 (molar).
* RATES[96] is d/dt c97 in component c97 (molar).
* RATES[97] is d/dt c98 in component c98 (molar).
* RATES[98] is d/dt c99 in component c99 (molar).
* RATES[99] is d/dt c100 in component c100 (molar).
* RATES[100] is d/dt c101 in component c101 (molar).
* RATES[101] is d/dt c102 in component c102 (molar).
* RATES[102] is d/dt c103 in component c103 (molar).
* There are a total of 0 condition variables.
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 5.0E-8;
STATES[1] = 5.0E4;
STATES[2] = 0.0;
STATES[3] = 0.0;
STATES[4] = 0.0;
STATES[5] = 0.0;
STATES[6] = 0.0;
STATES[7] = 0.0;
STATES[8] = 0.0;
STATES[9] = 0.0;
STATES[10] = 0.0;
STATES[11] = 8.1E4;
STATES[12] = 0.0;
STATES[13] = 1.2E4;
STATES[14] = 0.0;
STATES[15] = 0.0;
STATES[16] = 0.0;
STATES[17] = 0.0;
STATES[18] = 0.0;
STATES[19] = 0.0;
STATES[20] = 0.0;
STATES[21] = 1.1E4;
STATES[22] = 0.0;
STATES[23] = 2.63E4;
STATES[24] = 0.0;
STATES[25] = 7.2E4;
STATES[26] = 0.0;
STATES[27] = 0.0;
STATES[28] = 0.0;
STATES[29] = 4.0E4;
STATES[30] = 1.01E5;
STATES[31] = 0.0;
STATES[32] = 0.0;
STATES[33] = 0.0;
STATES[34] = 0.0;
STATES[35] = 0.0;
STATES[36] = 0.0;
STATES[37] = 0.0;
STATES[38] = 0.0;
STATES[39] = 0.0;
STATES[40] = 4.0E4;
STATES[41] = 0.0;
STATES[42] = 0.0;
STATES[43] = 4.0E4;
STATES[44] = 0.0;
STATES[45] = 0.0;
STATES[46] = 2.1E7;
STATES[47] = 0.0;
STATES[48] = 0.0;
STATES[49] = 0.0;
STATES[50] = 0.0;
STATES[51] = 0.0;
STATES[52] = 4.0E4;
STATES[53] = 0.0;
STATES[54] = 2.21E7;
STATES[55] = 0.0;
STATES[56] = 0.0;
STATES[57] = 0.0;
STATES[58] = 0.0;
STATES[59] = 1.0E7;
STATES[60] = 0.0;
STATES[61] = 0.0;
STATES[62] = 0.0;
STATES[63] = 0.0;
STATES[64] = 0.0;
STATES[65] = 0.0;
STATES[66] = 0.0;
STATES[67] = 0.0;
STATES[68] = 0.0;
STATES[69] = 0.0;
STATES[70] = 0.0;
STATES[71] = 0.0;
STATES[72] = 0.0;
STATES[73] = 0.0;
STATES[74] = 0.0;
STATES[75] = 0.0;
STATES[76] = 0.0;
STATES[77] = 0.0;
STATES[78] = 0.0;
STATES[79] = 0.0;
STATES[80] = 0.0;
STATES[81] = 0.0;
STATES[82] = 0.0;
STATES[83] = 0.0;
STATES[84] = 0.0;
STATES[85] = 0.0;
STATES[86] = 0.0;
STATES[87] = 0.0;
STATES[88] = 0.0;
STATES[89] = 0.0;
STATES[90] = 0.0;
STATES[91] = 0.0;
STATES[92] = 0.0;
STATES[93] = 0.0;
STATES[94] = 0.0;
STATES[95] = 0.0;
STATES[96] = 0.0;
STATES[97] = 0.0;
STATES[98] = 0.0;
STATES[99] = 0.0;
STATES[100] = 0.0;
STATES[101] = 0.0;
STATES[102] = 0.0;
CONSTANTS[0] = 3.0E7;
CONSTANTS[1] = 3.84E-3;
CONSTANTS[2] = 1.66E-5;
CONSTANTS[3] = 1E-1;
CONSTANTS[4] = 1E0;
CONSTANTS[5] = 1E-2;
CONSTANTS[6] = 1.73E-7;
CONSTANTS[7] = 1.66E-3;
CONSTANTS[8] = 0E0;
CONSTANTS[9] = 1.46E-2;
CONSTANTS[10] = 5E-4;
CONSTANTS[11] = 5E-3;
CONSTANTS[12] = 1.66E-6;
CONSTANTS[13] = 2E-1;
CONSTANTS[14] = 5.43E-2;
CONSTANTS[15] = 1.1E-2;
CONSTANTS[16] = 2.17E0;
CONSTANTS[17] = 0E0;
CONSTANTS[18] = 1E4;
CONSTANTS[19] = 0E0;
CONSTANTS[20] = 1.66E-5;
CONSTANTS[21] = 2.75E-1;
CONSTANTS[22] = 1.66E-5;
CONSTANTS[23] = 6E-2;
CONSTANTS[24] = 2.5E-5;
CONSTANTS[25] = 1.3E0;
CONSTANTS[26] = 1.66E-7;
CONSTANTS[27] = 5E-1;
CONSTANTS[28] = 3.5E-6;
CONSTANTS[29] = 4E-1;
CONSTANTS[30] = 3.66E-7;
CONSTANTS[31] = 2.3E-2;
CONSTANTS[32] = 3.5E-5;
CONSTANTS[33] = 1E-1;
CONSTANTS[34] = 6E0;
CONSTANTS[35] = 6E-2;
CONSTANTS[36] = 5.5E-1;
CONSTANTS[37] = 1.66E-5;
CONSTANTS[38] = 2.14E-2;
CONSTANTS[39] = 1.66E-6;
CONSTANTS[40] = 5.3E-3;
CONSTANTS[41] = 1.17E-6;
CONSTANTS[42] = 1E0;
CONSTANTS[43] = 4E-7;
CONSTANTS[44] = 1E-1;
CONSTANTS[45] = 3.5E-5;
CONSTANTS[46] = 2E-1;
CONSTANTS[47] = 7.5E-6;
CONSTANTS[48] = 3E-2;
CONSTANTS[49] = 7.5E-6;
CONSTANTS[50] = 1.5E-3;
CONSTANTS[51] = 5E-3;
CONSTANTS[52] = 0E0;
CONSTANTS[53] = 1.5E-6;
CONSTANTS[54] = 3E-1;
CONSTANTS[55] = 5E-5;
CONSTANTS[56] = 6.4E-2;
CONSTANTS[57] = 5E-5;
CONSTANTS[58] = 4.29E-2;
CONSTANTS[59] = 1.18E-4;
CONSTANTS[60] = 2E-1;
CONSTANTS[61] = 0E0;
CONSTANTS[62] = 1E0;
CONSTANTS[63] = 1.95E-5;
CONSTANTS[64] = 3.3E-2;
CONSTANTS[65] = 0E0;
CONSTANTS[66] = 3.5E0;
CONSTANTS[67] = 0E0;
CONSTANTS[68] = 2.9E0;
CONSTANTS[69] = 2.38E-5;
CONSTANTS[70] = 8E-1;
CONSTANTS[71] = 0E0;
CONSTANTS[72] = 5.68E-2;
CONSTANTS[73] = 4.5E-7;
CONSTANTS[74] = 5E-1;
CONSTANTS[75] = 8.91E-5;
CONSTANTS[76] = 1.83E-2;
CONSTANTS[77] = 0E0;
CONSTANTS[78] = 1.6E1;
CONSTANTS[79] = 0E0;
CONSTANTS[80] = 5.7E0;
CONSTANTS[81] = 2.35E-5;
CONSTANTS[82] = 6E-1;
CONSTANTS[83] = 0E0;
CONSTANTS[84] = 2.46E-1;
CONSTANTS[85] = 8.33E-6;
CONSTANTS[86] = 5E-1;
CONSTANTS[87] = 5.5E-3;
CONSTANTS[88] = 0E0;
CONSTANTS[89] = 6.7E-4;
CONSTANTS[90] = 0E0;
CONSTANTS[91] = 1.66E-7;
CONSTANTS[92] = 2E0;
CONSTANTS[93] = 0E0;
CONSTANTS[94] = 1E-4;
CONSTANTS[95] = 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;
RATES[94] = 0.1001;
RATES[95] = 0.1001;
RATES[96] = 0.1001;
RATES[97] = 0.1001;
RATES[98] = 0.1001;
RATES[99] = 0.1001;
RATES[100] = 0.1001;
RATES[101] = 0.1001;
RATES[102] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[1] - ALGEBRAIC[12] - (ALGEBRAIC[0]+ALGEBRAIC[5]);
resid[1] = RATES[2] - ALGEBRAIC[0] - 2.00000*ALGEBRAIC[1];
resid[2] = RATES[3] - ALGEBRAIC[1] - ALGEBRAIC[2];
resid[3] = RATES[4] - ALGEBRAIC[2] - (ALGEBRAIC[6]+ALGEBRAIC[7]);
resid[4] = RATES[5] - ALGEBRAIC[5] - (ALGEBRAIC[9]+ALGEBRAIC[59]);
resid[5] = RATES[6] - ALGEBRAIC[3]+ALGEBRAIC[4];
resid[6] = RATES[7] - (ALGEBRAIC[6]+ALGEBRAIC[11]) - (ALGEBRAIC[13]+ALGEBRAIC[61]);
resid[7] = RATES[8] - - (ALGEBRAIC[4]+ALGEBRAIC[14]+ALGEBRAIC[106]+ALGEBRAIC[109]+ALGEBRAIC[112]+ALGEBRAIC[115]+ALGEBRAIC[118]+ALGEBRAIC[121]+ALGEBRAIC[124]);
resid[8] = RATES[9] - ALGEBRAIC[9] - 2.00000*ALGEBRAIC[10];
resid[9] = RATES[10] - ALGEBRAIC[10] - ALGEBRAIC[11];
resid[10] = RATES[11] - ALGEBRAIC[14] - (ALGEBRAIC[3]+ALGEBRAIC[105]+ALGEBRAIC[108]+ALGEBRAIC[111]+ALGEBRAIC[114]+ALGEBRAIC[117]+ALGEBRAIC[120]+ALGEBRAIC[123]);
resid[11] = RATES[12] - ALGEBRAIC[60];
resid[12] = RATES[13] - - (ALGEBRAIC[7]+ALGEBRAIC[13]);
resid[13] = RATES[14] - ALGEBRAIC[7] - (ALGEBRAIC[15]+ALGEBRAIC[21]+ALGEBRAIC[31]+ALGEBRAIC[33]+ALGEBRAIC[36]+ALGEBRAIC[38]+ALGEBRAIC[101]);
resid[14] = RATES[15] - - (ALGEBRAIC[9]+ALGEBRAIC[60]);
resid[15] = RATES[16] - (ALGEBRAIC[13]+ALGEBRAIC[101]) - (ALGEBRAIC[62]+ALGEBRAIC[68]+ALGEBRAIC[78]+ALGEBRAIC[79]+ALGEBRAIC[80]+ALGEBRAIC[81]+ALGEBRAIC[131]);
resid[16] = RATES[17] - (ALGEBRAIC[8]+ALGEBRAIC[62]) - (ALGEBRAIC[4]+ALGEBRAIC[63]+ALGEBRAIC[132]);
resid[17] = RATES[18] - (ALGEBRAIC[79]+ALGEBRAIC[63]+ALGEBRAIC[104]) - (ALGEBRAIC[64]+ALGEBRAIC[65]+ALGEBRAIC[66]+ALGEBRAIC[67]+ALGEBRAIC[106]+ALGEBRAIC[126]+ALGEBRAIC[133]);
resid[18] = RATES[19] - (ALGEBRAIC[64]+ALGEBRAIC[65]+ALGEBRAIC[107]) - (ALGEBRAIC[109]+ALGEBRAIC[134]);
resid[19] = RATES[20] - (ALGEBRAIC[66]+ALGEBRAIC[67]+ALGEBRAIC[110]) - (ALGEBRAIC[112]+ALGEBRAIC[135]);
resid[20] = RATES[21] - - (ALGEBRAIC[15]+ALGEBRAIC[23]+ALGEBRAIC[34]+ALGEBRAIC[37]+ALGEBRAIC[62]+ALGEBRAIC[70]);
resid[21] = RATES[22] - ALGEBRAIC[15] - (ALGEBRAIC[3]+ALGEBRAIC[8]+ALGEBRAIC[16]);
resid[22] = RATES[23] - - (ALGEBRAIC[16]+ALGEBRAIC[24]+ALGEBRAIC[34]+ALGEBRAIC[39]+ALGEBRAIC[63]+ALGEBRAIC[71]+ALGEBRAIC[129]+ALGEBRAIC[130]);
resid[23] = RATES[24] - (ALGEBRAIC[16]+ALGEBRAIC[33]) - (ALGEBRAIC[17]+ALGEBRAIC[18]+ALGEBRAIC[19]+ALGEBRAIC[20]+ALGEBRAIC[104]+ALGEBRAIC[105]+ALGEBRAIC[125]);
resid[24] = RATES[25] - - (ALGEBRAIC[17]+ALGEBRAIC[20]+ALGEBRAIC[25]+ALGEBRAIC[30]+ALGEBRAIC[64]+ALGEBRAIC[67]+ALGEBRAIC[72]+ALGEBRAIC[77]);
resid[25] = RATES[26] - (ALGEBRAIC[17]+ALGEBRAIC[18]) - (ALGEBRAIC[107]+ALGEBRAIC[108]);
resid[26] = RATES[27] - - (ALGEBRAIC[18]+ALGEBRAIC[26]+ALGEBRAIC[27]);
resid[27] = RATES[28] - (ALGEBRAIC[19]+ALGEBRAIC[20]) - (ALGEBRAIC[110]+ALGEBRAIC[111]);
resid[28] = RATES[29] - ALGEBRAIC[34] - (ALGEBRAIC[32]+ALGEBRAIC[33]+ALGEBRAIC[40]+ALGEBRAIC[79]+ALGEBRAIC[82]);
resid[29] = RATES[30] - ALGEBRAIC[35] - (ALGEBRAIC[21]+ALGEBRAIC[68]);
resid[30] = RATES[31] - ALGEBRAIC[21] - (ALGEBRAIC[22]+ALGEBRAIC[102]);
resid[31] = RATES[32] - (ALGEBRAIC[22]+ALGEBRAIC[36]) - (ALGEBRAIC[23]+ALGEBRAIC[40]+ALGEBRAIC[103]);
resid[32] = RATES[33] - (ALGEBRAIC[38]+ALGEBRAIC[23]) - (ALGEBRAIC[24]+ALGEBRAIC[113]+ALGEBRAIC[114]);
resid[33] = RATES[34] - (ALGEBRAIC[24]+ALGEBRAIC[31]+ALGEBRAIC[40]) - (ALGEBRAIC[116]+ALGEBRAIC[117]+ALGEBRAIC[127]+ALGEBRAIC[25]+ALGEBRAIC[26]+ALGEBRAIC[29]+ALGEBRAIC[30]);
resid[34] = RATES[35] - (ALGEBRAIC[26]+ALGEBRAIC[25]) - (ALGEBRAIC[119]+ALGEBRAIC[120]);
resid[35] = RATES[36] - (ALGEBRAIC[29]+ALGEBRAIC[30]) - (ALGEBRAIC[122]+ALGEBRAIC[123]);
resid[36] = RATES[37] - (ALGEBRAIC[32]+ALGEBRAIC[39]) - (ALGEBRAIC[31]+ALGEBRAIC[78]);
resid[37] = RATES[38] - ALGEBRAIC[37] - (ALGEBRAIC[38]+ALGEBRAIC[39]+ALGEBRAIC[81]);
resid[38] = RATES[39] - - (ALGEBRAIC[32]+ALGEBRAIC[35]+ALGEBRAIC[36]+ALGEBRAIC[37]+ALGEBRAIC[80]);
resid[39] = RATES[40] - - (ALGEBRAIC[27]+ALGEBRAIC[42]+ALGEBRAIC[74]+ALGEBRAIC[84]);
resid[40] = RATES[41] - ALGEBRAIC[27]+ALGEBRAIC[28];
resid[41] = RATES[42] - - (ALGEBRAIC[28]+ALGEBRAIC[19]+ALGEBRAIC[29]);
resid[42] = RATES[43] - - (ALGEBRAIC[41]+ALGEBRAIC[42]+ALGEBRAIC[83]+ALGEBRAIC[84]);
resid[43] = RATES[44] - - (ALGEBRAIC[28]+ALGEBRAIC[41]+ALGEBRAIC[43]+ALGEBRAIC[44]+ALGEBRAIC[45]+ALGEBRAIC[46]);
resid[44] = RATES[45] - ALGEBRAIC[41]+ALGEBRAIC[42];
resid[45] = RATES[46] - - (ALGEBRAIC[43]+ALGEBRAIC[50]+ALGEBRAIC[85]+ALGEBRAIC[92]);
resid[46] = RATES[47] - ALGEBRAIC[43]+ALGEBRAIC[44];
resid[47] = RATES[48] - - (ALGEBRAIC[44]+ALGEBRAIC[45]+ALGEBRAIC[48]+ALGEBRAIC[49]);
resid[48] = RATES[49] - ALGEBRAIC[45]+ALGEBRAIC[46];
resid[49] = RATES[50] - - (ALGEBRAIC[46]+ALGEBRAIC[52]+ALGEBRAIC[54]+ALGEBRAIC[47]+ALGEBRAIC[51]+ALGEBRAIC[53]);
resid[50] = RATES[51] - ALGEBRAIC[47]+ALGEBRAIC[48];
resid[51] = RATES[52] - - (ALGEBRAIC[47]+ALGEBRAIC[48]+ALGEBRAIC[49]+ALGEBRAIC[50]+ALGEBRAIC[89]+ALGEBRAIC[90]+ALGEBRAIC[91]+ALGEBRAIC[92]);
resid[52] = RATES[53] - ALGEBRAIC[49]+ALGEBRAIC[50];
resid[53] = RATES[54] - - (ALGEBRAIC[51]+ALGEBRAIC[58]+ALGEBRAIC[93]+ALGEBRAIC[100]);
resid[54] = RATES[55] - ALGEBRAIC[51]+ALGEBRAIC[52];
resid[55] = RATES[56] - - (ALGEBRAIC[52]+ALGEBRAIC[53]+ALGEBRAIC[56]+ALGEBRAIC[57]);
resid[56] = RATES[57] - ALGEBRAIC[53]+ALGEBRAIC[54];
resid[57] = RATES[58] - - (ALGEBRAIC[54]+ALGEBRAIC[55]+ALGEBRAIC[125]+ALGEBRAIC[127]+ALGEBRAIC[129]+ALGEBRAIC[142]+ALGEBRAIC[143]+ALGEBRAIC[144]);
resid[58] = RATES[59] - - (ALGEBRAIC[55]+ALGEBRAIC[56]+ALGEBRAIC[57]+ALGEBRAIC[58]+ALGEBRAIC[97]+ALGEBRAIC[98]+ALGEBRAIC[99]+ALGEBRAIC[100]);
resid[59] = RATES[60] - ALGEBRAIC[55]+ALGEBRAIC[56];
resid[60] = RATES[61] - ALGEBRAIC[57]+ALGEBRAIC[58];
resid[61] = RATES[62] - (ALGEBRAIC[68]+ALGEBRAIC[102]) - (ALGEBRAIC[69]+ALGEBRAIC[136]);
resid[62] = RATES[63] - (ALGEBRAIC[69]+ALGEBRAIC[80]+ALGEBRAIC[103]) - (ALGEBRAIC[70]+ALGEBRAIC[82]+ALGEBRAIC[137]);
resid[63] = RATES[64] - (ALGEBRAIC[70]+ALGEBRAIC[81]+ALGEBRAIC[113]) - (ALGEBRAIC[71]+ALGEBRAIC[115]+ALGEBRAIC[138]);
resid[64] = RATES[65] - (ALGEBRAIC[71]+ALGEBRAIC[78]+ALGEBRAIC[82]+ALGEBRAIC[116]) - (ALGEBRAIC[72]+ALGEBRAIC[73]+ALGEBRAIC[76]+ALGEBRAIC[77]+ALGEBRAIC[118]+ALGEBRAIC[128]+ALGEBRAIC[139]);
resid[65] = RATES[66] - (ALGEBRAIC[72]+ALGEBRAIC[73]+ALGEBRAIC[119]) - (ALGEBRAIC[121]+ALGEBRAIC[140]);
resid[66] = RATES[67] - (ALGEBRAIC[76]+ALGEBRAIC[77]+ALGEBRAIC[122]) - (ALGEBRAIC[124]+ALGEBRAIC[141]);
resid[67] = RATES[68] - - (ALGEBRAIC[65]+ALGEBRAIC[73]+ALGEBRAIC[74]);
resid[68] = RATES[69] - ALGEBRAIC[74]+ALGEBRAIC[75];
resid[69] = RATES[70] - - (ALGEBRAIC[66]+ALGEBRAIC[75]+ALGEBRAIC[76]);
resid[70] = RATES[71] - - (ALGEBRAIC[75]+ALGEBRAIC[83]+ALGEBRAIC[85]+ALGEBRAIC[86]+ALGEBRAIC[87]+ALGEBRAIC[88]);
resid[71] = RATES[72] - ALGEBRAIC[83]+ALGEBRAIC[84];
resid[72] = RATES[73] - ALGEBRAIC[85]+ALGEBRAIC[86];
resid[73] = RATES[74] - - (ALGEBRAIC[86]+ALGEBRAIC[87]+ALGEBRAIC[90]+ALGEBRAIC[91]);
resid[74] = RATES[75] - ALGEBRAIC[87]+ALGEBRAIC[88];
resid[75] = RATES[76] - - (ALGEBRAIC[88]+ALGEBRAIC[89]+ALGEBRAIC[93]+ALGEBRAIC[94]+ALGEBRAIC[95]+ALGEBRAIC[96]);
resid[76] = RATES[77] - ALGEBRAIC[89]+ALGEBRAIC[90];
resid[77] = RATES[78] - ALGEBRAIC[91]+ALGEBRAIC[92];
resid[78] = RATES[79] - ALGEBRAIC[93]+ALGEBRAIC[94];
resid[79] = RATES[80] - - (ALGEBRAIC[94]+ALGEBRAIC[95]+ALGEBRAIC[98]+ALGEBRAIC[99]);
resid[80] = RATES[81] - ALGEBRAIC[95]+ALGEBRAIC[96];
resid[81] = RATES[82] - - (ALGEBRAIC[96]+ALGEBRAIC[97]+ALGEBRAIC[126]+ALGEBRAIC[128]+ALGEBRAIC[130]+ALGEBRAIC[145]+ALGEBRAIC[146]+ALGEBRAIC[147]);
resid[82] = RATES[83] - ALGEBRAIC[97]+ALGEBRAIC[98];
resid[83] = RATES[84] - ALGEBRAIC[99]+ALGEBRAIC[100];
resid[84] = RATES[85] - ALGEBRAIC[59];
resid[85] = RATES[86] - ALGEBRAIC[61]+ALGEBRAIC[131]+ALGEBRAIC[132]+ALGEBRAIC[133]+ALGEBRAIC[134]+ALGEBRAIC[135]+ALGEBRAIC[136]+ALGEBRAIC[137]+ALGEBRAIC[138]+ALGEBRAIC[139]+ALGEBRAIC[140]+ALGEBRAIC[141];
resid[86] = RATES[87] - ALGEBRAIC[105]+ALGEBRAIC[106];
resid[87] = RATES[88] - ALGEBRAIC[108]+ALGEBRAIC[109];
resid[88] = RATES[89] - ALGEBRAIC[111]+ALGEBRAIC[112];
resid[89] = RATES[90] - ALGEBRAIC[114]+ALGEBRAIC[115];
resid[90] = RATES[91] - ALGEBRAIC[117]+ALGEBRAIC[118];
resid[91] = RATES[92] - ALGEBRAIC[120]+ALGEBRAIC[121];
resid[92] = RATES[93] - ALGEBRAIC[123]+ALGEBRAIC[124];
resid[93] = RATES[94] - ALGEBRAIC[125]+ALGEBRAIC[142];
resid[94] = RATES[95] - ALGEBRAIC[126]+ALGEBRAIC[145];
resid[95] = RATES[96] - ALGEBRAIC[127]+ALGEBRAIC[143];
resid[96] = RATES[97] - ALGEBRAIC[128]+ALGEBRAIC[146];
resid[97] = RATES[98] - - (ALGEBRAIC[142]+ALGEBRAIC[143]);
resid[98] = RATES[99] - - (ALGEBRAIC[145]+ALGEBRAIC[146]);
resid[99] = RATES[100] - ALGEBRAIC[129]+ALGEBRAIC[144];
resid[100] = RATES[101] - ALGEBRAIC[130]+ALGEBRAIC[147];
resid[101] = RATES[102] - - (ALGEBRAIC[144]+ALGEBRAIC[147]);
resid[102] = RATES[0] - CONSTANTS[95];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = CONSTANTS[0]*STATES[0]*STATES[1] - CONSTANTS[1]*STATES[2];
ALGEBRAIC[1] = CONSTANTS[2]*STATES[2]*STATES[2] - CONSTANTS[3]*STATES[3];
ALGEBRAIC[2] = CONSTANTS[4]*STATES[3]*1.00000 - CONSTANTS[5]*STATES[4];
ALGEBRAIC[3] = CONSTANTS[6]*STATES[22]*STATES[11] - CONSTANTS[7]*STATES[6];
ALGEBRAIC[4] = CONSTANTS[8]*STATES[17]*STATES[8] - CONSTANTS[9]*STATES[6];
ALGEBRAIC[5] = CONSTANTS[10]*STATES[1] - CONSTANTS[11]*STATES[5];
ALGEBRAIC[6] = CONSTANTS[10]*STATES[4] - CONSTANTS[11]*STATES[7];
ALGEBRAIC[7] = CONSTANTS[12]*STATES[4]*STATES[13] - CONSTANTS[13]*STATES[14];
ALGEBRAIC[8] = CONSTANTS[10]*STATES[22] - CONSTANTS[11]*STATES[17];
ALGEBRAIC[9] = CONSTANTS[14]*STATES[5]*STATES[15] - CONSTANTS[15]*STATES[9];
ALGEBRAIC[10] = CONSTANTS[2]*STATES[9]*STATES[9] - CONSTANTS[3]*STATES[10];
ALGEBRAIC[11] = CONSTANTS[4]*STATES[10] - CONSTANTS[5]*STATES[7];
ALGEBRAIC[12] = CONSTANTS[16]*1.00000 - CONSTANTS[17]*STATES[1];
ALGEBRAIC[13] = CONSTANTS[12]*STATES[7]*STATES[13] - CONSTANTS[13]*STATES[16];
ALGEBRAIC[14] = CONSTANTS[18]*STATES[8] - CONSTANTS[19]*STATES[11];
ALGEBRAIC[15] = CONSTANTS[20]*STATES[21]*STATES[14] - CONSTANTS[21]*STATES[22];
ALGEBRAIC[16] = CONSTANTS[22]*STATES[23]*STATES[22] - CONSTANTS[23]*STATES[24];
ALGEBRAIC[17] = CONSTANTS[24]*STATES[25]*STATES[24] - CONSTANTS[25]*STATES[26];
ALGEBRAIC[18] = CONSTANTS[26]*STATES[27]*STATES[24] - CONSTANTS[27]*STATES[26];
ALGEBRAIC[19] = CONSTANTS[28]*STATES[24]*STATES[42] - CONSTANTS[29]*STATES[28];
ALGEBRAIC[20] = CONSTANTS[30]*STATES[24]*STATES[25] - CONSTANTS[31]*STATES[28];
ALGEBRAIC[21] = CONSTANTS[32]*STATES[30]*STATES[14] - CONSTANTS[33]*STATES[31];
ALGEBRAIC[22] = CONSTANTS[34]*STATES[31] - CONSTANTS[35]*STATES[32];
ALGEBRAIC[23] = CONSTANTS[20]*STATES[21]*STATES[32] - CONSTANTS[36]*STATES[33];
ALGEBRAIC[24] = CONSTANTS[37]*STATES[23]*STATES[33] - CONSTANTS[38]*STATES[34];
ALGEBRAIC[25] = CONSTANTS[24]*STATES[25]*STATES[34] - CONSTANTS[25]*STATES[35];
ALGEBRAIC[26] = CONSTANTS[26]*STATES[34]*STATES[27] - CONSTANTS[27]*STATES[35];
ALGEBRAIC[27] = CONSTANTS[39]*STATES[27]*STATES[40] - CONSTANTS[40]*STATES[41];
ALGEBRAIC[28] = CONSTANTS[41]*STATES[42]*STATES[44] - CONSTANTS[42]*STATES[41];
ALGEBRAIC[29] = CONSTANTS[28]*STATES[34]*STATES[42] - CONSTANTS[29]*STATES[36];
ALGEBRAIC[30] = CONSTANTS[30]*STATES[34]*STATES[25] - CONSTANTS[31]*STATES[36];
ALGEBRAIC[31] = CONSTANTS[43]*STATES[37]*STATES[14] - CONSTANTS[44]*STATES[34];
ALGEBRAIC[32] = CONSTANTS[45]*STATES[39]*STATES[29] - CONSTANTS[46]*STATES[37];
ALGEBRAIC[33] = CONSTANTS[47]*STATES[14]*STATES[29] - CONSTANTS[48]*STATES[24];
ALGEBRAIC[34] = CONSTANTS[49]*STATES[23]*STATES[21] - CONSTANTS[50]*STATES[29];
ALGEBRAIC[35] = CONSTANTS[51]*STATES[39] - CONSTANTS[52]*STATES[30];
ALGEBRAIC[36] = CONSTANTS[53]*STATES[14]*STATES[39] - CONSTANTS[54]*STATES[32];
ALGEBRAIC[37] = CONSTANTS[20]*STATES[21]*STATES[39] - CONSTANTS[36]*STATES[38];
ALGEBRAIC[38] = CONSTANTS[53]*STATES[14]*STATES[38] - CONSTANTS[54]*STATES[33];
ALGEBRAIC[39] = CONSTANTS[55]*STATES[23]*STATES[38] - CONSTANTS[56]*STATES[37];
ALGEBRAIC[40] = CONSTANTS[57]*STATES[29]*STATES[32] - CONSTANTS[58]*STATES[34];
ALGEBRAIC[41] = CONSTANTS[59]*STATES[43]*STATES[44] - CONSTANTS[60]*STATES[45];
ALGEBRAIC[42] = CONSTANTS[61]*STATES[40]*STATES[43] - CONSTANTS[62]*STATES[45];
ALGEBRAIC[43] = CONSTANTS[63]*STATES[46]*STATES[44] - CONSTANTS[64]*STATES[47];
ALGEBRAIC[44] = CONSTANTS[65]*STATES[48]*STATES[44] - CONSTANTS[66]*STATES[47];
ALGEBRAIC[45] = CONSTANTS[63]*STATES[48]*STATES[44] - CONSTANTS[64]*STATES[49];
ALGEBRAIC[46] = CONSTANTS[67]*STATES[50]*STATES[44] - CONSTANTS[68]*STATES[49];
ALGEBRAIC[47] = CONSTANTS[69]*STATES[50]*STATES[52] - CONSTANTS[70]*STATES[51];
ALGEBRAIC[48] = CONSTANTS[71]*STATES[48]*STATES[52] - CONSTANTS[72]*STATES[51];
ALGEBRAIC[49] = CONSTANTS[73]*STATES[52]*STATES[48] - CONSTANTS[74]*STATES[53];
ALGEBRAIC[50] = CONSTANTS[71]*STATES[46]*STATES[52] - CONSTANTS[72]*STATES[53];
ALGEBRAIC[51] = CONSTANTS[75]*STATES[54]*STATES[50] - CONSTANTS[76]*STATES[55];
ALGEBRAIC[52] = CONSTANTS[77]*STATES[50]*STATES[56] - CONSTANTS[78]*STATES[55];
ALGEBRAIC[53] = CONSTANTS[75]*STATES[50]*STATES[56] - CONSTANTS[76]*STATES[57];
ALGEBRAIC[54] = CONSTANTS[79]*STATES[58]*STATES[50] - CONSTANTS[80]*STATES[57];
ALGEBRAIC[55] = CONSTANTS[81]*STATES[58]*STATES[59] - CONSTANTS[82]*STATES[60];
ALGEBRAIC[56] = CONSTANTS[83]*STATES[56]*STATES[59] - CONSTANTS[84]*STATES[60];
ALGEBRAIC[57] = CONSTANTS[85]*STATES[59]*STATES[56] - CONSTANTS[86]*STATES[61];
ALGEBRAIC[58] = CONSTANTS[83]*STATES[54]*STATES[59] - CONSTANTS[84]*STATES[61];
ALGEBRAIC[59] = CONSTANTS[87]*STATES[5] - CONSTANTS[88]*STATES[85];
ALGEBRAIC[60] = CONSTANTS[89]*STATES[15] - CONSTANTS[90]*STATES[12];
ALGEBRAIC[61] = CONSTANTS[87]*STATES[7] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[62] = CONSTANTS[20]*STATES[16]*STATES[21] - CONSTANTS[21]*STATES[17];
ALGEBRAIC[63] = CONSTANTS[22]*STATES[23]*STATES[17] - CONSTANTS[23]*STATES[18];
ALGEBRAIC[64] = CONSTANTS[24]*STATES[25]*STATES[18] - CONSTANTS[25]*STATES[19];
ALGEBRAIC[65] = CONSTANTS[26]*STATES[68]*STATES[18] - CONSTANTS[27]*STATES[19];
ALGEBRAIC[66] = CONSTANTS[28]*STATES[70]*STATES[18] - CONSTANTS[29]*STATES[20];
ALGEBRAIC[67] = CONSTANTS[30]*STATES[18]*STATES[25] - CONSTANTS[31]*STATES[20];
ALGEBRAIC[68] = CONSTANTS[32]*STATES[30]*STATES[16] - CONSTANTS[33]*STATES[62];
ALGEBRAIC[69] = CONSTANTS[34]*STATES[62] - CONSTANTS[35]*STATES[63];
ALGEBRAIC[70] = CONSTANTS[20]*STATES[21]*STATES[63] - CONSTANTS[36]*STATES[64];
ALGEBRAIC[71] = CONSTANTS[37]*STATES[23]*STATES[64] - CONSTANTS[38]*STATES[65];
ALGEBRAIC[72] = CONSTANTS[24]*STATES[25]*STATES[65] - CONSTANTS[25]*STATES[66];
ALGEBRAIC[73] = CONSTANTS[26]*STATES[65]*STATES[68] - CONSTANTS[27]*STATES[66];
ALGEBRAIC[74] = CONSTANTS[39]*STATES[68]*STATES[40] - CONSTANTS[40]*STATES[69];
ALGEBRAIC[75] = CONSTANTS[41]*STATES[70]*STATES[71] - CONSTANTS[42]*STATES[69];
ALGEBRAIC[76] = CONSTANTS[28]*STATES[70]*STATES[65] - CONSTANTS[29]*STATES[67];
ALGEBRAIC[77] = CONSTANTS[30]*STATES[65]*STATES[25] - CONSTANTS[31]*STATES[67];
ALGEBRAIC[78] = CONSTANTS[43]*STATES[16]*STATES[37] - CONSTANTS[44]*STATES[65];
ALGEBRAIC[79] = CONSTANTS[47]*STATES[16]*STATES[29] - CONSTANTS[48]*STATES[18];
ALGEBRAIC[80] = CONSTANTS[53]*STATES[16]*STATES[39] - CONSTANTS[54]*STATES[63];
ALGEBRAIC[81] = CONSTANTS[53]*STATES[16]*STATES[38] - CONSTANTS[54]*STATES[64];
ALGEBRAIC[82] = CONSTANTS[57]*STATES[29]*STATES[63] - CONSTANTS[58]*STATES[65];
ALGEBRAIC[83] = CONSTANTS[59]*STATES[43]*STATES[71] - CONSTANTS[60]*STATES[72];
ALGEBRAIC[84] = CONSTANTS[61]*STATES[40]*STATES[43] - CONSTANTS[62]*STATES[72];
ALGEBRAIC[85] = CONSTANTS[63]*STATES[46]*STATES[71] - CONSTANTS[64]*STATES[73];
ALGEBRAIC[86] = CONSTANTS[65]*STATES[74]*STATES[71] - CONSTANTS[66]*STATES[73];
ALGEBRAIC[87] = CONSTANTS[63]*STATES[71]*STATES[74] - CONSTANTS[64]*STATES[75];
ALGEBRAIC[88] = CONSTANTS[67]*STATES[71]*STATES[76] - CONSTANTS[68]*STATES[75];
ALGEBRAIC[89] = CONSTANTS[69]*STATES[76]*STATES[52] - CONSTANTS[70]*STATES[77];
ALGEBRAIC[90] = CONSTANTS[71]*STATES[74]*STATES[52] - CONSTANTS[72]*STATES[77];
ALGEBRAIC[91] = CONSTANTS[73]*STATES[52]*STATES[74] - CONSTANTS[74]*STATES[78];
ALGEBRAIC[92] = CONSTANTS[71]*STATES[46]*STATES[52] - CONSTANTS[72]*STATES[78];
ALGEBRAIC[93] = CONSTANTS[75]*STATES[54]*STATES[76] - CONSTANTS[76]*STATES[79];
ALGEBRAIC[94] = CONSTANTS[77]*STATES[80]*STATES[76] - CONSTANTS[78]*STATES[79];
ALGEBRAIC[95] = CONSTANTS[75]*STATES[76]*STATES[80] - CONSTANTS[76]*STATES[81];
ALGEBRAIC[96] = CONSTANTS[79]*STATES[82]*STATES[76] - CONSTANTS[80]*STATES[81];
ALGEBRAIC[97] = CONSTANTS[81]*STATES[82]*STATES[59] - CONSTANTS[82]*STATES[83];
ALGEBRAIC[98] = CONSTANTS[83]*STATES[80]*STATES[59] - CONSTANTS[84]*STATES[83];
ALGEBRAIC[99] = CONSTANTS[85]*STATES[59]*STATES[80] - CONSTANTS[86]*STATES[84];
ALGEBRAIC[100] = CONSTANTS[83]*STATES[54]*STATES[59] - CONSTANTS[84]*STATES[84];
ALGEBRAIC[101] = CONSTANTS[10]*STATES[14] - CONSTANTS[11]*STATES[16];
ALGEBRAIC[102] = CONSTANTS[10]*STATES[31] - CONSTANTS[11]*STATES[62];
ALGEBRAIC[103] = CONSTANTS[10]*STATES[32] - CONSTANTS[11]*STATES[63];
ALGEBRAIC[104] = CONSTANTS[10]*STATES[24] - CONSTANTS[11]*STATES[18];
ALGEBRAIC[105] = CONSTANTS[6]*STATES[24]*STATES[11] - CONSTANTS[7]*STATES[87];
ALGEBRAIC[106] = CONSTANTS[8]*STATES[8]*STATES[18] - CONSTANTS[9]*STATES[87];
ALGEBRAIC[107] = CONSTANTS[10]*STATES[26] - CONSTANTS[11]*STATES[19];
ALGEBRAIC[108] = CONSTANTS[6]*STATES[26]*STATES[11] - CONSTANTS[7]*STATES[88];
ALGEBRAIC[109] = CONSTANTS[8]*STATES[8]*STATES[19] - CONSTANTS[9]*STATES[88];
ALGEBRAIC[110] = CONSTANTS[10]*STATES[28] - CONSTANTS[11]*STATES[20];
ALGEBRAIC[111] = CONSTANTS[6]*STATES[28]*STATES[11] - CONSTANTS[7]*STATES[89];
ALGEBRAIC[112] = CONSTANTS[8]*STATES[8]*STATES[20] - CONSTANTS[9]*STATES[89];
ALGEBRAIC[113] = CONSTANTS[10]*STATES[33] - CONSTANTS[11]*STATES[64];
ALGEBRAIC[114] = CONSTANTS[6]*STATES[33]*STATES[11] - CONSTANTS[7]*STATES[90];
ALGEBRAIC[115] = CONSTANTS[8]*STATES[8]*STATES[64] - CONSTANTS[9]*STATES[90];
ALGEBRAIC[116] = CONSTANTS[10]*STATES[34] - CONSTANTS[11]*STATES[65];
ALGEBRAIC[117] = CONSTANTS[6]*STATES[34]*STATES[11] - CONSTANTS[7]*STATES[91];
ALGEBRAIC[118] = CONSTANTS[8]*STATES[8]*STATES[65] - CONSTANTS[9]*STATES[91];
ALGEBRAIC[119] = CONSTANTS[10]*STATES[35] - CONSTANTS[11]*STATES[66];
ALGEBRAIC[120] = CONSTANTS[6]*STATES[35]*STATES[11] - CONSTANTS[7]*STATES[92];
ALGEBRAIC[121] = CONSTANTS[8]*STATES[8]*STATES[66] - CONSTANTS[9]*STATES[92];
ALGEBRAIC[122] = CONSTANTS[10]*STATES[36] - CONSTANTS[11]*STATES[67];
ALGEBRAIC[123] = CONSTANTS[6]*STATES[36]*STATES[11] - CONSTANTS[7]*STATES[93];
ALGEBRAIC[124] = CONSTANTS[8]*STATES[67]*STATES[8] - CONSTANTS[9]*STATES[93];
ALGEBRAIC[125] = CONSTANTS[91]*STATES[58]*STATES[24] - CONSTANTS[92]*STATES[94];
ALGEBRAIC[126] = CONSTANTS[91]*STATES[82]*STATES[18] - CONSTANTS[92]*STATES[95];
ALGEBRAIC[127] = CONSTANTS[91]*STATES[58]*STATES[34] - CONSTANTS[92]*STATES[96];
ALGEBRAIC[128] = CONSTANTS[91]*STATES[82]*STATES[65] - CONSTANTS[92]*STATES[97];
ALGEBRAIC[129] = CONSTANTS[91]*STATES[58]*STATES[23] - CONSTANTS[92]*STATES[100];
ALGEBRAIC[130] = CONSTANTS[91]*STATES[82]*STATES[23] - CONSTANTS[92]*STATES[101];
ALGEBRAIC[131] = CONSTANTS[87]*STATES[16] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[132] = CONSTANTS[87]*STATES[17] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[133] = CONSTANTS[87]*STATES[18] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[134] = CONSTANTS[87]*STATES[19] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[135] = CONSTANTS[87]*STATES[20] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[136] = CONSTANTS[87]*STATES[62] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[137] = CONSTANTS[87]*STATES[63] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[138] = CONSTANTS[87]*STATES[64] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[139] = CONSTANTS[87]*STATES[65] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[140] = CONSTANTS[87]*STATES[66] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[141] = CONSTANTS[87]*STATES[67] - CONSTANTS[88]*STATES[86];
ALGEBRAIC[142] = CONSTANTS[93]*STATES[58]*STATES[98] - CONSTANTS[94]*STATES[94];
ALGEBRAIC[143] = CONSTANTS[93]*STATES[58]*STATES[98] - CONSTANTS[94]*STATES[96];
ALGEBRAIC[144] = CONSTANTS[93]*STATES[58]*STATES[102] - CONSTANTS[94]*STATES[100];
ALGEBRAIC[145] = CONSTANTS[93]*STATES[82]*STATES[99] - CONSTANTS[94]*STATES[95];
ALGEBRAIC[146] = CONSTANTS[93]*STATES[82]*STATES[99] - CONSTANTS[94]*STATES[97];
ALGEBRAIC[147] = CONSTANTS[93]*STATES[82]*STATES[102] - CONSTANTS[94]*STATES[101];
}
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;
SI[94] = 1.0;
SI[95] = 1.0;
SI[96] = 1.0;
SI[97] = 1.0;
SI[98] = 1.0;
SI[99] = 1.0;
SI[100] = 1.0;
SI[101] = 1.0;
SI[102] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
}