/* There are a total of 10 entries in the algebraic variable array. There are a total of 31 entries in each of the rate and state variable arrays. There are a total of 50 entries in the constant variable array. */ /* * VOI is time in component environment (second). * ALGEBRAIC[1] is LRG in component beta_adrenergic_receptor_module (micromolar). * ALGEBRAIC[2] is RG in component beta_adrenergic_receptor_module (micromolar). * ALGEBRAIC[0] is LR in component beta_adrenergic_receptor_module (micromolar). * STATES[0] is beta1_AR in component beta_adrenergic_receptor_module (micromolar). * STATES[1] is Gs in component beta_adrenergic_receptor_module (micromolar). * STATES[2] is beta1_ARact in component beta_adrenergic_receptor_module (micromolar). * STATES[3] is beta1_AR_S464 in component beta_adrenergic_receptor_module (micromolar). * STATES[4] is beta1_AR_S301 in component beta_adrenergic_receptor_module (micromolar). * STATES[5] is L in component beta_adrenergic_receptor_module (micromolar). * CONSTANTS[0] is Ltotmax in component beta_adrenergic_receptor_module (micromolar). * CONSTANTS[1] is KL in component beta_adrenergic_receptor_module (micromolar). * CONSTANTS[2] is KR in component beta_adrenergic_receptor_module (micromolar). * CONSTANTS[3] is KC in component beta_adrenergic_receptor_module (micromolar). * CONSTANTS[4] is Gstot in component beta_adrenergic_receptor_module (micromolar). * CONSTANTS[5] is k_betaARK_plus in component beta_adrenergic_receptor_module (first_order_rate_constant). * CONSTANTS[6] is k_betaARK_minus in component beta_adrenergic_receptor_module (first_order_rate_constant). * CONSTANTS[7] is k_PKA_plus in component beta_adrenergic_receptor_module (second_order_rate_constant). * CONSTANTS[8] is k_PKA_minus in component beta_adrenergic_receptor_module (first_order_rate_constant). * STATES[6] is Gs_beta_gamma in component Gs_activation_module (micromolar). * STATES[20] is PKACI in component PKA_activation_module (micromolar). * STATES[7] is Gs_alpha_GTPtot in component Gs_activation_module (micromolar). * STATES[8] is Gs_alpha_GDP in component Gs_activation_module (micromolar). * CONSTANTS[9] is k_gact in component Gs_activation_module (first_order_rate_constant). * CONSTANTS[10] is k_hyd in component Gs_activation_module (first_order_rate_constant). * CONSTANTS[11] is k_reassoc in component Gs_activation_module (second_order_rate_constant). * STATES[9] is cAMPtot in component cyclic_AMP_metabolism_module (micromolar). * STATES[10] is AC in component cyclic_AMP_metabolism_module (micromolar). * STATES[11] is Gs_alpha_GTP in component cyclic_AMP_metabolism_module (micromolar). * ALGEBRAIC[3] is Gs_alpha_GTP_AC in component cyclic_AMP_metabolism_module (micromolar). * STATES[12] is PDE in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[12] is ACtot in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[13] is ATP in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[14] is PDEtot in component cyclic_AMP_metabolism_module (micromolar). * ALGEBRAIC[5] is PDEinhib in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[15] is IBMXtot in component cyclic_AMP_metabolism_module (micromolar). * STATES[13] is IBMX in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[16] is fsktot in component cyclic_AMP_metabolism_module (micromolar). * STATES[14] is fsk in component cyclic_AMP_metabolism_module (micromolar). * ALGEBRAIC[4] is fsk_AC in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[17] is kAC_basal in component cyclic_AMP_metabolism_module (first_order_rate_constant). * CONSTANTS[18] is kAC_fsk in component cyclic_AMP_metabolism_module (first_order_rate_constant). * CONSTANTS[19] is k_PDE in component cyclic_AMP_metabolism_module (first_order_rate_constant). * CONSTANTS[20] is kAC_Gs_alpha in component cyclic_AMP_metabolism_module (first_order_rate_constant). * CONSTANTS[21] is Km_basal in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[22] is Km_PDE in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[23] is Km_fsk in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[24] is K_fsk in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[25] is Km_Gs_alpha_GTP in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[26] is K_Gs_alpha in component cyclic_AMP_metabolism_module (micromolar). * CONSTANTS[27] is K_IBMX in component cyclic_AMP_metabolism_module (micromolar). * STATES[21] is cAMP in component PKA_activation_module (micromolar). * STATES[22] is PKACII in component PKA_activation_module (micromolar). * STATES[23] is PKAtemp in component PKA_activation_module (micromolar). * STATES[24] is ARCI in component PKA_activation_module (micromolar). * STATES[25] is A2RCI in component PKA_activation_module (micromolar). * STATES[26] is A2RI in component PKA_activation_module (micromolar). * STATES[27] is A2RCII in component PKA_activation_module (micromolar). * STATES[28] is A2RII in component PKA_activation_module (micromolar). * STATES[29] is ARCII in component PKA_activation_module (micromolar). * CONSTANTS[28] is Ki_pki in component PKA_activation_module (micromolar). * CONSTANTS[29] is PKAItot in component PKA_activation_module (micromolar). * CONSTANTS[30] is PKAIItot in component PKA_activation_module (micromolar). * CONSTANTS[31] is PKItot in component PKA_activation_module (micromolar). * CONSTANTS[32] is KA in component PKA_activation_module (micromolar). * CONSTANTS[33] is KB in component PKA_activation_module (micromolar). * CONSTANTS[34] is KD in component PKA_activation_module (micromolar). * CONSTANTS[35] is KPKI in component PKA_activation_module (micromolar). * STATES[30] is PKI in component PKA_activation_module (micromolar). * ALGEBRAIC[8] is fracPLBp in component phospholamban_regulation_module (dimensionless). * ALGEBRAIC[9] is fracPLB in component phospholamban_regulation_module (dimensionless). * CONSTANTS[36] is fracPLBo in component phospholamban_regulation_module (dimensionless). * STATES[15] is PLBp in component phospholamban_regulation_module (micromolar). * ALGEBRAIC[7] is PLB in component phospholamban_regulation_module (micromolar). * STATES[16] is Inhib1ptot in component phospholamban_regulation_module (micromolar). * STATES[17] is Inhib1 in component phospholamban_regulation_module (micromolar). * STATES[18] is Inhib1p in component phospholamban_regulation_module (micromolar). * ALGEBRAIC[6] is PP1_Inhib1p in component phospholamban_regulation_module (micromolar). * STATES[19] is PP1 in component phospholamban_regulation_module (micromolar). * CONSTANTS[37] is PP1tot in component phospholamban_regulation_module (micromolar). * CONSTANTS[38] is PLBtot in component phospholamban_regulation_module (micromolar). * CONSTANTS[39] is Inhib1tot in component phospholamban_regulation_module (micromolar). * CONSTANTS[40] is epsilon in component phospholamban_regulation_module (dimensionless). * CONSTANTS[41] is kPKA_PLB in component phospholamban_regulation_module (first_order_rate_constant). * CONSTANTS[42] is KmPKA_PLB in component phospholamban_regulation_module (micromolar). * CONSTANTS[43] is kPKA_Inhib1 in component phospholamban_regulation_module (first_order_rate_constant). * CONSTANTS[44] is kPP1_PLB in component phospholamban_regulation_module (first_order_rate_constant). * CONSTANTS[45] is KmPP1_PLB in component phospholamban_regulation_module (micromolar). * CONSTANTS[46] is KmPKA_Inhib1 in component phospholamban_regulation_module (micromolar). * CONSTANTS[47] is VmaxPP2A_Inhib1 in component phospholamban_regulation_module (flux). * CONSTANTS[48] is KmPP2A_Inhib1 in component phospholamban_regulation_module (micromolar). * CONSTANTS[49] is KInhib1 in component phospholamban_regulation_module (micromolar). * RATES[5] is d/dt L in component beta_adrenergic_receptor_module (micromolar). * RATES[1] is d/dt Gs in component beta_adrenergic_receptor_module (micromolar). * RATES[0] is d/dt beta1_AR in component beta_adrenergic_receptor_module (micromolar). * RATES[2] is d/dt beta1_ARact in component beta_adrenergic_receptor_module (micromolar). * RATES[3] is d/dt beta1_AR_S464 in component beta_adrenergic_receptor_module (micromolar). * RATES[4] is d/dt beta1_AR_S301 in component beta_adrenergic_receptor_module (micromolar). * RATES[7] is d/dt Gs_alpha_GTPtot in component Gs_activation_module (micromolar). * RATES[6] is d/dt Gs_beta_gamma in component Gs_activation_module (micromolar). * RATES[8] is d/dt Gs_alpha_GDP in component Gs_activation_module (micromolar). * RATES[11] is d/dt Gs_alpha_GTP in component cyclic_AMP_metabolism_module (micromolar). * RATES[14] is d/dt fsk in component cyclic_AMP_metabolism_module (micromolar). * RATES[10] is d/dt AC in component cyclic_AMP_metabolism_module (micromolar). * RATES[12] is d/dt PDE in component cyclic_AMP_metabolism_module (micromolar). * RATES[13] is d/dt IBMX in component cyclic_AMP_metabolism_module (micromolar). * RATES[9] is d/dt cAMPtot in component cyclic_AMP_metabolism_module (micromolar). * RATES[19] is d/dt PP1 in component phospholamban_regulation_module (micromolar). * RATES[18] is d/dt Inhib1p in component phospholamban_regulation_module (micromolar). * RATES[15] is d/dt PLBp in component phospholamban_regulation_module (micromolar). * RATES[16] is d/dt Inhib1ptot in component phospholamban_regulation_module (micromolar). * RATES[17] is d/dt Inhib1 in component phospholamban_regulation_module (micromolar). * There are a total of 0 condition variables. */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { STATES[0] = 0.0001; STATES[1] = 3.182; STATES[2] = 0.01205; STATES[3] = 0; STATES[4] = 1.154e-3; STATES[5] = 0.988; CONSTANTS[0] = 1; CONSTANTS[1] = 0.285; CONSTANTS[2] = 0.062; CONSTANTS[3] = 33; CONSTANTS[4] = 3.83; CONSTANTS[5] = 1.1e-3; CONSTANTS[6] = 2.2e-3; CONSTANTS[7] = 3.6e-3; CONSTANTS[8] = 2.2e-3; STATES[6] = 0.02569; STATES[7] = 0.02505; STATES[8] = 6.44e-4; CONSTANTS[9] = 16; CONSTANTS[10] = 0.8; CONSTANTS[11] = 1.21e3; STATES[9] = 0.8453; STATES[10] = 0.04706295; STATES[11] = 0.02241295; STATES[12] = 0.0389; CONSTANTS[12] = 49.7e-3; CONSTANTS[13] = 5e3; CONSTANTS[14] = 38.9e-3; CONSTANTS[15] = 0; STATES[13] = 0; CONSTANTS[16] = 0; STATES[14] = 0; CONSTANTS[17] = 0.2; CONSTANTS[18] = 7.3; CONSTANTS[19] = 5; CONSTANTS[20] = 8.5; CONSTANTS[21] = 1.03e3; CONSTANTS[22] = 1.3; CONSTANTS[23] = 860; CONSTANTS[24] = 44; CONSTANTS[25] = 0.4; CONSTANTS[26] = 0.4; CONSTANTS[27] = 30; CONSTANTS[28] = 0.2e-3; CONSTANTS[29] = 0.59; CONSTANTS[30] = 0.025; CONSTANTS[31] = 0.18; CONSTANTS[32] = 9.14e-3; CONSTANTS[33] = 1.64e-3; CONSTANTS[34] = 4.375e-3; CONSTANTS[35] = 2e-4; CONSTANTS[36] = 0.9613; STATES[15] = 4.105; STATES[16] = 0.0526; STATES[17] = 0.2474; STATES[18] = 6.27339e-5; STATES[19] = 0.8374627; CONSTANTS[37] = 0.89; CONSTANTS[38] = 106; CONSTANTS[39] = 0.3; CONSTANTS[40] = 10; CONSTANTS[41] = 54; CONSTANTS[42] = 21; CONSTANTS[43] = 60; CONSTANTS[44] = 8.5; CONSTANTS[45] = 7; CONSTANTS[46] = 1; CONSTANTS[47] = 14; CONSTANTS[48] = 1; CONSTANTS[49] = 1e-3; STATES[20] = 0.1001; STATES[21] = 0.1001; STATES[22] = 0.1001; STATES[23] = 0.1001; STATES[24] = 0.1001; STATES[25] = 0.1001; STATES[26] = 0.1001; STATES[27] = 0.1001; STATES[28] = 0.1001; STATES[29] = 0.1001; STATES[30] = 0.1001; RATES[5] = 0.1001; RATES[1] = 0.1001; RATES[0] = 0.1001; RATES[2] = 0.1001; RATES[3] = 0.1001; RATES[4] = 0.1001; RATES[7] = 0.1001; RATES[6] = 0.1001; RATES[8] = 0.1001; RATES[11] = 0.1001; RATES[14] = 0.1001; RATES[10] = 0.1001; RATES[12] = 0.1001; RATES[13] = 0.1001; RATES[9] = 0.1001; RATES[19] = 0.1001; RATES[18] = 0.1001; RATES[15] = 0.1001; RATES[16] = 0.1001; RATES[17] = 0.1001; } void computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { resid[0] = RATES[5] - CONSTANTS[0] - (ALGEBRAIC[0]+ALGEBRAIC[1]+STATES[5]); resid[1] = RATES[1] - CONSTANTS[4] - (ALGEBRAIC[2]+ALGEBRAIC[1]+STATES[1]); resid[2] = RATES[0] - STATES[2] - (ALGEBRAIC[0]+ALGEBRAIC[1]+ALGEBRAIC[2]+STATES[0]); resid[3] = RATES[2] - ( CONSTANTS[6]*STATES[3] - CONSTANTS[5]*(ALGEBRAIC[0]+ALGEBRAIC[1]))+( CONSTANTS[8]*STATES[4] - CONSTANTS[7]*STATES[20]*STATES[2]); resid[4] = RATES[3] - CONSTANTS[5]*(ALGEBRAIC[0]+ALGEBRAIC[1]) - CONSTANTS[6]*STATES[3]; resid[5] = RATES[4] - CONSTANTS[7]*STATES[20]*STATES[2] - CONSTANTS[8]*STATES[4]; resid[6] = RATES[7] - CONSTANTS[9]*(ALGEBRAIC[2]+ALGEBRAIC[1]) - CONSTANTS[10]*STATES[7]; resid[7] = RATES[6] - CONSTANTS[9]*(ALGEBRAIC[2]+ALGEBRAIC[1]) - CONSTANTS[11]*STATES[8]*STATES[6]; resid[8] = RATES[8] - CONSTANTS[10]*STATES[7] - CONSTANTS[11]*STATES[8]*STATES[6]; resid[9] = RATES[11] - STATES[7] - (ALGEBRAIC[3]+STATES[11]); resid[10] = RATES[14] - CONSTANTS[16] - (ALGEBRAIC[4]+STATES[14]); resid[11] = RATES[10] - CONSTANTS[12] - (ALGEBRAIC[3]+STATES[10]); resid[12] = RATES[12] - CONSTANTS[14] - (ALGEBRAIC[5]+STATES[12]); resid[13] = RATES[13] - CONSTANTS[15] - (ALGEBRAIC[5]+STATES[13]); resid[14] = RATES[9] - (( CONSTANTS[17]*STATES[10]*CONSTANTS[13])/(CONSTANTS[21]+CONSTANTS[13])+( CONSTANTS[20]*ALGEBRAIC[3]*CONSTANTS[13])/(CONSTANTS[25]+CONSTANTS[13])+( CONSTANTS[18]*ALGEBRAIC[4]*CONSTANTS[13])/(CONSTANTS[23]+CONSTANTS[13])) - ( CONSTANTS[19]*STATES[12]*STATES[9])/(CONSTANTS[22]+STATES[9]); resid[15] = STATES[30] - ( CONSTANTS[31]*CONSTANTS[28])/(CONSTANTS[28]+STATES[20]+STATES[22]); resid[16] = STATES[25] - (STATES[20]/CONSTANTS[34])*STATES[20]*(1.00000+STATES[30]/CONSTANTS[28]); resid[17] = STATES[26] - STATES[20]*(1.00000+STATES[30]/CONSTANTS[28]); resid[18] = STATES[27] - (STATES[22]/CONSTANTS[34])*STATES[22]*(1.00000+STATES[30]/CONSTANTS[28]); resid[19] = STATES[28] - STATES[22]*(1.00000+STATES[30]/CONSTANTS[28]); resid[20] = STATES[24] - (CONSTANTS[32]/STATES[21])*STATES[25]; resid[21] = STATES[29] - (CONSTANTS[32]/STATES[21])*STATES[27]; resid[22] = STATES[21] - (STATES[9] - (STATES[24]+ 2.00000*STATES[25]+ 2.00000*STATES[26])) - (STATES[29]+ 2.00000*STATES[27]+ 2.00000*STATES[28]); resid[23] = STATES[23] - ( CONSTANTS[32]*CONSTANTS[33])/CONSTANTS[34]+( CONSTANTS[32]*STATES[21])/CONSTANTS[34]+( STATES[21]*STATES[21])/CONSTANTS[34]; resid[24] = STATES[20] - 2.00000*CONSTANTS[29]*STATES[21]*STATES[21] - STATES[20]*(1.00000+STATES[30]/CONSTANTS[28])*( STATES[23]*STATES[20]+ STATES[21]*STATES[21]); resid[25] = STATES[22] - 2.00000*CONSTANTS[30]*STATES[21]*STATES[21] - STATES[22]*(1.00000+STATES[30]/CONSTANTS[28])*( STATES[23]*STATES[22]+ STATES[21]*STATES[21]); resid[26] = RATES[19] - CONSTANTS[37] - (ALGEBRAIC[6]+STATES[19]); resid[27] = RATES[18] - STATES[16] - (ALGEBRAIC[6]+STATES[18]); resid[28] = RATES[15] - ( CONSTANTS[41]*STATES[20]*ALGEBRAIC[7])/(CONSTANTS[42]+ALGEBRAIC[7]) - ( CONSTANTS[44]*STATES[19]*STATES[15])/(CONSTANTS[45]+STATES[15]); resid[29] = RATES[16] - ( CONSTANTS[43]*STATES[20]*STATES[16])/(CONSTANTS[46]+STATES[16]) - ( CONSTANTS[47]*STATES[17])/(CONSTANTS[48]+STATES[17]); resid[30] = RATES[17] - ( CONSTANTS[47]*STATES[17])/(CONSTANTS[48]+STATES[17]) - ( CONSTANTS[43]*STATES[20]*STATES[16])/(CONSTANTS[46]+STATES[16]); } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[8] = STATES[15]/CONSTANTS[38]; ALGEBRAIC[9] = ALGEBRAIC[7]/CONSTANTS[38]; } void computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[0] = ( STATES[5]*STATES[0])/CONSTANTS[1]; ALGEBRAIC[1] = ( STATES[5]*STATES[0]*STATES[1])/( CONSTANTS[1]*CONSTANTS[2]); ALGEBRAIC[2] = ( STATES[0]*STATES[1])/CONSTANTS[3]; ALGEBRAIC[3] = ( STATES[11]*STATES[10])/CONSTANTS[26]; ALGEBRAIC[4] = ( STATES[14]*STATES[10])/CONSTANTS[24]; ALGEBRAIC[5] = ( STATES[12]*STATES[13])/CONSTANTS[27]; ALGEBRAIC[6] = ( STATES[19]*STATES[18])/CONSTANTS[49]; ALGEBRAIC[7] = CONSTANTS[38] - STATES[15]; } 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[20] = 0.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[21] = 0.0; SI[22] = 0.0; SI[23] = 0.0; SI[24] = 0.0; SI[25] = 0.0; SI[26] = 0.0; SI[27] = 0.0; SI[28] = 0.0; SI[29] = 0.0; SI[30] = 0.0; SI[15] = 1.0; SI[16] = 1.0; SI[17] = 1.0; SI[18] = 1.0; SI[19] = 1.0; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { }