/* There are a total of 0 entries in the algebraic variable array. There are a total of 9 entries in each of the rate and state variable arrays. There are a total of 20 entries in the constant variable array. */ /* * VOI is t in component environment (s). * CONSTANTS[0] is kf1 in component PKC (per_s). * CONSTANTS[1] is kb1 in component PKC (per_s). * CONSTANTS[2] is kf2 in component PKC (per_um_s). * CONSTANTS[3] is kb2 in component PKC (per_s). * CONSTANTS[4] is kf3 in component PKC (per_s). * CONSTANTS[5] is kb3 in component PKC (per_s). * CONSTANTS[6] is kf4 in component PKC (per_um_s). * CONSTANTS[7] is kb4 in component PKC (per_s). * CONSTANTS[8] is kf5 in component PKC (per_s). * CONSTANTS[9] is kb5 in component PKC (per_s). * CONSTANTS[10] is kf6 in component PKC (per_s). * CONSTANTS[11] is kb6 in component PKC (per_s). * CONSTANTS[12] is kf7 in component PKC (per_um_s). * CONSTANTS[13] is kb7 in component PKC (per_s). * CONSTANTS[14] is kf8 in component PKC (per_um_s). * CONSTANTS[15] is kb8 in component PKC (per_s). * CONSTANTS[16] is kf9 in component PKC (per_um_s). * CONSTANTS[17] is kb9 in component PKC (per_s). * CONSTANTS[18] is kf10 in component PKC (per_um_s). * CONSTANTS[19] is kb10 in component PKC (per_s). * STATES[0] is PKC_i in component PKC (um). * STATES[1] is AA in component PKC (um). * STATES[2] is DAG in component PKC (um). * STATES[3] is Ca in component PKC (um). * STATES[4] is DAGPKC in component PKC (um). * STATES[5] is AADAGPKC in component PKC (um). * STATES[6] is DAGCaPKC in component PKC (um). * STATES[7] is CaPKC in component PKC (um). * STATES[8] is PKC_a in component PKC (um). * RATES[0] is d/dt PKC_i in component PKC (um). * RATES[4] is d/dt DAGPKC in component PKC (um). * RATES[5] is d/dt AADAGPKC in component PKC (um). * RATES[7] is d/dt CaPKC in component PKC (um). * RATES[6] is d/dt DAGCaPKC in component PKC (um). * RATES[8] is d/dt PKC_a in component PKC (um). * RATES[1] is d/dt AA in component PKC (um). * RATES[2] is d/dt DAG in component PKC (um). * RATES[3] is d/dt Ca in component PKC (um). */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { CONSTANTS[0] = 1; CONSTANTS[1] = 50; CONSTANTS[2] = 0.0000000002; CONSTANTS[3] = 0.1; CONSTANTS[4] = 1.2705; CONSTANTS[5] = 3.5026; CONSTANTS[6] = 0.000000002; CONSTANTS[7] = 0.1; CONSTANTS[8] = 1; CONSTANTS[9] = 0.1; CONSTANTS[10] = 2; CONSTANTS[11] = 0.2; CONSTANTS[12] = 0.000001; CONSTANTS[13] = 0.5; CONSTANTS[14] = 0.000000013333; CONSTANTS[15] = 8.6348; CONSTANTS[16] = 0.000000001; CONSTANTS[17] = 0.1; CONSTANTS[18] = 0.00000003; CONSTANTS[19] = 2; STATES[0] = 1; STATES[1] = 1; STATES[2] = 0.1; STATES[3] = 0.1; STATES[4] = 0; STATES[5] = 0; STATES[6] = 0; STATES[7] = 0; STATES[8] = 0; } void computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { RATES[0] = (((((( - STATES[0]*CONSTANTS[0]+ STATES[8]*CONSTANTS[1]) - STATES[0]*STATES[1]*CONSTANTS[2])+ STATES[8]*CONSTANTS[3]) - STATES[0]*STATES[2]*CONSTANTS[16])+ STATES[4]*CONSTANTS[17]) - STATES[0]*STATES[3]*CONSTANTS[12])+ STATES[7]*CONSTANTS[13]; RATES[4] = (( STATES[0]*STATES[2]*CONSTANTS[16] - STATES[4]*CONSTANTS[17]) - STATES[4]*STATES[1]*CONSTANTS[18])+ STATES[5]*CONSTANTS[19]; RATES[5] = (( STATES[4]*STATES[1]*CONSTANTS[18] - STATES[5]*CONSTANTS[19]) - STATES[5]*CONSTANTS[10])+ STATES[8]*CONSTANTS[11]; RATES[7] = (((((( STATES[0]*STATES[3]*CONSTANTS[12] - STATES[7]*CONSTANTS[13]) - STATES[7]*CONSTANTS[4])+ STATES[8]*CONSTANTS[5]) - STATES[7]*STATES[1]*CONSTANTS[6])+ STATES[8]*CONSTANTS[7]) - STATES[2]*STATES[7]*CONSTANTS[14])+ STATES[6]*CONSTANTS[15]; RATES[6] = (( STATES[2]*STATES[7]*CONSTANTS[14] - STATES[6]*CONSTANTS[15]) - STATES[6]*CONSTANTS[8])+ STATES[8]*CONSTANTS[9]; RATES[8] = (((((((((( STATES[0]*CONSTANTS[0] - STATES[8]*CONSTANTS[1])+ STATES[0]*STATES[1]*CONSTANTS[2]) - STATES[8]*CONSTANTS[3])+ STATES[7]*CONSTANTS[4]) - STATES[8]*CONSTANTS[5])+ STATES[7]*STATES[1]*CONSTANTS[6]) - STATES[8]*CONSTANTS[7])+ STATES[6]*CONSTANTS[8]) - STATES[8]*CONSTANTS[9])+ STATES[5]*CONSTANTS[10]) - STATES[8]*CONSTANTS[11]; RATES[1] = (((( - STATES[0]*STATES[1]*CONSTANTS[2]+ STATES[8]*CONSTANTS[3]) - STATES[7]*STATES[1]*CONSTANTS[6])+ STATES[8]*CONSTANTS[7]) - STATES[4]*STATES[1]*CONSTANTS[18])+ STATES[5]*CONSTANTS[19]; RATES[2] = (( - STATES[0]*STATES[2]*CONSTANTS[16]+ STATES[4]*CONSTANTS[17]) - STATES[2]*STATES[7]*CONSTANTS[14])+ STATES[6]*CONSTANTS[15]; RATES[3] = - STATES[0]*STATES[3]*CONSTANTS[12]+ STATES[7]*CONSTANTS[13]; } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { }