Generated Code

/*
   There are a total of 30 entries in the algebraic variable array.
   There are a total of 36 entries in each of the rate and state variable arrays.
   There are a total of 136 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * STATES[0] is mass in component mass (dimensionless).
 * CONSTANTS[0] is kg in component model_parameters (first_order_rate_constant).
 * STATES[1] is Cln2 in component Cln2 (dimensionless).
 * CONSTANTS[1] is ks_n2_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[2] is ks_n2__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[3] is kd_n2 in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[29] is SBF in component SBF_MBF (dimensionless).
 * STATES[2] is Clb5 in component Clb5 (dimensionless).
 * CONSTANTS[4] is ks_b5_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[5] is ks_b5__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[6] is kdi_f5 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[7] is kdi_b5 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[8] is kas_b5 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[9] is kas_f5 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[10] is kd3_f6 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[11] is kd3_c1 in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[12] is Vd_b5 in component Vd_b5 (first_order_rate_constant).
 * ALGEBRAIC[28] is MBF in component SBF_MBF (dimensionless).
 * STATES[3] is C5P in component C5P (dimensionless).
 * STATES[4] is C5 in component C5 (dimensionless).
 * STATES[5] is F5P in component F5P (dimensionless).
 * STATES[6] is F5 in component F5 (dimensionless).
 * STATES[7] is Sic1 in component Sic1 (dimensionless).
 * STATES[8] is Cdc6 in component Cdc6 (dimensionless).
 * STATES[9] is Clb2 in component Clb2 (dimensionless).
 * CONSTANTS[12] is ks_b2_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[13] is ks_b2__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[14] is kdi_b2 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[15] is kdi_f2 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[16] is kas_b2 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[17] is kas_f2 in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[13] is Vd_b2 in component Vd_b2 (first_order_rate_constant).
 * ALGEBRAIC[1] is Mcm1 in component Mcm1 (dimensionless).
 * STATES[10] is C2P in component C2P (dimensionless).
 * STATES[11] is C2 in component C2 (dimensionless).
 * STATES[12] is F2P in component F2P (dimensionless).
 * STATES[13] is F2 in component F2 (dimensionless).
 * CONSTANTS[18] is ks_c1_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[19] is ks_c1__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[20] is kpp_c1 in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[16] is Vkp_c1 in component Vkp_c1 (first_order_rate_constant).
 * STATES[14] is Swi5 in component Swi5 (dimensionless).
 * STATES[15] is Cdc14 in component Cdc14 (dimensionless).
 * STATES[16] is Sic1P in component Sic1P (dimensionless).
 * CONSTANTS[21] is ks_f6_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[22] is ks_f6__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[23] is ks_f6___ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[24] is kpp_f6 in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[17] is Vkp_f6 in component Vkp_f6 (first_order_rate_constant).
 * STATES[17] is Cdc6P in component Cdc6P (dimensionless).
 * STATES[18] is Pds1 in component Pds1 (dimensionless).
 * CONSTANTS[25] is ks_pds_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[26] is ks1_pds__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[27] is ks2_pds__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[28] is kdi_esp in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[23] is Vd_pds in component Vd_pds (first_order_rate_constant).
 * CONSTANTS[29] is kas_esp in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[11] is PE in component PE (dimensionless).
 * STATES[19] is Esp1 in component Esp1 (dimensionless).
 * STATES[20] is ORI in component ORI (dimensionless).
 * CONSTANTS[30] is ks_ori in component model_parameters (first_order_rate_constant).
 * CONSTANTS[31] is kd_ori in component model_parameters (first_order_rate_constant).
 * CONSTANTS[32] is epsilon_ori_b5 in component model_parameters (dimensionless).
 * CONSTANTS[33] is epsilon_ori_b2 in component model_parameters (dimensionless).
 * STATES[21] is BUD in component BUD (dimensionless).
 * CONSTANTS[34] is ks_bud in component model_parameters (first_order_rate_constant).
 * CONSTANTS[35] is kd_bud in component model_parameters (first_order_rate_constant).
 * CONSTANTS[36] is epsilon_bud_n2 in component model_parameters (dimensionless).
 * CONSTANTS[37] is epsilon_bud_n3 in component model_parameters (dimensionless).
 * CONSTANTS[38] is epsilon_bud_b5 in component model_parameters (dimensionless).
 * ALGEBRAIC[2] is Cln3 in component Cln3 (dimensionless).
 * STATES[22] is SPN in component SPN (dimensionless).
 * CONSTANTS[39] is ks_spn in component model_parameters (first_order_rate_constant).
 * CONSTANTS[40] is kd_spn in component model_parameters (first_order_rate_constant).
 * CONSTANTS[41] is Jspn in component model_parameters (dimensionless).
 * ALGEBRAIC[27] is G_sbf in component G_sbf (dimensionless).
 * CONSTANTS[42] is Ji_sbf in component model_parameters (dimensionless).
 * CONSTANTS[43] is Ja_sbf in component model_parameters (dimensionless).
 * ALGEBRAIC[15] is Vi_sbf in component Vi_sbf (first_order_rate_constant).
 * ALGEBRAIC[14] is Va_sbf in component Va_sbf (first_order_rate_constant).
 * ALGEBRAIC[0] is G_mcm in component G_mcm (dimensionless).
 * CONSTANTS[44] is Ji_mcm in component model_parameters (dimensionless).
 * CONSTANTS[45] is Ja_mcm in component model_parameters (dimensionless).
 * CONSTANTS[46] is ki_mcm in component model_parameters (first_order_rate_constant).
 * CONSTANTS[47] is ka_mcm in component model_parameters (first_order_rate_constant).
 * CONSTANTS[48] is C0 in component model_parameters (dimensionless).
 * CONSTANTS[49] is Dn3 in component model_parameters (dimensionless).
 * CONSTANTS[50] is Jn3 in component model_parameters (dimensionless).
 * ALGEBRAIC[3] is Bck2 in component Bck2 (dimensionless).
 * CONSTANTS[51] is B0 in component model_parameters (dimensionless).
 * ALGEBRAIC[4] is Clb5_T in component Clb5_T (dimensionless).
 * ALGEBRAIC[5] is Clb2_T in component Clb2_T (dimensionless).
 * ALGEBRAIC[6] is Sic1_T in component Sic1_T (dimensionless).
 * STATES[23] is Swi5_T in component Swi5_T (dimensionless).
 * CONSTANTS[52] is ks_swi_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[53] is ks_swi__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[54] is kd_swi in component model_parameters (first_order_rate_constant).
 * CONSTANTS[55] is ka_swi in component model_parameters (first_order_rate_constant).
 * CONSTANTS[56] is ki_swi in component model_parameters (first_order_rate_constant).
 * STATES[24] is APC_P in component APC_P (dimensionless).
 * CONSTANTS[57] is ka_apc in component model_parameters (first_order_rate_constant).
 * CONSTANTS[58] is ki_apc in component model_parameters (first_order_rate_constant).
 * CONSTANTS[59] is Ja_apc in component model_parameters (dimensionless).
 * CONSTANTS[60] is Ji_apc in component model_parameters (dimensionless).
 * STATES[25] is Cdc20_T in component Cdc20_T (dimensionless).
 * CONSTANTS[61] is ks_20_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[62] is ks_20__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[63] is kd_20 in component model_parameters (first_order_rate_constant).
 * STATES[26] is Cdc20_A in component Cdc20_A (dimensionless).
 * CONSTANTS[64] is ka_20_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[65] is ka_20__ in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[24] is kmad2 in component model_parameters (first_order_rate_constant).
 * STATES[27] is Cdh1_T in component Cdh1_T (dimensionless).
 * CONSTANTS[66] is ks_cdh in component model_parameters (first_order_rate_constant).
 * CONSTANTS[67] is kd_cdh in component model_parameters (first_order_rate_constant).
 * STATES[28] is Cdh1 in component Cdh1 (dimensionless).
 * CONSTANTS[68] is Ja_cdh in component model_parameters (dimensionless).
 * CONSTANTS[69] is Ji_cdh in component model_parameters (dimensionless).
 * ALGEBRAIC[18] is Va_cdh in component Va_cdh (first_order_rate_constant).
 * ALGEBRAIC[19] is Vi_cdh in component Vi_cdh (first_order_rate_constant).
 * STATES[29] is Tem1 in component Tem1 (dimensionless).
 * ALGEBRAIC[25] is kbub2 in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[26] is klte1 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[70] is Ja_tem in component model_parameters (dimensionless).
 * CONSTANTS[71] is Ji_tem in component model_parameters (dimensionless).
 * CONSTANTS[72] is Tem1_T in component model_parameters (dimensionless).
 * STATES[30] is Cdc15 in component Cdc15 (dimensionless).
 * CONSTANTS[73] is ka_15_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[74] is ka_15__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[75] is ka_15___ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[76] is ki_15 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[77] is Cdc15_T in component model_parameters (dimensionless).
 * STATES[31] is Cdc14_T in component Cdc14_T (dimensionless).
 * CONSTANTS[78] is ks_14 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[79] is kd_14 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[80] is kd_net in component model_parameters (first_order_rate_constant).
 * CONSTANTS[81] is kdi_rent in component model_parameters (first_order_rate_constant).
 * CONSTANTS[82] is kdi_rentp in component model_parameters (first_order_rate_constant).
 * CONSTANTS[83] is kas_rent in component model_parameters (first_order_rate_constant).
 * CONSTANTS[84] is kas_rentp in component model_parameters (first_order_rate_constant).
 * STATES[32] is RENT in component RENT (dimensionless).
 * ALGEBRAIC[9] is RENTP in component RENTP (dimensionless).
 * STATES[33] is Net1 in component Net1 (dimensionless).
 * ALGEBRAIC[10] is Net1P in component Net1P (dimensionless).
 * STATES[34] is Net1_T in component Net1_T (dimensionless).
 * CONSTANTS[85] is ks_net in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[7] is Cdc6_T in component Cdc6_T (dimensionless).
 * ALGEBRAIC[8] is CKI_T in component CKI_T (dimensionless).
 * CONSTANTS[86] is Esp1_T in component model_parameters (dimensionless).
 * CONSTANTS[87] is kd_b5_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[88] is kd_b5__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[89] is kd_b2_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[90] is kd_b2__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[91] is kd_b2p in component model_parameters (first_order_rate_constant).
 * CONSTANTS[92] is ka_sbf in component model_parameters (first_order_rate_constant).
 * CONSTANTS[93] is epsilon_sbf_n2 in component model_parameters (dimensionless).
 * CONSTANTS[94] is epsilon_sbf_n3 in component model_parameters (dimensionless).
 * CONSTANTS[95] is epsilon_sbf_b5 in component model_parameters (dimensionless).
 * CONSTANTS[96] is ki_sbf_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[97] is ki_sbf__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[98] is kd1_c1 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[99] is kd2_c1 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[100] is Jd2_c1 in component model_parameters (dimensionless).
 * CONSTANTS[101] is epsilon_c1_n2 in component model_parameters (dimensionless).
 * CONSTANTS[102] is epsilon_c1_n3 in component model_parameters (dimensionless).
 * CONSTANTS[103] is epsilon_c1_k2 in component model_parameters (dimensionless).
 * CONSTANTS[104] is epsilon_c1_b5 in component model_parameters (dimensionless).
 * CONSTANTS[105] is epsilon_c1_b2 in component model_parameters (dimensionless).
 * CONSTANTS[106] is Jd2_f6 in component model_parameters (dimensionless).
 * CONSTANTS[107] is kd1_f6 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[108] is kd2_f6 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[109] is epsilon_f6_n2 in component model_parameters (dimensionless).
 * CONSTANTS[110] is epsilon_f6_n3 in component model_parameters (dimensionless).
 * CONSTANTS[111] is epsilon_f6_k2 in component model_parameters (dimensionless).
 * CONSTANTS[112] is epsilon_f6_b5 in component model_parameters (dimensionless).
 * CONSTANTS[113] is epsilon_f6_b2 in component model_parameters (dimensionless).
 * CONSTANTS[114] is ka_cdh_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[115] is ka_cdh__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[116] is ki_cdh_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[117] is ki_cdh__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[118] is epsilon_cdh_n2 in component model_parameters (dimensionless).
 * CONSTANTS[119] is epsilon_cdh_n3 in component model_parameters (dimensionless).
 * CONSTANTS[120] is epsilon_cdh_b5 in component model_parameters (dimensionless).
 * CONSTANTS[121] is epsilon_cdh_b2 in component model_parameters (dimensionless).
 * ALGEBRAIC[20] is Vpp_net in component Vpp_net (first_order_rate_constant).
 * CONSTANTS[122] is kpp_net_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[123] is kpp_net__ in component model_parameters (first_order_rate_constant).
 * STATES[35] is PPX in component PPX (dimensionless).
 * ALGEBRAIC[21] is Vkp_net in component Vkp_net (first_order_rate_constant).
 * CONSTANTS[124] is kkp_net_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[125] is kkp_net__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[126] is ks_ppx in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[22] is Vd_ppx in component Vd_ppx (first_order_rate_constant).
 * CONSTANTS[127] is kd_ppx_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[128] is kd_ppx__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[129] is Jpds in component model_parameters (dimensionless).
 * CONSTANTS[130] is J20_ppx in component model_parameters (dimensionless).
 * CONSTANTS[131] is kd1_pds_ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[132] is kd2_pds__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[133] is kd3_pds__ in component model_parameters (first_order_rate_constant).
 * CONSTANTS[134] is Kez in component model_parameters (dimensionless).
 * CONSTANTS[135] is Kez2 in component model_parameters (dimensionless).
 * RATES[0] is d/dt mass in component mass (dimensionless).
 * RATES[1] is d/dt Cln2 in component Cln2 (dimensionless).
 * RATES[2] is d/dt Clb5 in component Clb5 (dimensionless).
 * RATES[9] is d/dt Clb2 in component Clb2 (dimensionless).
 * RATES[7] is d/dt Sic1 in component Sic1 (dimensionless).
 * RATES[16] is d/dt Sic1P in component Sic1P (dimensionless).
 * RATES[11] is d/dt C2 in component C2 (dimensionless).
 * RATES[4] is d/dt C5 in component C5 (dimensionless).
 * RATES[10] is d/dt C2P in component C2P (dimensionless).
 * RATES[3] is d/dt C5P in component C5P (dimensionless).
 * RATES[8] is d/dt Cdc6 in component Cdc6 (dimensionless).
 * RATES[17] is d/dt Cdc6P in component Cdc6P (dimensionless).
 * RATES[13] is d/dt F2 in component F2 (dimensionless).
 * RATES[18] is d/dt Pds1 in component Pds1 (dimensionless).
 * RATES[19] is d/dt Esp1 in component Esp1 (dimensionless).
 * RATES[20] is d/dt ORI in component ORI (dimensionless).
 * RATES[21] is d/dt BUD in component BUD (dimensionless).
 * RATES[22] is d/dt SPN in component SPN (dimensionless).
 * RATES[6] is d/dt F5 in component F5 (dimensionless).
 * RATES[12] is d/dt F2P in component F2P (dimensionless).
 * RATES[5] is d/dt F5P in component F5P (dimensionless).
 * RATES[23] is d/dt Swi5_T in component Swi5_T (dimensionless).
 * RATES[14] is d/dt Swi5 in component Swi5 (dimensionless).
 * RATES[24] is d/dt APC_P in component APC_P (dimensionless).
 * RATES[25] is d/dt Cdc20_T in component Cdc20_T (dimensionless).
 * RATES[26] is d/dt Cdc20_A in component Cdc20_A (dimensionless).
 * RATES[27] is d/dt Cdh1_T in component Cdh1_T (dimensionless).
 * RATES[28] is d/dt Cdh1 in component Cdh1 (dimensionless).
 * RATES[29] is d/dt Tem1 in component Tem1 (dimensionless).
 * RATES[30] is d/dt Cdc15 in component Cdc15 (dimensionless).
 * RATES[31] is d/dt Cdc14_T in component Cdc14_T (dimensionless).
 * RATES[15] is d/dt Cdc14 in component Cdc14 (dimensionless).
 * RATES[34] is d/dt Net1_T in component Net1_T (dimensionless).
 * RATES[33] is d/dt Net1 in component Net1 (dimensionless).
 * RATES[32] is d/dt RENT in component RENT (dimensionless).
 * RATES[35] is d/dt PPX in component PPX (dimensionless).
 * There are a total of 6 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 1.206;
CONSTANTS[0] = 0.007702;
STATES[1] = 0.0652;
CONSTANTS[1] = 0;
CONSTANTS[2] = 0.15;
CONSTANTS[3] = 0.12;
STATES[2] = 0.0518;
CONSTANTS[4] = 0.0008;
CONSTANTS[5] = 0.005;
CONSTANTS[6] = 0.01;
CONSTANTS[7] = 0.06;
CONSTANTS[8] = 50;
CONSTANTS[9] = 0.01;
CONSTANTS[10] = 1;
CONSTANTS[11] = 1;
STATES[3] = 0.0069;
STATES[4] = 0.0701;
STATES[5] = 7.9e-6;
STATES[6] = 7.2e-5;
STATES[7] = 0.0229;
STATES[8] = 0.1076;
STATES[9] = 0.1469;
CONSTANTS[12] = 0.001;
CONSTANTS[13] = 0.04;
CONSTANTS[14] = 0.05;
CONSTANTS[15] = 0.5;
CONSTANTS[16] = 50;
CONSTANTS[17] = 15;
STATES[10] = 0.024;
STATES[11] = 0.2384;
STATES[12] = 0.0274;
STATES[13] = 0.2361;
CONSTANTS[18] = 0.012;
CONSTANTS[19] = 0.12;
CONSTANTS[20] = 4;
STATES[14] = 0.9562;
STATES[15] = 0.4683;
STATES[16] = 0.0064;
CONSTANTS[21] = 0.024;
CONSTANTS[22] = 0.12;
CONSTANTS[23] = 0.004;
CONSTANTS[24] = 4;
STATES[17] = 0.0155;
STATES[18] = 0.0256;
CONSTANTS[25] = 0;
CONSTANTS[26] = 0.03;
CONSTANTS[27] = 0.055;
CONSTANTS[28] = 0.5;
CONSTANTS[29] = 50;
STATES[19] = 0.3013;
STATES[20] = 0.0009;
CONSTANTS[30] = 2;
CONSTANTS[31] = 0.06;
CONSTANTS[32] = 0.9;
CONSTANTS[33] = 0.45;
STATES[21] = 0.0085;
CONSTANTS[34] = 0.2;
CONSTANTS[35] = 0.06;
CONSTANTS[36] = 0.25;
CONSTANTS[37] = 0.05;
CONSTANTS[38] = 1;
STATES[22] = 0.0305;
CONSTANTS[39] = 0.1;
CONSTANTS[40] = 0.06;
CONSTANTS[41] = 0.14;
CONSTANTS[42] = 0.01;
CONSTANTS[43] = 0.01;
CONSTANTS[44] = 0.1;
CONSTANTS[45] = 0.1;
CONSTANTS[46] = 0.15;
CONSTANTS[47] = 1;
CONSTANTS[48] = 0.4;
CONSTANTS[49] = 1;
CONSTANTS[50] = 6;
CONSTANTS[51] = 0.054;
STATES[23] = 0.9765;
CONSTANTS[52] = 0.005;
CONSTANTS[53] = 0.08;
CONSTANTS[54] = 0.08;
CONSTANTS[55] = 2;
CONSTANTS[56] = 0.05;
STATES[24] = 0.1015;
CONSTANTS[57] = 0.1;
CONSTANTS[58] = 0.15;
CONSTANTS[59] = 0.1;
CONSTANTS[60] = 0.1;
STATES[25] = 1.9163;
CONSTANTS[61] = 0.006;
CONSTANTS[62] = 0.6;
CONSTANTS[63] = 0.3;
STATES[26] = 0.4443;
CONSTANTS[64] = 0.05;
CONSTANTS[65] = 0.2;
STATES[27] = 1;
CONSTANTS[66] = 0.01;
CONSTANTS[67] = 0.01;
STATES[28] = 0.9305;
CONSTANTS[68] = 0.03;
CONSTANTS[69] = 0.03;
STATES[29] = 0.9039;
CONSTANTS[70] = 0.1;
CONSTANTS[71] = 0.1;
CONSTANTS[72] = 1;
STATES[30] = 0.6565;
CONSTANTS[73] = 0.002;
CONSTANTS[74] = 1;
CONSTANTS[75] = 0.001;
CONSTANTS[76] = 0.5;
CONSTANTS[77] = 1;
STATES[31] = 2;
CONSTANTS[78] = 0.2;
CONSTANTS[79] = 0.1;
CONSTANTS[80] = 0.03;
CONSTANTS[81] = 1;
CONSTANTS[82] = 2;
CONSTANTS[83] = 200;
CONSTANTS[84] = 1;
STATES[32] = 1.0495;
STATES[33] = 0.0186;
STATES[34] = 2.8;
CONSTANTS[85] = 0.084;
CONSTANTS[86] = 1;
CONSTANTS[87] = 0.01;
CONSTANTS[88] = 0.16;
CONSTANTS[89] = 0.003;
CONSTANTS[90] = 0.4;
CONSTANTS[91] = 0.15;
CONSTANTS[92] = 0.38;
CONSTANTS[93] = 2;
CONSTANTS[94] = 10;
CONSTANTS[95] = 2;
CONSTANTS[96] = 0.6;
CONSTANTS[97] = 8;
CONSTANTS[98] = 0.01;
CONSTANTS[99] = 1;
CONSTANTS[100] = 0.05;
CONSTANTS[101] = 0.06;
CONSTANTS[102] = 0.3;
CONSTANTS[103] = 0.03;
CONSTANTS[104] = 0.1;
CONSTANTS[105] = 0.45;
CONSTANTS[106] = 0.05;
CONSTANTS[107] = 0.01;
CONSTANTS[108] = 1;
CONSTANTS[109] = 0.06;
CONSTANTS[110] = 0.3;
CONSTANTS[111] = 0.03;
CONSTANTS[112] = 0.1;
CONSTANTS[113] = 0.55;
CONSTANTS[114] = 0.01;
CONSTANTS[115] = 0.8;
CONSTANTS[116] = 0.001;
CONSTANTS[117] = 0.08;
CONSTANTS[118] = 0.4;
CONSTANTS[119] = 0.25;
CONSTANTS[120] = 8;
CONSTANTS[121] = 1.2;
CONSTANTS[122] = 0.05;
CONSTANTS[123] = 3;
STATES[35] = 0.1232;
CONSTANTS[124] = 0.01;
CONSTANTS[125] = 0.6;
CONSTANTS[126] = 0.1;
CONSTANTS[127] = 0.17;
CONSTANTS[128] = 2;
CONSTANTS[129] = 0.04;
CONSTANTS[130] = 0.15;
CONSTANTS[131] = 0.01;
CONSTANTS[132] = 0.2;
CONSTANTS[133] = 0.04;
CONSTANTS[134] = 0.3;
CONSTANTS[135] = 0.2;
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[2] = 0.1001;
RATES[9] = 0.1001;
RATES[7] = 0.1001;
RATES[16] = 0.1001;
RATES[11] = 0.1001;
RATES[4] = 0.1001;
RATES[10] = 0.1001;
RATES[3] = 0.1001;
RATES[8] = 0.1001;
RATES[17] = 0.1001;
RATES[13] = 0.1001;
RATES[18] = 0.1001;
RATES[19] = 0.1001;
RATES[20] = 0.1001;
RATES[21] = 0.1001;
RATES[22] = 0.1001;
RATES[6] = 0.1001;
RATES[12] = 0.1001;
RATES[5] = 0.1001;
RATES[23] = 0.1001;
RATES[14] = 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[15] = 0.1001;
RATES[34] = 0.1001;
RATES[33] = 0.1001;
RATES[32] = 0.1001;
RATES[35] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] -  STATES[0]*CONSTANTS[0];
resid[1] = RATES[1] -  (CONSTANTS[1]+ CONSTANTS[2]*ALGEBRAIC[29])*STATES[0] -  CONSTANTS[3]*STATES[1];
resid[2] = RATES[2] - ( (CONSTANTS[4]+ CONSTANTS[5]*ALGEBRAIC[28])*STATES[0]+ CONSTANTS[11]*STATES[3]+ CONSTANTS[7]*STATES[4]+ CONSTANTS[10]*STATES[5]+ CONSTANTS[6]*STATES[6]) -  (ALGEBRAIC[12]+ CONSTANTS[8]*STATES[7]+ CONSTANTS[9]*STATES[8])*STATES[2];
resid[3] = RATES[9] - ( (CONSTANTS[12]+ CONSTANTS[13]*ALGEBRAIC[1])*STATES[0]+ CONSTANTS[11]*STATES[10]+ CONSTANTS[14]*STATES[11]+ CONSTANTS[10]*STATES[12]+ CONSTANTS[15]*STATES[13]) -  (ALGEBRAIC[13]+ CONSTANTS[16]*STATES[7]+ CONSTANTS[17]*STATES[8])*STATES[9];
resid[4] = RATES[7] - (CONSTANTS[18]+ CONSTANTS[19]*STATES[14]+ (ALGEBRAIC[13]+CONSTANTS[14])*STATES[11]+ (ALGEBRAIC[12]+CONSTANTS[7])*STATES[4]+ CONSTANTS[20]*STATES[15]*STATES[16]) -  ( CONSTANTS[16]*STATES[9]+ CONSTANTS[8]*STATES[2]+ALGEBRAIC[16])*STATES[7];
resid[5] = RATES[16] - ( ALGEBRAIC[16]*STATES[7] -  ( CONSTANTS[20]*STATES[15]+CONSTANTS[11])*STATES[16])+ ALGEBRAIC[13]*STATES[10]+ ALGEBRAIC[12]*STATES[3];
resid[6] = RATES[11] - ( CONSTANTS[16]*STATES[9]*STATES[7]+ CONSTANTS[20]*STATES[15]*STATES[10]) -  (CONSTANTS[14]+ALGEBRAIC[13]+ALGEBRAIC[16])*STATES[11];
resid[7] = RATES[4] - ( CONSTANTS[8]*STATES[2]*STATES[7]+ CONSTANTS[20]*STATES[15]*STATES[3]) -  (CONSTANTS[7]+ALGEBRAIC[12]+ALGEBRAIC[16])*STATES[4];
resid[8] = RATES[10] -  ALGEBRAIC[16]*STATES[11] -  ( CONSTANTS[20]*STATES[15]+CONSTANTS[11]+ALGEBRAIC[13])*STATES[10];
resid[9] = RATES[3] -  ALGEBRAIC[16]*STATES[4] -  ( CONSTANTS[20]*STATES[15]+CONSTANTS[11]+ALGEBRAIC[12])*STATES[3];
resid[10] = RATES[8] - (CONSTANTS[21]+ CONSTANTS[22]*STATES[14]+ CONSTANTS[23]*ALGEBRAIC[29]+ (ALGEBRAIC[13]+CONSTANTS[15])*STATES[13]+ (ALGEBRAIC[12]+CONSTANTS[6])*STATES[6]+ CONSTANTS[24]*STATES[15]*STATES[17]) -  ( CONSTANTS[17]*STATES[9]+ CONSTANTS[9]*STATES[2]+ALGEBRAIC[17])*STATES[8];
resid[11] = RATES[17] - ( ALGEBRAIC[17]*STATES[8] -  ( CONSTANTS[24]*STATES[15]+CONSTANTS[10])*STATES[17])+ ALGEBRAIC[13]*STATES[12]+ ALGEBRAIC[12]*STATES[5];
resid[12] = RATES[13] - ( CONSTANTS[17]*STATES[9]*STATES[8]+ CONSTANTS[24]*STATES[15]*STATES[12]) -  (CONSTANTS[15]+ALGEBRAIC[13]+ALGEBRAIC[17])*STATES[13];
resid[13] = RATES[18] - (CONSTANTS[25]+ CONSTANTS[26]*ALGEBRAIC[29]+ CONSTANTS[27]*ALGEBRAIC[1]+ CONSTANTS[28]*ALGEBRAIC[11]) -  (ALGEBRAIC[23]+ CONSTANTS[29]*STATES[19])*STATES[18];
resid[14] = RATES[19] -  - CONSTANTS[29]*STATES[18]*STATES[19]+ (CONSTANTS[28]+ALGEBRAIC[23])*ALGEBRAIC[11];
resid[15] = RATES[20] -  CONSTANTS[30]*( CONSTANTS[32]*STATES[2]+ CONSTANTS[33]*STATES[9]) -  CONSTANTS[31]*STATES[20];
resid[16] = RATES[21] -  CONSTANTS[34]*( CONSTANTS[36]*STATES[1]+ CONSTANTS[37]*ALGEBRAIC[2]+ CONSTANTS[38]*STATES[2]) -  CONSTANTS[35]*STATES[21];
resid[17] = RATES[22] - ( CONSTANTS[39]*STATES[9])/(CONSTANTS[41]+STATES[9]) -  CONSTANTS[40]*STATES[22];
resid[18] = RATES[6] - ( CONSTANTS[9]*STATES[2]*STATES[8]+ CONSTANTS[24]*STATES[15]*STATES[5]) -  (CONSTANTS[6]+ALGEBRAIC[12]+ALGEBRAIC[17])*STATES[6];
resid[19] = RATES[12] -  ALGEBRAIC[17]*STATES[13] -  ( CONSTANTS[24]*STATES[15]+CONSTANTS[10]+ALGEBRAIC[13])*STATES[12];
resid[20] = RATES[5] -  ALGEBRAIC[17]*STATES[6] -  ( CONSTANTS[24]*STATES[15]+CONSTANTS[10]+ALGEBRAIC[12])*STATES[5];
resid[21] = RATES[23] - (CONSTANTS[52]+ CONSTANTS[53]*ALGEBRAIC[1]) -  CONSTANTS[54]*STATES[23];
resid[22] = RATES[14] - (CONSTANTS[52]+ CONSTANTS[53]*ALGEBRAIC[1]+ CONSTANTS[55]*STATES[15]*(STATES[23] - STATES[14])) -  (CONSTANTS[54]+ CONSTANTS[56]*STATES[9])*STATES[14];
resid[23] = RATES[24] - ( CONSTANTS[57]*STATES[9]*(1.00000 - STATES[24]))/((CONSTANTS[59]+1.00000) - STATES[24]) - ( CONSTANTS[58]*STATES[24])/(CONSTANTS[60]+STATES[24]);
resid[24] = RATES[25] - (CONSTANTS[61]+ CONSTANTS[62]*ALGEBRAIC[1]) -  CONSTANTS[63]*STATES[25];
resid[25] = RATES[26] -  (CONSTANTS[64]+ CONSTANTS[65]*STATES[24])*(STATES[25] - STATES[26]) -  (ALGEBRAIC[24]+CONSTANTS[63])*STATES[26];
resid[26] = RATES[27] - CONSTANTS[66] -  CONSTANTS[67]*STATES[27];
resid[27] = RATES[28] - (CONSTANTS[66]+( ALGEBRAIC[18]*(STATES[27] - STATES[28]))/((CONSTANTS[68]+STATES[27]) - STATES[28])) - ( CONSTANTS[67]*STATES[28]+( ALGEBRAIC[19]*STATES[28])/(CONSTANTS[69]+STATES[28]));
resid[28] = RATES[29] - ( ALGEBRAIC[26]*(CONSTANTS[72] - STATES[29]))/((CONSTANTS[70]+CONSTANTS[72]) - STATES[29]) - ( ALGEBRAIC[25]*STATES[29])/(CONSTANTS[71]+STATES[29]);
resid[29] = RATES[30] -  ( CONSTANTS[73]*(CONSTANTS[72] - STATES[29])+ CONSTANTS[74]*STATES[29]+ CONSTANTS[75]*STATES[15])*(CONSTANTS[77] - STATES[30]) -  CONSTANTS[76]*STATES[30];
resid[30] = RATES[31] - CONSTANTS[78] -  CONSTANTS[79]*STATES[31];
resid[31] = RATES[15] - (CONSTANTS[78]+ CONSTANTS[80]*(STATES[32]+ALGEBRAIC[9])+ CONSTANTS[81]*STATES[32]+ CONSTANTS[82]*ALGEBRAIC[9]) - ( CONSTANTS[79]*STATES[15]+ ( CONSTANTS[83]*STATES[33]+ CONSTANTS[84]*ALGEBRAIC[10])*STATES[15]);
resid[32] = RATES[34] - CONSTANTS[85] -  CONSTANTS[80]*STATES[34];
resid[33] = RATES[33] - (CONSTANTS[85]+ CONSTANTS[79]*STATES[32]+ CONSTANTS[81]*STATES[32]+ ALGEBRAIC[20]*ALGEBRAIC[10]) - ( CONSTANTS[80]*STATES[33]+ CONSTANTS[83]*STATES[15]*STATES[33]+ ALGEBRAIC[21]*STATES[33]);
resid[34] = RATES[32] - ( CONSTANTS[83]*STATES[15]*STATES[33]+ ALGEBRAIC[20]*ALGEBRAIC[9]) - ( (CONSTANTS[79]+CONSTANTS[80])*STATES[32]+ CONSTANTS[81]*STATES[32]+ ALGEBRAIC[21]*STATES[32]);
resid[35] = RATES[35] - CONSTANTS[126] -  ALGEBRAIC[22]*STATES[35];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[4] = STATES[2]+STATES[4]+STATES[3]+STATES[6]+STATES[5];
ALGEBRAIC[5] = STATES[9]+STATES[11]+STATES[10]+STATES[13]+STATES[12];
ALGEBRAIC[8] = ALGEBRAIC[6]+ALGEBRAIC[7];
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = ( 2.00000*CONSTANTS[44]*CONSTANTS[47]*STATES[9])/((CONSTANTS[46]+ CONSTANTS[45]*CONSTANTS[46]+ CONSTANTS[44]*CONSTANTS[47]*STATES[9]+ pow((pow((CONSTANTS[46]+ CONSTANTS[45]*CONSTANTS[46]+ CONSTANTS[44]*CONSTANTS[47]*STATES[9]) -  CONSTANTS[47]*STATES[9], 2.00000) -  4.00000*(CONSTANTS[46] -  CONSTANTS[47]*STATES[9])*CONSTANTS[44]*CONSTANTS[47]*STATES[9]), 1.0 / 2)) -  CONSTANTS[47]*STATES[9]);
ALGEBRAIC[1] = ALGEBRAIC[0];
ALGEBRAIC[2] = ( CONSTANTS[48]*CONSTANTS[49]*STATES[0])/(CONSTANTS[50]+ CONSTANTS[49]*STATES[0]);
ALGEBRAIC[9] = STATES[31] - (STATES[32]+STATES[15]);
ALGEBRAIC[10] = (STATES[34]+STATES[15]) - (STATES[33]+STATES[31]);
ALGEBRAIC[11] = CONSTANTS[86] - STATES[19];
ALGEBRAIC[12] = CONSTANTS[87]+ CONSTANTS[88]*STATES[26];
ALGEBRAIC[13] = CONSTANTS[89]+ CONSTANTS[90]*STATES[28]+ CONSTANTS[91]*STATES[26];
ALGEBRAIC[3] =  CONSTANTS[51]*STATES[0];
ALGEBRAIC[6] = STATES[7]+STATES[16]+STATES[11]+STATES[10]+STATES[4]+STATES[3];
ALGEBRAIC[16] = CONSTANTS[98]+( CONSTANTS[99]*( CONSTANTS[102]*ALGEBRAIC[2]+ CONSTANTS[103]*ALGEBRAIC[3]+ CONSTANTS[101]*STATES[1]+ CONSTANTS[104]*STATES[2]+ CONSTANTS[105]*STATES[9]))/(CONSTANTS[100]+ALGEBRAIC[6]);
ALGEBRAIC[7] = STATES[8]+STATES[17]+STATES[13]+STATES[12]+STATES[6]+STATES[5];
ALGEBRAIC[17] = CONSTANTS[107]+( CONSTANTS[108]*( CONSTANTS[110]*ALGEBRAIC[2]+ CONSTANTS[111]*ALGEBRAIC[3]+ CONSTANTS[109]*STATES[1]+ CONSTANTS[112]*STATES[2]+ CONSTANTS[113]*STATES[9]))/(CONSTANTS[106]+ALGEBRAIC[7]);
ALGEBRAIC[18] = CONSTANTS[114]+ CONSTANTS[115]*STATES[15];
ALGEBRAIC[19] = CONSTANTS[116]+ CONSTANTS[117]*( CONSTANTS[119]*ALGEBRAIC[2]+ CONSTANTS[118]*STATES[1]+ CONSTANTS[120]*STATES[2]+ CONSTANTS[121]*STATES[9]);
ALGEBRAIC[20] = CONSTANTS[122]+ CONSTANTS[123]*STATES[35];
ALGEBRAIC[21] =  (CONSTANTS[124]+ CONSTANTS[125]*STATES[30])*STATES[0];
ALGEBRAIC[22] = CONSTANTS[127]+( CONSTANTS[128]*(CONSTANTS[130]+STATES[26])*CONSTANTS[129])/(CONSTANTS[129]+STATES[18]);
ALGEBRAIC[23] = CONSTANTS[131]+ CONSTANTS[132]*STATES[26]+ CONSTANTS[133]*STATES[28];
ALGEBRAIC[24] = (CONDVAR[0]>0.00000&&CONDVAR[1]<0.00000 ? 8.00000 : 0.0100000);
ALGEBRAIC[25] = (CONDVAR[2]>0.00000&&CONDVAR[3]<0.00000 ? 1.00000 : 0.200000);
ALGEBRAIC[26] = (CONDVAR[4]>0.00000&&CONDVAR[5]>0.00000 ? 1.00000 : 0.100000);
ALGEBRAIC[15] = CONSTANTS[96]+ CONSTANTS[97]*STATES[9];
ALGEBRAIC[14] =  CONSTANTS[92]*( CONSTANTS[93]*STATES[1]+ CONSTANTS[94]*(ALGEBRAIC[2]+ALGEBRAIC[3])+ CONSTANTS[95]*STATES[2]);
ALGEBRAIC[27] = ( 2.00000*CONSTANTS[42]*ALGEBRAIC[14])/((ALGEBRAIC[15]+ CONSTANTS[43]*ALGEBRAIC[15]+ CONSTANTS[42]*ALGEBRAIC[14]+ pow((pow((ALGEBRAIC[15]+ CONSTANTS[43]*ALGEBRAIC[15]+ CONSTANTS[42]*ALGEBRAIC[14]) - ALGEBRAIC[14], 2.00000) -  4.00000*(ALGEBRAIC[15] - ALGEBRAIC[14])*CONSTANTS[42]*ALGEBRAIC[14]), 1.0 / 2)) - ALGEBRAIC[14]);
ALGEBRAIC[28] = ALGEBRAIC[27];
ALGEBRAIC[29] = ALGEBRAIC[28];
}
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;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = STATES[20] - 1.00000;
CONDVAR[1] = STATES[22] - 1.00000;
CONDVAR[2] = STATES[20] - 1.00000;
CONDVAR[3] = STATES[22] - 1.00000;
CONDVAR[4] = STATES[22] - 1.00000;
CONDVAR[5] = STATES[9] - CONSTANTS[134];
}