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]; }