A model for human ventricular tissue

Generated code

• The required array sizes are:
  • VARIABLES[85]
  • CONSTANTS[46]
  • BOUND[1]
  • RATES[17]
• The following variables are defined:
  • Variable xr1_inf(VARIABLES[17]; recomputed after rate changes)
  • Variable alpha_xr1(VARIABLES[18]; recomputed after rate changes)
  • Variable beta_xr1(VARIABLES[19]; recomputed after rate changes)
  • Variable tau_xr1(VARIABLES[20]; recomputed after rate changes)
  • Variable Xr2(RATES[0] and VARIABLES[0])
  • Variable xr2_inf(VARIABLES[21]; recomputed after rate changes)
  • Variable alpha_xr2(VARIABLES[22]; recomputed after rate changes)
  • Variable beta_xr2(VARIABLES[23]; recomputed after rate changes)
  • Variable tau_xr2(VARIABLES[24]; recomputed after rate changes)
  • Variable i_Ks(VARIABLES[25]; recomputed after rate changes)
  • Variable g_Ks(CONSTANTS[0])
  • Variable Xs(RATES[1] and VARIABLES[1])
  • Variable xs_inf(VARIABLES[26]; recomputed after rate changes)
  • Variable alpha_xs(VARIABLES[27]; recomputed after rate changes)
  • Variable beta_xs(VARIABLES[28]; recomputed after rate changes)
  • Variable tau_xs(VARIABLES[29]; recomputed after rate changes)
  • Variable i_Na(VARIABLES[30]; recomputed after rate changes)
  • Variable g_Na(CONSTANTS[1])
  • Variable i_K1(VARIABLES[31]; recomputed after rate changes)
  • Variable g_K1(CONSTANTS[2])
  • Variable xK1_inf(VARIABLES[32]; recomputed after rate changes)
  • Variable alpha_K1(VARIABLES[33]; recomputed after rate changes)
  • Variable beta_K1(VARIABLES[34]; recomputed after rate changes)
  • Variable i_Kr(VARIABLES[35]; recomputed after rate changes)
  • Variable g_Kr(CONSTANTS[3])
  • Variable Xr1(RATES[2] and VARIABLES[2])
  • Variable time(BOUND[0])
  • Variable V(RATES[3] and VARIABLES[3])
  • Variable R(CONSTANTS[4])
  • Variable g_bca(CONSTANTS[5])
  • Variable i_to(VARIABLES[36]; recomputed after rate changes)
  • Variable g_to(CONSTANTS[6])
  • Variable s(RATES[4] and VARIABLES[4])
  • Variable s_inf(VARIABLES[37]; recomputed after rate changes)
  • Variable tau_s(VARIABLES[38]; recomputed after rate changes)
  • Variable r(RATES[5] and VARIABLES[5])
  • Variable r_inf(VARIABLES[39]; recomputed after rate changes)
  • Variable tau_r(VARIABLES[40]; recomputed after rate changes)
  • Variable i_NaK(VARIABLES[41]; recomputed after rate changes)
  • Variable P_NaK(CONSTANTS[7])
  • Variable K_mk(CONSTANTS[8])
  • Variable K_mNa(CONSTANTS[9])
  • Variable i_NaCa(VARIABLES[42]; recomputed after rate changes)
  • Variable K_NaCa(CONSTANTS[10])
  • Variable K_sat(CONSTANTS[11])
  • Variable alpha(CONSTANTS[12])
  • Variable gamma(CONSTANTS[13])
  • Variable Km_Ca(CONSTANTS[14])
  • Variable m(RATES[6] and VARIABLES[6])
  • Variable m_inf(VARIABLES[43]; recomputed after rate changes)
  • Variable alpha_m(VARIABLES[44]; recomputed after rate changes)
  • Variable beta_m(VARIABLES[45]; recomputed after rate changes)
  • Variable tau_m(VARIABLES[46]; recomputed after rate changes)
  • Variable h(RATES[7] and VARIABLES[7])
  • Variable h_inf(VARIABLES[47]; recomputed after rate changes)
  • Variable alpha_h(VARIABLES[48]; recomputed after rate changes)
  • Variable beta_h(VARIABLES[49]; recomputed after rate changes)
  • Variable tau_h(VARIABLES[50]; recomputed after rate changes)
  • Variable j(RATES[8] and VARIABLES[8])
  • Variable j_inf(VARIABLES[51]; recomputed after rate changes)
  • Variable alpha_j(VARIABLES[52]; recomputed after rate changes)
  • Variable beta_j(VARIABLES[53]; recomputed after rate changes)
  • Variable tau_j(VARIABLES[54]; recomputed after rate changes)
  • Variable i_b_Na(VARIABLES[55]; recomputed after rate changes)
  • Variable g_bna(CONSTANTS[15])
  • Variable i_CaL(VARIABLES[56]; recomputed after rate changes)
  • Variable g_CaL(CONSTANTS[16])
  • Variable d(RATES[9] and VARIABLES[9])
  • Variable d_inf(VARIABLES[57]; recomputed after rate changes)
  • Variable alpha_d(VARIABLES[58]; recomputed after rate changes)
  • Variable beta_d(VARIABLES[59]; recomputed after rate changes)
  • Variable gamma_d(VARIABLES[60]; recomputed after rate changes)
  • Variable tau_d(VARIABLES[61]; recomputed after rate changes)
  • Variable f(RATES[10] and VARIABLES[10])
  • Variable f_inf(VARIABLES[62]; recomputed after rate changes)
  • Variable tau_f(VARIABLES[63]; recomputed after rate changes)
  • Variable fCa(RATES[11] and VARIABLES[11])
  • Variable alpha_fCa(VARIABLES[64]; recomputed after rate changes)
  • Variable beta_fCa(VARIABLES[65]; recomputed after rate changes)
  • Variable gama_fCa(VARIABLES[66]; recomputed after rate changes)
  • Variable fCa_inf(VARIABLES[67]; recomputed after rate changes)
  • Variable tau_fCa(VARIABLES[68]; computed once)
  • Variable d_fCa(VARIABLES[69]; recomputed after rate changes)
  • Variable i_b_Ca(VARIABLES[70]; recomputed after rate changes)
  • Variable d_g(VARIABLES[71]; recomputed after rate changes)
  • Variable Na_i(RATES[12] and VARIABLES[12])
  • Variable Na_o(CONSTANTS[17])
  • Variable K_i(RATES[13] and VARIABLES[13])
  • Variable K_o(CONSTANTS[18])
  • Variable Istim(VARIABLES[72]; recomputed after rate changes)
  • Variable IstimStart(CONSTANTS[19])
  • Variable IstimEnd(CONSTANTS[20])
  • Variable IstimAmplitude(CONSTANTS[21])
  • Variable IstimPeriod(CONSTANTS[22])
  • Variable IstimPulseDuration(CONSTANTS[23])
  • Variable T(CONSTANTS[24])
  • Variable F(CONSTANTS[25])
  • Variable Cm(CONSTANTS[26])
  • Variable V_c(CONSTANTS[27])
  • Variable E_Na(VARIABLES[73]; recomputed after rate changes)
  • Variable E_K(VARIABLES[74]; recomputed after rate changes)
  • Variable E_Ks(VARIABLES[75]; recomputed after rate changes)
  • Variable E_Ca(VARIABLES[76]; recomputed after rate changes)
  • Variable P_kna(CONSTANTS[28])
  • Variable Km_Nai(CONSTANTS[29])
  • Variable i_p_Ca(VARIABLES[77]; recomputed after rate changes)
  • Variable g_pCa(CONSTANTS[30])
  • Variable K_pCa(CONSTANTS[31])
  • Variable i_p_K(VARIABLES[78]; recomputed after rate changes)
  • Variable g_pK(CONSTANTS[32])
  • Variable Ca_i(RATES[14] and VARIABLES[14])
  • Variable Ca_o(CONSTANTS[33])
  • Variable Ca_SR(RATES[15] and VARIABLES[15])
  • Variable i_rel(VARIABLES[79]; recomputed after rate changes)
  • Variable i_up(VARIABLES[80]; recomputed after rate changes)
  • Variable i_leak(VARIABLES[81]; recomputed after rate changes)
  • Variable g(RATES[16] and VARIABLES[16])
  • Variable tau_g(CONSTANTS[34])
  • Variable g_inf(VARIABLES[82]; recomputed after rate changes)
  • Variable a_rel(CONSTANTS[35])
  • Variable b_rel(CONSTANTS[36])
  • Variable c_rel(CONSTANTS[37])
  • Variable K_up(CONSTANTS[38])
  • Variable V_leak(CONSTANTS[39])
  • Variable Vmax_up(CONSTANTS[40])
  • Variable Ca_i_bufc(VARIABLES[83]; recomputed after rate changes)
  • Variable Ca_sr_bufsr(VARIABLES[84]; recomputed after rate changes)
  • Variable Buf_c(CONSTANTS[41])
  • Variable K_buf_c(CONSTANTS[42])
  • Variable Buf_sr(CONSTANTS[43])
  • Variable K_buf_sr(CONSTANTS[44])
  • Variable V_sr(CONSTANTS[45])
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 0.245;
CONSTANTS[1] = 14.838;
CONSTANTS[2] = 5.405;
CONSTANTS[3] = 0.096;
CONSTANTS[4] = 8314.47;
CONSTANTS[5] = 0.000592;
CONSTANTS[6] = 0.073;
CONSTANTS[7] = 1.362;
CONSTANTS[8] = 1;
CONSTANTS[9] = 40;
CONSTANTS[10] = 1000;
CONSTANTS[11] = 0.1;
CONSTANTS[12] = 2.5;
CONSTANTS[13] = 0.35;
CONSTANTS[14] = 1.38;
CONSTANTS[15] = 0.00029;
CONSTANTS[16] = 0.000175;
CONSTANTS[17] = 140;
CONSTANTS[18] = 5.4;
CONSTANTS[19] = 100;
CONSTANTS[20] = 50000;
CONSTANTS[21] = -52;
CONSTANTS[22] = 1000;
CONSTANTS[23] = 1;
CONSTANTS[24] = 310;
CONSTANTS[25] = 96485.3;
CONSTANTS[26] = 0.185;
CONSTANTS[27] = 0.016404;
CONSTANTS[28] = 0.03;
CONSTANTS[29] = 87.5;
CONSTANTS[30] = 0.825;
CONSTANTS[31] = 0.0005;
CONSTANTS[32] = 0.0146;
CONSTANTS[33] = 2;
CONSTANTS[34] = 2;
CONSTANTS[35] = 0.016464;
CONSTANTS[36] = 0.25;
CONSTANTS[37] = 0.008232;
CONSTANTS[38] = 0.00025;
CONSTANTS[39] = 8e-05;
CONSTANTS[40] = 0.000425;
CONSTANTS[41] = 0.15;
CONSTANTS[42] = 0.001;
CONSTANTS[43] = 10;
CONSTANTS[44] = 0.3;
CONSTANTS[45] = 0.001094;
}

