Generated Code

The following is c code generated by the CellML API from this CellML file. (Back to language selection)

The raw code is available.

/*
   There are a total of 29 entries in the algebraic variable array.
   There are a total of 29 entries in each of the rate and state variable arrays.
   There are a total of 67 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * STATES[0] is Rs in component Rs (dimensionless).
 * ALGEBRAIC[0] is v1 in component v1 (first_order_rate_constant).
 * ALGEBRAIC[7] is v3 in component v3 (first_order_rate_constant).
 * STATES[1] is RL in component RL (dimensionless).
 * ALGEBRAIC[4] is v2 in component v2 (first_order_rate_constant).
 * ALGEBRAIC[10] is v4 in component v4 (first_order_rate_constant).
 * STATES[2] is Ri in component Ri (dimensionless).
 * ALGEBRAIC[15] is v6 in component v6 (first_order_rate_constant).
 * STATES[3] is L in component L (dimensionless).
 * STATES[4] is R2L2 in component R2L2 (dimensionless).
 * ALGEBRAIC[13] is v5 in component v5 (first_order_rate_constant).
 * ALGEBRAIC[18] is v7 in component v7 (first_order_rate_constant).
 * STATES[5] is R2_CPP in component R2_CPP (dimensionless).
 * ALGEBRAIC[22] is v8 in component v8 (first_order_rate_constant).
 * STATES[6] is Li in component Li (dimensionless).
 * STATES[7] is R2i in component R2i (dimensionless).
 * STATES[8] is Shc in component Shc (dimensionless).
 * ALGEBRAIC[1] is v9 in component v9 (first_order_rate_constant).
 * ALGEBRAIC[5] is v10 in component v10 (first_order_rate_constant).
 * STATES[9] is ShcP in component ShcP (dimensionless).
 * ALGEBRAIC[8] is v11 in component v11 (first_order_rate_constant).
 * ALGEBRAIC[16] is v27 in component v27 (first_order_rate_constant).
 * STATES[10] is ShcGS in component ShcGS (dimensionless).
 * ALGEBRAIC[14] is v13 in component v13 (first_order_rate_constant).
 * ALGEBRAIC[11] is v12 in component v12 (first_order_rate_constant).
 * STATES[11] is GS in component GS (dimensionless).
 * ALGEBRAIC[19] is v28 in component v28 (first_order_rate_constant).
 * STATES[12] is GSP in component GSP (dimensionless).
 * STATES[13] is RasGDP in component RasGDP (dimensionless).
 * ALGEBRAIC[20] is v15 in component v15 (first_order_rate_constant).
 * STATES[14] is Ras_ShcGS in component Ras_ShcGS (dimensionless).
 * STATES[15] is RasGTP in component RasGTP (dimensionless).
 * ALGEBRAIC[23] is v17 in component v17 (first_order_rate_constant).
 * ALGEBRAIC[17] is v14 in component v14 (first_order_rate_constant).
 * ALGEBRAIC[21] is v16 in component v16 (first_order_rate_constant).
 * STATES[16] is GAP in component GAP (dimensionless).
 * STATES[17] is Ras_GAP in component Ras_GAP (dimensionless).
 * STATES[18] is Raf in component Raf (dimensionless).
 * ALGEBRAIC[24] is v18 in component v18 (first_order_rate_constant).
 * STATES[19] is Ras_Raf in component Ras_Raf (dimensionless).
 * STATES[20] is Rafa in component Rafa (dimensionless).
 * ALGEBRAIC[25] is v19 in component v19 (first_order_rate_constant).
 * ALGEBRAIC[27] is v21 in component v21 (first_order_rate_constant).
 * STATES[21] is MEK in component MEK (dimensionless).
 * ALGEBRAIC[26] is v20 in component v20 (first_order_rate_constant).
 * STATES[22] is MEKP in component MEKP (dimensionless).
 * ALGEBRAIC[28] is v22 in component v22 (first_order_rate_constant).
 * STATES[23] is MEKPP in component MEKPP (dimensionless).
 * STATES[24] is ERK in component ERK (dimensionless).
 * ALGEBRAIC[6] is v24 in component v24 (first_order_rate_constant).
 * ALGEBRAIC[2] is v23 in component v23 (first_order_rate_constant).
 * STATES[25] is ERKP in component ERKP (dimensionless).
 * ALGEBRAIC[12] is v26 in component v26 (first_order_rate_constant).
 * ALGEBRAIC[9] is v25 in component v25 (first_order_rate_constant).
 * STATES[26] is ERKPP in component ERKPP (dimensionless).
 * STATES[27] is t in component t (dimensionless).
 * CONSTANTS[65] is v29 in component v29 (first_order_rate_constant).
 * STATES[28] is X in component X (dimensionless).
 * ALGEBRAIC[3] is v30 in component v30 (first_order_rate_constant).
 * CONSTANTS[0] is k1 in component v1 (first_order_rate_constant).
 * CONSTANTS[1] is kn1 in component v1 (first_order_rate_constant).
 * CONSTANTS[2] is DT in component v2 (dimensionless).
 * CONSTANTS[3] is E in component v2 (dimensionless).
 * CONSTANTS[4] is k2 in component v2 (first_order_rate_constant).
 * CONSTANTS[5] is f in component v2 (dimensionless).
 * CONSTANTS[6] is DT in component v3 (dimensionless).
 * CONSTANTS[7] is E in component v3 (dimensionless).
 * CONSTANTS[8] is kn3 in component v3 (first_order_rate_constant).
 * CONSTANTS[9] is f in component v3 (dimensionless).
 * CONSTANTS[10] is k3 in component v3 (first_order_rate_constant).
 * CONSTANTS[11] is k2_4 in component v4 (first_order_rate_constant).
 * CONSTANTS[12] is k4 in component v4 (first_order_rate_constant).
 * CONSTANTS[13] is DT in component v5 (dimensionless).
 * CONSTANTS[14] is E in component v5 (dimensionless).
 * CONSTANTS[15] is k5 in component v5 (first_order_rate_constant).
 * CONSTANTS[16] is f in component v5 (dimensionless).
 * CONSTANTS[17] is DT in component v6 (dimensionless).
 * CONSTANTS[18] is E in component v6 (dimensionless).
 * CONSTANTS[19] is k6 in component v6 (first_order_rate_constant).
 * CONSTANTS[20] is k7 in component v7 (first_order_rate_constant).
 * CONSTANTS[21] is f in component v7 (dimensionless).
 * CONSTANTS[22] is kn7 in component v7 (first_order_rate_constant).
 * CONSTANTS[23] is DT in component v8 (dimensionless).
 * CONSTANTS[24] is E in component v8 (dimensionless).
 * CONSTANTS[25] is k8 in component v8 (first_order_rate_constant).
 * CONSTANTS[26] is k9 in component v9 (first_order_rate_constant).
 * CONSTANTS[27] is K_9 in component v9 (dimensionless).
 * CONSTANTS[28] is V_10 in component v10 (first_order_rate_constant).
 * CONSTANTS[29] is K_10 in component v10 (dimensionless).
 * CONSTANTS[30] is k11 in component v11 (first_order_rate_constant).
 * CONSTANTS[31] is kn11 in component v11 (first_order_rate_constant).
 * CONSTANTS[32] is k12 in component v12 (first_order_rate_constant).
 * CONSTANTS[33] is kn12 in component v12 (first_order_rate_constant).
 * CONSTANTS[34] is k_13 in component v13 (first_order_rate_constant).
 * CONSTANTS[35] is k14 in component v14 (first_order_rate_constant).
 * CONSTANTS[36] is kn14 in component v14 (first_order_rate_constant).
 * CONSTANTS[37] is k15 in component v15 (first_order_rate_constant).
 * CONSTANTS[38] is k16 in component v16 (first_order_rate_constant).
 * CONSTANTS[39] is kn16 in component v16 (first_order_rate_constant).
 * CONSTANTS[40] is k17 in component v17 (first_order_rate_constant).
 * CONSTANTS[41] is V_18 in component v18 (first_order_rate_constant).
 * CONSTANTS[42] is K_18 in component v18 (dimensionless).
 * CONSTANTS[43] is k19 in component v19 (first_order_rate_constant).
 * CONSTANTS[44] is K_19 in component v19 (dimensionless).
 * CONSTANTS[45] is V_20 in component v20 (first_order_rate_constant).
 * CONSTANTS[46] is K_20 in component v20 (dimensionless).
 * CONSTANTS[47] is k21 in component v21 (first_order_rate_constant).
 * CONSTANTS[48] is K_21 in component v21 (dimensionless).
 * CONSTANTS[49] is V_22 in component v22 (first_order_rate_constant).
 * CONSTANTS[50] is K_22 in component v22 (dimensionless).
 * CONSTANTS[51] is k23 in component v23 (first_order_rate_constant).
 * CONSTANTS[52] is K_23 in component v23 (dimensionless).
 * CONSTANTS[53] is V_24 in component v24 (first_order_rate_constant).
 * CONSTANTS[54] is K_24 in component v24 (dimensionless).
 * CONSTANTS[55] is k25 in component v25 (first_order_rate_constant).
 * CONSTANTS[56] is K_25 in component v25 (dimensionless).
 * CONSTANTS[57] is V_26 in component v26 (first_order_rate_constant).
 * CONSTANTS[58] is K_26 in component v26 (dimensionless).
 * CONSTANTS[59] is k27 in component v27 (first_order_rate_constant).
 * CONSTANTS[60] is K_27 in component v27 (dimensionless).
 * CONSTANTS[61] is V_28 in component v28 (first_order_rate_constant).
 * CONSTANTS[62] is K_28 in component v28 (dimensionless).
 * CONSTANTS[63] is v_1 in component v29 (first_order_rate_constant).
 * CONSTANTS[64] is k_11 in component v30 (first_order_rate_constant).
 * RATES[0] is d/dt Rs in component Rs (dimensionless).
 * RATES[1] is d/dt RL in component RL (dimensionless).
 * RATES[2] is d/dt Ri in component Ri (dimensionless).
 * RATES[3] is d/dt L in component L (dimensionless).
 * RATES[4] is d/dt R2L2 in component R2L2 (dimensionless).
 * RATES[5] is d/dt R2_CPP in component R2_CPP (dimensionless).
 * RATES[6] is d/dt Li in component Li (dimensionless).
 * RATES[7] is d/dt R2i in component R2i (dimensionless).
 * RATES[8] is d/dt Shc in component Shc (dimensionless).
 * RATES[9] is d/dt ShcP in component ShcP (dimensionless).
 * RATES[10] is d/dt ShcGS in component ShcGS (dimensionless).
 * RATES[11] is d/dt GS in component GS (dimensionless).
 * RATES[12] is d/dt GSP in component GSP (dimensionless).
 * RATES[13] is d/dt RasGDP in component RasGDP (dimensionless).
 * RATES[14] is d/dt Ras_ShcGS in component Ras_ShcGS (dimensionless).
 * RATES[15] is d/dt RasGTP in component RasGTP (dimensionless).
 * RATES[16] is d/dt GAP in component GAP (dimensionless).
 * RATES[17] is d/dt Ras_GAP in component Ras_GAP (dimensionless).
 * RATES[18] is d/dt Raf in component Raf (dimensionless).
 * RATES[19] is d/dt Ras_Raf in component Ras_Raf (dimensionless).
 * RATES[20] is d/dt Rafa in component Rafa (dimensionless).
 * RATES[21] is d/dt MEK in component MEK (dimensionless).
 * RATES[22] is d/dt MEKP in component MEKP (dimensionless).
 * RATES[23] is d/dt MEKPP in component MEKPP (dimensionless).
 * RATES[24] is d/dt ERK in component ERK (dimensionless).
 * RATES[25] is d/dt ERKP in component ERKP (dimensionless).
 * RATES[26] is d/dt ERKPP in component ERKPP (dimensionless).
 * RATES[27] is d/dt t in component t (dimensionless).
 * RATES[28] is d/dt X in component X (dimensionless).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 11100;
STATES[1] = 0;
STATES[2] = 3900;
STATES[3] = 0.0000001;
STATES[4] = 0;
STATES[5] = 0;
STATES[6] = 0;
STATES[7] = 0;
STATES[8] = 30000;
STATES[9] = 0;
STATES[10] = 0;
STATES[11] = 20000;
STATES[12] = 0;
STATES[13] = 19800;
STATES[14] = 0;
STATES[15] = 200;
STATES[16] = 15000;
STATES[17] = 0;
STATES[18] = 10000;
STATES[19] = 0;
STATES[20] = 0;
STATES[21] = 360000;
STATES[22] = 0;
STATES[23] = 0;
STATES[24] = 750000;
STATES[25] = 0;
STATES[26] = 0;
STATES[27] = 0;
STATES[28] = 0;
CONSTANTS[0] = 384210000;
CONSTANTS[1] = 0.73;
CONSTANTS[2] = 6.5;
CONSTANTS[3] = 0.12;
CONSTANTS[4] = 0.7;
CONSTANTS[5] = 0.2;
CONSTANTS[6] = 6.5;
CONSTANTS[7] = 0.12;
CONSTANTS[8] = 0.7;
CONSTANTS[9] = 0.2;
CONSTANTS[10] = 0.0484;
CONSTANTS[11] = 0.0000001;
CONSTANTS[12] = 0.001383;
CONSTANTS[13] = 6.5;
CONSTANTS[14] = 0.12;
CONSTANTS[15] = 0.35;
CONSTANTS[16] = 0.2;
CONSTANTS[17] = 6.5;
CONSTANTS[18] = 0.12;
CONSTANTS[19] = 0.35;
CONSTANTS[20] = 1;
CONSTANTS[21] = 0.2;
CONSTANTS[22] = 0.000347;
CONSTANTS[23] = 6.5;
CONSTANTS[24] = 0.12;
CONSTANTS[25] = 0.35;
CONSTANTS[26] = 12;
CONSTANTS[27] = 6000;
CONSTANTS[28] = 300000;
CONSTANTS[29] = 6000;
CONSTANTS[30] = 0.002;
CONSTANTS[31] = 3.81;
CONSTANTS[32] = 0.0163;
CONSTANTS[33] = 10;
CONSTANTS[34] = 15;
CONSTANTS[35] = 0.005;
CONSTANTS[36] = 60;
CONSTANTS[37] = 720;
CONSTANTS[38] = 0.0012;
CONSTANTS[39] = 3;
CONSTANTS[40] = 27;
CONSTANTS[41] = 97000;
CONSTANTS[42] = 6000;
CONSTANTS[43] = 50;
CONSTANTS[44] = 9000;
CONSTANTS[45] = 920000;
CONSTANTS[46] = 600000;
CONSTANTS[47] = 50;
CONSTANTS[48] = 9000;
CONSTANTS[49] = 920000;
CONSTANTS[50] = 600000;
CONSTANTS[51] = 8.3;
CONSTANTS[52] = 90000;
CONSTANTS[53] = 200000;
CONSTANTS[54] = 600000;
CONSTANTS[55] = 8.3;
CONSTANTS[56] = 90000;
CONSTANTS[57] = 400000;
CONSTANTS[58] = 600000;
CONSTANTS[59] = 1.6;
CONSTANTS[60] = 600000;
CONSTANTS[61] = 75;
CONSTANTS[62] = 20000;
CONSTANTS[63] = 1;
CONSTANTS[64] = 0;
CONSTANTS[65] = CONSTANTS[63];
CONSTANTS[66] = CONSTANTS[65];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[27] = CONSTANTS[66];
ALGEBRAIC[0] =  CONSTANTS[0]*STATES[0]*STATES[3] -  CONSTANTS[1]*STATES[1];
RATES[3] = - ALGEBRAIC[0];
ALGEBRAIC[3] =  CONSTANTS[64]*STATES[27];
RATES[28] = ALGEBRAIC[3];
ALGEBRAIC[1] = ( CONSTANTS[26]*2.00000*(STATES[4]+STATES[7]+STATES[5])*STATES[8])/(CONSTANTS[27]+STATES[8]);
ALGEBRAIC[5] = ( CONSTANTS[28]*STATES[9])/(CONSTANTS[29]+STATES[9]);
RATES[8] = ALGEBRAIC[5] - ALGEBRAIC[1];
ALGEBRAIC[6] = ( CONSTANTS[53]*STATES[25])/(CONSTANTS[54]+STATES[25]);
ALGEBRAIC[2] = ( CONSTANTS[51]*STATES[24]*(STATES[22]+STATES[23]))/(CONSTANTS[52]+STATES[24]);
RATES[24] = ALGEBRAIC[6] - ALGEBRAIC[2];
ALGEBRAIC[7] =  CONSTANTS[10]*STATES[2] -  CONSTANTS[9]*CONSTANTS[8]*(CONSTANTS[7]+ (1.00000 - CONSTANTS[7])*(1.00000 - exp(- pow(STATES[27]/CONSTANTS[6], 3.00000))))*STATES[0];
RATES[0] = ALGEBRAIC[7] - ALGEBRAIC[0];
ALGEBRAIC[4] =  CONSTANTS[5]*CONSTANTS[4]*(CONSTANTS[3]+ (1.00000 - CONSTANTS[3])*(1.00000 - exp(- pow(STATES[27]/CONSTANTS[2], 3.00000))))*STATES[1];
ALGEBRAIC[10] =  CONSTANTS[12]*STATES[1]*STATES[1] -  CONSTANTS[11]*STATES[4];
RATES[1] = (ALGEBRAIC[0] - ALGEBRAIC[4]) - ALGEBRAIC[10];
ALGEBRAIC[12] = ( CONSTANTS[57]*STATES[26])/(CONSTANTS[58]+STATES[26]);
ALGEBRAIC[9] = ( CONSTANTS[55]*STATES[25]*(STATES[22]+STATES[23]))/(CONSTANTS[56]+STATES[25]);
RATES[25] = ((ALGEBRAIC[2]+ALGEBRAIC[12]) - ALGEBRAIC[6]) - ALGEBRAIC[9];
RATES[26] = ALGEBRAIC[9] - ALGEBRAIC[12];
ALGEBRAIC[14] =  CONSTANTS[34]*STATES[14];
ALGEBRAIC[11] =  CONSTANTS[32]*STATES[13]*STATES[10] -  CONSTANTS[33]*STATES[14];
RATES[14] = ALGEBRAIC[11] - ALGEBRAIC[14];
ALGEBRAIC[15] =  CONSTANTS[19]*(CONSTANTS[18]+ (1.00000 - CONSTANTS[18])*(1.00000 - exp(- pow(STATES[27]/CONSTANTS[17], 3.00000))))*STATES[7];
RATES[2] = ALGEBRAIC[4]+ALGEBRAIC[7]+ALGEBRAIC[15];
RATES[6] = ALGEBRAIC[4]+ALGEBRAIC[15];
ALGEBRAIC[8] =  CONSTANTS[30]*STATES[9]*STATES[11] -  CONSTANTS[31]*STATES[10];
ALGEBRAIC[16] = ( STATES[26]*STATES[10]*CONSTANTS[59])/(CONSTANTS[60]+STATES[10]);
RATES[9] = (ALGEBRAIC[16] - ALGEBRAIC[5]) - ALGEBRAIC[8];
RATES[10] = ((ALGEBRAIC[14]+ALGEBRAIC[8]) - ALGEBRAIC[16]) - ALGEBRAIC[11];
ALGEBRAIC[13] =  CONSTANTS[16]*CONSTANTS[15]*(CONSTANTS[14]+ (1.00000 - CONSTANTS[14])*(1.00000 - exp(- pow(STATES[27]/CONSTANTS[13], 3.00000))))*STATES[4];
ALGEBRAIC[18] =  CONSTANTS[20]*CONSTANTS[21]*STATES[4] -  CONSTANTS[22]*STATES[5];
RATES[4] = (ALGEBRAIC[10] - ALGEBRAIC[13]) - ALGEBRAIC[18];
ALGEBRAIC[19] = ( CONSTANTS[61]*STATES[12])/(CONSTANTS[62]+STATES[12]);
RATES[11] = ALGEBRAIC[19] - ALGEBRAIC[8];
RATES[12] = ALGEBRAIC[16] - ALGEBRAIC[19];
ALGEBRAIC[20] =  CONSTANTS[37]*STATES[17];
RATES[13] = ALGEBRAIC[20] - ALGEBRAIC[11];
ALGEBRAIC[17] =  CONSTANTS[35]*STATES[15]*STATES[16] -  CONSTANTS[36]*STATES[17];
RATES[16] = ALGEBRAIC[20] - ALGEBRAIC[17];
RATES[17] = ALGEBRAIC[17] - ALGEBRAIC[20];
ALGEBRAIC[22] =  CONSTANTS[25]*(CONSTANTS[24]+ (1.00000 - CONSTANTS[24])*(1.00000 - exp(- pow(STATES[27]/CONSTANTS[23], 3.00000))))*STATES[5];
RATES[5] = ALGEBRAIC[18] - ALGEBRAIC[22];
RATES[7] = (ALGEBRAIC[13]+ALGEBRAIC[22]) - ALGEBRAIC[15];
ALGEBRAIC[23] =  CONSTANTS[40]*STATES[19];
ALGEBRAIC[21] =  CONSTANTS[38]*STATES[15]*STATES[18] -  CONSTANTS[39]*STATES[19];
RATES[15] = ((ALGEBRAIC[14]+ALGEBRAIC[23]) - ALGEBRAIC[17]) - ALGEBRAIC[21];
RATES[19] = ALGEBRAIC[21] - ALGEBRAIC[23];
ALGEBRAIC[24] = ( CONSTANTS[41]*STATES[20])/(CONSTANTS[42]+STATES[20]);
RATES[18] = ALGEBRAIC[24] - ALGEBRAIC[21];
ALGEBRAIC[25] = ( STATES[20]*STATES[21]*CONSTANTS[43])/(CONSTANTS[44]+STATES[21]);
ALGEBRAIC[26] = ( CONSTANTS[45]*STATES[22])/(CONSTANTS[46]+STATES[22]);
RATES[21] = ALGEBRAIC[26] - ALGEBRAIC[25];
ALGEBRAIC[27] = ( STATES[20]*STATES[22]*CONSTANTS[47])/(CONSTANTS[48]+STATES[22]);
RATES[20] = ((ALGEBRAIC[23] - ALGEBRAIC[25]) - ALGEBRAIC[27]) - ALGEBRAIC[24];
ALGEBRAIC[28] = ( CONSTANTS[49]*STATES[23])/(CONSTANTS[50]+STATES[23]);
RATES[22] = ((ALGEBRAIC[25]+ALGEBRAIC[28]) - ALGEBRAIC[26]) - ALGEBRAIC[27];
RATES[23] = ALGEBRAIC[27] - ALGEBRAIC[28];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] =  CONSTANTS[0]*STATES[0]*STATES[3] -  CONSTANTS[1]*STATES[1];
ALGEBRAIC[3] =  CONSTANTS[64]*STATES[27];
ALGEBRAIC[1] = ( CONSTANTS[26]*2.00000*(STATES[4]+STATES[7]+STATES[5])*STATES[8])/(CONSTANTS[27]+STATES[8]);
ALGEBRAIC[5] = ( CONSTANTS[28]*STATES[9])/(CONSTANTS[29]+STATES[9]);
ALGEBRAIC[6] = ( CONSTANTS[53]*STATES[25])/(CONSTANTS[54]+STATES[25]);
ALGEBRAIC[2] = ( CONSTANTS[51]*STATES[24]*(STATES[22]+STATES[23]))/(CONSTANTS[52]+STATES[24]);
ALGEBRAIC[7] =  CONSTANTS[10]*STATES[2] -  CONSTANTS[9]*CONSTANTS[8]*(CONSTANTS[7]+ (1.00000 - CONSTANTS[7])*(1.00000 - exp(- pow(STATES[27]/CONSTANTS[6], 3.00000))))*STATES[0];
ALGEBRAIC[4] =  CONSTANTS[5]*CONSTANTS[4]*(CONSTANTS[3]+ (1.00000 - CONSTANTS[3])*(1.00000 - exp(- pow(STATES[27]/CONSTANTS[2], 3.00000))))*STATES[1];
ALGEBRAIC[10] =  CONSTANTS[12]*STATES[1]*STATES[1] -  CONSTANTS[11]*STATES[4];
ALGEBRAIC[12] = ( CONSTANTS[57]*STATES[26])/(CONSTANTS[58]+STATES[26]);
ALGEBRAIC[9] = ( CONSTANTS[55]*STATES[25]*(STATES[22]+STATES[23]))/(CONSTANTS[56]+STATES[25]);
ALGEBRAIC[14] =  CONSTANTS[34]*STATES[14];
ALGEBRAIC[11] =  CONSTANTS[32]*STATES[13]*STATES[10] -  CONSTANTS[33]*STATES[14];
ALGEBRAIC[15] =  CONSTANTS[19]*(CONSTANTS[18]+ (1.00000 - CONSTANTS[18])*(1.00000 - exp(- pow(STATES[27]/CONSTANTS[17], 3.00000))))*STATES[7];
ALGEBRAIC[8] =  CONSTANTS[30]*STATES[9]*STATES[11] -  CONSTANTS[31]*STATES[10];
ALGEBRAIC[16] = ( STATES[26]*STATES[10]*CONSTANTS[59])/(CONSTANTS[60]+STATES[10]);
ALGEBRAIC[13] =  CONSTANTS[16]*CONSTANTS[15]*(CONSTANTS[14]+ (1.00000 - CONSTANTS[14])*(1.00000 - exp(- pow(STATES[27]/CONSTANTS[13], 3.00000))))*STATES[4];
ALGEBRAIC[18] =  CONSTANTS[20]*CONSTANTS[21]*STATES[4] -  CONSTANTS[22]*STATES[5];
ALGEBRAIC[19] = ( CONSTANTS[61]*STATES[12])/(CONSTANTS[62]+STATES[12]);
ALGEBRAIC[20] =  CONSTANTS[37]*STATES[17];
ALGEBRAIC[17] =  CONSTANTS[35]*STATES[15]*STATES[16] -  CONSTANTS[36]*STATES[17];
ALGEBRAIC[22] =  CONSTANTS[25]*(CONSTANTS[24]+ (1.00000 - CONSTANTS[24])*(1.00000 - exp(- pow(STATES[27]/CONSTANTS[23], 3.00000))))*STATES[5];
ALGEBRAIC[23] =  CONSTANTS[40]*STATES[19];
ALGEBRAIC[21] =  CONSTANTS[38]*STATES[15]*STATES[18] -  CONSTANTS[39]*STATES[19];
ALGEBRAIC[24] = ( CONSTANTS[41]*STATES[20])/(CONSTANTS[42]+STATES[20]);
ALGEBRAIC[25] = ( STATES[20]*STATES[21]*CONSTANTS[43])/(CONSTANTS[44]+STATES[21]);
ALGEBRAIC[26] = ( CONSTANTS[45]*STATES[22])/(CONSTANTS[46]+STATES[22]);
ALGEBRAIC[27] = ( STATES[20]*STATES[22]*CONSTANTS[47])/(CONSTANTS[48]+STATES[22]);
ALGEBRAIC[28] = ( CONSTANTS[49]*STATES[23])/(CONSTANTS[50]+STATES[23]);
}