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 41 entries in the algebraic variable array.
There are a total of 27 entries in each of the rate and state variable arrays.
There are a total of 65 entries in the constant variable array.
*/
/*
* VOI is time in component Environment (second).
* CONSTANTS[0] is RT in component fixed_parameters (kilojoule_per_mole).
* CONSTANTS[1] is F in component fixed_parameters (kilojoule_per_mole_per_millivolt).
* CONSTANTS[2] is dG_C1o in component fixed_parameters (kilojoule_per_mole).
* CONSTANTS[3] is dG_C3o in component fixed_parameters (kilojoule_per_mole).
* CONSTANTS[4] is dG_C4o in component fixed_parameters (kilojoule_per_mole).
* ALGEBRAIC[3] is dG_C1op in component fixed_parameters (kilojoule_per_mole).
* ALGEBRAIC[6] is dG_C3op in component fixed_parameters (kilojoule_per_mole).
* ALGEBRAIC[7] is dG_C4op in component fixed_parameters (kilojoule_per_mole).
* CONSTANTS[5] is dG_F1o in component fixed_parameters (kilojoule_per_mole).
* ALGEBRAIC[0] is dG_H in component fixed_parameters (kilojoule_per_mole).
* CONSTANTS[6] is k_dHPI in component fixed_parameters (molar).
* CONSTANTS[7] is k_dHatp in component fixed_parameters (molar).
* CONSTANTS[8] is k_dHadp in component fixed_parameters (molar).
* CONSTANTS[9] is K_DT in component fixed_parameters (molar).
* CONSTANTS[10] is K_DD in component fixed_parameters (molar).
* CONSTANTS[11] is K_AK in component fixed_parameters (dimensionless).
* CONSTANTS[12] is K_CK in component fixed_parameters (per_molar).
* CONSTANTS[13] is W_c in component fixed_parameters (l_water_per_l_cyto).
* CONSTANTS[14] is W_x in component fixed_parameters (l_water_per_l_mito).
* CONSTANTS[15] is W_i in component fixed_parameters (l_water_per_l_mito).
* CONSTANTS[16] is V_cyto in component fixed_parameters (cyto_per_cell).
* CONSTANTS[17] is V_mito in component fixed_parameters (mito_per_cell).
* CONSTANTS[61] is Rm_cyto in component fixed_parameters (mito_per_cyto).
* CONSTANTS[58] is Rm_cell in component fixed_parameters (mito_per_cell).
* CONSTANTS[63] is Rc_cell in component fixed_parameters (cyto_per_cell).
* CONSTANTS[18] is n_A in component fixed_parameters (dimensionless).
* CONSTANTS[19] is C_tot in component fixed_parameters (molar).
* CONSTANTS[20] is CR_tot in component fixed_parameters (molar).
* CONSTANTS[21] is Q_tot in component fixed_parameters (molar).
* CONSTANTS[22] is NAD_tot in component fixed_parameters (molar).
* CONSTANTS[23] is pH_C in component fixed_parameters (dimensionless).
* CONSTANTS[59] is H_c in component fixed_parameters (molar).
* CONSTANTS[24] is K_c in component fixed_parameters (molar).
* CONSTANTS[25] is ATP_c in component fixed_parameters (molar).
* CONSTANTS[26] is AMP_c in component fixed_parameters (molar).
* CONSTANTS[62] is H_i in component fixed_parameters (molar).
* CONSTANTS[60] is K_i in component fixed_parameters (molar).
* STATES[0] is H_x in component dH_x_dt (molar).
* STATES[1] is deltaPsi in component ddeltaPsi_dt (millivolt).
* CONSTANTS[27] is k_PI1 in component adjustable_parameters (molar).
* CONSTANTS[28] is k_PI2 in component adjustable_parameters (molar).
* CONSTANTS[29] is k_PI3 in component adjustable_parameters (molar).
* CONSTANTS[30] is k_PI4 in component adjustable_parameters (molar).
* CONSTANTS[31] is k_PIH in component adjustable_parameters (molar).
* CONSTANTS[32] is r in component adjustable_parameters (dimensionless).
* CONSTANTS[33] is X_DH in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
* CONSTANTS[34] is X_C1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[35] is X_C3 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_half_molar).
* CONSTANTS[36] is X_C4 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
* CONSTANTS[37] is X_F1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[38] is X_ANT in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
* CONSTANTS[39] is X_PI1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
* CONSTANTS[40] is X_KH in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[41] is X_Hle in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_millivolt).
* CONSTANTS[42] is X_K in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_millivolt).
* CONSTANTS[43] is k_mADP in component adjustable_parameters (molar).
* CONSTANTS[44] is X_AK in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[45] is X_CK in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[46] is X_MgA in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[47] is k_O2 in component adjustable_parameters (molar).
* CONSTANTS[48] is beta in component adjustable_parameters (molar).
* CONSTANTS[49] is C_IM in component adjustable_parameters (mole_per_l_mito_per_millivolt).
* CONSTANTS[50] is X_A in component adjustable_parameters (micron_per_second).
* CONSTANTS[51] is X_PI2 in component adjustable_parameters (micron_per_second).
* CONSTANTS[52] is gamma in component adjustable_parameters (per_micron).
* ALGEBRAIC[8] is NAD_x in component dNAD_x_dt (molar).
* STATES[2] is NADH_x in component dNADH_x_dt (molar).
* ALGEBRAIC[9] is QH2_x in component dQH2_x_dt (molar).
* STATES[3] is Q_x in component dQ_x_dt (molar).
* ALGEBRAIC[10] is Cred_x in component dCred_x_dt (molar).
* STATES[4] is Cox_x in component dCox_x_dt (molar).
* ALGEBRAIC[11] is ATP_fx in component dATP_fx_dt (molar).
* STATES[5] is ATP_x in component dATP_x_dt (molar).
* STATES[6] is ATP_mx in component dATP_mx_dt (molar).
* ALGEBRAIC[12] is ADP_fx in component dADP_fx_dt (molar).
* STATES[7] is ADP_x in component dADP_x_dt (molar).
* STATES[8] is ADP_mx in component dADP_mx_dt (molar).
* ALGEBRAIC[19] is J_DH in component J_DH (mole_per_second_per_l_mito).
* ALGEBRAIC[21] is J_C1 in component J_C1 (mole_per_second_per_l_mito).
* ALGEBRAIC[22] is J_C3 in component J_C3 (mole_per_second_per_l_mito).
* ALGEBRAIC[23] is J_C4 in component J_C4 (mole_per_second_per_l_mito).
* ALGEBRAIC[24] is J_F1 in component J_F1 (mole_per_second_per_l_mito).
* ALGEBRAIC[32] is J_PI1 in component J_PI1 (mole_per_second_per_l_mito).
* ALGEBRAIC[38] is J_KH in component J_KH (mole_per_second_per_l_mito).
* ALGEBRAIC[34] is J_Hle in component J_Hle (mole_per_second_per_l_mito).
* STATES[9] is K_x in component dK_x_dt (molar).
* ALGEBRAIC[40] is J_K in component J_K (mole_per_second_per_l_mito).
* STATES[10] is Mg_x in component dMg_x_dt (molar).
* ALGEBRAIC[14] is J_MgATP_x in component J_MgATP_x (mole_per_second_per_l_mito).
* ALGEBRAIC[16] is J_MgADP_x in component J_MgADP_x (mole_per_second_per_l_mito).
* ALGEBRAIC[27] is J_ANT in component J_ANT (mole_per_second_per_l_mito).
* STATES[11] is PI_x in component dPI_x_dt (molar).
* ALGEBRAIC[13] is ATP_fi in component dATP_fi_dt (molar).
* STATES[12] is ATP_i in component dATP_i_dt (molar).
* STATES[13] is ATP_mi in component dATP_mi_dt (molar).
* ALGEBRAIC[15] is ADP_fi in component dADP_fi_dt (molar).
* STATES[14] is ADP_i in component dADP_i_dt (molar).
* STATES[15] is ADP_mi in component dADP_mi_dt (molar).
* ALGEBRAIC[33] is J_ATP in component J_ATP (mole_per_second_per_l_mito).
* ALGEBRAIC[35] is J_AKi in component J_AKi (mole_per_second_per_l_mito).
* ALGEBRAIC[31] is J_ADP in component J_ADP (mole_per_second_per_l_mito).
* STATES[16] is AMP_i in component dAMP_i_dt (molar).
* ALGEBRAIC[29] is J_AMP in component J_AMP (mole_per_second_per_l_mito).
* ALGEBRAIC[18] is J_MgATP_i in component J_MgATP_i (mole_per_second_per_l_mito).
* ALGEBRAIC[20] is J_MgADP_i in component J_MgADP_i (mole_per_second_per_l_mito).
* STATES[17] is PI_i in component dPI_i_dt (molar).
* ALGEBRAIC[36] is J_PI2 in component J_PI2 (mole_per_second_per_l_mito).
* STATES[18] is Mg_i in component dMg_i_dt (molar).
* STATES[19] is ATP_c in component dATP_c_dt (molar).
* ALGEBRAIC[37] is J_AKc in component J_AKc (mole_per_second_per_l_cyto).
* CONSTANTS[64] is J_AtC in component J_AtC (mole_per_second_per_l_cell).
* ALGEBRAIC[39] is J_CKc in component J_CKc (mole_per_second_per_l_cyto).
* STATES[20] is ADP_c in component dADP_c_dt (molar).
* STATES[21] is AMP_c in component dAMP_c_dt (molar).
* STATES[22] is ATP_mc in component dATP_mc_dt (molar).
* ALGEBRAIC[4] is J_MgATP_c in component J_MgATP_c (mole_per_second_per_l_cyto).
* STATES[23] is ADP_mc in component dADP_mc_dt (molar).
* ALGEBRAIC[5] is J_MgADP_c in component J_MgADP_c (mole_per_second_per_l_cyto).
* STATES[24] is PI_c in component dPI_c_dt (molar).
* STATES[25] is Mg_c in component dMg_c_dt (molar).
* STATES[26] is PCr_c in component dPCr_c_dt (molar).
* ALGEBRAIC[1] is ATP_fc in component dATP_fc_dt (molar).
* ALGEBRAIC[2] is ADP_fc in component dADP_fc_dt (molar).
* ALGEBRAIC[17] is Cr_c in component dCr_c_dt (molar).
* CONSTANTS[53] is O2_x in component J_C4 (molar).
* CONSTANTS[54] is mincon in component J_ANT (molar).
* ALGEBRAIC[25] is Psi_x in component J_ANT (millivolt).
* ALGEBRAIC[26] is Psi_i in component J_ANT (millivolt).
* ALGEBRAIC[28] is H2PIi in component J_PI1 (molar).
* ALGEBRAIC[30] is H2PIx in component J_PI1 (molar).
* CONSTANTS[55] is mincond in component J_Hle (millivolt).
* CONSTANTS[56] is mincond in component J_K (millivolt).
* CONSTANTS[57] is X_AtC in component J_AtC (mole_per_second_per_l_cell).
* RATES[0] is d/dt H_x in component dH_x_dt (molar).
* RATES[9] is d/dt K_x in component dK_x_dt (molar).
* RATES[10] is d/dt Mg_x in component dMg_x_dt (molar).
* RATES[2] is d/dt NADH_x in component dNADH_x_dt (molar).
* RATES[3] is d/dt Q_x in component dQ_x_dt (molar).
* RATES[4] is d/dt Cox_x in component dCox_x_dt (molar).
* RATES[5] is d/dt ATP_x in component dATP_x_dt (molar).
* RATES[7] is d/dt ADP_x in component dADP_x_dt (molar).
* RATES[6] is d/dt ATP_mx in component dATP_mx_dt (molar).
* RATES[8] is d/dt ADP_mx in component dADP_mx_dt (molar).
* RATES[11] is d/dt PI_x in component dPI_x_dt (molar).
* RATES[12] is d/dt ATP_i in component dATP_i_dt (molar).
* RATES[14] is d/dt ADP_i in component dADP_i_dt (molar).
* RATES[16] is d/dt AMP_i in component dAMP_i_dt (molar).
* RATES[13] is d/dt ATP_mi in component dATP_mi_dt (molar).
* RATES[15] is d/dt ADP_mi in component dADP_mi_dt (molar).
* RATES[17] is d/dt PI_i in component dPI_i_dt (molar).
* RATES[18] is d/dt Mg_i in component dMg_i_dt (molar).
* RATES[19] is d/dt ATP_c in component dATP_c_dt (molar).
* RATES[20] is d/dt ADP_c in component dADP_c_dt (molar).
* RATES[21] is d/dt AMP_c in component dAMP_c_dt (molar).
* RATES[22] is d/dt ATP_mc in component dATP_mc_dt (molar).
* RATES[23] is d/dt ADP_mc in component dADP_mc_dt (molar).
* RATES[24] is d/dt PI_c in component dPI_c_dt (molar).
* RATES[25] is d/dt Mg_c in component dMg_c_dt (molar).
* RATES[26] is d/dt PCr_c in component dPCr_c_dt (molar).
* RATES[1] is d/dt deltaPsi in component ddeltaPsi_dt (millivolt).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 2.5775;
CONSTANTS[1] = 0.096484;
CONSTANTS[2] = -69.37;
CONSTANTS[3] = -32.53;
CONSTANTS[4] = -122.94;
CONSTANTS[5] = 36.03;
CONSTANTS[6] = 1.7782754100389227e-7;
CONSTANTS[7] = 3.311311214825911e-7;
CONSTANTS[8] = 5.128613839913649e-7;
CONSTANTS[9] = 24e-6;
CONSTANTS[10] = 347e-6;
CONSTANTS[11] = 0.4331;
CONSTANTS[12] = 1.66e9;
CONSTANTS[13] = 0.8425;
CONSTANTS[14] = 0.6514;
CONSTANTS[15] = 0.0724;
CONSTANTS[16] = 0.894;
CONSTANTS[17] = 0.056;
CONSTANTS[18] = 3;
CONSTANTS[19] = 2.7e-3;
CONSTANTS[20] = 42.7e-3;
CONSTANTS[21] = 1.35e-3;
CONSTANTS[22] = 2.97e-3;
CONSTANTS[23] = 7;
CONSTANTS[24] = 0.15;
CONSTANTS[25] = 0;
CONSTANTS[26] = 0;
STATES[0] = 1e-7;
STATES[1] = 166.67;
CONSTANTS[27] = 0.1553e-3;
CONSTANTS[28] = 0.8222e-3;
CONSTANTS[29] = 0.3601e-3;
CONSTANTS[30] = 5.924e-3;
CONSTANTS[31] = 2.542e-4;
CONSTANTS[32] = 4.599;
CONSTANTS[33] = 0.0866;
CONSTANTS[34] = 4.405e3;
CONSTANTS[35] = 4.887;
CONSTANTS[36] = 6.766e-5;
CONSTANTS[37] = 1e3;
CONSTANTS[38] = 8.123e-3;
CONSTANTS[39] = 3.85e5;
CONSTANTS[40] = 5.651e7;
CONSTANTS[41] = 200;
CONSTANTS[42] = 0;
CONSTANTS[43] = 3.5e-6;
CONSTANTS[44] = 1e7;
CONSTANTS[45] = 1e7;
CONSTANTS[46] = 1e7;
CONSTANTS[47] = 1.2e-4;
CONSTANTS[48] = 0.01;
CONSTANTS[49] = 6.75e-6;
CONSTANTS[50] = 85;
CONSTANTS[51] = 327;
CONSTANTS[52] = 5.99;
STATES[2] = 0.0015723;
STATES[3] = 6.75e-4;
STATES[4] = 1.35e-3;
STATES[5] = 0.0026657;
STATES[6] = 0.0026046;
STATES[7] = 0.0073343;
STATES[8] = 0.0054765;
STATES[9] = 0.14661;
STATES[10] = 1.0229e-3;
STATES[11] = 1.72e-4;
STATES[12] = 0.0065339;
STATES[13] = 0.0063812;
STATES[14] = 6.5773e-5;
STATES[15] = 4.8866e-5;
STATES[16] = 2.8837e-7;
STATES[17] = 1.72e-4;
STATES[18] = 0.0010029;
STATES[19] = 0.0081312;
STATES[20] = 6.85e-5;
STATES[21] = 3.0911e-7;
STATES[22] = 0.0079786;
STATES[23] = 5.1958e-5;
STATES[24] = 8.7702e-3;
STATES[25] = 0.001003;
STATES[26] = 23.589e-3;
CONSTANTS[53] = 0;
CONSTANTS[54] = 1e-9;
CONSTANTS[55] = 1e-6;
CONSTANTS[56] = 1e-6;
CONSTANTS[57] = 0;
CONSTANTS[58] = CONSTANTS[17];
CONSTANTS[59] = 1.00000*pow(10.0000, - CONSTANTS[23]);
CONSTANTS[60] = CONSTANTS[24];
CONSTANTS[61] = CONSTANTS[17]/CONSTANTS[16];
CONSTANTS[62] = CONSTANTS[59];
CONSTANTS[63] = CONSTANTS[16];
CONSTANTS[64] = CONSTANTS[57];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = STATES[19] - STATES[22];
ALGEBRAIC[4] = CONSTANTS[46]*( ALGEBRAIC[1]*STATES[25] - CONSTANTS[9]*STATES[22]);
RATES[22] = ALGEBRAIC[4]/CONSTANTS[13];
ALGEBRAIC[2] = STATES[20] - STATES[23];
ALGEBRAIC[5] = CONSTANTS[46]*( ALGEBRAIC[2]*STATES[25] - CONSTANTS[10]*STATES[23]);
RATES[23] = ALGEBRAIC[5]/CONSTANTS[13];
RATES[25] = (- ALGEBRAIC[4] - ALGEBRAIC[4])/CONSTANTS[13];
ALGEBRAIC[11] = STATES[5] - STATES[6];
ALGEBRAIC[14] = CONSTANTS[46]*( ALGEBRAIC[11]*STATES[10] - CONSTANTS[9]*STATES[6]);
RATES[6] = ALGEBRAIC[14]/CONSTANTS[14];
ALGEBRAIC[12] = STATES[7] - STATES[8];
ALGEBRAIC[16] = CONSTANTS[46]*( ALGEBRAIC[12]*STATES[10] - CONSTANTS[10]*STATES[8]);
RATES[10] = (- ALGEBRAIC[16] - ALGEBRAIC[14])/CONSTANTS[14];
RATES[8] = ALGEBRAIC[16]/CONSTANTS[14];
ALGEBRAIC[13] = STATES[12] - STATES[13];
ALGEBRAIC[18] = CONSTANTS[46]*( ALGEBRAIC[13]*STATES[18] - CONSTANTS[9]*STATES[13]);
RATES[13] = ALGEBRAIC[18]/CONSTANTS[15];
ALGEBRAIC[15] = STATES[14] - STATES[15];
ALGEBRAIC[20] = CONSTANTS[46]*( ALGEBRAIC[15]*STATES[18] - CONSTANTS[10]*STATES[15]);
RATES[15] = ALGEBRAIC[20]/CONSTANTS[15];
RATES[18] = (- ALGEBRAIC[20] - ALGEBRAIC[18])/CONSTANTS[15];
ALGEBRAIC[8] = CONSTANTS[22] - STATES[2];
ALGEBRAIC[19] = ( CONSTANTS[33]*( CONSTANTS[32]*ALGEBRAIC[8] - STATES[2])*(1.00000+STATES[11]/CONSTANTS[27]))/(1.00000+STATES[11]/CONSTANTS[28]);
ALGEBRAIC[3] = CONSTANTS[2] - 1.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[0] = CONSTANTS[1]*STATES[1]+ CONSTANTS[0]*arbitrary_log(CONSTANTS[62]/STATES[0], 10);
ALGEBRAIC[9] = CONSTANTS[21] - STATES[3];
ALGEBRAIC[21] = CONSTANTS[34]*( exp(- (ALGEBRAIC[3]+ 4.00000*ALGEBRAIC[0])/CONSTANTS[0])*STATES[2]*STATES[3] - ALGEBRAIC[8]*ALGEBRAIC[9]);
RATES[2] = (ALGEBRAIC[19] - ALGEBRAIC[21])/CONSTANTS[14];
ALGEBRAIC[6] = CONSTANTS[3]+ 2.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[10] = CONSTANTS[19] - STATES[4];
ALGEBRAIC[22] = (( CONSTANTS[35]*(1.00000+STATES[11]/CONSTANTS[29]))/(1.00000+STATES[11]/CONSTANTS[30]))*( exp(- ((ALGEBRAIC[6]+ 4.00000*ALGEBRAIC[0]) - 2.00000*CONSTANTS[1]*STATES[1])/( 2.00000*CONSTANTS[0]))*STATES[4]* pow(ALGEBRAIC[9], 1.0 / 2) - ALGEBRAIC[10]* pow(STATES[3], 1.0 / 2));
RATES[3] = (ALGEBRAIC[21] - ALGEBRAIC[22])/CONSTANTS[14];
ALGEBRAIC[7] = CONSTANTS[4] - 2.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[23] = (CONSTANTS[53]<1.00000e-12 ? 0.00000 : (( (( CONSTANTS[36]*CONSTANTS[53])/(CONSTANTS[53]+CONSTANTS[47]))*ALGEBRAIC[10])/CONSTANTS[19])*( exp(- (ALGEBRAIC[7]+ 2.00000*ALGEBRAIC[0])/( 2.00000*CONSTANTS[0]))*ALGEBRAIC[10]*pow(CONSTANTS[53], 0.250000) - STATES[4]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0])));
RATES[4] = ( 2.00000*ALGEBRAIC[22] - 2.00000*ALGEBRAIC[23])/CONSTANTS[14];
ALGEBRAIC[24] = CONSTANTS[37]*( (( exp(- (CONSTANTS[5] - CONSTANTS[18]*ALGEBRAIC[0])/CONSTANTS[0])*CONSTANTS[10])/CONSTANTS[9])*STATES[8]*STATES[11] - STATES[6]*1.00000);
ALGEBRAIC[25] = - 0.650000*STATES[1];
ALGEBRAIC[26] = 0.350000*STATES[1];
ALGEBRAIC[27] = (ALGEBRAIC[15]>CONSTANTS[54]&&ALGEBRAIC[13]>CONSTANTS[54] ? ( CONSTANTS[38]*(ALGEBRAIC[15]/(ALGEBRAIC[15]+ ALGEBRAIC[13]*exp(( - CONSTANTS[1]*ALGEBRAIC[26])/CONSTANTS[0])) - ALGEBRAIC[12]/(ALGEBRAIC[12]+ ALGEBRAIC[11]*exp(( - CONSTANTS[1]*ALGEBRAIC[25])/CONSTANTS[0])))*ALGEBRAIC[15])/(ALGEBRAIC[15]+CONSTANTS[43]) : 0.00000);
RATES[5] = (ALGEBRAIC[24] - ALGEBRAIC[27])/CONSTANTS[14];
RATES[7] = (- ALGEBRAIC[24]+ALGEBRAIC[27])/CONSTANTS[14];
ALGEBRAIC[28] = ( STATES[17]*CONSTANTS[62])/(CONSTANTS[62]+CONSTANTS[6]);
ALGEBRAIC[30] = ( STATES[11]*STATES[0])/(STATES[0]+CONSTANTS[6]);
ALGEBRAIC[32] = ( CONSTANTS[39]*( CONSTANTS[62]*ALGEBRAIC[28] - STATES[0]*ALGEBRAIC[30]))/(ALGEBRAIC[28]+CONSTANTS[31]);
RATES[11] = (- ALGEBRAIC[24] - ALGEBRAIC[32])/CONSTANTS[14];
ALGEBRAIC[33] = CONSTANTS[52]*CONSTANTS[50]*(STATES[19] - STATES[12]);
ALGEBRAIC[35] = CONSTANTS[44]*( CONSTANTS[11]*STATES[14]*STATES[14] - STATES[16]*STATES[12]);
RATES[12] = (ALGEBRAIC[33]+ALGEBRAIC[27]+ALGEBRAIC[35])/CONSTANTS[15];
ALGEBRAIC[31] = CONSTANTS[52]*CONSTANTS[50]*(STATES[20] - STATES[14]);
RATES[14] = ((ALGEBRAIC[31] - ALGEBRAIC[27]) - 2.00000*ALGEBRAIC[35])/CONSTANTS[15];
ALGEBRAIC[29] = CONSTANTS[52]*CONSTANTS[50]*(STATES[21] - STATES[16]);
RATES[16] = (ALGEBRAIC[29]+ALGEBRAIC[35])/CONSTANTS[15];
ALGEBRAIC[36] = CONSTANTS[52]*CONSTANTS[51]*(STATES[24] - STATES[17]);
RATES[17] = (- ALGEBRAIC[32]+ALGEBRAIC[36])/CONSTANTS[15];
ALGEBRAIC[37] = CONSTANTS[44]*( CONSTANTS[11]*STATES[20]*STATES[20] - STATES[21]*STATES[19]);
RATES[21] = ( - CONSTANTS[61]*ALGEBRAIC[29]+ALGEBRAIC[37])/CONSTANTS[13];
RATES[24] = ( - ALGEBRAIC[36]*CONSTANTS[61]+CONSTANTS[64]/CONSTANTS[63])/CONSTANTS[13];
ALGEBRAIC[38] = CONSTANTS[40]*( CONSTANTS[60]*STATES[0] - STATES[9]*CONSTANTS[62]);
ALGEBRAIC[34] = (fabs(STATES[1])>CONSTANTS[55] ? ( CONSTANTS[41]*STATES[1]*( CONSTANTS[62]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - STATES[0]))/(exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - 1.00000) : ( CONSTANTS[41]*CONSTANTS[0]*(CONSTANTS[62] - STATES[0]))/CONSTANTS[1]);
RATES[0] = ( (STATES[0]/CONSTANTS[48])*((((((ALGEBRAIC[19] - 5.00000*ALGEBRAIC[21]) - 2.00000*ALGEBRAIC[22]) - 4.00000*ALGEBRAIC[23])+ (CONSTANTS[18] - 1.00000)*ALGEBRAIC[24]+ 2.00000*ALGEBRAIC[32]) - ALGEBRAIC[38])+ALGEBRAIC[34]))/CONSTANTS[14];
ALGEBRAIC[17] = CONSTANTS[20] - STATES[26];
ALGEBRAIC[39] = CONSTANTS[45]*( CONSTANTS[12]*STATES[20]*STATES[26]*CONSTANTS[59] - STATES[19]*ALGEBRAIC[17]);
RATES[19] = ((( - CONSTANTS[61]*ALGEBRAIC[33]+ALGEBRAIC[37]) - CONSTANTS[64]/CONSTANTS[63])+ALGEBRAIC[39])/CONSTANTS[13];
RATES[20] = ((( - CONSTANTS[61]*ALGEBRAIC[31] - 2.00000*ALGEBRAIC[37])+CONSTANTS[64]/CONSTANTS[63]) - ALGEBRAIC[39])/CONSTANTS[13];
RATES[26] = - ALGEBRAIC[39]/CONSTANTS[13];
ALGEBRAIC[40] = (fabs(STATES[1])>CONSTANTS[56] ? ( CONSTANTS[42]*STATES[1]*( CONSTANTS[60]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - STATES[9]))/(exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - 1.00000) : ( CONSTANTS[42]*CONSTANTS[0]*(CONSTANTS[60] - STATES[9]))/CONSTANTS[1]);
RATES[9] = (ALGEBRAIC[38]+ALGEBRAIC[40])/CONSTANTS[14];
RATES[1] = ((((( 4.00000*ALGEBRAIC[21]+ 2.00000*ALGEBRAIC[22]+ 4.00000*ALGEBRAIC[23]) - CONSTANTS[18]*ALGEBRAIC[24]) - ALGEBRAIC[27]) - ALGEBRAIC[34]) - ALGEBRAIC[40])/CONSTANTS[49];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = STATES[19] - STATES[22];
ALGEBRAIC[4] = CONSTANTS[46]*( ALGEBRAIC[1]*STATES[25] - CONSTANTS[9]*STATES[22]);
ALGEBRAIC[2] = STATES[20] - STATES[23];
ALGEBRAIC[5] = CONSTANTS[46]*( ALGEBRAIC[2]*STATES[25] - CONSTANTS[10]*STATES[23]);
ALGEBRAIC[11] = STATES[5] - STATES[6];
ALGEBRAIC[14] = CONSTANTS[46]*( ALGEBRAIC[11]*STATES[10] - CONSTANTS[9]*STATES[6]);
ALGEBRAIC[12] = STATES[7] - STATES[8];
ALGEBRAIC[16] = CONSTANTS[46]*( ALGEBRAIC[12]*STATES[10] - CONSTANTS[10]*STATES[8]);
ALGEBRAIC[13] = STATES[12] - STATES[13];
ALGEBRAIC[18] = CONSTANTS[46]*( ALGEBRAIC[13]*STATES[18] - CONSTANTS[9]*STATES[13]);
ALGEBRAIC[15] = STATES[14] - STATES[15];
ALGEBRAIC[20] = CONSTANTS[46]*( ALGEBRAIC[15]*STATES[18] - CONSTANTS[10]*STATES[15]);
ALGEBRAIC[8] = CONSTANTS[22] - STATES[2];
ALGEBRAIC[19] = ( CONSTANTS[33]*( CONSTANTS[32]*ALGEBRAIC[8] - STATES[2])*(1.00000+STATES[11]/CONSTANTS[27]))/(1.00000+STATES[11]/CONSTANTS[28]);
ALGEBRAIC[3] = CONSTANTS[2] - 1.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[0] = CONSTANTS[1]*STATES[1]+ CONSTANTS[0]*arbitrary_log(CONSTANTS[62]/STATES[0], 10);
ALGEBRAIC[9] = CONSTANTS[21] - STATES[3];
ALGEBRAIC[21] = CONSTANTS[34]*( exp(- (ALGEBRAIC[3]+ 4.00000*ALGEBRAIC[0])/CONSTANTS[0])*STATES[2]*STATES[3] - ALGEBRAIC[8]*ALGEBRAIC[9]);
ALGEBRAIC[6] = CONSTANTS[3]+ 2.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[10] = CONSTANTS[19] - STATES[4];
ALGEBRAIC[22] = (( CONSTANTS[35]*(1.00000+STATES[11]/CONSTANTS[29]))/(1.00000+STATES[11]/CONSTANTS[30]))*( exp(- ((ALGEBRAIC[6]+ 4.00000*ALGEBRAIC[0]) - 2.00000*CONSTANTS[1]*STATES[1])/( 2.00000*CONSTANTS[0]))*STATES[4]* pow(ALGEBRAIC[9], 1.0 / 2) - ALGEBRAIC[10]* pow(STATES[3], 1.0 / 2));
ALGEBRAIC[7] = CONSTANTS[4] - 2.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[23] = (CONSTANTS[53]<1.00000e-12 ? 0.00000 : (( (( CONSTANTS[36]*CONSTANTS[53])/(CONSTANTS[53]+CONSTANTS[47]))*ALGEBRAIC[10])/CONSTANTS[19])*( exp(- (ALGEBRAIC[7]+ 2.00000*ALGEBRAIC[0])/( 2.00000*CONSTANTS[0]))*ALGEBRAIC[10]*pow(CONSTANTS[53], 0.250000) - STATES[4]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0])));
ALGEBRAIC[24] = CONSTANTS[37]*( (( exp(- (CONSTANTS[5] - CONSTANTS[18]*ALGEBRAIC[0])/CONSTANTS[0])*CONSTANTS[10])/CONSTANTS[9])*STATES[8]*STATES[11] - STATES[6]*1.00000);
ALGEBRAIC[25] = - 0.650000*STATES[1];
ALGEBRAIC[26] = 0.350000*STATES[1];
ALGEBRAIC[27] = (ALGEBRAIC[15]>CONSTANTS[54]&&ALGEBRAIC[13]>CONSTANTS[54] ? ( CONSTANTS[38]*(ALGEBRAIC[15]/(ALGEBRAIC[15]+ ALGEBRAIC[13]*exp(( - CONSTANTS[1]*ALGEBRAIC[26])/CONSTANTS[0])) - ALGEBRAIC[12]/(ALGEBRAIC[12]+ ALGEBRAIC[11]*exp(( - CONSTANTS[1]*ALGEBRAIC[25])/CONSTANTS[0])))*ALGEBRAIC[15])/(ALGEBRAIC[15]+CONSTANTS[43]) : 0.00000);
ALGEBRAIC[28] = ( STATES[17]*CONSTANTS[62])/(CONSTANTS[62]+CONSTANTS[6]);
ALGEBRAIC[30] = ( STATES[11]*STATES[0])/(STATES[0]+CONSTANTS[6]);
ALGEBRAIC[32] = ( CONSTANTS[39]*( CONSTANTS[62]*ALGEBRAIC[28] - STATES[0]*ALGEBRAIC[30]))/(ALGEBRAIC[28]+CONSTANTS[31]);
ALGEBRAIC[33] = CONSTANTS[52]*CONSTANTS[50]*(STATES[19] - STATES[12]);
ALGEBRAIC[35] = CONSTANTS[44]*( CONSTANTS[11]*STATES[14]*STATES[14] - STATES[16]*STATES[12]);
ALGEBRAIC[31] = CONSTANTS[52]*CONSTANTS[50]*(STATES[20] - STATES[14]);
ALGEBRAIC[29] = CONSTANTS[52]*CONSTANTS[50]*(STATES[21] - STATES[16]);
ALGEBRAIC[36] = CONSTANTS[52]*CONSTANTS[51]*(STATES[24] - STATES[17]);
ALGEBRAIC[37] = CONSTANTS[44]*( CONSTANTS[11]*STATES[20]*STATES[20] - STATES[21]*STATES[19]);
ALGEBRAIC[38] = CONSTANTS[40]*( CONSTANTS[60]*STATES[0] - STATES[9]*CONSTANTS[62]);
ALGEBRAIC[34] = (fabs(STATES[1])>CONSTANTS[55] ? ( CONSTANTS[41]*STATES[1]*( CONSTANTS[62]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - STATES[0]))/(exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - 1.00000) : ( CONSTANTS[41]*CONSTANTS[0]*(CONSTANTS[62] - STATES[0]))/CONSTANTS[1]);
ALGEBRAIC[17] = CONSTANTS[20] - STATES[26];
ALGEBRAIC[39] = CONSTANTS[45]*( CONSTANTS[12]*STATES[20]*STATES[26]*CONSTANTS[59] - STATES[19]*ALGEBRAIC[17]);
ALGEBRAIC[40] = (fabs(STATES[1])>CONSTANTS[56] ? ( CONSTANTS[42]*STATES[1]*( CONSTANTS[60]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - STATES[9]))/(exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - 1.00000) : ( CONSTANTS[42]*CONSTANTS[0]*(CONSTANTS[60] - STATES[9]))/CONSTANTS[1]);
}