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 53 entries in the algebraic variable array.
There are a total of 22 entries in each of the rate and state variable arrays.
There are a total of 83 entries in the constant variable array.
*/
/*
* VOI is time in component environment (hour).
* CONSTANTS[0] is eps in component E (dimensionless).
* CONSTANTS[1] is K1 in component A (per_hour).
* CONSTANTS[2] is K1a in component A (per_hour).
* CONSTANTS[3] is K2 in component A (per_hour).
* CONSTANTS[4] is K2a in component A (per_hour).
* CONSTANTS[5] is K2aa in component A (per_hour).
* CONSTANTS[6] is K3 in component A (per_hour).
* CONSTANTS[7] is K3a in component A (per_hour).
* CONSTANTS[8] is K4 in component A (per_hour).
* CONSTANTS[9] is J1 in component A (dimensionless).
* CONSTANTS[10] is J3 in component A (dimensionless).
* CONSTANTS[11] is J4 in component A (dimensionless).
* ALGEBRAIC[1] is v_20 in component A (per_hour).
* ALGEBRAIC[21] is v_21 in component A (per_hour).
* ALGEBRAIC[0] is v_42 in component A (per_hour).
* ALGEBRAIC[22] is v_19 in component A (per_hour).
* CONSTANTS[12] is GA in component A (dimensionless).
* CONSTANTS[13] is GB in component A (dimensionless).
* CONSTANTS[14] is GE in component A (dimensionless).
* ALGEBRAIC[12] is V_2 in component A (per_hour).
* ALGEBRAIC[13] is V_4 in component A (per_hour).
* STATES[0] is Cdh1 in component A (dimensionless).
* CONSTANTS[15] is Cdc20 in component D (dimensionless).
* STATES[1] is CYCA in component D (dimensionless).
* STATES[2] is CYCB in component A (dimensionless).
* STATES[3] is CYCE in component D (dimensionless).
* CONSTANTS[16] is k15 in component B (per_hour).
* CONSTANTS[17] is k16 in component B (per_hour).
* CONSTANTS[18] is k17 in component B (per_hour).
* CONSTANTS[19] is k17a in component B (per_hour).
* CONSTANTS[20] is k18 in component B (per_hour).
* CONSTANTS[21] is K9 in component B (per_hour).
* CONSTANTS[22] is J17 in component B (dimensionless).
* CONSTANTS[23] is J15 in component B (dimensionless).
* ALGEBRAIC[14] is v_1 in component B (per_hour).
* ALGEBRAIC[15] is v_2 in component B (per_hour).
* ALGEBRAIC[2] is v_34 in component B (per_hour).
* ALGEBRAIC[5] is v_39 in component B (per_hour).
* ALGEBRAIC[3] is v_41 in component B (per_hour).
* STATES[4] is DRG in component B (dimensionless).
* STATES[5] is ERG in component B (dimensionless).
* CONSTANTS[24] is K19 in component C (per_hour).
* CONSTANTS[25] is K19a in component C (per_hour).
* CONSTANTS[26] is K20 in component C (per_hour).
* CONSTANTS[27] is K21 in component C (dimensionless).
* CONSTANTS[28] is K22 in component C (per_hour).
* CONSTANTS[29] is K23a in component C (per_hour).
* CONSTANTS[30] is K23 in component C (per_hour).
* CONSTANTS[31] is K26 in component C (per_hour).
* CONSTANTS[32] is K26R in component C (per_hour).
* ALGEBRAIC[16] is v_29 in component C (per_hour).
* ALGEBRAIC[23] is v_30 in component C (per_hour).
* ALGEBRAIC[25] is v_43 in component C (per_hour).
* ALGEBRAIC[27] is v_44 in component C (per_hour).
* ALGEBRAIC[29] is v_45 in component C (per_hour).
* ALGEBRAIC[31] is v_46 in component C (per_hour).
* ALGEBRAIC[33] is v_47 in component C (per_hour).
* ALGEBRAIC[35] is v_48 in component C (per_hour).
* ALGEBRAIC[37] is v_49 in component C (per_hour).
* ALGEBRAIC[40] is v_50 in component C (per_hour).
* ALGEBRAIC[42] is v_51 in component C (per_hour).
* ALGEBRAIC[44] is v_52 in component C (per_hour).
* CONSTANTS[33] is LA in component C (dimensionless).
* CONSTANTS[34] is LB in component C (dimensionless).
* CONSTANTS[35] is LD in component C (dimensionless).
* CONSTANTS[36] is LE in component C (dimensionless).
* CONSTANTS[37] is FE in component C (dimensionless).
* CONSTANTS[38] is FB in component C (dimensionless).
* STATES[6] is Rb in component C (dimensionless).
* STATES[7] is PPRb in component C (dimensionless).
* STATES[8] is E2F in component C (dimensionless).
* STATES[9] is PE2F in component C (dimensionless).
* STATES[10] is E2FRb in component C (dimensionless).
* STATES[11] is PE2FRb in component C (dimensionless).
* ALGEBRAIC[4] is PP1A in component C (dimensionless).
* CONSTANTS[39] is PP1T in component C (dimensionless).
* CONSTANTS[40] is CYCDT in component C (dimensionless).
* CONSTANTS[41] is K5 in component D (per_hour).
* CONSTANTS[42] is K6a in component D (per_hour).
* CONSTANTS[43] is K6 in component D (per_hour).
* CONSTANTS[44] is K7a in component D (per_hour).
* CONSTANTS[45] is K7 in component D (per_hour).
* CONSTANTS[46] is K8a in component D (per_hour).
* CONSTANTS[47] is K8 in component D (per_hour).
* CONSTANTS[48] is K10 in component D (per_hour).
* CONSTANTS[49] is k24 in component D (per_hour).
* CONSTANTS[50] is k24r in component D (per_hour).
* CONSTANTS[51] is K25 in component D (per_hour).
* CONSTANTS[52] is K25R in component D (per_hour).
* CONSTANTS[53] is K29 in component D (per_hour).
* CONSTANTS[54] is K30 in component D (per_hour).
* CONSTANTS[55] is J8 in component D (dimensionless).
* ALGEBRAIC[26] is v_3 in component D (per_hour).
* ALGEBRAIC[28] is v_4 in component D (per_hour).
* ALGEBRAIC[30] is v_5 in component D (per_hour).
* ALGEBRAIC[32] is v_6 in component D (per_hour).
* ALGEBRAIC[34] is v_7 in component D (per_hour).
* ALGEBRAIC[36] is v_8 in component D (per_hour).
* ALGEBRAIC[38] is v_9 in component D (per_hour).
* ALGEBRAIC[41] is v_10 in component D (per_hour).
* ALGEBRAIC[43] is v_11 in component D (per_hour).
* ALGEBRAIC[45] is v_12 in component D (per_hour).
* ALGEBRAIC[46] is v_13 in component D (per_hour).
* ALGEBRAIC[49] is v_14 in component D (per_hour).
* ALGEBRAIC[48] is v_15 in component D (per_hour).
* ALGEBRAIC[51] is v_16 in component D (per_hour).
* ALGEBRAIC[39] is v_17 in component D (per_hour).
* ALGEBRAIC[47] is v_18 in component D (per_hour).
* ALGEBRAIC[50] is v_36 in component D (per_hour).
* ALGEBRAIC[52] is v_38 in component D (per_hour).
* CONSTANTS[79] is v_40 in component D (per_hour).
* CONSTANTS[56] is HA in component D (dimensionless).
* CONSTANTS[57] is HB in component D (dimensionless).
* CONSTANTS[58] is HE in component D (dimensionless).
* CONSTANTS[59] is YE in component D (dimensionless).
* CONSTANTS[60] is YB in component D (dimensionless).
* STATES[12] is p27 in component D (dimensionless).
* ALGEBRAIC[17] is V_6 in component D (per_hour).
* ALGEBRAIC[24] is V_8 in component D (per_hour).
* STATES[13] is MASS in component E (dimensionless).
* STATES[14] is CA in component D (dimensionless).
* STATES[15] is CD in component D (dimensionless).
* STATES[16] is CE in component D (dimensionless).
* STATES[17] is CYCD in component D (dimensionless).
* CONSTANTS[61] is CYCET in component D (dimensionless).
* CONSTANTS[62] is K27 in component E (per_hour).
* CONSTANTS[63] is K28 in component E (per_hour).
* ALGEBRAIC[6] is v_31 in component E (per_hour).
* ALGEBRAIC[18] is v_32 in component E (per_hour).
* ALGEBRAIC[7] is v_33 in component E (per_hour).
* CONSTANTS[64] is r31switch in component E (dimensionless).
* CONSTANTS[65] is MU in component E (per_hour).
* STATES[18] is GM in component E (dimensionless).
* CONSTANTS[66] is K11a in component F (per_hour).
* CONSTANTS[67] is K11 in component F (per_hour).
* CONSTANTS[68] is K12 in component F (per_hour).
* CONSTANTS[69] is K13 in component F (per_hour).
* CONSTANTS[70] is K14 in component F (per_hour).
* CONSTANTS[71] is K31 in component F (per_hour).
* CONSTANTS[72] is K32 in component F (per_hour).
* CONSTANTS[73] is K33 in component F (per_hour).
* CONSTANTS[74] is K34 in component F (per_hour).
* CONSTANTS[75] is J13 in component F (dimensionless).
* CONSTANTS[76] is J14 in component F (dimensionless).
* CONSTANTS[77] is J31 in component F (dimensionless).
* CONSTANTS[78] is J32 in component F (dimensionless).
* ALGEBRAIC[8] is v_22 in component F (per_hour).
* ALGEBRAIC[9] is v_23 in component F (per_hour).
* ALGEBRAIC[19] is v_24 in component F (per_hour).
* ALGEBRAIC[11] is v_25 in component F (per_hour).
* ALGEBRAIC[10] is v_26 in component F (per_hour).
* CONSTANTS[80] is v_27 in component F (per_hour).
* CONSTANTS[81] is v_28 in component F (per_hour).
* ALGEBRAIC[20] is v_35 in component F (per_hour).
* CONSTANTS[82] is v_37 in component F (per_hour).
* STATES[19] is IEP in component F (dimensionless).
* STATES[20] is PPX in component F (dimensionless).
* STATES[21] is Cdc20T in component F (dimensionless).
* RATES[0] is d/dt Cdh1 in component A (dimensionless).
* RATES[2] is d/dt CYCB in component A (dimensionless).
* RATES[4] is d/dt DRG in component B (dimensionless).
* RATES[5] is d/dt ERG in component B (dimensionless).
* RATES[7] is d/dt PPRb in component C (dimensionless).
* RATES[8] is d/dt E2F in component C (dimensionless).
* RATES[9] is d/dt PE2F in component C (dimensionless).
* RATES[6] is d/dt Rb in component C (dimensionless).
* RATES[10] is d/dt E2FRb in component C (dimensionless).
* RATES[11] is d/dt PE2FRb in component C (dimensionless).
* RATES[14] is d/dt CA in component D (dimensionless).
* RATES[15] is d/dt CD in component D (dimensionless).
* RATES[16] is d/dt CE in component D (dimensionless).
* RATES[1] is d/dt CYCA in component D (dimensionless).
* RATES[17] is d/dt CYCD in component D (dimensionless).
* RATES[3] is d/dt CYCE in component D (dimensionless).
* RATES[12] is d/dt p27 in component D (dimensionless).
* RATES[18] is d/dt GM in component E (dimensionless).
* RATES[13] is d/dt MASS in component E (dimensionless).
* RATES[21] is d/dt Cdc20T in component F (dimensionless).
* RATES[19] is d/dt IEP in component F (dimensionless).
* RATES[20] is d/dt PPX in component F (dimensionless).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1;
CONSTANTS[1] = 0.6;
CONSTANTS[2] = 0.1;
CONSTANTS[3] = 20;
CONSTANTS[4] = 0.05;
CONSTANTS[5] = 1;
CONSTANTS[6] = 140;
CONSTANTS[7] = 7.5;
CONSTANTS[8] = 40;
CONSTANTS[9] = 0.1;
CONSTANTS[10] = 0.01;
CONSTANTS[11] = 0.04;
CONSTANTS[12] = 0.3;
CONSTANTS[13] = 1;
CONSTANTS[14] = 0;
STATES[0] = 0.000653278;
CONSTANTS[15] = 0.00220177;
STATES[1] = 1.4094;
STATES[2] = 2.72898;
STATES[3] = 0.0229112;
CONSTANTS[16] = 0.25;
CONSTANTS[17] = 0.25;
CONSTANTS[18] = 10;
CONSTANTS[19] = 0.35;
CONSTANTS[20] = 10;
CONSTANTS[21] = 2.5;
CONSTANTS[22] = 0.3;
CONSTANTS[23] = 0.1;
STATES[4] = 0.900533;
STATES[5] = 0.0121809;
CONSTANTS[24] = 20;
CONSTANTS[25] = 0;
CONSTANTS[26] = 10;
CONSTANTS[27] = 1;
CONSTANTS[28] = 1;
CONSTANTS[29] = 0.005;
CONSTANTS[30] = 1;
CONSTANTS[31] = 10000;
CONSTANTS[32] = 200;
CONSTANTS[33] = 3;
CONSTANTS[34] = 5;
CONSTANTS[35] = 3.3;
CONSTANTS[36] = 5;
CONSTANTS[37] = 25;
CONSTANTS[38] = 2;
STATES[6] = 0.000190871;
STATES[7] = 9.97574;
STATES[8] = 0.989986;
STATES[9] = 3.98594;
STATES[10] = 0.00478911;
STATES[11] = 0.0192822;
CONSTANTS[39] = 1;
CONSTANTS[40] = 0.010976;
CONSTANTS[41] = 20;
CONSTANTS[42] = 10;
CONSTANTS[43] = 100;
CONSTANTS[44] = 0;
CONSTANTS[45] = 0.6;
CONSTANTS[46] = 0.1;
CONSTANTS[47] = 2;
CONSTANTS[48] = 5;
CONSTANTS[49] = 1000;
CONSTANTS[50] = 10;
CONSTANTS[51] = 1000;
CONSTANTS[52] = 10;
CONSTANTS[53] = 0.05;
CONSTANTS[54] = 20;
CONSTANTS[55] = 0.1;
CONSTANTS[56] = 0.5;
CONSTANTS[57] = 1;
CONSTANTS[58] = 0.5;
CONSTANTS[59] = 1;
CONSTANTS[60] = 0.05;
STATES[12] = 0.00922806;
STATES[13] = 1.68776;
STATES[14] = 0.0356927;
STATES[15] = 0.010976;
STATES[16] = 0.000542587;
STATES[17] = 0.43929;
CONSTANTS[61] = 0.000542587;
CONSTANTS[62] = 0.2;
CONSTANTS[63] = 0.2;
CONSTANTS[64] = 1;
CONSTANTS[65] = 0.061;
STATES[18] = 1.35565;
CONSTANTS[66] = 0;
CONSTANTS[67] = 1.5;
CONSTANTS[68] = 1.5;
CONSTANTS[69] = 5;
CONSTANTS[70] = 2.5;
CONSTANTS[71] = 0.7;
CONSTANTS[72] = 1.8;
CONSTANTS[73] = 0.05;
CONSTANTS[74] = 0.05;
CONSTANTS[75] = 0.005;
CONSTANTS[76] = 0.005;
CONSTANTS[77] = 0.01;
CONSTANTS[78] = 0.01;
STATES[19] = 0.154655;
STATES[20] = 1;
STATES[21] = 2.36733;
CONSTANTS[79] = CONSTANTS[0]*CONSTANTS[41];
CONSTANTS[80] = ( CONSTANTS[70]*CONSTANTS[15])/(CONSTANTS[76]+CONSTANTS[15]);
CONSTANTS[81] = CONSTANTS[68]*CONSTANTS[15];
CONSTANTS[82] = CONSTANTS[0]*CONSTANTS[73];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[7] = CONSTANTS[0]*CONSTANTS[65]*STATES[18];
RATES[13] = ALGEBRAIC[7];
ALGEBRAIC[8] = CONSTANTS[74]*STATES[20];
RATES[20] = CONSTANTS[82] - ALGEBRAIC[8];
ALGEBRAIC[15] = CONSTANTS[20]*STATES[4];
ALGEBRAIC[3] = CONSTANTS[0]*(( CONSTANTS[18]*pow(STATES[4]/CONSTANTS[22], 2.00000))/(pow(STATES[4]/CONSTANTS[22], 2.00000)+1.00000)+ CONSTANTS[19]*STATES[5]);
RATES[4] = ALGEBRAIC[3] - ALGEBRAIC[15];
ALGEBRAIC[14] = CONSTANTS[17]*STATES[5];
ALGEBRAIC[2] = ( CONSTANTS[0]*CONSTANTS[16])/(pow(STATES[4]/CONSTANTS[23], 2.00000)+1.00000);
RATES[5] = ALGEBRAIC[2] - ALGEBRAIC[14];
ALGEBRAIC[6] = CONSTANTS[62]*STATES[13]*CONSTANTS[64];
ALGEBRAIC[18] = CONSTANTS[63]*STATES[18];
RATES[18] = ALGEBRAIC[6] - ALGEBRAIC[18];
ALGEBRAIC[11] = CONSTANTS[68]*STATES[21];
ALGEBRAIC[20] = CONSTANTS[0]*CONSTANTS[66]+ CONSTANTS[67]*STATES[2];
RATES[21] = ALGEBRAIC[20] - ALGEBRAIC[11];
ALGEBRAIC[9] = ( CONSTANTS[71]*STATES[2]*(1.00000 - STATES[19]))/((CONSTANTS[77] - STATES[19])+1.00000);
ALGEBRAIC[19] = ( CONSTANTS[72]*STATES[20]*STATES[19])/(CONSTANTS[78]+STATES[19]);
RATES[19] = ALGEBRAIC[9] - ALGEBRAIC[19];
ALGEBRAIC[1] = ( (CONSTANTS[7]+ CONSTANTS[6]*CONSTANTS[15])*(1.00000 - STATES[0]))/((CONSTANTS[10] - STATES[0])+1.00000);
ALGEBRAIC[13] = CONSTANTS[8]*( CONSTANTS[12]*STATES[1]+ CONSTANTS[13]*STATES[2]+ CONSTANTS[14]*STATES[3]);
ALGEBRAIC[21] = ( ALGEBRAIC[13]*STATES[0])/(CONSTANTS[11]+STATES[0]);
RATES[0] = ALGEBRAIC[1] - ALGEBRAIC[21];
ALGEBRAIC[0] = CONSTANTS[0]*(( CONSTANTS[1]*pow(STATES[2]/CONSTANTS[9], 2.00000))/(pow(STATES[2]/CONSTANTS[9], 2.00000)+1.00000)+CONSTANTS[2]);
ALGEBRAIC[12] = CONSTANTS[5]*CONSTANTS[15]+ CONSTANTS[4]*(1.00000 - STATES[0])+ CONSTANTS[3]*STATES[0];
ALGEBRAIC[22] = ALGEBRAIC[12]*STATES[2];
RATES[2] = ALGEBRAIC[0] - ALGEBRAIC[22];
ALGEBRAIC[16] = STATES[10]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[23] = STATES[11]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[25] = STATES[6]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[4] = CONSTANTS[39]/( CONSTANTS[27]*CONSTANTS[37]*(STATES[1]+STATES[3]+ CONSTANTS[38]*STATES[2])+1.00000);
ALGEBRAIC[27] = STATES[7]*( CONSTANTS[25]*(CONSTANTS[39] - ALGEBRAIC[4])+ CONSTANTS[24]*ALGEBRAIC[4]);
RATES[7] = (ALGEBRAIC[16]+ALGEBRAIC[23]+ALGEBRAIC[25]) - ALGEBRAIC[27];
ALGEBRAIC[29] = STATES[10]*CONSTANTS[32];
ALGEBRAIC[31] = STATES[8]*(CONSTANTS[29]+ CONSTANTS[30]*(STATES[1]+STATES[2]));
ALGEBRAIC[33] = STATES[9]*CONSTANTS[28];
ALGEBRAIC[35] = STATES[8]*STATES[6]*CONSTANTS[31];
RATES[8] = ((ALGEBRAIC[16]+ALGEBRAIC[29]+ALGEBRAIC[33]) - ALGEBRAIC[31]) - ALGEBRAIC[35];
ALGEBRAIC[26] = CONSTANTS[48]*STATES[15];
ALGEBRAIC[34] = CONSTANTS[49]*STATES[12]*STATES[17];
ALGEBRAIC[36] = CONSTANTS[50]*STATES[15];
ALGEBRAIC[17] = CONSTANTS[42]+ CONSTANTS[43]*( CONSTANTS[56]*STATES[1]+ CONSTANTS[57]*STATES[2]+ CONSTANTS[58]*STATES[3]);
ALGEBRAIC[39] = ALGEBRAIC[17]*STATES[15];
RATES[15] = ((ALGEBRAIC[34] - ALGEBRAIC[36]) - ALGEBRAIC[39]) - ALGEBRAIC[26];
ALGEBRAIC[5] = CONSTANTS[0]*CONSTANTS[21]*STATES[4];
ALGEBRAIC[28] = CONSTANTS[48]*STATES[17];
RATES[17] = ((ALGEBRAIC[5]+ALGEBRAIC[39]+ALGEBRAIC[36]) - ALGEBRAIC[34]) - ALGEBRAIC[28];
ALGEBRAIC[37] = STATES[11]*CONSTANTS[32];
ALGEBRAIC[40] = STATES[6]*STATES[9]*CONSTANTS[31];
RATES[9] = ((ALGEBRAIC[23]+ALGEBRAIC[37]+ALGEBRAIC[31]) - ALGEBRAIC[33]) - ALGEBRAIC[40];
RATES[6] = (((ALGEBRAIC[27]+ALGEBRAIC[29]+ALGEBRAIC[37]) - ALGEBRAIC[35]) - ALGEBRAIC[40]) - ALGEBRAIC[25];
ALGEBRAIC[42] = STATES[11]*CONSTANTS[28];
ALGEBRAIC[44] = STATES[10]*(CONSTANTS[29]+ CONSTANTS[30]*(STATES[1]+STATES[2]));
RATES[10] = (((ALGEBRAIC[42]+ALGEBRAIC[35]) - ALGEBRAIC[44]) - ALGEBRAIC[16]) - ALGEBRAIC[29];
RATES[11] = (((ALGEBRAIC[44]+ALGEBRAIC[40]) - ALGEBRAIC[42]) - ALGEBRAIC[23]) - ALGEBRAIC[37];
ALGEBRAIC[32] = CONSTANTS[51]*STATES[12]*STATES[1];
ALGEBRAIC[41] = CONSTANTS[54]*CONSTANTS[15]*STATES[14];
ALGEBRAIC[45] = CONSTANTS[52]*STATES[14];
ALGEBRAIC[47] = ALGEBRAIC[17]*STATES[14];
RATES[14] = ((ALGEBRAIC[32] - ALGEBRAIC[45]) - ALGEBRAIC[47]) - ALGEBRAIC[41];
ALGEBRAIC[38] = CONSTANTS[54]*CONSTANTS[15]*STATES[1];
ALGEBRAIC[50] = CONSTANTS[0]*CONSTANTS[53]*STATES[8]*STATES[13];
RATES[1] = ((ALGEBRAIC[50] - ALGEBRAIC[38]) - ALGEBRAIC[32])+ALGEBRAIC[45]+ALGEBRAIC[47];
ALGEBRAIC[30] = CONSTANTS[51]*STATES[12]*STATES[3];
ALGEBRAIC[43] = CONSTANTS[52]*STATES[16];
ALGEBRAIC[24] = ( CONSTANTS[47]*( CONSTANTS[59]*(STATES[1]+STATES[3])+ CONSTANTS[60]*STATES[2]))/(CONSTANTS[61]+CONSTANTS[55])+CONSTANTS[46];
ALGEBRAIC[46] = ALGEBRAIC[24]*STATES[16];
ALGEBRAIC[48] = ALGEBRAIC[17]*STATES[12];
RATES[12] = (((((CONSTANTS[79]+ALGEBRAIC[26]+ALGEBRAIC[36]) - ALGEBRAIC[48]) - ALGEBRAIC[30]) - ALGEBRAIC[32]) - ALGEBRAIC[34])+ALGEBRAIC[43]+ALGEBRAIC[45]+ALGEBRAIC[46]+ALGEBRAIC[41];
ALGEBRAIC[51] = ALGEBRAIC[17]*STATES[16];
RATES[16] = ((ALGEBRAIC[30] - ALGEBRAIC[43]) - ALGEBRAIC[46]) - ALGEBRAIC[51];
ALGEBRAIC[49] = ALGEBRAIC[24]*STATES[3];
ALGEBRAIC[52] = CONSTANTS[0]*(CONSTANTS[44]+ CONSTANTS[45]*STATES[8]);
RATES[3] = ((ALGEBRAIC[52] - ALGEBRAIC[49]) - ALGEBRAIC[30])+ALGEBRAIC[43]+ALGEBRAIC[51];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[7] = CONSTANTS[0]*CONSTANTS[65]*STATES[18];
ALGEBRAIC[8] = CONSTANTS[74]*STATES[20];
ALGEBRAIC[15] = CONSTANTS[20]*STATES[4];
ALGEBRAIC[3] = CONSTANTS[0]*(( CONSTANTS[18]*pow(STATES[4]/CONSTANTS[22], 2.00000))/(pow(STATES[4]/CONSTANTS[22], 2.00000)+1.00000)+ CONSTANTS[19]*STATES[5]);
ALGEBRAIC[14] = CONSTANTS[17]*STATES[5];
ALGEBRAIC[2] = ( CONSTANTS[0]*CONSTANTS[16])/(pow(STATES[4]/CONSTANTS[23], 2.00000)+1.00000);
ALGEBRAIC[6] = CONSTANTS[62]*STATES[13]*CONSTANTS[64];
ALGEBRAIC[18] = CONSTANTS[63]*STATES[18];
ALGEBRAIC[11] = CONSTANTS[68]*STATES[21];
ALGEBRAIC[20] = CONSTANTS[0]*CONSTANTS[66]+ CONSTANTS[67]*STATES[2];
ALGEBRAIC[9] = ( CONSTANTS[71]*STATES[2]*(1.00000 - STATES[19]))/((CONSTANTS[77] - STATES[19])+1.00000);
ALGEBRAIC[19] = ( CONSTANTS[72]*STATES[20]*STATES[19])/(CONSTANTS[78]+STATES[19]);
ALGEBRAIC[1] = ( (CONSTANTS[7]+ CONSTANTS[6]*CONSTANTS[15])*(1.00000 - STATES[0]))/((CONSTANTS[10] - STATES[0])+1.00000);
ALGEBRAIC[13] = CONSTANTS[8]*( CONSTANTS[12]*STATES[1]+ CONSTANTS[13]*STATES[2]+ CONSTANTS[14]*STATES[3]);
ALGEBRAIC[21] = ( ALGEBRAIC[13]*STATES[0])/(CONSTANTS[11]+STATES[0]);
ALGEBRAIC[0] = CONSTANTS[0]*(( CONSTANTS[1]*pow(STATES[2]/CONSTANTS[9], 2.00000))/(pow(STATES[2]/CONSTANTS[9], 2.00000)+1.00000)+CONSTANTS[2]);
ALGEBRAIC[12] = CONSTANTS[5]*CONSTANTS[15]+ CONSTANTS[4]*(1.00000 - STATES[0])+ CONSTANTS[3]*STATES[0];
ALGEBRAIC[22] = ALGEBRAIC[12]*STATES[2];
ALGEBRAIC[16] = STATES[10]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[23] = STATES[11]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[25] = STATES[6]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[4] = CONSTANTS[39]/( CONSTANTS[27]*CONSTANTS[37]*(STATES[1]+STATES[3]+ CONSTANTS[38]*STATES[2])+1.00000);
ALGEBRAIC[27] = STATES[7]*( CONSTANTS[25]*(CONSTANTS[39] - ALGEBRAIC[4])+ CONSTANTS[24]*ALGEBRAIC[4]);
ALGEBRAIC[29] = STATES[10]*CONSTANTS[32];
ALGEBRAIC[31] = STATES[8]*(CONSTANTS[29]+ CONSTANTS[30]*(STATES[1]+STATES[2]));
ALGEBRAIC[33] = STATES[9]*CONSTANTS[28];
ALGEBRAIC[35] = STATES[8]*STATES[6]*CONSTANTS[31];
ALGEBRAIC[26] = CONSTANTS[48]*STATES[15];
ALGEBRAIC[34] = CONSTANTS[49]*STATES[12]*STATES[17];
ALGEBRAIC[36] = CONSTANTS[50]*STATES[15];
ALGEBRAIC[17] = CONSTANTS[42]+ CONSTANTS[43]*( CONSTANTS[56]*STATES[1]+ CONSTANTS[57]*STATES[2]+ CONSTANTS[58]*STATES[3]);
ALGEBRAIC[39] = ALGEBRAIC[17]*STATES[15];
ALGEBRAIC[5] = CONSTANTS[0]*CONSTANTS[21]*STATES[4];
ALGEBRAIC[28] = CONSTANTS[48]*STATES[17];
ALGEBRAIC[37] = STATES[11]*CONSTANTS[32];
ALGEBRAIC[40] = STATES[6]*STATES[9]*CONSTANTS[31];
ALGEBRAIC[42] = STATES[11]*CONSTANTS[28];
ALGEBRAIC[44] = STATES[10]*(CONSTANTS[29]+ CONSTANTS[30]*(STATES[1]+STATES[2]));
ALGEBRAIC[32] = CONSTANTS[51]*STATES[12]*STATES[1];
ALGEBRAIC[41] = CONSTANTS[54]*CONSTANTS[15]*STATES[14];
ALGEBRAIC[45] = CONSTANTS[52]*STATES[14];
ALGEBRAIC[47] = ALGEBRAIC[17]*STATES[14];
ALGEBRAIC[38] = CONSTANTS[54]*CONSTANTS[15]*STATES[1];
ALGEBRAIC[50] = CONSTANTS[0]*CONSTANTS[53]*STATES[8]*STATES[13];
ALGEBRAIC[30] = CONSTANTS[51]*STATES[12]*STATES[3];
ALGEBRAIC[43] = CONSTANTS[52]*STATES[16];
ALGEBRAIC[24] = ( CONSTANTS[47]*( CONSTANTS[59]*(STATES[1]+STATES[3])+ CONSTANTS[60]*STATES[2]))/(CONSTANTS[61]+CONSTANTS[55])+CONSTANTS[46];
ALGEBRAIC[46] = ALGEBRAIC[24]*STATES[16];
ALGEBRAIC[48] = ALGEBRAIC[17]*STATES[12];
ALGEBRAIC[51] = ALGEBRAIC[17]*STATES[16];
ALGEBRAIC[49] = ALGEBRAIC[24]*STATES[3];
ALGEBRAIC[52] = CONSTANTS[0]*(CONSTANTS[44]+ CONSTANTS[45]*STATES[8]);
ALGEBRAIC[10] = ( CONSTANTS[69]*STATES[19]*(STATES[21] - CONSTANTS[15]))/((CONSTANTS[75] - CONSTANTS[15])+STATES[21]);
}