/* There are a total of 28 entries in the algebraic variable array. There are a total of 14 entries in each of the rate and state variable arrays. There are a total of 98 entries in the constant variable array. */ /* * VOI is time in component environment (minute). * STATES[0] is actCycA in component actCycA (dimensionless). * CONSTANTS[0] is k_asa in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[1] is k_sap in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[2] is k_sapp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[3] is k_dia in component kinetic_parameters (first_order_rate_constant). * ALGEBRAIC[12] is V_di in component V_di (first_order_rate_constant). * ALGEBRAIC[9] is V_da in component V_da (first_order_rate_constant). * STATES[1] is mass in component mass (dimensionless). * ALGEBRAIC[1] is Tri_A in component Tri_A (dimensionless). * ALGEBRAIC[7] is TF_E in component TF_E (dimensionless). * ALGEBRAIC[3] is freeCKI in component freeCKI (dimensionless). * STATES[2] is actCycB in component actCycB (dimensionless). * CONSTANTS[4] is k_asb in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[5] is k_dib in component kinetic_parameters (first_order_rate_constant). * ALGEBRAIC[15] is V_sb in component V_sb (first_order_rate_constant). * ALGEBRAIC[16] is V_db in component V_db (first_order_rate_constant). * ALGEBRAIC[24] is V_25 in component V_25 (first_order_rate_constant). * ALGEBRAIC[20] is V_wee in component V_wee (first_order_rate_constant). * STATES[3] is CycB in component CycB (dimensionless). * STATES[4] is Tri_B in component Tri_B (dimensionless). * STATES[5] is preMPF in component preMPF (dimensionless). * STATES[6] is actCycE in component actCycE (dimensionless). * CONSTANTS[6] is k_die in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[7] is k_ase in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[8] is k_sep in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[9] is k_sepp in component kinetic_parameters (first_order_rate_constant). * ALGEBRAIC[8] is V_de in component V_de (first_order_rate_constant). * ALGEBRAIC[2] is Tri_E in component Tri_E (dimensionless). * STATES[7] is CycA in component CycA (dimensionless). * STATES[8] is CycE in component CycE (dimensionless). * ALGEBRAIC[0] is CycD in component CycD (dimensionless). * CONSTANTS[10] is CycD_0 in component CycD (dimensionless). * STATES[9] is Cdc20_A in component Cdc20_A (dimensionless). * CONSTANTS[11] is k_a20 in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[12] is k_i20 in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[13] is k_d20 in component kinetic_parameters (first_order_rate_constant). * STATES[10] is APCP in component APCP (dimensionless). * STATES[11] is Cdc20_T in component Cdc20_T (dimensionless). * CONSTANTS[14] is J_a20 in component Cdc20_A (dimensionless). * CONSTANTS[15] is J_i20 in component Cdc20_A (dimensionless). * CONSTANTS[16] is k_s20p in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[17] is k_s20pp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[18] is n in component Cdc20_T (dimensionless). * CONSTANTS[19] is J_20 in component Cdc20_T (dimensionless). * CONSTANTS[20] is k_aAPC in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[21] is k_iAPC in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[22] is J_aAPC in component APCP (dimensionless). * CONSTANTS[23] is J_iAPC in component APCP (dimensionless). * STATES[12] is CKI in component CKI (dimensionless). * STATES[13] is Cdh1 in component Cdh1 (dimensionless). * CONSTANTS[24] is k_ah1p in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[25] is k_ah1pp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[26] is k_ih1p in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[27] is k_ih1pp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[28] is k_ih1ppp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[29] is k_ih1pppp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[30] is k_ih1ppppp in component kinetic_parameters (first_order_rate_constant). * ALGEBRAIC[11] is Cdc14 in component TF_I (dimensionless). * CONSTANTS[31] is J_ah1 in component Cdh1 (dimensionless). * CONSTANTS[32] is J_ih1 in component Cdh1 (dimensionless). * CONSTANTS[33] is mu in component mass (first_order_rate_constant). * CONSTANTS[34] is maxmass in component mass (dimensionless). * ALGEBRAIC[4] is V_atf in component V_atf (first_order_rate_constant). * CONSTANTS[35] is k_atfp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[36] is k_atfpp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[37] is k_atfppp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[38] is k_atfpppp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[39] is k_itfp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[40] is k_itfpp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[41] is k_itfppp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[42] is J_itf in component TF_E (dimensionless). * CONSTANTS[43] is J_atf in component TF_E (dimensionless). * ALGEBRAIC[5] is A1 in component TF_E (first_order_rate_constant). * ALGEBRAIC[6] is A2 in component TF_E (first_order_rate_constant). * CONSTANTS[88] is A3 in component TF_E (dimensionless). * CONSTANTS[91] is A4 in component TF_E (dimensionless). * CONSTANTS[44] is k_dep in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[45] is k_depp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[46] is k_deppp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[47] is k_depppp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[48] is k_dap in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[49] is k_dapp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[50] is k_dappp in component kinetic_parameters (first_order_rate_constant). * ALGEBRAIC[26] is TF_I in component TF_I (dimensionless). * CONSTANTS[51] is k_afi in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[52] is k_ifip in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[53] is k_ifipp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[54] is J_ifi in component TF_I (dimensionless). * CONSTANTS[55] is J_afi in component TF_I (dimensionless). * ALGEBRAIC[25] is A1 in component TF_I (first_order_rate_constant). * ALGEBRAIC[10] is A2 in component TF_I (first_order_rate_constant). * CONSTANTS[87] is A3 in component TF_I (dimensionless). * CONSTANTS[90] is A4 in component TF_I (dimensionless). * ALGEBRAIC[27] is V_si in component V_si (first_order_rate_constant). * CONSTANTS[56] is k_sip in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[57] is k_sipp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[58] is k_dip in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[59] is k_dipp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[60] is k_dippp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[61] is k_dipppp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[62] is k_dippppp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[63] is J_14di in component kinetic_parameters (dimensionless). * ALGEBRAIC[14] is TF_B in component TF_B (dimensionless). * CONSTANTS[64] is k_afb in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[65] is k_ifb in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[66] is J_ifb in component TF_B (dimensionless). * CONSTANTS[67] is J_afb in component TF_B (dimensionless). * ALGEBRAIC[13] is A1 in component TF_B (first_order_rate_constant). * CONSTANTS[89] is A2 in component TF_B (first_order_rate_constant). * CONSTANTS[92] is A3 in component TF_B (dimensionless). * CONSTANTS[94] is A4 in component TF_B (dimensionless). * CONSTANTS[68] is k_sbp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[69] is k_sbpp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[70] is k_dbp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[71] is k_dbpp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[72] is k_dbppp in component kinetic_parameters (first_order_rate_constant). * ALGEBRAIC[19] is Wee1 in component Wee1 (dimensionless). * CONSTANTS[73] is k_aweep in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[74] is k_aweepp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[75] is k_iwee in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[76] is J_iwee in component Wee1 (dimensionless). * CONSTANTS[77] is J_awee in component Wee1 (dimensionless). * ALGEBRAIC[17] is A1 in component Wee1 (first_order_rate_constant). * ALGEBRAIC[18] is A2 in component Wee1 (first_order_rate_constant). * CONSTANTS[93] is A3 in component Wee1 (dimensionless). * CONSTANTS[95] is A4 in component Wee1 (dimensionless). * CONSTANTS[78] is k_weep in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[79] is k_weepp in component kinetic_parameters (first_order_rate_constant). * ALGEBRAIC[23] is Cdc25 in component Cdc25 (dimensionless). * CONSTANTS[80] is k_a25 in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[81] is k_i25p in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[82] is k_i25pp in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[83] is J_i25 in component Cdc25 (dimensionless). * CONSTANTS[84] is J_a25 in component Cdc25 (dimensionless). * ALGEBRAIC[21] is A1 in component Cdc25 (first_order_rate_constant). * ALGEBRAIC[22] is A2 in component Cdc25 (first_order_rate_constant). * CONSTANTS[96] is A3 in component Cdc25 (dimensionless). * CONSTANTS[97] is A4 in component Cdc25 (dimensionless). * CONSTANTS[85] is k_25p in component kinetic_parameters (first_order_rate_constant). * CONSTANTS[86] is k_25pp in component kinetic_parameters (first_order_rate_constant). * RATES[0] is d/dt actCycA in component actCycA (dimensionless). * RATES[2] is d/dt actCycB in component actCycB (dimensionless). * RATES[6] is d/dt actCycE in component actCycE (dimensionless). * RATES[7] is d/dt CycA in component CycA (dimensionless). * RATES[3] is d/dt CycB in component CycB (dimensionless). * RATES[8] is d/dt CycE in component CycE (dimensionless). * RATES[9] is d/dt Cdc20_A in component Cdc20_A (dimensionless). * RATES[11] is d/dt Cdc20_T in component Cdc20_T (dimensionless). * RATES[10] is d/dt APCP in component APCP (dimensionless). * RATES[13] is d/dt Cdh1 in component Cdh1 (dimensionless). * RATES[1] is d/dt mass in component mass (dimensionless). * RATES[5] is d/dt preMPF in component preMPF (dimensionless). * RATES[4] is d/dt Tri_B in component Tri_B (dimensionless). * RATES[12] is d/dt CKI in component CKI (dimensionless). * There are a total of 0 condition variables. */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { STATES[0] = 0.09450308233499527; CONSTANTS[0] = 50; CONSTANTS[1] = 0.0015; CONSTANTS[2] = 0.01; CONSTANTS[3] = 0.06; STATES[1] = 1.338267803192139; STATES[2] = 0.1903585940599442; CONSTANTS[4] = 60; CONSTANTS[5] = 0.05; STATES[3] = 0.09450308233499527; STATES[4] = 0.01; STATES[5] = 0.01; STATES[6] = 0.2092024385929108; CONSTANTS[6] = 0; CONSTANTS[7] = 0; CONSTANTS[8] = 0; CONSTANTS[9] = 0.15; STATES[7] = 0.09450308233499527; STATES[8] = 0.09450308233499527; CONSTANTS[10] = 0.108; STATES[9] = 0.3572723865509033; CONSTANTS[11] = 1; CONSTANTS[12] = 0.16; CONSTANTS[13] = 0.05; STATES[10] = 0.07591402530670166; STATES[11] = 0.7702387571334839; CONSTANTS[14] = 1; CONSTANTS[15] = 1; CONSTANTS[16] = 0.001; CONSTANTS[17] = 1; CONSTANTS[18] = 1; CONSTANTS[19] = 10; CONSTANTS[20] = 0.1; CONSTANTS[21] = 0.15; CONSTANTS[22] = 0.1; CONSTANTS[23] = 0.1; STATES[12] = 0.01; STATES[13] = 0.7189393639564514; CONSTANTS[24] = 0.02; CONSTANTS[25] = 0.8; CONSTANTS[26] = 0.001; CONSTANTS[27] = 0.35; CONSTANTS[28] = 0.1; CONSTANTS[29] = 0.06; CONSTANTS[30] = 0.005; CONSTANTS[31] = 0.03; CONSTANTS[32] = 0.03; CONSTANTS[33] = 0.005776; CONSTANTS[34] = 10000; CONSTANTS[35] = 0; CONSTANTS[36] = 1.5; CONSTANTS[37] = 0.38; CONSTANTS[38] = 3; CONSTANTS[39] = 0.75; CONSTANTS[40] = 8; CONSTANTS[41] = 0; CONSTANTS[42] = 0.01; CONSTANTS[43] = 0.01; CONSTANTS[44] = 0.12; CONSTANTS[45] = 0; CONSTANTS[46] = 0; CONSTANTS[47] = 0; CONSTANTS[48] = 0.01; CONSTANTS[49] = 0.16; CONSTANTS[50] = 0; CONSTANTS[51] = 6; CONSTANTS[52] = 0.008; CONSTANTS[53] = 0.05; CONSTANTS[54] = 2; CONSTANTS[55] = 1; CONSTANTS[56] = 0.018; CONSTANTS[57] = 0.18; CONSTANTS[58] = 0.02; CONSTANTS[59] = 0.1; CONSTANTS[60] = 0.8; CONSTANTS[61] = 0.12; CONSTANTS[62] = 0.1; CONSTANTS[63] = 12; CONSTANTS[64] = 1; CONSTANTS[65] = 0.15; CONSTANTS[66] = 0.1; CONSTANTS[67] = 0.1; CONSTANTS[68] = 0.004; CONSTANTS[69] = 0.04; CONSTANTS[70] = 0.003; CONSTANTS[71] = 0.4; CONSTANTS[72] = 0.15; CONSTANTS[73] = 0.3; CONSTANTS[74] = 0; CONSTANTS[75] = 1; CONSTANTS[76] = 0.05; CONSTANTS[77] = 0.05; CONSTANTS[78] = 0.02; CONSTANTS[79] = 0.2; CONSTANTS[80] = 0; CONSTANTS[81] = 0.3; CONSTANTS[82] = 0; CONSTANTS[83] = 0.1; CONSTANTS[84] = 0.1; CONSTANTS[85] = 0.01; CONSTANTS[86] = 5; CONSTANTS[87] = CONSTANTS[55]; CONSTANTS[88] = CONSTANTS[43]; CONSTANTS[89] = CONSTANTS[65]; CONSTANTS[90] = CONSTANTS[54]; CONSTANTS[91] = CONSTANTS[42]; CONSTANTS[92] = CONSTANTS[67]; CONSTANTS[93] = CONSTANTS[77]; CONSTANTS[94] = CONSTANTS[66]; CONSTANTS[95] = CONSTANTS[76]; CONSTANTS[96] = CONSTANTS[84]; CONSTANTS[97] = CONSTANTS[83]; RATES[0] = 0.1001; RATES[2] = 0.1001; RATES[6] = 0.1001; RATES[7] = 0.1001; RATES[3] = 0.1001; RATES[8] = 0.1001; RATES[9] = 0.1001; RATES[11] = 0.1001; RATES[10] = 0.1001; RATES[13] = 0.1001; RATES[1] = 0.1001; RATES[5] = 0.1001; RATES[4] = 0.1001; RATES[12] = 0.1001; } void computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { resid[0] = RATES[0] - ( (CONSTANTS[1]+ CONSTANTS[2]*ALGEBRAIC[7])*STATES[1]+ (ALGEBRAIC[12]+CONSTANTS[3])*ALGEBRAIC[1]) - (ALGEBRAIC[9]+ CONSTANTS[0]*ALGEBRAIC[3])*STATES[0]; resid[1] = RATES[2] - ( ALGEBRAIC[15]*STATES[1]+ ALGEBRAIC[24]*(STATES[3] - (STATES[4]+STATES[2]))+ (CONSTANTS[5]+ALGEBRAIC[12])*(STATES[3] - (STATES[5]+STATES[2]))) - (ALGEBRAIC[16]+ALGEBRAIC[20]+ CONSTANTS[4]*ALGEBRAIC[3])*STATES[2]; resid[2] = RATES[6] - ( (CONSTANTS[8]+ CONSTANTS[9]*ALGEBRAIC[7])*STATES[1]+ (ALGEBRAIC[12]+CONSTANTS[6])*ALGEBRAIC[2]) - (ALGEBRAIC[8]+ CONSTANTS[7]*ALGEBRAIC[3])*STATES[6]; resid[3] = RATES[7] - (CONSTANTS[1]+ CONSTANTS[2]*ALGEBRAIC[7])*STATES[1] - ALGEBRAIC[9]*STATES[7]; resid[4] = RATES[3] - ALGEBRAIC[15]*STATES[1] - ALGEBRAIC[16]*STATES[3]; resid[5] = RATES[8] - (CONSTANTS[8]+ CONSTANTS[9]*ALGEBRAIC[7])*STATES[1] - ALGEBRAIC[8]*STATES[8]; resid[6] = RATES[9] - ( CONSTANTS[11]*STATES[10]*(STATES[11] - STATES[9]))/((CONSTANTS[14]+STATES[11]) - STATES[9]) - (CONSTANTS[12]/(CONSTANTS[15]+STATES[9])+CONSTANTS[13])*STATES[9]; resid[7] = RATES[11] - (CONSTANTS[16]+ CONSTANTS[17]*pow(STATES[2], CONSTANTS[18]))/(pow(CONSTANTS[19], CONSTANTS[18])+pow(STATES[2], CONSTANTS[18])) - CONSTANTS[13]*STATES[11]; resid[8] = RATES[10] - ( CONSTANTS[20]*STATES[2]*(1.00000 - STATES[10]))/((CONSTANTS[22]+1.00000) - STATES[10]) - ( CONSTANTS[21]*STATES[10])/(CONSTANTS[23]+STATES[10]); resid[9] = RATES[13] - ( (CONSTANTS[24]+ CONSTANTS[25]*ALGEBRAIC[11])*(1.00000 - STATES[13]))/((CONSTANTS[31]+1.00000) - STATES[13]) - ( (CONSTANTS[26]+ CONSTANTS[27]*STATES[0]+ CONSTANTS[28]*STATES[2]+ CONSTANTS[29]*STATES[6]+ CONSTANTS[30]*ALGEBRAIC[0])*STATES[13])/(CONSTANTS[32]+STATES[13]); resid[10] = RATES[1] - CONSTANTS[33]*STATES[1]*(1.00000 - STATES[1]/CONSTANTS[34]); resid[11] = RATES[5] - ALGEBRAIC[20]*(STATES[3] - STATES[5]) - (ALGEBRAIC[24]+ALGEBRAIC[16])*STATES[5]; resid[12] = RATES[4] - CONSTANTS[4]*(STATES[3] - STATES[4])*ALGEBRAIC[3] - (CONSTANTS[5]+ALGEBRAIC[16]+ALGEBRAIC[12])*STATES[4]; resid[13] = RATES[12] - ALGEBRAIC[27] - ALGEBRAIC[12]*STATES[12]; } 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[10]*STATES[1]; ALGEBRAIC[1] = STATES[7] - STATES[0]; ALGEBRAIC[2] = STATES[8] - STATES[6]; ALGEBRAIC[3] = STATES[12] - (STATES[4]+ALGEBRAIC[1]+ALGEBRAIC[2]); ALGEBRAIC[4] = CONSTANTS[35]+ CONSTANTS[36]*STATES[0]+ CONSTANTS[37]*STATES[6]+ CONSTANTS[38]*ALGEBRAIC[0]; ALGEBRAIC[5] = ALGEBRAIC[4]; ALGEBRAIC[6] = CONSTANTS[39]+ CONSTANTS[40]*STATES[2]+ CONSTANTS[41]*STATES[0]; ALGEBRAIC[7] = ( 2.00000*CONSTANTS[91]*ALGEBRAIC[5])/((ALGEBRAIC[6]+- ALGEBRAIC[5]+ CONSTANTS[88]*ALGEBRAIC[6]+ CONSTANTS[91]*ALGEBRAIC[5])+ pow((pow(ALGEBRAIC[6]+- ALGEBRAIC[5]+ CONSTANTS[88]*ALGEBRAIC[6]+ CONSTANTS[91]*ALGEBRAIC[5], 2.00000) - 4.00000*(ALGEBRAIC[6] - ALGEBRAIC[5])*CONSTANTS[91]*ALGEBRAIC[5]), 1.0 / 2)); ALGEBRAIC[8] = CONSTANTS[44]+ CONSTANTS[45]*STATES[6]+ CONSTANTS[46]*STATES[0]+ CONSTANTS[47]*STATES[2]; ALGEBRAIC[9] = CONSTANTS[48]+ CONSTANTS[49]*STATES[9]+ CONSTANTS[50]*STATES[11]; ALGEBRAIC[11] = STATES[9]; ALGEBRAIC[12] = (CONSTANTS[58]+ CONSTANTS[59]*STATES[0]+ CONSTANTS[60]*STATES[2]+ CONSTANTS[61]*STATES[6]+ CONSTANTS[62]*ALGEBRAIC[0])/(1.00000+ALGEBRAIC[11]/CONSTANTS[63]); ALGEBRAIC[13] = CONSTANTS[64]*STATES[2]; ALGEBRAIC[14] = ( 2.00000*CONSTANTS[94]*ALGEBRAIC[13])/((CONSTANTS[89]+- ALGEBRAIC[13]+ CONSTANTS[92]*CONSTANTS[89]+ CONSTANTS[94]*ALGEBRAIC[13])+ pow((pow(CONSTANTS[89]+- ALGEBRAIC[13]+ CONSTANTS[92]*CONSTANTS[89]+ CONSTANTS[94]*ALGEBRAIC[13], 2.00000) - 4.00000*(CONSTANTS[89] - ALGEBRAIC[13])*CONSTANTS[94]*ALGEBRAIC[13]), 1.0 / 2)); ALGEBRAIC[15] = CONSTANTS[68]+ CONSTANTS[69]*ALGEBRAIC[14]; ALGEBRAIC[16] = CONSTANTS[70]+ CONSTANTS[71]*STATES[13]+ CONSTANTS[72]*STATES[9]; ALGEBRAIC[17] = CONSTANTS[73]+ CONSTANTS[74]*ALGEBRAIC[11]; ALGEBRAIC[18] = CONSTANTS[75]*STATES[2]; ALGEBRAIC[19] = ( 2.00000*CONSTANTS[95]*ALGEBRAIC[17])/((ALGEBRAIC[18]+- ALGEBRAIC[17]+ CONSTANTS[93]*ALGEBRAIC[18]+ CONSTANTS[95]*ALGEBRAIC[17])+ pow((pow(ALGEBRAIC[18]+- ALGEBRAIC[17]+ CONSTANTS[93]*ALGEBRAIC[18]+ CONSTANTS[95]*ALGEBRAIC[17], 2.00000) - 4.00000*(ALGEBRAIC[18] - ALGEBRAIC[17])*CONSTANTS[95]*ALGEBRAIC[17]), 1.0 / 2)); ALGEBRAIC[20] = CONSTANTS[78]+ CONSTANTS[79]*ALGEBRAIC[19]; ALGEBRAIC[21] = CONSTANTS[80]*STATES[2]; ALGEBRAIC[22] = CONSTANTS[81]+ CONSTANTS[82]*ALGEBRAIC[11]; ALGEBRAIC[23] = ( 2.00000*CONSTANTS[97]*ALGEBRAIC[21])/((ALGEBRAIC[22]+- ALGEBRAIC[21]+ CONSTANTS[96]*ALGEBRAIC[22]+ CONSTANTS[97]*ALGEBRAIC[21])+ pow((pow(ALGEBRAIC[22]+- ALGEBRAIC[21]+ CONSTANTS[96]*ALGEBRAIC[22]+ CONSTANTS[97]*ALGEBRAIC[21], 2.00000) - 4.00000*(ALGEBRAIC[22] - ALGEBRAIC[21])*CONSTANTS[97]*ALGEBRAIC[21]), 1.0 / 2)); ALGEBRAIC[24] = CONSTANTS[85]+ CONSTANTS[86]*ALGEBRAIC[23]; ALGEBRAIC[25] = CONSTANTS[51]*ALGEBRAIC[11]; ALGEBRAIC[10] = CONSTANTS[52]+ CONSTANTS[53]*STATES[2]; ALGEBRAIC[26] = ( 2.00000*CONSTANTS[90]*ALGEBRAIC[25])/((ALGEBRAIC[10]+- ALGEBRAIC[25]+ CONSTANTS[87]*ALGEBRAIC[10]+ CONSTANTS[90]*ALGEBRAIC[25])+ pow((pow(ALGEBRAIC[10]+- ALGEBRAIC[25]+ CONSTANTS[87]*ALGEBRAIC[10]+ CONSTANTS[90]*ALGEBRAIC[25], 2.00000) - 4.00000*(ALGEBRAIC[10] - ALGEBRAIC[25])*CONSTANTS[90]*ALGEBRAIC[25]), 1.0 / 2)); ALGEBRAIC[27] = CONSTANTS[56]+ CONSTANTS[57]*ALGEBRAIC[26]; } 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; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { }