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 34 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 54 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (second).
 * CONSTANTS[0] is isotonic in component parameters (dimensionless).
 * CONSTANTS[1] is alpha_1 in component parameters (per_micrometre).
 * CONSTANTS[2] is beta_1 in component parameters (millinewton).
 * CONSTANTS[3] is alpha_2 in component parameters (per_micrometre).
 * CONSTANTS[4] is beta_2 in component parameters (millinewton).
 * CONSTANTS[5] is alpha_3 in component parameters (per_micrometre).
 * CONSTANTS[6] is beta_3 in component parameters (millinewton).
 * CONSTANTS[7] is lambda in component parameters (millinewton).
 * CONSTANTS[8] is A_half in component parameters (micromole).
 * CONSTANTS[9] is mu in component parameters (dimensionless).
 * CONSTANTS[10] is chi in component parameters (dimensionless).
 * CONSTANTS[11] is chi_0 in component parameters (dimensionless).
 * CONSTANTS[12] is m_0 in component parameters (dimensionless).
 * CONSTANTS[13] is v_max in component parameters (micrometre_per_second).
 * CONSTANTS[14] is a in component parameters (dimensionless).
 * CONSTANTS[15] is d_h in component parameters (dimensionless).
 * CONSTANTS[16] is alpha_P in component parameters (dimensionless).
 * ALGEBRAIC[9] is l in component length (micrometre).
 * ALGEBRAIC[5] is F_muscle in component force (millinewton).
 * ALGEBRAIC[11] is flag in component isotonic (dimensionless).
 * CONSTANTS[17] is F_afterload in component isotonic (millinewton).
 * ALGEBRAIC[13] is isotonic_mode in component isotonic (dimensionless).
 * CONSTANTS[18] is l_0 in component isotonic (micrometre).
 * CONSTANTS[19] is S_0 in component parameters_izakov_et_al_1991 (micrometre).
 * ALGEBRAIC[0] is q_v in component parameters_izakov_et_al_1991 (per_second).
 * CONSTANTS[20] is q_1 in component parameters_izakov_et_al_1991 (per_second).
 * CONSTANTS[21] is q_2 in component parameters_izakov_et_al_1991 (per_second).
 * CONSTANTS[22] is q_3 in component parameters_izakov_et_al_1991 (per_second).
 * CONSTANTS[23] is q_4 in component parameters_izakov_et_al_1991 (per_second).
 * CONSTANTS[24] is v_star in component parameters_izakov_et_al_1991 (micrometre_per_second).
 * CONSTANTS[50] is v_1 in component parameters_izakov_et_al_1991 (micrometre_per_second).
 * CONSTANTS[25] is alpha_G in component parameters_izakov_et_al_1991 (dimensionless).
 * CONSTANTS[26] is a_on in component parameters_izakov_et_al_1991 (per_second_micromole).
 * CONSTANTS[27] is a_off in component parameters_izakov_et_al_1991 (per_second).
 * CONSTANTS[28] is k_A in component parameters_izakov_et_al_1991 (per_micromole).
 * STATES[0] is v in component CE_velocity (micrometre_per_second).
 * CONSTANTS[29] is alpha_Q in component parameters_izakov_et_al_1991 (dimensionless).
 * CONSTANTS[30] is beta_Q in component parameters_izakov_et_al_1991 (dimensionless).
 * ALGEBRAIC[31] is F_CE in component force (millinewton).
 * ALGEBRAIC[4] is F_XSE in component force (millinewton).
 * ALGEBRAIC[1] is F_SE in component force (millinewton).
 * ALGEBRAIC[2] is F_PE in component force (millinewton).
 * STATES[1] is N in component crossbridge_kinetics (dimensionless).
 * ALGEBRAIC[17] is k_P_vis in component CE_velocity (millinewton_second_per_micrometre).
 * ALGEBRAIC[20] is k_S_vis in component PE_velocity (millinewton_second_per_micrometre).
 * STATES[2] is w in component PE_velocity (micrometre_per_second).
 * STATES[3] is l_1 in component length (micrometre).
 * STATES[4] is l_2 in component length (micrometre).
 * STATES[5] is l_3 in component length (micrometre).
 * ALGEBRAIC[30] is p_v in component average_crossbridge_force (dimensionless).
 * ALGEBRAIC[27] is K_chi in component crossbridge_kinetics (per_second).
 * ALGEBRAIC[6] is M_A in component crossbridge_kinetics (dimensionless).
 * ALGEBRAIC[7] is n_1 in component crossbridge_kinetics (dimensionless).
 * ALGEBRAIC[8] is L_oz in component crossbridge_kinetics (dimensionless).
 * ALGEBRAIC[25] is k_p_v in component crossbridge_kinetics (per_second).
 * ALGEBRAIC[26] is k_m_v in component crossbridge_kinetics (per_second).
 * STATES[6] is A in component calcium_handling (micromole).
 * ALGEBRAIC[24] is G_star in component average_crossbridge_force (dimensionless).
 * ALGEBRAIC[21] is P_star in component average_crossbridge_force (dimensionless).
 * CONSTANTS[31] is v_0 in component crossbridge_kinetics (micrometre_per_second).
 * CONSTANTS[32] is q_star in component crossbridge_kinetics (per_second).
 * ALGEBRAIC[3] is dl_1_dt in component length (micrometre_per_second).
 * ALGEBRAIC[22] is dl_2_dt in component length (micrometre_per_second).
 * ALGEBRAIC[23] is dl_3_dt in component length (micrometre_per_second).
 * ALGEBRAIC[18] is phi_chi_2 in component CE_velocity (micrometre_per_second).
 * ALGEBRAIC[33] is phi_chi in component CE_velocity (micrometre_per_second2).
 * ALGEBRAIC[32] is p_prime_v in component average_crossbridge_force (second_per_micrometre).
 * CONSTANTS[33] is alpha_P_lengthening in component CE_velocity (per_micrometre).
 * CONSTANTS[34] is beta_P_lengthening in component CE_velocity (millinewton_second_per_micrometre).
 * CONSTANTS[35] is alpha_P_shortening in component CE_velocity (per_micrometre).
 * CONSTANTS[36] is beta_P_shortening in component CE_velocity (millinewton_second_per_micrometre).
 * ALGEBRAIC[15] is alp_p in component CE_velocity (per_micrometre).
 * CONSTANTS[37] is alpha_S_lengthening in component PE_velocity (per_micrometre).
 * CONSTANTS[38] is beta_S_lengthening in component PE_velocity (millinewton_second_per_micrometre).
 * CONSTANTS[39] is alpha_S_shortening in component PE_velocity (per_micrometre).
 * CONSTANTS[40] is beta_S_shortening in component PE_velocity (millinewton_second_per_micrometre).
 * ALGEBRAIC[19] is alp_s in component PE_velocity (per_micrometre).
 * CONSTANTS[53] is gamma in component average_crossbridge_force (dimensionless).
 * CONSTANTS[51] is case_1 in component average_crossbridge_force (second_per_micrometre).
 * ALGEBRAIC[28] is case_2 in component average_crossbridge_force (second_per_micrometre).
 * CONSTANTS[52] is case_3 in component average_crossbridge_force (second_per_micrometre).
 * ALGEBRAIC[29] is case_4 in component average_crossbridge_force (second_per_micrometre).
 * ALGEBRAIC[14] is dA_dt in component calcium_handling (micromole_per_second).
 * ALGEBRAIC[10] is N_A in component calcium_handling (dimensionless).
 * ALGEBRAIC[12] is pi_N_A in component calcium_handling (dimensionless).
 * STATES[7] is B in component calcium_handling (micromole).
 * ALGEBRAIC[16] is dB_dt in component calcium_handling (micromole_per_second).
 * STATES[8] is Ca_C in component calcium_handling (micromole).
 * CONSTANTS[41] is A_tot in component calcium_handling (micromole).
 * CONSTANTS[42] is B_tot in component calcium_handling (micromole).
 * CONSTANTS[43] is b_on in component calcium_handling (per_second_micromole).
 * CONSTANTS[44] is b_off in component calcium_handling (per_second).
 * CONSTANTS[45] is a_c in component calcium_handling (per_second2).
 * CONSTANTS[46] is r_Ca in component calcium_handling (per_second).
 * CONSTANTS[47] is q_Ca in component calcium_handling (per_micromole).
 * CONSTANTS[48] is t_d in component calcium_handling (second).
 * CONSTANTS[49] is Ca_m in component calcium_handling (micromole).
 * RATES[1] is d/dt N in component crossbridge_kinetics (dimensionless).
 * RATES[3] is d/dt l_1 in component length (micrometre).
 * RATES[4] is d/dt l_2 in component length (micrometre).
 * RATES[5] is d/dt l_3 in component length (micrometre).
 * RATES[0] is d/dt v in component CE_velocity (micrometre_per_second).
 * RATES[2] is d/dt w in component PE_velocity (micrometre_per_second).
 * RATES[6] is d/dt A in component calcium_handling (micromole).
 * RATES[7] is d/dt B in component calcium_handling (micromole).
 * RATES[8] is d/dt Ca_C in component calcium_handling (micromole).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 0;
CONSTANTS[1] = 19;
CONSTANTS[2] = 0.29;
CONSTANTS[3] = 14.6;
CONSTANTS[4] = 0.000924;
CONSTANTS[5] = 48;
CONSTANTS[6] = 0.01;
CONSTANTS[7] = 96;
CONSTANTS[8] = 42;
CONSTANTS[9] = 3;
CONSTANTS[10] = 0.705;
CONSTANTS[11] = 3;
CONSTANTS[12] = 0.9;
CONSTANTS[13] = 5.6;
CONSTANTS[14] = 0.25;
CONSTANTS[15] = 0.5;
CONSTANTS[16] = 4;
CONSTANTS[17] = 2;
CONSTANTS[18] = 0.527;
CONSTANTS[19] = 1.14;
CONSTANTS[20] = 17.3;
CONSTANTS[21] = 259;
CONSTANTS[22] = 17.3;
CONSTANTS[23] = 15;
CONSTANTS[24] = 5.3035675;
CONSTANTS[25] = 1;
CONSTANTS[26] = 32.85;
CONSTANTS[27] = 290;
CONSTANTS[28] = 0.04;
STATES[0] = 0;
CONSTANTS[29] = 10;
CONSTANTS[30] = 5000;
STATES[1] = 0.0001;
STATES[2] = 0;
STATES[3] = 0.437;
STATES[4] = 0.439;
STATES[5] = 0.089;
STATES[6] = 0.7;
CONSTANTS[31] = 560;
CONSTANTS[32] = 1000;
CONSTANTS[33] = 16;
CONSTANTS[34] = 0.0015;
CONSTANTS[35] = 16;
CONSTANTS[36] = 0.0015;
CONSTANTS[37] = 39;
CONSTANTS[38] = 0.008;
CONSTANTS[39] = 46;
CONSTANTS[40] = 0.006;
STATES[7] = 0;
STATES[8] = 0;
CONSTANTS[41] = 70;
CONSTANTS[42] = 28;
CONSTANTS[43] = 37.143;
CONSTANTS[44] = 182;
CONSTANTS[45] = 5200;
CONSTANTS[46] = 650;
CONSTANTS[47] = 0.714;
CONSTANTS[48] = 0.033;
CONSTANTS[49] = 2.1;
CONSTANTS[50] = CONSTANTS[13]/10.0000;
CONSTANTS[51] = ( CONSTANTS[14]*(0.400000+ 0.400000*CONSTANTS[14]))/( CONSTANTS[13]*pow( (CONSTANTS[14]+1.00000)*0.400000, 2.00000));
CONSTANTS[52] = ( 0.400000*CONSTANTS[14]+1.00000)/( CONSTANTS[14]*CONSTANTS[13]);
CONSTANTS[53] = ( CONSTANTS[14]*CONSTANTS[15]*pow(CONSTANTS[50]/CONSTANTS[13], 2.00000))/( 3.00000*CONSTANTS[14]*CONSTANTS[15] - ( (CONSTANTS[14]+1.00000)*CONSTANTS[50])/CONSTANTS[13]);
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[3] = STATES[0];
RATES[3] = ALGEBRAIC[3];
ALGEBRAIC[8] = (STATES[3]+CONSTANTS[19])/(0.460000+CONSTANTS[19]);
ALGEBRAIC[10] = ( CONSTANTS[41]*STATES[1])/( ALGEBRAIC[8]*STATES[6]);
ALGEBRAIC[12] = (ALGEBRAIC[10]>=1.00000 ? 1.00000 : pow(0.0200000, ALGEBRAIC[10]));
ALGEBRAIC[14] =  CONSTANTS[26]*(CONSTANTS[41] - STATES[6])*STATES[8] -  CONSTANTS[27]*exp( - CONSTANTS[28]*STATES[6])*ALGEBRAIC[12]*STATES[6];
RATES[6] = ALGEBRAIC[14];
ALGEBRAIC[16] =  CONSTANTS[43]*(CONSTANTS[42] - STATES[7])*STATES[8] -  CONSTANTS[44]*STATES[7];
RATES[7] = ALGEBRAIC[16];
RATES[8] = (VOI<CONSTANTS[48] ?  4.00000*CONSTANTS[45]*CONSTANTS[49]*VOI*(1.00000 - exp( - CONSTANTS[45]*pow(VOI, 2.00000)))*exp( - CONSTANTS[45]*pow(VOI, 2.00000)) : (- ALGEBRAIC[14] - ALGEBRAIC[16]) -  CONSTANTS[46]*exp( - CONSTANTS[47]*STATES[8])*STATES[8]);
ALGEBRAIC[20] = (STATES[2]<=STATES[0] ?  CONSTANTS[38]*exp( CONSTANTS[37]*(STATES[4] - STATES[3])) :  CONSTANTS[40]*exp( CONSTANTS[39]*(STATES[4] - STATES[3])));
ALGEBRAIC[9] = STATES[4]+STATES[5];
ALGEBRAIC[4] =  CONSTANTS[6]*(exp( CONSTANTS[5]*STATES[5]) - 1.00000);
ALGEBRAIC[5] = ALGEBRAIC[4];
ALGEBRAIC[11] = (ALGEBRAIC[9]>=CONSTANTS[18]&&VOI<0.150000 ? 0.00000 : 1.00000);
ALGEBRAIC[13] = (CONSTANTS[0]==0.00000 ? 0.00000 : CONSTANTS[0]==1.00000&&ALGEBRAIC[5]>=CONSTANTS[17] ? 1.00000 : CONSTANTS[0]==1.00000&&ALGEBRAIC[9]>=CONSTANTS[18]&&ALGEBRAIC[11]==1.00000 ? 0.00000 : 0.0/0.0);
ALGEBRAIC[18] = (ALGEBRAIC[13]==1.00000 ? ( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*STATES[0])/( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))+ CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])) : ( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*STATES[0])/( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))+ CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])+ CONSTANTS[5]*CONSTANTS[6]*exp( CONSTANTS[5]*STATES[5])));
ALGEBRAIC[22] = (ALGEBRAIC[20]==0.00000 ? ALGEBRAIC[18] : STATES[2]);
RATES[4] = ALGEBRAIC[22];
ALGEBRAIC[23] = (ALGEBRAIC[13]==1.00000 ? 0.00000 : ALGEBRAIC[13]==0.00000&&ALGEBRAIC[20]==0.00000 ? - ALGEBRAIC[18] : - STATES[2]);
RATES[5] = ALGEBRAIC[23];
ALGEBRAIC[6] = pow(STATES[6], CONSTANTS[9])/(pow(STATES[6], CONSTANTS[9])+pow(CONSTANTS[8], CONSTANTS[9]));
ALGEBRAIC[7] =  0.600000*STATES[3]+0.500000;
ALGEBRAIC[0] = (STATES[0]<=0.00000 ? CONSTANTS[20] - ( CONSTANTS[21]*STATES[0])/CONSTANTS[13] : STATES[0]<=CONSTANTS[24]&&0.00000<STATES[0] ? ( (CONSTANTS[23] - CONSTANTS[22])*STATES[0])/CONSTANTS[24]+CONSTANTS[22] : CONSTANTS[23]/pow(1.00000+( CONSTANTS[30]*(STATES[0] - CONSTANTS[24]))/CONSTANTS[13], CONSTANTS[29]));
ALGEBRAIC[21] = (STATES[0]<=0.00000 ? ( CONSTANTS[14]*(1.00000+STATES[0]/CONSTANTS[13]))/(CONSTANTS[14] - STATES[0]/CONSTANTS[13]) : (1.00000+CONSTANTS[15]) - ( pow(CONSTANTS[15], 2.00000)*CONSTANTS[14])/( (( CONSTANTS[14]*CONSTANTS[15])/CONSTANTS[53])*pow(STATES[0]/CONSTANTS[13], 2.00000)+( (CONSTANTS[14]+1.00000)*STATES[0])/CONSTANTS[13]+ CONSTANTS[14]*CONSTANTS[15]));
ALGEBRAIC[24] = (- CONSTANTS[13]<=STATES[0]&&STATES[0]<=0.00000 ? 1.00000+( 0.600000*STATES[0])/CONSTANTS[13] : 0.00000<STATES[0]&&STATES[0]<=CONSTANTS[50] ? ALGEBRAIC[21]/(( ( 0.400000*CONSTANTS[14]+1.00000)*STATES[0])/( CONSTANTS[14]*CONSTANTS[13])+1.00000) : ( ALGEBRAIC[21]*exp( - CONSTANTS[25]*pow((STATES[0] - CONSTANTS[50])/CONSTANTS[13], CONSTANTS[16])))/(( ( 0.400000*CONSTANTS[14]+1.00000)*STATES[0])/( CONSTANTS[14]*CONSTANTS[13])+1.00000));
ALGEBRAIC[25] =  CONSTANTS[10]*CONSTANTS[11]*ALGEBRAIC[0]*CONSTANTS[12]*ALGEBRAIC[24];
ALGEBRAIC[26] = (STATES[0]<=CONSTANTS[24] ?  CONSTANTS[11]*ALGEBRAIC[0]*(1.00000 -  CONSTANTS[10]*CONSTANTS[12]*ALGEBRAIC[24]) :  CONSTANTS[11]*( CONSTANTS[23]*(1.00000 -  CONSTANTS[10]*CONSTANTS[12]*ALGEBRAIC[24])+( CONSTANTS[32]*(STATES[0] - CONSTANTS[24]))/(CONSTANTS[31] - CONSTANTS[24])));
ALGEBRAIC[27] =  ALGEBRAIC[25]*ALGEBRAIC[6]*ALGEBRAIC[7]*ALGEBRAIC[8]*(1.00000 - STATES[1]) -  ALGEBRAIC[26]*STATES[1];
RATES[1] = ALGEBRAIC[27];
ALGEBRAIC[17] = (STATES[0]<=0.00000 ?  CONSTANTS[34]*exp( CONSTANTS[33]*STATES[3]) :  CONSTANTS[36]*exp( CONSTANTS[35]*STATES[3]));
ALGEBRAIC[30] = ALGEBRAIC[21]/ALGEBRAIC[24];
ALGEBRAIC[28] = ( CONSTANTS[14]*1.00000*(1.00000+ 0.400000*CONSTANTS[14]+( 1.20000*STATES[0])/CONSTANTS[13]+ 0.600000*pow(STATES[0]/CONSTANTS[13], 2.00000)))/( CONSTANTS[13]*pow( (CONSTANTS[14] - STATES[0]/CONSTANTS[13])*(1.00000+( 0.600000*STATES[0])/CONSTANTS[13]), 2.00000));
ALGEBRAIC[29] =  (1.00000/CONSTANTS[13])*exp( - CONSTANTS[25]*pow(STATES[0]/CONSTANTS[13] - CONSTANTS[50]/CONSTANTS[13], CONSTANTS[16]))*(( 0.400000*CONSTANTS[14]+1.00000)/CONSTANTS[14]+ CONSTANTS[25]*CONSTANTS[16]*(1.00000+( ( 0.400000*CONSTANTS[14]+1.00000)*STATES[0])/( CONSTANTS[14]*CONSTANTS[13]))*pow(STATES[0]/CONSTANTS[13] - CONSTANTS[50]/CONSTANTS[13], CONSTANTS[16] - 1.00000));
ALGEBRAIC[32] = (STATES[0]<=- CONSTANTS[13] ? CONSTANTS[51] : - CONSTANTS[13]<STATES[0]&&STATES[0]<=0.00000 ? ALGEBRAIC[28] : 0.00000<STATES[0]&&STATES[0]<=CONSTANTS[50] ? CONSTANTS[52] : ALGEBRAIC[29]);
ALGEBRAIC[15] = (STATES[0]<=0.00000 ? CONSTANTS[33] : CONSTANTS[35]);
ALGEBRAIC[33] = (ALGEBRAIC[13]==1.00000 ? - ( CONSTANTS[7]*ALGEBRAIC[27]*ALGEBRAIC[30]+ ALGEBRAIC[15]*ALGEBRAIC[17]*pow(STATES[0], 2.00000)+ CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])*STATES[2])/( CONSTANTS[7]*STATES[1]*ALGEBRAIC[32]+ALGEBRAIC[17]) : - ( CONSTANTS[7]*ALGEBRAIC[27]*ALGEBRAIC[30]+ ALGEBRAIC[15]*ALGEBRAIC[17]*pow(STATES[0], 2.00000)+ ( CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])+ CONSTANTS[5]*CONSTANTS[6]*exp( CONSTANTS[5]*STATES[5]))*STATES[2])/( CONSTANTS[7]*STATES[1]*ALGEBRAIC[32]+ALGEBRAIC[17]));
RATES[0] = (ALGEBRAIC[20]==0.00000 ? ( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*(ALGEBRAIC[18] - STATES[0]) - ( CONSTANTS[7]*ALGEBRAIC[27]*ALGEBRAIC[30]+ ALGEBRAIC[15]*ALGEBRAIC[17]*pow(STATES[0], 2.00000)))/( CONSTANTS[7]*STATES[1]*ALGEBRAIC[32]+ALGEBRAIC[17]) : ALGEBRAIC[33]);
ALGEBRAIC[19] = (STATES[2]<=STATES[0] ? CONSTANTS[37] : CONSTANTS[39]);
RATES[2] = (ALGEBRAIC[13]==1.00000 ? (( ALGEBRAIC[20]*(ALGEBRAIC[33] -  ALGEBRAIC[19]*pow(STATES[2] - STATES[0], 2.00000)) -  CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*(STATES[2] - STATES[0])) -  CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])*STATES[2])/ALGEBRAIC[20] : (ALGEBRAIC[33] -  ALGEBRAIC[19]*pow(STATES[2] - STATES[0], 2.00000)) - ( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*(STATES[2] - STATES[0])+ ( CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])+ CONSTANTS[5]*CONSTANTS[6]*exp( CONSTANTS[5]*STATES[5]))*STATES[2])/ALGEBRAIC[20]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[3] = STATES[0];
ALGEBRAIC[8] = (STATES[3]+CONSTANTS[19])/(0.460000+CONSTANTS[19]);
ALGEBRAIC[10] = ( CONSTANTS[41]*STATES[1])/( ALGEBRAIC[8]*STATES[6]);
ALGEBRAIC[12] = (ALGEBRAIC[10]>=1.00000 ? 1.00000 : pow(0.0200000, ALGEBRAIC[10]));
ALGEBRAIC[14] =  CONSTANTS[26]*(CONSTANTS[41] - STATES[6])*STATES[8] -  CONSTANTS[27]*exp( - CONSTANTS[28]*STATES[6])*ALGEBRAIC[12]*STATES[6];
ALGEBRAIC[16] =  CONSTANTS[43]*(CONSTANTS[42] - STATES[7])*STATES[8] -  CONSTANTS[44]*STATES[7];
ALGEBRAIC[20] = (STATES[2]<=STATES[0] ?  CONSTANTS[38]*exp( CONSTANTS[37]*(STATES[4] - STATES[3])) :  CONSTANTS[40]*exp( CONSTANTS[39]*(STATES[4] - STATES[3])));
ALGEBRAIC[9] = STATES[4]+STATES[5];
ALGEBRAIC[4] =  CONSTANTS[6]*(exp( CONSTANTS[5]*STATES[5]) - 1.00000);
ALGEBRAIC[5] = ALGEBRAIC[4];
ALGEBRAIC[11] = (ALGEBRAIC[9]>=CONSTANTS[18]&&VOI<0.150000 ? 0.00000 : 1.00000);
ALGEBRAIC[13] = (CONSTANTS[0]==0.00000 ? 0.00000 : CONSTANTS[0]==1.00000&&ALGEBRAIC[5]>=CONSTANTS[17] ? 1.00000 : CONSTANTS[0]==1.00000&&ALGEBRAIC[9]>=CONSTANTS[18]&&ALGEBRAIC[11]==1.00000 ? 0.00000 : 0.0/0.0);
ALGEBRAIC[18] = (ALGEBRAIC[13]==1.00000 ? ( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*STATES[0])/( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))+ CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])) : ( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*STATES[0])/( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))+ CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])+ CONSTANTS[5]*CONSTANTS[6]*exp( CONSTANTS[5]*STATES[5])));
ALGEBRAIC[22] = (ALGEBRAIC[20]==0.00000 ? ALGEBRAIC[18] : STATES[2]);
ALGEBRAIC[23] = (ALGEBRAIC[13]==1.00000 ? 0.00000 : ALGEBRAIC[13]==0.00000&&ALGEBRAIC[20]==0.00000 ? - ALGEBRAIC[18] : - STATES[2]);
ALGEBRAIC[6] = pow(STATES[6], CONSTANTS[9])/(pow(STATES[6], CONSTANTS[9])+pow(CONSTANTS[8], CONSTANTS[9]));
ALGEBRAIC[7] =  0.600000*STATES[3]+0.500000;
ALGEBRAIC[0] = (STATES[0]<=0.00000 ? CONSTANTS[20] - ( CONSTANTS[21]*STATES[0])/CONSTANTS[13] : STATES[0]<=CONSTANTS[24]&&0.00000<STATES[0] ? ( (CONSTANTS[23] - CONSTANTS[22])*STATES[0])/CONSTANTS[24]+CONSTANTS[22] : CONSTANTS[23]/pow(1.00000+( CONSTANTS[30]*(STATES[0] - CONSTANTS[24]))/CONSTANTS[13], CONSTANTS[29]));
ALGEBRAIC[21] = (STATES[0]<=0.00000 ? ( CONSTANTS[14]*(1.00000+STATES[0]/CONSTANTS[13]))/(CONSTANTS[14] - STATES[0]/CONSTANTS[13]) : (1.00000+CONSTANTS[15]) - ( pow(CONSTANTS[15], 2.00000)*CONSTANTS[14])/( (( CONSTANTS[14]*CONSTANTS[15])/CONSTANTS[53])*pow(STATES[0]/CONSTANTS[13], 2.00000)+( (CONSTANTS[14]+1.00000)*STATES[0])/CONSTANTS[13]+ CONSTANTS[14]*CONSTANTS[15]));
ALGEBRAIC[24] = (- CONSTANTS[13]<=STATES[0]&&STATES[0]<=0.00000 ? 1.00000+( 0.600000*STATES[0])/CONSTANTS[13] : 0.00000<STATES[0]&&STATES[0]<=CONSTANTS[50] ? ALGEBRAIC[21]/(( ( 0.400000*CONSTANTS[14]+1.00000)*STATES[0])/( CONSTANTS[14]*CONSTANTS[13])+1.00000) : ( ALGEBRAIC[21]*exp( - CONSTANTS[25]*pow((STATES[0] - CONSTANTS[50])/CONSTANTS[13], CONSTANTS[16])))/(( ( 0.400000*CONSTANTS[14]+1.00000)*STATES[0])/( CONSTANTS[14]*CONSTANTS[13])+1.00000));
ALGEBRAIC[25] =  CONSTANTS[10]*CONSTANTS[11]*ALGEBRAIC[0]*CONSTANTS[12]*ALGEBRAIC[24];
ALGEBRAIC[26] = (STATES[0]<=CONSTANTS[24] ?  CONSTANTS[11]*ALGEBRAIC[0]*(1.00000 -  CONSTANTS[10]*CONSTANTS[12]*ALGEBRAIC[24]) :  CONSTANTS[11]*( CONSTANTS[23]*(1.00000 -  CONSTANTS[10]*CONSTANTS[12]*ALGEBRAIC[24])+( CONSTANTS[32]*(STATES[0] - CONSTANTS[24]))/(CONSTANTS[31] - CONSTANTS[24])));
ALGEBRAIC[27] =  ALGEBRAIC[25]*ALGEBRAIC[6]*ALGEBRAIC[7]*ALGEBRAIC[8]*(1.00000 - STATES[1]) -  ALGEBRAIC[26]*STATES[1];
ALGEBRAIC[17] = (STATES[0]<=0.00000 ?  CONSTANTS[34]*exp( CONSTANTS[33]*STATES[3]) :  CONSTANTS[36]*exp( CONSTANTS[35]*STATES[3]));
ALGEBRAIC[30] = ALGEBRAIC[21]/ALGEBRAIC[24];
ALGEBRAIC[28] = ( CONSTANTS[14]*1.00000*(1.00000+ 0.400000*CONSTANTS[14]+( 1.20000*STATES[0])/CONSTANTS[13]+ 0.600000*pow(STATES[0]/CONSTANTS[13], 2.00000)))/( CONSTANTS[13]*pow( (CONSTANTS[14] - STATES[0]/CONSTANTS[13])*(1.00000+( 0.600000*STATES[0])/CONSTANTS[13]), 2.00000));
ALGEBRAIC[29] =  (1.00000/CONSTANTS[13])*exp( - CONSTANTS[25]*pow(STATES[0]/CONSTANTS[13] - CONSTANTS[50]/CONSTANTS[13], CONSTANTS[16]))*(( 0.400000*CONSTANTS[14]+1.00000)/CONSTANTS[14]+ CONSTANTS[25]*CONSTANTS[16]*(1.00000+( ( 0.400000*CONSTANTS[14]+1.00000)*STATES[0])/( CONSTANTS[14]*CONSTANTS[13]))*pow(STATES[0]/CONSTANTS[13] - CONSTANTS[50]/CONSTANTS[13], CONSTANTS[16] - 1.00000));
ALGEBRAIC[32] = (STATES[0]<=- CONSTANTS[13] ? CONSTANTS[51] : - CONSTANTS[13]<STATES[0]&&STATES[0]<=0.00000 ? ALGEBRAIC[28] : 0.00000<STATES[0]&&STATES[0]<=CONSTANTS[50] ? CONSTANTS[52] : ALGEBRAIC[29]);
ALGEBRAIC[15] = (STATES[0]<=0.00000 ? CONSTANTS[33] : CONSTANTS[35]);
ALGEBRAIC[33] = (ALGEBRAIC[13]==1.00000 ? - ( CONSTANTS[7]*ALGEBRAIC[27]*ALGEBRAIC[30]+ ALGEBRAIC[15]*ALGEBRAIC[17]*pow(STATES[0], 2.00000)+ CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])*STATES[2])/( CONSTANTS[7]*STATES[1]*ALGEBRAIC[32]+ALGEBRAIC[17]) : - ( CONSTANTS[7]*ALGEBRAIC[27]*ALGEBRAIC[30]+ ALGEBRAIC[15]*ALGEBRAIC[17]*pow(STATES[0], 2.00000)+ ( CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])+ CONSTANTS[5]*CONSTANTS[6]*exp( CONSTANTS[5]*STATES[5]))*STATES[2])/( CONSTANTS[7]*STATES[1]*ALGEBRAIC[32]+ALGEBRAIC[17]));
ALGEBRAIC[19] = (STATES[2]<=STATES[0] ? CONSTANTS[37] : CONSTANTS[39]);
ALGEBRAIC[1] =  CONSTANTS[2]*(exp( CONSTANTS[1]*(STATES[4] - STATES[3])) - 1.00000);
ALGEBRAIC[2] =  CONSTANTS[4]*(exp( CONSTANTS[3]*STATES[4]) - 1.00000);
ALGEBRAIC[31] =  CONSTANTS[7]*ALGEBRAIC[30]*STATES[1];
}