/* There are a total of 4 entries in the algebraic variable array. There are a total of 17 entries in each of the rate and state variable arrays. There are a total of 43 entries in the constant variable array. */ /* * VOI is time in component environment (minute). * STATES[0] is MPF_n in component MPF_n (dimensionless). * CONSTANTS[0] is k_in in component parameters (first_order_rate_constant). * CONSTANTS[1] is k_out in component parameters (first_order_rate_constant). * CONSTANTS[2] is k_dn_ in component parameters (first_order_rate_constant). * CONSTANTS[3] is k_dn__ in component parameters (first_order_rate_constant). * CONSTANTS[4] is k_wee_ in component parameters (first_order_rate_constant). * CONSTANTS[5] is k_wee__ in component parameters (first_order_rate_constant). * CONSTANTS[6] is k_stg_ in component parameters (first_order_rate_constant). * CONSTANTS[7] is k_stg__ in component parameters (first_order_rate_constant). * STATES[1] is MPF_c in component MPF_c (dimensionless). * STATES[2] is FZY in component FZY (dimensionless). * STATES[3] is Wee1_n in component Wee1_n (dimensionless). * STATES[4] is StgP_n in component StgP_n (dimensionless). * STATES[5] is preMPF_n in component preMPF_n (dimensionless). * STATES[6] is preMPF_c in component preMPF_c (dimensionless). * CONSTANTS[8] is k_sc in component parameters (first_order_rate_constant). * CONSTANTS[9] is k_dc_ in component parameters (first_order_rate_constant). * CONSTANTS[10] is epsilon in component parameters (dimensionless). * STATES[7] is N in component N (dimensionless). * STATES[8] is StgP_c in component StgP_c (dimensionless). * STATES[9] is Wee1_c in component Wee1_c (dimensionless). * STATES[10] is IE in component IE (dimensionless). * CONSTANTS[11] is j_aie in component parameters (dimensionless). * CONSTANTS[12] is j_iie in component parameters (dimensionless). * CONSTANTS[13] is k_aie in component parameters (first_order_rate_constant). * CONSTANTS[14] is k_iie in component parameters (first_order_rate_constant). * CONSTANTS[15] is j_afz in component parameters (dimensionless). * CONSTANTS[16] is j_ifz in component parameters (dimensionless). * CONSTANTS[17] is k_afz in component parameters (first_order_rate_constant). * CONSTANTS[18] is k_ifz in component parameters (first_order_rate_constant). * STATES[11] is Stg_m in component Stg_m (dimensionless). * CONSTANTS[19] is k_dm_ in component parameters (first_order_rate_constant). * CONSTANTS[20] is k_dm__ in component parameters (first_order_rate_constant). * CONSTANTS[21] is j_m in component parameters (dimensionless). * STATES[12] is Xp in component Xp (dimensionless). * STATES[13] is Xm in component Xm (dimensionless). * CONSTANTS[22] is k_sxm in component parameters (first_order_rate_constant). * CONSTANTS[23] is k_sxp in component parameters (first_order_rate_constant). * CONSTANTS[24] is k_ins in component parameters (first_order_rate_constant). * CONSTANTS[25] is k_outs in component parameters (first_order_rate_constant). * CONSTANTS[26] is k_astg_ in component parameters (first_order_rate_constant). * CONSTANTS[27] is k_astg__ in component parameters (first_order_rate_constant). * CONSTANTS[28] is k_istg in component parameters (first_order_rate_constant). * CONSTANTS[29] is k_dstg in component parameters (first_order_rate_constant). * CONSTANTS[30] is j_astg in component parameters (dimensionless). * CONSTANTS[31] is j_istg in component parameters (dimensionless). * STATES[14] is Stg_n in component Stg_n (dimensionless). * STATES[15] is Stg_c in component Stg_c (dimensionless). * CONSTANTS[32] is k_sstg in component parameters (first_order_rate_constant). * CONSTANTS[33] is k_inw in component parameters (first_order_rate_constant). * CONSTANTS[34] is k_outw in component parameters (first_order_rate_constant). * CONSTANTS[35] is k_awee in component parameters (first_order_rate_constant). * CONSTANTS[36] is k_iwee_ in component parameters (first_order_rate_constant). * CONSTANTS[37] is k_iwee__ in component parameters (first_order_rate_constant). * CONSTANTS[38] is j_awee in component parameters (dimensionless). * CONSTANTS[39] is j_iwee in component parameters (dimensionless). * STATES[16] is Wee1P_n in component Wee1P_n (dimensionless). * ALGEBRAIC[0] is Wee1P_c in component Wee1P_c (dimensionless). * CONSTANTS[40] is Wee1_T in component Wee1P_c (dimensionless). * ALGEBRAIC[1] is CycB_T in component CycB_T (dimensionless). * ALGEBRAIC[2] is Stg_T in component Stg_T (dimensionless). * ALGEBRAIC[3] is StgP_T in component StgP_T (dimensionless). * CONSTANTS[41] is k_ez in component parameters (first_order_rate_constant). * RATES[0] is d/dt MPF_n in component MPF_n (dimensionless). * RATES[5] is d/dt preMPF_n in component preMPF_n (dimensionless). * RATES[1] is d/dt MPF_c in component MPF_c (dimensionless). * RATES[6] is d/dt preMPF_c in component preMPF_c (dimensionless). * RATES[10] is d/dt IE in component IE (dimensionless). * RATES[2] is d/dt FZY in component FZY (dimensionless). * RATES[11] is d/dt Stg_m in component Stg_m (dimensionless). * RATES[13] is d/dt Xm in component Xm (dimensionless). * RATES[12] is d/dt Xp in component Xp (dimensionless). * RATES[4] is d/dt StgP_n in component StgP_n (dimensionless). * RATES[14] is d/dt Stg_n in component Stg_n (dimensionless). * RATES[8] is d/dt StgP_c in component StgP_c (dimensionless). * RATES[15] is d/dt Stg_c in component Stg_c (dimensionless). * RATES[3] is d/dt Wee1_n in component Wee1_n (dimensionless). * RATES[9] is d/dt Wee1_c in component Wee1_c (dimensionless). * RATES[16] is d/dt Wee1P_n in component Wee1P_n (dimensionless). * RATES[7] is d/dt N in component N (dimensionless). * There are a total of 0 condition variables. */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { STATES[0] = 0; CONSTANTS[0] = 0.15; CONSTANTS[1] = 0; CONSTANTS[2] = 0.01; CONSTANTS[3] = 1.5; CONSTANTS[4] = 0.005; CONSTANTS[5] = 1; CONSTANTS[6] = 0.2; CONSTANTS[7] = 2; STATES[1] = 1; STATES[2] = 0; STATES[3] = 0; STATES[4] = 0; STATES[5] = 0; STATES[6] = 0; CONSTANTS[8] = 0.01; CONSTANTS[9] = 0.01; CONSTANTS[10] = 0.00007; STATES[7] = 1; STATES[8] = 0; STATES[9] = 0; STATES[10] = 0; CONSTANTS[11] = 0.01; CONSTANTS[12] = 0.01; CONSTANTS[13] = 1; CONSTANTS[14] = 0.4; CONSTANTS[15] = 0.01; CONSTANTS[16] = 0.01; CONSTANTS[17] = 1; CONSTANTS[18] = 0.2; STATES[11] = 1; CONSTANTS[19] = 0.002; CONSTANTS[20] = 0.2; CONSTANTS[21] = 0.05; STATES[12] = 0; STATES[13] = 0; CONSTANTS[22] = 0.0005; CONSTANTS[23] = 0.001; CONSTANTS[24] = 0.08; CONSTANTS[25] = 0.02; CONSTANTS[26] = 0; CONSTANTS[27] = 1; CONSTANTS[28] = 0.3; CONSTANTS[29] = 0.015; CONSTANTS[30] = 0.05; CONSTANTS[31] = 0.05; STATES[14] = 0; STATES[15] = 1; CONSTANTS[32] = 0.02; CONSTANTS[33] = 0.04; CONSTANTS[34] = 0.01; CONSTANTS[35] = 0.3; CONSTANTS[36] = 0.01; CONSTANTS[37] = 1; CONSTANTS[38] = 0.05; CONSTANTS[39] = 0.05; STATES[16] = 0; CONSTANTS[40] = 0.8; CONSTANTS[41] = 0.5; CONSTANTS[42] = 0.00000; RATES[0] = 0.1001; RATES[5] = 0.1001; RATES[1] = 0.1001; RATES[6] = 0.1001; RATES[10] = 0.1001; RATES[2] = 0.1001; RATES[11] = 0.1001; RATES[13] = 0.1001; RATES[12] = 0.1001; RATES[4] = 0.1001; RATES[14] = 0.1001; RATES[8] = 0.1001; RATES[15] = 0.1001; RATES[3] = 0.1001; RATES[9] = 0.1001; RATES[16] = 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[0]*STATES[1] - CONSTANTS[1]*STATES[0]) - (CONSTANTS[2]+ CONSTANTS[3]*STATES[2])*STATES[0]) - (CONSTANTS[4]+ CONSTANTS[5]*STATES[3])*STATES[0])+ (CONSTANTS[6]+ CONSTANTS[7]*STATES[4])*STATES[5]; resid[1] = RATES[5] - ((( CONSTANTS[0]*STATES[6] - CONSTANTS[1]*STATES[5]) - (CONSTANTS[2]+ CONSTANTS[3]*STATES[2])*STATES[5])+ (CONSTANTS[4]+ CONSTANTS[5]*STATES[3])*STATES[0]) - (CONSTANTS[4]+ CONSTANTS[7]*STATES[4])*STATES[5]; resid[2] = RATES[1] - (((CONSTANTS[8] - ( CONSTANTS[10]*STATES[7]*( CONSTANTS[0]*STATES[1] - CONSTANTS[1]*STATES[0]))/(1.00000 - STATES[7]*CONSTANTS[10])) - CONSTANTS[9]*STATES[1])+ (CONSTANTS[6]+ CONSTANTS[7]*STATES[8])*STATES[6]) - (CONSTANTS[4]+ CONSTANTS[5]*STATES[9])*STATES[1]; resid[3] = RATES[6] - ((( - CONSTANTS[10]*STATES[7]*( CONSTANTS[0]*STATES[6] - CONSTANTS[1]*STATES[5]))/(1.00000 - STATES[7]*CONSTANTS[10]) - CONSTANTS[9]*STATES[6]) - (CONSTANTS[6]+ CONSTANTS[7]*STATES[8])*STATES[6])+ (CONSTANTS[4]+ CONSTANTS[5]*STATES[9])*STATES[1]; resid[4] = RATES[10] - ( CONSTANTS[13]*(1.00000 - STATES[10])*STATES[0])/((CONSTANTS[11]+1.00000) - STATES[10]) - ( CONSTANTS[14]*STATES[10])/(CONSTANTS[12]+STATES[10]); resid[5] = RATES[2] - ( CONSTANTS[17]*STATES[10]*(1.00000 - STATES[2]))/((CONSTANTS[15]+1.00000) - STATES[2]) - ( CONSTANTS[18]*STATES[2])/(CONSTANTS[16]+STATES[2]); resid[6] = RATES[11] - - (CONSTANTS[19]/(CONSTANTS[21]+STATES[11])+ CONSTANTS[20]*STATES[12])*STATES[11]; resid[7] = RATES[13] - CONSTANTS[22]*STATES[7]; resid[8] = RATES[12] - CONSTANTS[23]*STATES[13]; resid[9] = RATES[4] - ((( CONSTANTS[24]*STATES[8] - CONSTANTS[25]*STATES[4])+( (CONSTANTS[26]+ CONSTANTS[27]*STATES[0])*STATES[14])/(CONSTANTS[30]+STATES[14])) - ( CONSTANTS[28]*STATES[4])/(CONSTANTS[31]+STATES[4])) - CONSTANTS[29]*STATES[4]; resid[10] = RATES[14] - ((( CONSTANTS[24]*STATES[15] - CONSTANTS[25]*STATES[14]) - ( (CONSTANTS[26]+ CONSTANTS[27]*STATES[0])*STATES[14])/(CONSTANTS[30]+STATES[14]))+( CONSTANTS[28]*STATES[4])/(CONSTANTS[31]+STATES[4])) - CONSTANTS[29]*STATES[14]; resid[11] = RATES[8] - (( - CONSTANTS[29]*STATES[8] - ( CONSTANTS[10]*STATES[7]*( CONSTANTS[24]*STATES[8] - CONSTANTS[25]*STATES[4]))/(1.00000 - STATES[7]*CONSTANTS[10]))+( (CONSTANTS[26]+ CONSTANTS[27]*STATES[1])*STATES[15])/(CONSTANTS[30]+STATES[15])) - ( CONSTANTS[28]*STATES[8])/(CONSTANTS[31]+STATES[8]); resid[12] = RATES[15] - ((( CONSTANTS[32]*STATES[11] - CONSTANTS[29]*STATES[15]) - ( CONSTANTS[10]*STATES[7]*( CONSTANTS[24]*STATES[15] - CONSTANTS[25]*STATES[14]))/(1.00000 - STATES[7]*CONSTANTS[10])) - ( (CONSTANTS[26]+ CONSTANTS[27]*STATES[1])*STATES[15])/(CONSTANTS[30]+STATES[15]))+( CONSTANTS[28]*STATES[8])/(CONSTANTS[31]+STATES[8]); resid[13] = RATES[3] - (( CONSTANTS[33]*STATES[9] - CONSTANTS[34]*STATES[3])+( CONSTANTS[35]*STATES[16])/(CONSTANTS[38]+STATES[16])) - ( (CONSTANTS[36]+ CONSTANTS[37]*STATES[0])*STATES[3])/(CONSTANTS[39]+STATES[3]); resid[14] = RATES[9] - (( - ( CONSTANTS[33]*STATES[9] - CONSTANTS[34]*STATES[3])*STATES[7]*CONSTANTS[10])/(1.00000 - STATES[7]*CONSTANTS[10])+( CONSTANTS[35]*ALGEBRAIC[0])/(CONSTANTS[38]+ALGEBRAIC[0])) - ( (CONSTANTS[36]+ CONSTANTS[37]*STATES[1])*STATES[9])/(CONSTANTS[39]+STATES[9]); resid[15] = RATES[16] - (( CONSTANTS[33]*ALGEBRAIC[0] - CONSTANTS[34]*STATES[16]) - ( CONSTANTS[35]*STATES[16])/(CONSTANTS[38]+STATES[16]))+( (CONSTANTS[36]+ CONSTANTS[37]*STATES[0])*STATES[3])/(CONSTANTS[39]+STATES[3]); resid[16] = RATES[7] - CONSTANTS[42]; } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[1] = (1.00000 - STATES[7]*CONSTANTS[10])*(STATES[1]+STATES[6])+ STATES[7]*CONSTANTS[10]*(STATES[0]+STATES[5]); ALGEBRAIC[2] = (1.00000 - STATES[7]*CONSTANTS[10])*(STATES[15]+STATES[8])+ STATES[7]*CONSTANTS[10]*(STATES[14]+STATES[4]); ALGEBRAIC[3] = (1.00000 - STATES[7]*CONSTANTS[10])*STATES[8]+ STATES[7]*CONSTANTS[10]*STATES[4]; } void computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[0] = (CONSTANTS[40] - STATES[7]*CONSTANTS[10]*(STATES[3]+STATES[16]))/(1.00000 - STATES[7]*CONSTANTS[10]) - STATES[9]; } 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; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { }