Intra- and inter-islet synchronization of metabolically driven insulin secretion

Generated code

• The required array sizes are:
  • VARIABLES[94]
  • CONSTANTS[72]
  • BOUND[1]
  • RATES[11]
• The following variables are defined:
  • Variable time(BOUND[0])
  • Variable Vm(RATES[0] and VARIABLES[0])
  • Variable cm(CONSTANTS[0])
  • Variable Ik(VARIABLES[11]; recomputed after rate changes)
  • Variable gK(CONSTANTS[1])
  • Variable n(RATES[1] and VARIABLES[1])
  • Variable n_infinity(VARIABLES[12]; recomputed after rate changes)
  • Variable tau_n(CONSTANTS[2])
  • Variable vn(CONSTANTS[3])
  • Variable sn(CONSTANTS[4])
  • Variable Ica(VARIABLES[13]; recomputed after rate changes)
  • Variable gCa(CONSTANTS[5])
  • Variable m_infinity(VARIABLES[14]; recomputed after rate changes)
  • Variable v(CONSTANTS[6])
  • Variable sm(CONSTANTS[7])
  • Variable Ikca(VARIABLES[15]; recomputed after rate changes)
  • Variable gkCa(CONSTANTS[8])
  • Variable kd(CONSTANTS[9])
  • Variable Ikatp(VARIABLES[16]; recomputed after rate changes)
  • Variable gkATP_(CONSTANTS[10])
  • Variable katpo(VARIABLES[17]; recomputed after rate changes)
  • Variable topo(VARIABLES[18]; recomputed after rate changes)
  • Variable bottomo(VARIABLES[19]; recomputed after rate changes)
  • Variable mgadp(VARIABLES[20]; recomputed after rate changes)
  • Variable adp3m(VARIABLES[21]; recomputed after rate changes)
  • Variable atp4m(VARIABLES[22]; recomputed after rate changes)
  • Variable JGPDH(VARIABLES[23]; recomputed after rate changes)
  • Variable kGPDH(CONSTANTS[11])
  • Variable G6P(RATES[2] and VARIABLES[2])
  • Variable FBP(RATES[3] and VARIABLES[3])
  • Variable F6P(VARIABLES[24]; recomputed after rate changes)
  • Variable Juni(VARIABLES[25]; recomputed after rate changes)
  • Variable p21(CONSTANTS[12])
  • Variable p22(CONSTANTS[13])
  • Variable JNaCa(VARIABLES[26]; recomputed after rate changes)
  • Variable p23(CONSTANTS[14])
  • Variable p24(CONSTANTS[15])
  • Variable Cam(RATES[4] and VARIABLES[4])
  • Variable fmito(CONSTANTS[16])
  • Variable Jmito(VARIABLES[27]; recomputed after rate changes)
  • Variable ADPm(RATES[5] and VARIABLES[5])
  • Variable ATPm(VARIABLES[28]; recomputed after rate changes)
  • Variable Amtot(CONSTANTS[17])
  • Variable adp(RATES[6] and VARIABLES[6])
  • Variable Jhyd(VARIABLES[29]; recomputed after rate changes)
  • Variable khyd(CONSTANTS[18])
  • Variable khydbas(CONSTANTS[19])
  • Variable atp(VARIABLES[30]; recomputed after rate changes)
  • Variable atot(CONSTANTS[20])
  • Variable c(RATES[7] and VARIABLES[7])
  • Variable fcyt(CONSTANTS[21])
  • Variable Jmem(VARIABLES[31]; recomputed after rate changes)
  • Variable kPMCA(CONSTANTS[22])
  • Variable alpha(CONSTANTS[23])
  • Variable Cbas(CONSTANTS[24])
  • Variable Jleak(VARIABLES[32]; recomputed after rate changes)
  • Variable pleak(CONSTANTS[25])
  • Variable JSERCA(VARIABLES[33]; recomputed after rate changes)
  • Variable kSERCA(CONSTANTS[26])
  • Variable Jer(VARIABLES[34]; recomputed after rate changes)
  • Variable Caer(RATES[8] and VARIABLES[8])
  • Variable fer(CONSTANTS[27])
  • Variable Vc_Ver(CONSTANTS[28])
  • Variable gamma(CONSTANTS[29])
  • Variable kappa(CONSTANTS[30])
  • Variable delta(VARIABLES[35]; computed once)
  • Variable VK(CONSTANTS[31])
  • Variable VCa(CONSTANTS[32])
  • Variable AMP(CONSTANTS[33])
  • Variable JPFK(VARIABLES[36]; recomputed after rate changes)
  • Variable JPFK_ms(VARIABLES[37]; recomputed after rate changes)
  • Variable bottom1(CONSTANTS[34])
  • Variable weight1(CONSTANTS[35])
  • Variable topa1(CONSTANTS[36])
  • Variable k1(CONSTANTS[37])
  • Variable k2(CONSTANTS[38])
  • Variable k3(CONSTANTS[39])
  • Variable k4(CONSTANTS[40])
  • Variable VmaxPFK(CONSTANTS[41])
  • Variable weight2(VARIABLES[38]; recomputed after rate changes)
  • Variable topa2(VARIABLES[39]; computed once)
  • Variable bottom2(VARIABLES[40]; recomputed after rate changes)
  • Variable topa3(VARIABLES[41]; recomputed after rate changes)
  • Variable weight3(VARIABLES[42]; recomputed after rate changes)
  • Variable bottom3(VARIABLES[43]; recomputed after rate changes)
  • Variable f13(CONSTANTS[42])
  • Variable f43(CONSTANTS[43])
  • Variable f23(CONSTANTS[44])
  • Variable f42(CONSTANTS[45])
  • Variable f41(CONSTANTS[46])
  • Variable weight4(VARIABLES[44]; recomputed after rate changes)
  • Variable topa4(VARIABLES[45]; recomputed after rate changes)
  • Variable bottom4(VARIABLES[46]; recomputed after rate changes)
  • Variable weight5(VARIABLES[47]; recomputed after rate changes)
  • Variable topa5(VARIABLES[48]; recomputed after rate changes)
  • Variable bottom5(VARIABLES[49]; recomputed after rate changes)
  • Variable weight6(VARIABLES[50]; recomputed after rate changes)
  • Variable topa6(VARIABLES[51]; recomputed after rate changes)
  • Variable bottom6(VARIABLES[52]; recomputed after rate changes)
  • Variable weight7(VARIABLES[53]; recomputed after rate changes)
  • Variable topa7(VARIABLES[54]; recomputed after rate changes)
  • Variable bottom7(VARIABLES[55]; recomputed after rate changes)
  • Variable weight8(VARIABLES[56]; recomputed after rate changes)
  • Variable topa8(VARIABLES[57]; recomputed after rate changes)
  • Variable bottom8(VARIABLES[58]; recomputed after rate changes)
  • Variable weight9(VARIABLES[59]; computed once)
  • Variable topa9(VARIABLES[60]; recomputed after rate changes)
  • Variable bottom9(VARIABLES[61]; recomputed after rate changes)
  • Variable weight10(VARIABLES[62]; recomputed after rate changes)
  • Variable topa10(VARIABLES[63]; recomputed after rate changes)
  • Variable bottom10(VARIABLES[64]; recomputed after rate changes)
  • Variable weight11(VARIABLES[65]; recomputed after rate changes)
  • Variable topa11(VARIABLES[66]; recomputed after rate changes)
  • Variable bottom11(VARIABLES[67]; recomputed after rate changes)
  • Variable weight12(VARIABLES[68]; recomputed after rate changes)
  • Variable topa12(VARIABLES[69]; recomputed after rate changes)
  • Variable bottom12(VARIABLES[70]; recomputed after rate changes)
  • Variable weight13(VARIABLES[71]; recomputed after rate changes)
  • Variable topa13(VARIABLES[72]; recomputed after rate changes)
  • Variable bottom13(VARIABLES[73]; recomputed after rate changes)
  • Variable weight14(VARIABLES[74]; recomputed after rate changes)
  • Variable topa14(VARIABLES[75]; recomputed after rate changes)
  • Variable bottom14(VARIABLES[76]; recomputed after rate changes)
  • Variable weight15(VARIABLES[77]; recomputed after rate changes)
  • Variable topa15(VARIABLES[78]; recomputed after rate changes)
  • Variable bottom15(VARIABLES[79]; recomputed after rate changes)
  • Variable weight16(VARIABLES[80]; recomputed after rate changes)
  • Variable topa16(VARIABLES[81]; recomputed after rate changes)
  • Variable bottom16(VARIABLES[82]; recomputed after rate changes)
  • Variable topb(VARIABLES[83]; recomputed after rate changes)
  • Variable lambda(CONSTANTS[47])
  • Variable JPDH(VARIABLES[84]; recomputed after rate changes)
  • Variable p1(CONSTANTS[48])
  • Variable p2(CONSTANTS[49])
  • Variable p3(CONSTANTS[50])
  • Variable JGPDHbas(CONSTANTS[51])
  • Variable NADHm(RATES[9] and VARIABLES[9])
  • Variable JO(VARIABLES[85]; recomputed after rate changes)
  • Variable p4(CONSTANTS[52])
  • Variable p5(CONSTANTS[53])
  • Variable p6(CONSTANTS[54])
  • Variable p7(CONSTANTS[55])
  • Variable NADm(VARIABLES[86]; recomputed after rate changes)
  • Variable NADmtot(CONSTANTS[56])
  • Variable delta_psi(RATES[10] and VARIABLES[10])
  • Variable Cmito(CONSTANTS[57])
  • Variable JHres(VARIABLES[87]; recomputed after rate changes)
  • Variable p8(CONSTANTS[58])
  • Variable p9(CONSTANTS[59])
  • Variable p10(CONSTANTS[60])
  • Variable p11(CONSTANTS[61])
  • Variable JF1F0(VARIABLES[88]; recomputed after rate changes)
  • Variable p13(CONSTANTS[62])
  • Variable p14(CONSTANTS[63])
  • Variable p15(CONSTANTS[64])
  • Variable p16(CONSTANTS[65])
  • Variable JHatp(VARIABLES[89]; recomputed after rate changes)
  • Variable JGK_ms(VARIABLES[90]; computed once)
  • Variable JGK(CONSTANTS[66])
  • Variable JHleak(VARIABLES[91]; recomputed after rate changes)
  • Variable p17(CONSTANTS[67])
  • Variable p18(CONSTANTS[68])
  • Variable JANT(VARIABLES[92]; recomputed after rate changes)
  • Variable p19(CONSTANTS[69])
  • Variable p20(CONSTANTS[70])
  • Variable FRT(CONSTANTS[71])
  • Variable RATm(VARIABLES[93]; recomputed after rate changes)
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 5300;
CONSTANTS[1] = 2700;
CONSTANTS[2] = 20;
CONSTANTS[3] = -16;
CONSTANTS[4] = 5;
CONSTANTS[5] = 1000;
CONSTANTS[6] = -20;
CONSTANTS[7] = 12;
CONSTANTS[8] = 300;
CONSTANTS[9] = 0.5;
CONSTANTS[10] = 16000;
CONSTANTS[11] = 0.0005;
CONSTANTS[12] = 0.04;
CONSTANTS[13] = 1.1;
CONSTANTS[14] = 0.01;
CONSTANTS[15] = 0.016;
CONSTANTS[16] = 0.01;
CONSTANTS[17] = 15;
CONSTANTS[18] = 5e-05;
CONSTANTS[19] = 5e-05;
CONSTANTS[20] = 2500;
CONSTANTS[21] = 0.01;
CONSTANTS[22] = 0.1;
CONSTANTS[23] = 4.5e-06;
CONSTANTS[24] = 0.05;
CONSTANTS[25] = 0.0002;
CONSTANTS[26] = 0.4;
CONSTANTS[27] = 0.01;
CONSTANTS[28] = 31;
CONSTANTS[29] = 0.001;
CONSTANTS[30] = 0.001;
CONSTANTS[31] = -75;
CONSTANTS[32] = 25;
CONSTANTS[33] = 500;
CONSTANTS[34] = 1;
CONSTANTS[35] = 1;
CONSTANTS[36] = 0;
CONSTANTS[37] = 30;
CONSTANTS[38] = 1;
CONSTANTS[39] = 50000;
CONSTANTS[40] = 1000;
CONSTANTS[41] = 5;
CONSTANTS[42] = 0.02;
CONSTANTS[43] = 20;
CONSTANTS[44] = 0.2;
CONSTANTS[45] = 20;
CONSTANTS[46] = 20;
CONSTANTS[47] = 0.06;
CONSTANTS[48] = 400;
CONSTANTS[49] = 1;
CONSTANTS[50] = 0.01;
CONSTANTS[51] = 0.0005;
CONSTANTS[52] = 0.6;
CONSTANTS[53] = 0.1;
CONSTANTS[54] = 177;
CONSTANTS[55] = 5;
CONSTANTS[56] = 10;
CONSTANTS[57] = 1.8;
CONSTANTS[58] = 7;
CONSTANTS[59] = 0.1;
CONSTANTS[60] = 177;
CONSTANTS[61] = 5;
CONSTANTS[62] = 10;
CONSTANTS[63] = 190;
CONSTANTS[64] = 8.5;
CONSTANTS[65] = 35;
CONSTANTS[66] = 0.4;
CONSTANTS[67] = 0.002;
CONSTANTS[68] = -0.03;
CONSTANTS[69] = 0.35;
CONSTANTS[70] = 2;
CONSTANTS[71] = 0.0374101;
}

void SetupComputedConstants(double* CONSTANTS, double* VARIABLES)
{
VARIABLES[35] = (3.90000/53.2000);
VARIABLES[39] = CONSTANTS[36];
VARIABLES[59] = (CONSTANTS[33]/CONSTANTS[37]);
VARIABLES[90] = (CONSTANTS[30]*CONSTANTS[66]);
VARIABLES[0] = -60.0000;
VARIABLES[1] = 0.00000;
VARIABLES[2] = 301.000;
VARIABLES[3] = 2.16000;
VARIABLES[4] = 0.200000;
VARIABLES[5] = 11.1000;
VARIABLES[6] = 1137.00;
VARIABLES[7] = 0.170000;
VARIABLES[8] = 345.000;
VARIABLES[9] = 0.400000;
VARIABLES[10] = 164.000;
}

void ComputeRates(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
RATES[0] = ((-(VARIABLES[11]+VARIABLES[13]+VARIABLES[15]+VARIABLES[16]))/CONSTANTS[0]);
RATES[1] = ((VARIABLES[12]-VARIABLES[1])/CONSTANTS[2]);
RATES[2] = (VARIABLES[90]-VARIABLES[37]);
RATES[3] = (VARIABLES[37]-(0.500000*VARIABLES[23]));
RATES[9] = (CONSTANTS[29]*(VARIABLES[84]-VARIABLES[85]));
RATES[10] = ((VARIABLES[87]-(VARIABLES[89]+VARIABLES[92]+VARIABLES[91]+VARIABLES[26]+(2.00000*VARIABLES[25])))/CONSTANTS[57]);
RATES[4] = ((-CONSTANTS[16])*VARIABLES[27]);
RATES[5] = (CONSTANTS[29]*(VARIABLES[92]-VARIABLES[88]));
RATES[6] = (((-VARIABLES[35])*VARIABLES[92])+VARIABLES[29]);
RATES[7] = (CONSTANTS[21]*(VARIABLES[31]+VARIABLES[34]+(VARIABLES[35]*VARIABLES[27])));
RATES[8] = ((-CONSTANTS[27])*CONSTANTS[28]*VARIABLES[34]);
}