void SetupComputedConstants(double* CONSTANTS, double* VARIABLES)
{
VARIABLES[68] = 2.00000;
VARIABLES[0] = 1.00000;
VARIABLES[1] = 0.00000;
VARIABLES[2] = 0.00000;
VARIABLES[3] = -86.2000;
VARIABLES[4] = 1.00000;
VARIABLES[5] = 0.00000;
VARIABLES[6] = 0.00000;
VARIABLES[7] = 0.750000;
VARIABLES[8] = 0.750000;
VARIABLES[9] = 0.00000;
VARIABLES[10] = 1.00000;
VARIABLES[11] = 1.00000;
VARIABLES[12] = 11.6000;
VARIABLES[13] = 138.300;
VARIABLES[14] = 0.000200000;
VARIABLES[15] = 0.200000;
VARIABLES[16] = 1.00000;
}

void ComputeRates(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
RATES[3] = (((-1.00000)/1.00000)*(VARIABLES[31]+VARIABLES[36]+VARIABLES[35]+VARIABLES[25]+VARIABLES[56]+VARIABLES[41]+VARIABLES[30]+VARIABLES[55]+VARIABLES[42]+VARIABLES[70]+VARIABLES[78]+VARIABLES[77]+VARIABLES[72]));
RATES[2] = ((VARIABLES[17]-VARIABLES[2])/VARIABLES[20]);
RATES[0] = ((VARIABLES[21]-VARIABLES[0])/VARIABLES[24]);
RATES[1] = ((VARIABLES[26]-VARIABLES[1])/VARIABLES[29]);
RATES[6] = ((VARIABLES[43]-VARIABLES[6])/VARIABLES[46]);
RATES[7] = ((VARIABLES[47]-VARIABLES[7])/VARIABLES[50]);
RATES[8] = ((VARIABLES[51]-VARIABLES[8])/VARIABLES[54]);
RATES[9] = ((VARIABLES[57]-VARIABLES[9])/VARIABLES[61]);
RATES[10] = ((VARIABLES[62]-VARIABLES[10])/VARIABLES[63]);
RATES[11] = (((((0.0100000*VARIABLES[69])>0.00000))&&((VARIABLES[3]>(-60.0000))))) ? (0.00000) : (VARIABLES[69]);
RATES[4] = ((VARIABLES[37]-VARIABLES[4])/VARIABLES[38]);
RATES[5] = ((VARIABLES[39]-VARIABLES[5])/VARIABLES[40]);
RATES[16] = (((((0.0100000*VARIABLES[71])>0.00000))&&((VARIABLES[3]>(-60.0000))))) ? (0.00000) : (VARIABLES[71]);
RATES[14] = (VARIABLES[83]*(((VARIABLES[81]-VARIABLES[80])+VARIABLES[79])-(((1.00000*((VARIABLES[56]+VARIABLES[70]+VARIABLES[77])-(2.00000*VARIABLES[42])))/(2.00000*1.00000*CONSTANTS[27]*CONSTANTS[25]))*CONSTANTS[26])));
RATES[15] = (((VARIABLES[84]*CONSTANTS[27])/CONSTANTS[45])*(VARIABLES[80]-(VARIABLES[79]+VARIABLES[81])));
RATES[12] = ((((-1.00000)*(VARIABLES[30]+VARIABLES[55]+(3.00000*VARIABLES[41])+(3.00000*VARIABLES[42])))/(1.00000*CONSTANTS[27]*CONSTANTS[25]))*CONSTANTS[26]);
RATES[13] = ((((-1.00000)*((VARIABLES[31]+VARIABLES[36]+VARIABLES[35]+VARIABLES[25]+VARIABLES[78]+VARIABLES[72])-(2.00000*VARIABLES[41])))/(1.00000*CONSTANTS[27]*CONSTANTS[25]))*CONSTANTS[26]);
}

void ComputeVariables(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
VARIABLES[18] = (450.000/(1.00000+exp((((-45.0000)-VARIABLES[3])/10.0000))));
VARIABLES[19] = (6.00000/(1.00000+exp(((VARIABLES[3]+30.0000)/11.5000))));
VARIABLES[20] = (1.00000*VARIABLES[18]*VARIABLES[19]);
VARIABLES[21] = (1.00000/(1.00000+exp(((VARIABLES[3]+88.0000)/24.0000))));
VARIABLES[22] = (3.00000/(1.00000+exp((((-60.0000)-VARIABLES[3])/20.0000))));
VARIABLES[23] = (1.12000/(1.00000+exp(((VARIABLES[3]-60.0000)/20.0000))));
VARIABLES[24] = (1.00000*VARIABLES[22]*VARIABLES[23]);
VARIABLES[75] = (((CONSTANTS[4]*CONSTANTS[24])/CONSTANTS[25])*log(((CONSTANTS[18]+(CONSTANTS[28]*CONSTANTS[17]))/(VARIABLES[13]+(CONSTANTS[28]*VARIABLES[12])))));
VARIABLES[25] = (CONSTANTS[0]*pow(VARIABLES[1],2.00000)*(VARIABLES[3]-VARIABLES[75]));
VARIABLES[26] = (1.00000/(1.00000+exp((((-5.00000)-VARIABLES[3])/14.0000))));
VARIABLES[27] = (1100.00/pow((1.00000+exp((((-10.0000)-VARIABLES[3])/6.00000))), 1.0 / (2.0)));
VARIABLES[28] = (1.00000/(1.00000+exp(((VARIABLES[3]-60.0000)/20.0000))));
VARIABLES[29] = (1.00000*VARIABLES[27]*VARIABLES[28]);
VARIABLES[73] = (((CONSTANTS[4]*CONSTANTS[24])/CONSTANTS[25])*log((CONSTANTS[17]/VARIABLES[12])));
VARIABLES[30] = (CONSTANTS[1]*pow(VARIABLES[6],3.00000)*VARIABLES[7]*VARIABLES[8]*(VARIABLES[3]-VARIABLES[73]));
VARIABLES[74] = (((CONSTANTS[4]*CONSTANTS[24])/CONSTANTS[25])*log((CONSTANTS[18]/VARIABLES[13])));
VARIABLES[33] = (0.100000/(1.00000+exp((0.0600000*((VARIABLES[3]-VARIABLES[74])-200.000)))));
VARIABLES[34] = (((3.00000*exp((0.000200000*((VARIABLES[3]-VARIABLES[74])+100.000))))+exp((0.100000*((VARIABLES[3]-VARIABLES[74])-10.0000))))/(1.00000+exp(((-0.500000)*(VARIABLES[3]-VARIABLES[74])))));
VARIABLES[32] = (VARIABLES[33]/(VARIABLES[33]+VARIABLES[34]));
VARIABLES[31] = (CONSTANTS[2]*VARIABLES[32]*pow((CONSTANTS[18]/5.40000), 1.0 / (2.0))*(VARIABLES[3]-VARIABLES[74]));
VARIABLES[35] = (CONSTANTS[3]*pow((CONSTANTS[18]/5.40000), 1.0 / (2.0))*VARIABLES[2]*VARIABLES[0]*(VARIABLES[3]-VARIABLES[74]));
VARIABLES[36] = (CONSTANTS[6]*VARIABLES[5]*VARIABLES[4]*(VARIABLES[3]-VARIABLES[74]));
VARIABLES[37] = (1.00000/(1.00000+exp(((VARIABLES[3]+28.0000)/5.00000))));
VARIABLES[38] = ((1000.00*exp(((-pow((VARIABLES[3]+67.0000),2.00000))/1000.00)))+8.00000);
VARIABLES[39] = (1.00000/(1.00000+exp(((20.0000-VARIABLES[3])/6.00000))));
VARIABLES[40] = ((9.50000*exp(((-pow((VARIABLES[3]+40.0000),2.00000))/1800.00)))+0.800000);
VARIABLES[76] = (((0.500000*CONSTANTS[4]*CONSTANTS[24])/CONSTANTS[25])*log((CONSTANTS[33]/VARIABLES[14])));
VARIABLES[77] = ((CONSTANTS[30]*VARIABLES[14])/(VARIABLES[14]+CONSTANTS[31]));
VARIABLES[78] = ((CONSTANTS[32]*(VARIABLES[3]-VARIABLES[74]))/(1.00000+exp(((25.0000-VARIABLES[3])/5.98000))));
VARIABLES[79] = ((((CONSTANTS[35]*pow(VARIABLES[15],2.00000))/(pow(CONSTANTS[36],2.00000)+pow(VARIABLES[15],2.00000)))+CONSTANTS[37])*VARIABLES[9]*VARIABLES[16]);
VARIABLES[80] = (CONSTANTS[40]/(1.00000+(pow(CONSTANTS[38],2.00000)/pow(VARIABLES[14],2.00000))));
VARIABLES[81] = (CONSTANTS[39]*(VARIABLES[15]-VARIABLES[14]));
VARIABLES[82] = (((VARIABLES[14]<0.000350000))) ? ((1.00000/(1.00000+pow((VARIABLES[14]/0.000350000),6.00000)))) : ((1.00000/(1.00000+pow((VARIABLES[14]/0.000350000),16.0000))));
VARIABLES[83] = (1.00000/(1.00000+((CONSTANTS[41]*CONSTANTS[42])/pow((VARIABLES[14]+CONSTANTS[42]),2.00000))));
VARIABLES[84] = (1.00000/(1.00000+((CONSTANTS[43]*CONSTANTS[44])/pow((VARIABLES[15]+CONSTANTS[44]),2.00000))));
VARIABLES[41] = (((((CONSTANTS[7]*CONSTANTS[18])/(CONSTANTS[18]+CONSTANTS[8]))*VARIABLES[12])/(VARIABLES[12]+CONSTANTS[9]))/(1.00000+(0.124500*exp((((-0.100000)*VARIABLES[3]*CONSTANTS[25])/(CONSTANTS[4]*CONSTANTS[24]))))+(0.0353000*exp((((-VARIABLES[3])*CONSTANTS[25])/(CONSTANTS[4]*CONSTANTS[24]))))));
VARIABLES[42] = ((CONSTANTS[10]*((exp(((CONSTANTS[13]*VARIABLES[3]*CONSTANTS[25])/(CONSTANTS[4]*CONSTANTS[24])))*pow(VARIABLES[12],3.00000)*CONSTANTS[33])-(exp((((CONSTANTS[13]-1.00000)*VARIABLES[3]*CONSTANTS[25])/(CONSTANTS[4]*CONSTANTS[24])))*pow(CONSTANTS[17],3.00000)*VARIABLES[14]*CONSTANTS[12])))/((pow(CONSTANTS[29],3.00000)+pow(CONSTANTS[17],3.00000))*(CONSTANTS[14]+CONSTANTS[33])*(1.00000+(CONSTANTS[11]*exp((((CONSTANTS[13]-1.00000)*VARIABLES[3]*CONSTANTS[25])/(CONSTANTS[4]*CONSTANTS[24])))))));
VARIABLES[43] = (1.00000/pow((1.00000+exp((((-56.8600)-VARIABLES[3])/9.03000))),2.00000));
VARIABLES[44] = (1.00000/(1.00000+exp((((-60.0000)-VARIABLES[3])/5.00000))));
VARIABLES[45] = ((0.100000/(1.00000+exp(((VARIABLES[3]+35.0000)/5.00000))))+(0.100000/(1.00000+exp(((VARIABLES[3]-50.0000)/200.000)))));
VARIABLES[46] = (1.00000*VARIABLES[44]*VARIABLES[45]);
VARIABLES[47] = (1.00000/pow((1.00000+exp(((VARIABLES[3]+71.5500)/7.43000))),2.00000));
VARIABLES[48] = (((VARIABLES[3]<(-40.0000)))) ? ((0.0570000*exp(((-(VARIABLES[3]+80.0000))/6.80000)))) : (0.00000);
VARIABLES[49] = (((VARIABLES[3]<(-40.0000)))) ? (((2.70000*exp((0.0790000*VARIABLES[3])))+(310000.*exp((0.348500*VARIABLES[3]))))) : ((0.770000/(0.130000*(1.00000+exp(((VARIABLES[3]+10.6600)/(-11.1000)))))));
VARIABLES[50] = (1.00000/(VARIABLES[48]+VARIABLES[49]));
VARIABLES[51] = (1.00000/pow((1.00000+exp(((VARIABLES[3]+71.5500)/7.43000))),2.00000));
VARIABLES[52] = (((VARIABLES[3]<(-40.0000)))) ? (((((((-25428.0)*exp((0.244400*VARIABLES[3])))-(6.94800e-06*exp(((-0.0439100)*VARIABLES[3]))))*(VARIABLES[3]+37.7800))/1.00000)/(1.00000+exp((0.311000*(VARIABLES[3]+79.2300)))))) : (0.00000);
VARIABLES[53] = (((VARIABLES[3]<(-40.0000)))) ? (((0.0242400*exp(((-0.0105200)*VARIABLES[3])))/(1.00000+exp(((-0.137800)*(VARIABLES[3]+40.1400)))))) : (((0.600000*exp((0.0570000*VARIABLES[3])))/(1.00000+exp(((-0.100000)*(VARIABLES[3]+32.0000))))));
VARIABLES[54] = (1.00000/(VARIABLES[52]+VARIABLES[53]));
VARIABLES[55] = (CONSTANTS[15]*(VARIABLES[3]-VARIABLES[73]));
VARIABLES[56] = ((((CONSTANTS[16]*VARIABLES[9]*VARIABLES[10]*VARIABLES[11]*4.00000*VARIABLES[3]*pow(CONSTANTS[25],2.00000))/(CONSTANTS[4]*CONSTANTS[24]))*((VARIABLES[14]*exp(((2.00000*VARIABLES[3]*CONSTANTS[25])/(CONSTANTS[4]*CONSTANTS[24]))))-(0.341000*CONSTANTS[33])))/(exp(((2.00000*VARIABLES[3]*CONSTANTS[25])/(CONSTANTS[4]*CONSTANTS[24])))-1.00000));
VARIABLES[57] = (1.00000/(1.00000+exp((((-5.00000)-VARIABLES[3])/7.50000))));
VARIABLES[58] = ((1.40000/(1.00000+exp((((-35.0000)-VARIABLES[3])/13.0000))))+0.250000);
VARIABLES[59] = (1.40000/(1.00000+exp(((VARIABLES[3]+5.00000)/5.00000))));
VARIABLES[60] = (1.00000/(1.00000+exp(((50.0000-VARIABLES[3])/20.0000))));
VARIABLES[61] = ((1.00000*VARIABLES[58]*VARIABLES[59])+VARIABLES[60]);
VARIABLES[62] = (1.00000/(1.00000+exp(((VARIABLES[3]+20.0000)/7.00000))));
VARIABLES[63] = ((1125.00*exp(((-pow((VARIABLES[3]+27.0000),2.00000))/240.000)))+80.0000+(165.000/(1.00000+exp(((25.0000-VARIABLES[3])/10.0000)))));
VARIABLES[64] = (1.00000/(1.00000+pow((VARIABLES[14]/0.000325000),8.00000)));
VARIABLES[65] = (0.100000/(1.00000+exp(((VARIABLES[14]-0.000500000)/0.000100000))));
VARIABLES[66] = (0.200000/(1.00000+exp(((VARIABLES[14]-0.000750000)/0.000800000))));
VARIABLES[67] = ((VARIABLES[64]+VARIABLES[65]+VARIABLES[66]+0.230000)/1.46000);
VARIABLES[69] = ((VARIABLES[67]-VARIABLES[11])/VARIABLES[68]);
VARIABLES[70] = (CONSTANTS[5]*(VARIABLES[3]-VARIABLES[76]));
VARIABLES[71] = ((VARIABLES[82]-VARIABLES[16])/CONSTANTS[34]);
VARIABLES[72] = ((((BOUND[0]>=CONSTANTS[19]))&&((BOUND[0]<=CONSTANTS[20]))&&((((BOUND[0]-CONSTANTS[19])-(floor(((BOUND[0]-CONSTANTS[19])/CONSTANTS[22]))*CONSTANTS[22]))<=CONSTANTS[23])))) ? (CONSTANTS[21]) : (0.00000);
VARIABLES[17] = (1.00000/(1.00000+exp((((-26.0000)-VARIABLES[3])/7.00000))));
#ifndef VARIABLES_FOR_RATES_ONLY
#endif
}