void ComputeVariables(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
VARIABLES[11] = (CONSTANTS[1]*VARIABLES[1]*(VARIABLES[0]-CONSTANTS[31]));
VARIABLES[12] = (1.00000/(1.00000+exp(((CONSTANTS[3]-VARIABLES[0])/CONSTANTS[4]))));
VARIABLES[14] = (1.00000/(1.00000+exp(((CONSTANTS[6]-VARIABLES[0])/CONSTANTS[7]))));
VARIABLES[13] = (CONSTANTS[5]*VARIABLES[14]*(VARIABLES[0]-CONSTANTS[32]));
VARIABLES[15] = ((CONSTANTS[8]/(1.00000+pow((CONSTANTS[9]/VARIABLES[7]),2.00000)))*(VARIABLES[0]-CONSTANTS[31]));
VARIABLES[20] = (0.165000*VARIABLES[6]);
VARIABLES[21] = (0.135000*VARIABLES[6]);
VARIABLES[30] = (CONSTANTS[20]-VARIABLES[6]);
VARIABLES[22] = (0.0500000*VARIABLES[30]);
VARIABLES[19] = (pow((1.00000+(VARIABLES[20]/17.0000)),2.00000)*(1.00000+(VARIABLES[21]/26.0000)+(VARIABLES[22]/1.00000)));
VARIABLES[18] = ((0.0800000*(1.00000+((2.00000*VARIABLES[20])/17.0000)))+(0.890000*pow((VARIABLES[20]/17.0000),2.00000)));
VARIABLES[17] = (VARIABLES[18]/VARIABLES[19]);
VARIABLES[16] = (CONSTANTS[10]*VARIABLES[17]*(VARIABLES[0]-CONSTANTS[31]));
VARIABLES[23] = (CONSTANTS[11]*pow(fabs(VARIABLES[3]), 1.0 / (2.0)));
VARIABLES[24] = (0.300000*VARIABLES[2]);
VARIABLES[25] = (((CONSTANTS[12]*VARIABLES[10])-CONSTANTS[13])*pow(VARIABLES[7],2.00000));
VARIABLES[26] = ((CONSTANTS[14]/VARIABLES[7])*VARIABLES[4]*exp((CONSTANTS[15]*VARIABLES[10])));
VARIABLES[27] = (VARIABLES[26]-VARIABLES[25]);
VARIABLES[28] = (CONSTANTS[17]-VARIABLES[5]);
VARIABLES[29] = (((CONSTANTS[18]*VARIABLES[7])+CONSTANTS[19])*VARIABLES[30]);
VARIABLES[31] = (-((CONSTANTS[23]*VARIABLES[13])+(CONSTANTS[22]*(VARIABLES[7]-CONSTANTS[24]))));
VARIABLES[32] = (CONSTANTS[25]*(VARIABLES[8]-VARIABLES[7]));
VARIABLES[33] = (CONSTANTS[26]*VARIABLES[7]);
VARIABLES[34] = (VARIABLES[32]-VARIABLES[33]);
VARIABLES[80] = ((CONSTANTS[33]*VARIABLES[3]*pow(VARIABLES[24],2.00000)*pow(VARIABLES[30],2.00000))/(CONSTANTS[37]*CONSTANTS[38]*CONSTANTS[39]*CONSTANTS[40]*CONSTANTS[44]*CONSTANTS[42]*CONSTANTS[45]*CONSTANTS[46]*CONSTANTS[43]));
VARIABLES[74] = ((CONSTANTS[33]*VARIABLES[3]*pow(VARIABLES[30],2.00000))/(CONSTANTS[37]*CONSTANTS[38]*CONSTANTS[40]*CONSTANTS[45]*CONSTANTS[46]));
VARIABLES[62] = ((CONSTANTS[33]*pow(VARIABLES[30],2.00000))/(CONSTANTS[37]*CONSTANTS[40]*CONSTANTS[46]));
VARIABLES[47] = (VARIABLES[3]/CONSTANTS[38]);
VARIABLES[44] = (pow((VARIABLES[24]*VARIABLES[30]),2.00000)/(CONSTANTS[43]*CONSTANTS[39]*CONSTANTS[40]));
VARIABLES[42] = (pow(VARIABLES[24],2.00000)/CONSTANTS[39]);
VARIABLES[38] = (pow(VARIABLES[30],2.00000)/CONSTANTS[40]);
VARIABLES[40] = (CONSTANTS[34]+VARIABLES[38]);
VARIABLES[43] = (VARIABLES[40]+VARIABLES[42]);
VARIABLES[46] = (VARIABLES[43]+VARIABLES[44]);
VARIABLES[49] = (VARIABLES[46]+VARIABLES[47]);
VARIABLES[50] = ((VARIABLES[3]*pow(VARIABLES[30],2.00000))/(CONSTANTS[38]*CONSTANTS[40]*CONSTANTS[45]));
VARIABLES[52] = (VARIABLES[49]+VARIABLES[50]);
VARIABLES[53] = ((VARIABLES[3]*pow(VARIABLES[24],2.00000))/(CONSTANTS[38]*CONSTANTS[39]*CONSTANTS[44]));
VARIABLES[55] = (VARIABLES[52]+VARIABLES[53]);
VARIABLES[56] = ((VARIABLES[3]*pow(VARIABLES[24],2.00000)*pow(VARIABLES[30],2.00000))/(CONSTANTS[38]*CONSTANTS[39]*CONSTANTS[40]*CONSTANTS[44]*CONSTANTS[45]*CONSTANTS[43]));
VARIABLES[58] = (VARIABLES[55]+VARIABLES[56]);
VARIABLES[61] = (VARIABLES[58]+VARIABLES[59]);
VARIABLES[64] = (VARIABLES[61]+VARIABLES[62]);
VARIABLES[65] = ((CONSTANTS[33]*pow(VARIABLES[24],2.00000))/(CONSTANTS[37]*CONSTANTS[39]*CONSTANTS[42]));
VARIABLES[67] = (VARIABLES[64]+VARIABLES[65]);
VARIABLES[68] = ((CONSTANTS[33]*pow(VARIABLES[24],2.00000)*pow(VARIABLES[30],2.00000))/(CONSTANTS[37]*CONSTANTS[39]*CONSTANTS[40]*CONSTANTS[42]*CONSTANTS[46]*CONSTANTS[43]));
VARIABLES[70] = (VARIABLES[67]+VARIABLES[68]);
VARIABLES[71] = ((CONSTANTS[33]*VARIABLES[3])/(CONSTANTS[37]*CONSTANTS[38]));
VARIABLES[73] = (VARIABLES[70]+VARIABLES[71]);
VARIABLES[76] = (VARIABLES[73]+VARIABLES[74]);
VARIABLES[77] = ((CONSTANTS[33]*VARIABLES[3]*pow(VARIABLES[24],2.00000))/(CONSTANTS[37]*CONSTANTS[38]*CONSTANTS[39]*CONSTANTS[44]*CONSTANTS[42]));
VARIABLES[79] = (VARIABLES[76]+VARIABLES[77]);
VARIABLES[82] = (VARIABLES[79]+VARIABLES[80]);
VARIABLES[83] = VARIABLES[77];
VARIABLES[41] = (VARIABLES[39]+VARIABLES[42]);
VARIABLES[45] = (VARIABLES[41]+VARIABLES[44]);
VARIABLES[48] = VARIABLES[45];
VARIABLES[51] = VARIABLES[48];
VARIABLES[54] = (VARIABLES[51]+VARIABLES[53]);
VARIABLES[57] = (VARIABLES[54]+VARIABLES[56]);
VARIABLES[60] = VARIABLES[57];
VARIABLES[63] = VARIABLES[60];
VARIABLES[66] = (VARIABLES[63]+VARIABLES[65]);
VARIABLES[69] = (VARIABLES[66]+VARIABLES[68]);
VARIABLES[72] = VARIABLES[69];
VARIABLES[75] = VARIABLES[72];
VARIABLES[78] = VARIABLES[75];
VARIABLES[81] = (VARIABLES[78]+VARIABLES[80]);
VARIABLES[36] = (((CONSTANTS[47]*CONSTANTS[41]*VARIABLES[81])+(CONSTANTS[41]*VARIABLES[83]))/VARIABLES[82]);
VARIABLES[37] = (CONSTANTS[30]*VARIABLES[36]);
VARIABLES[86] = (CONSTANTS[56]-VARIABLES[9]);
VARIABLES[84] = ((CONSTANTS[48]/(CONSTANTS[49]+(VARIABLES[9]/VARIABLES[86])))*(VARIABLES[4]/(CONSTANTS[50]+VARIABLES[4]))*(VARIABLES[23]+CONSTANTS[51]));
VARIABLES[85] = ((CONSTANTS[52]*(VARIABLES[9]/(CONSTANTS[53]+VARIABLES[9])))/(1.00000+exp(((VARIABLES[10]-CONSTANTS[54])/CONSTANTS[55]))));
VARIABLES[87] = ((CONSTANTS[58]*(VARIABLES[9]/(CONSTANTS[59]+VARIABLES[9])))/(1.00000+exp(((VARIABLES[10]-CONSTANTS[60])/CONSTANTS[61]))));
VARIABLES[88] = (((CONSTANTS[65]*CONSTANTS[62])/(CONSTANTS[62]+VARIABLES[28]))/(1.00000+exp(((CONSTANTS[63]-VARIABLES[10])/CONSTANTS[64]))));
VARIABLES[89] = (3.00000*VARIABLES[88]);
VARIABLES[91] = ((CONSTANTS[67]*VARIABLES[10])+CONSTANTS[68]);
VARIABLES[93] = (VARIABLES[28]/VARIABLES[5]);
VARIABLES[92] = (CONSTANTS[69]*((VARIABLES[93]/(VARIABLES[93]+CONSTANTS[70]))/exp((-0.500000*CONSTANTS[71]*VARIABLES[10]))));
#ifndef VARIABLES_FOR_RATES_ONLY
#endif
}