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 108 entries in the algebraic variable array.
There are a total of 26 entries in each of the rate and state variable arrays.
There are a total of 80 entries in the constant variable array.
*/
/*
* VOI is time in component Environment (ms).
* CONSTANTS[0] is R in component Environment (J_per_moleK).
* CONSTANTS[1] is T in component Environment (kelvin).
* CONSTANTS[2] is F in component Environment (coulomb_per_mmole).
* CONSTANTS[3] is K_o in component Environment (mM).
* CONSTANTS[4] is Ca_o in component Environment (mM).
* CONSTANTS[5] is Na_o in component Environment (mM).
* CONSTANTS[75] is FonRT in component Environment (per_mV).
* STATES[0] is V in component cell (mV).
* ALGEBRAIC[49] is xik1 in component IK1 (nA_per_nF).
* ALGEBRAIC[60] is xito in component Ito (nA_per_nF).
* ALGEBRAIC[65] is xiNaK in component INaK (nA_per_nF).
* CONSTANTS[6] is wca in component cell (mV_per_uM).
* ALGEBRAIC[102] is xiNaCa in component INaCa (nA_per_nF).
* ALGEBRAIC[93] is xica in component ICaL (nA_per_nF).
* ALGEBRAIC[107] is Itotal in component cell (nA_per_nF).
* ALGEBRAIC[106] is xina in component INa (nA_per_nF).
* ALGEBRAIC[51] is xikr in component IKr (nA_per_nF).
* ALGEBRAIC[104] is xiks in component IKs (nA_per_nF).
* ALGEBRAIC[10] is i_Stim in component cell (nA_per_nF).
* CONSTANTS[7] is stim_offset in component cell (ms).
* CONSTANTS[8] is stim_period in component cell (ms).
* CONSTANTS[9] is stim_duration in component cell (ms).
* CONSTANTS[10] is stim_amplitude in component cell (nA_per_nF).
* ALGEBRAIC[3] is past in component cell (ms).
* ALGEBRAIC[105] is ena in component reversal_potentials (mV).
* STATES[1] is xm in component INa (dimensionless).
* STATES[2] is xh in component INa (dimensionless).
* STATES[3] is xj in component INa (dimensionless).
* CONSTANTS[11] is gna in component INa (uS_per_nF).
* ALGEBRAIC[0] is am in component INa (per_ms).
* ALGEBRAIC[7] is bm in component INa (per_ms).
* ALGEBRAIC[1] is ah in component INa (per_ms).
* ALGEBRAIC[8] is bh in component INa (per_ms).
* ALGEBRAIC[2] is aj in component INa (per_ms).
* ALGEBRAIC[9] is bj in component INa (per_ms).
* STATES[4] is Ca_dyad in component Ca (uM).
* ALGEBRAIC[88] is csm in component Ca (mM).
* STATES[5] is c1 in component ICaL (dimensionless).
* STATES[6] is c2 in component ICaL (dimensionless).
* STATES[7] is xi1ca in component ICaL (dimensionless).
* STATES[8] is xi1ba in component ICaL (dimensionless).
* STATES[9] is xi2ca in component ICaL (dimensionless).
* STATES[10] is xi2ba in component ICaL (dimensionless).
* CONSTANTS[12] is gca in component ICaL (mmole_per_coulomb_cm).
* CONSTANTS[13] is pca in component ICaL (cm_per_s).
* ALGEBRAIC[14] is za in component ICaL (dimensionless).
* ALGEBRAIC[18] is poinf in component ICaL (dimensionless).
* ALGEBRAIC[26] is fca in component ICaL (dimensionless).
* CONSTANTS[14] is vth in component ICaL (mV).
* CONSTANTS[15] is s6 in component ICaL (mV).
* CONSTANTS[16] is vx in component ICaL (mV).
* CONSTANTS[17] is sx in component ICaL (mV).
* CONSTANTS[18] is vy in component ICaL (mV).
* CONSTANTS[19] is sy in component ICaL (mV).
* CONSTANTS[20] is vyr in component ICaL (mV).
* CONSTANTS[21] is syr in component ICaL (mV).
* CONSTANTS[22] is cat in component ICaL (uM).
* CONSTANTS[23] is cpt in component ICaL (uM).
* ALGEBRAIC[23] is alpha in component ICaL (per_ms).
* ALGEBRAIC[25] is beta in component ICaL (per_ms).
* ALGEBRAIC[28] is k1 in component ICaL (per_ms).
* CONSTANTS[24] is k2 in component ICaL (per_ms).
* CONSTANTS[25] is k1t in component ICaL (per_ms).
* CONSTANTS[26] is k2t in component ICaL (per_ms).
* ALGEBRAIC[31] is k3 in component ICaL (per_ms).
* ALGEBRAIC[32] is k3t in component ICaL (per_ms).
* ALGEBRAIC[39] is k6 in component ICaL (per_ms).
* ALGEBRAIC[40] is k5 in component ICaL (per_ms).
* ALGEBRAIC[41] is k6t in component ICaL (per_ms).
* ALGEBRAIC[42] is k5t in component ICaL (per_ms).
* ALGEBRAIC[43] is k4 in component ICaL (per_ms).
* ALGEBRAIC[44] is k4t in component ICaL (per_ms).
* CONSTANTS[27] is r1 in component ICaL (per_ms).
* CONSTANTS[28] is r2 in component ICaL (per_ms).
* ALGEBRAIC[27] is s1 in component ICaL (per_ms).
* CONSTANTS[29] is s1t in component ICaL (per_ms).
* ALGEBRAIC[29] is s2 in component ICaL (per_ms).
* CONSTANTS[76] is s2t in component ICaL (per_ms).
* ALGEBRAIC[34] is recov in component ICaL (ms).
* CONSTANTS[30] is tca in component ICaL (ms).
* ALGEBRAIC[35] is tau_ca in component ICaL (ms).
* ALGEBRAIC[36] is tauca in component ICaL (ms).
* ALGEBRAIC[37] is tauba in component ICaL (ms).
* CONSTANTS[31] is taupo in component ICaL (ms).
* CONSTANTS[32] is tau3 in component ICaL (ms).
* ALGEBRAIC[33] is Pr in component ICaL (dimensionless).
* ALGEBRAIC[38] is Ps in component ICaL (dimensionless).
* ALGEBRAIC[30] is poi in component ICaL (dimensionless).
* ALGEBRAIC[45] is po in component ICaL (dimensionless).
* ALGEBRAIC[89] is rxa in component ICaL (mA_per_cm2).
* ALGEBRAIC[90] is jca in component ICaL (uM_per_ms).
* CONSTANTS[79] is ek in component reversal_potentials (mV).
* CONSTANTS[33] is gkix in component IK1 (uS_per_nF).
* ALGEBRAIC[46] is aki in component IK1 (per_ms).
* ALGEBRAIC[47] is bki in component IK1 (per_ms).
* ALGEBRAIC[48] is xkin in component IK1 (dimensionless).
* STATES[11] is xr in component IKr (dimensionless).
* CONSTANTS[34] is gkr in component IKr (uS_per_nF).
* ALGEBRAIC[4] is xkrv1 in component IKr (per_ms).
* ALGEBRAIC[11] is xkrv2 in component IKr (per_ms).
* ALGEBRAIC[15] is taukr in component IKr (ms).
* ALGEBRAIC[19] is xkrinf in component IKr (dimensionless).
* ALGEBRAIC[50] is rg in component IKr (dimensionless).
* STATES[12] is Ca_i in component Ca (uM).
* CONSTANTS[35] is gks in component IKs (uS_per_nF).
* STATES[13] is xs1 in component IKs (dimensionless).
* STATES[14] is xs2 in component IKs (dimensionless).
* ALGEBRAIC[103] is eks in component reversal_potentials (mV).
* ALGEBRAIC[5] is xs1ss in component IKs (dimensionless).
* ALGEBRAIC[12] is xs2ss in component IKs (dimensionless).
* ALGEBRAIC[16] is tauxs1 in component IKs (ms).
* ALGEBRAIC[20] is tauxs2 in component IKs (ms).
* ALGEBRAIC[52] is gksx in component IKs (dimensionless).
* ALGEBRAIC[56] is xitos in component Ito (nA_per_nF).
* ALGEBRAIC[58] is xitof in component Ito (nA_per_nF).
* STATES[15] is xtos in component Ito (dimensionless).
* STATES[16] is ytos in component Ito (dimensionless).
* STATES[17] is xtof in component Ito (dimensionless).
* STATES[18] is ytof in component Ito (dimensionless).
* CONSTANTS[36] is gtos in component Ito (uS_per_nF).
* CONSTANTS[37] is gtof in component Ito (uS_per_nF).
* ALGEBRAIC[6] is rt1 in component Ito (dimensionless).
* ALGEBRAIC[53] is rt2 in component Ito (dimensionless).
* ALGEBRAIC[55] is rt3 in component Ito (dimensionless).
* ALGEBRAIC[13] is rt4 in component Ito (dimensionless).
* ALGEBRAIC[57] is rt5 in component Ito (dimensionless).
* ALGEBRAIC[17] is xtos_inf in component Ito (dimensionless).
* ALGEBRAIC[59] is ytos_inf in component Ito (dimensionless).
* ALGEBRAIC[21] is xtof_inf in component Ito (dimensionless).
* ALGEBRAIC[61] is ytof_inf in component Ito (dimensionless).
* ALGEBRAIC[54] is rs_inf in component Ito (dimensionless).
* ALGEBRAIC[22] is txs in component Ito (ms).
* ALGEBRAIC[62] is tys in component Ito (ms).
* ALGEBRAIC[24] is txf in component Ito (ms).
* ALGEBRAIC[64] is tyf in component Ito (ms).
* STATES[19] is Na_i in component Na (mM).
* CONSTANTS[38] is gNaK in component INaK (nA_per_nF).
* CONSTANTS[39] is xkmko in component INaK (mM).
* CONSTANTS[40] is xkmnai in component INaK (mM).
* CONSTANTS[77] is sigma in component INaK (dimensionless).
* ALGEBRAIC[63] is fNaK in component INaK (dimensionless).
* STATES[20] is Ca_submem in component Ca (uM).
* CONSTANTS[41] is gNaCa in component INaCa (uM_per_ms).
* ALGEBRAIC[67] is aloss in component INaCa (dimensionless).
* ALGEBRAIC[97] is yz1 in component INaCa (mM4).
* ALGEBRAIC[98] is yz2 in component INaCa (mM4).
* ALGEBRAIC[68] is yz3 in component INaCa (mM4).
* ALGEBRAIC[99] is yz4 in component INaCa (mM4).
* ALGEBRAIC[95] is zw3 in component INaCa (mM4).
* ALGEBRAIC[66] is zw4 in component INaCa (dimensionless).
* ALGEBRAIC[100] is zw8 in component INaCa (mM4).
* ALGEBRAIC[101] is jNaCa in component INaCa (uM_per_ms).
* CONSTANTS[42] is xkdna in component INaCa (uM).
* CONSTANTS[43] is xmcao in component INaCa (mM).
* CONSTANTS[44] is xmnao in component INaCa (mM).
* CONSTANTS[45] is xmnai in component INaCa (mM).
* CONSTANTS[46] is xmcai in component INaCa (mM).
* STATES[21] is Ca_NSR in component Ca (uM).
* ALGEBRAIC[87] is dCa_JSR in component Ca (uM_per_ms).
* CONSTANTS[47] is cstar in component Irel (uM).
* STATES[22] is Ca_JSR in component Irel (uM).
* CONSTANTS[48] is gryr in component Irel (per_ms).
* CONSTANTS[49] is gbarsr in component Irel (dimensionless).
* CONSTANTS[50] is gdyad in component Irel (mmole_per_coulomb_cm).
* CONSTANTS[51] is ax in component Irel (per_mV).
* CONSTANTS[52] is ay in component Irel (per_mV).
* CONSTANTS[53] is av in component Irel (per_ms).
* CONSTANTS[78] is bv in component Irel (uM_per_ms).
* ALGEBRAIC[69] is Qr0 in component Irel (uM_per_ms).
* ALGEBRAIC[70] is Qr in component Irel (uM_per_ms).
* ALGEBRAIC[71] is sparkV in component Irel (dimensionless).
* ALGEBRAIC[91] is spark_rate in component Irel (per_ms).
* CONSTANTS[54] is taua in component Irel (ms).
* CONSTANTS[55] is taur in component Irel (ms).
* ALGEBRAIC[92] is xirp in component Irel (uM_per_ms).
* ALGEBRAIC[94] is xicap in component Irel (uM_per_ms).
* ALGEBRAIC[96] is xiryr in component Irel (uM_per_ms).
* STATES[23] is xir in component Irel (uM_per_ms).
* ALGEBRAIC[72] is jup in component Ileak_Iup_Ixfer (uM_per_ms).
* ALGEBRAIC[73] is jleak in component Ileak_Iup_Ixfer (uM_per_ms).
* CONSTANTS[56] is cup in component Ileak_Iup_Ixfer (uM).
* CONSTANTS[57] is kj in component Ileak_Iup_Ixfer (uM).
* CONSTANTS[58] is vup in component Ileak_Iup_Ixfer (uM_per_ms).
* CONSTANTS[59] is gleak in component Ileak_Iup_Ixfer (per_ms).
* CONSTANTS[60] is bcal in component Ca (uM).
* CONSTANTS[61] is xkcal in component Ca (uM).
* CONSTANTS[62] is srmax in component Ca (uM).
* CONSTANTS[63] is srkd in component Ca (uM).
* CONSTANTS[64] is bmem in component Ca (uM).
* CONSTANTS[65] is kmem in component Ca (uM).
* CONSTANTS[66] is bsar in component Ca (uM).
* CONSTANTS[67] is ksar in component Ca (uM).
* ALGEBRAIC[74] is bpxs in component Ca (dimensionless).
* ALGEBRAIC[75] is spxs in component Ca (dimensionless).
* ALGEBRAIC[76] is mempxs in component Ca (dimensionless).
* ALGEBRAIC[77] is sarpxs in component Ca (dimensionless).
* ALGEBRAIC[78] is dcsib in component Ca (dimensionless).
* ALGEBRAIC[79] is bpxi in component Ca (dimensionless).
* ALGEBRAIC[80] is spxi in component Ca (dimensionless).
* ALGEBRAIC[81] is mempxi in component Ca (dimensionless).
* ALGEBRAIC[82] is sarpxi in component Ca (dimensionless).
* ALGEBRAIC[83] is dciib in component Ca (dimensionless).
* CONSTANTS[68] is xkon in component Ca (per_uM_per_ms).
* CONSTANTS[69] is xkoff in component Ca (per_ms).
* CONSTANTS[70] is btrop in component Ca (uM).
* ALGEBRAIC[86] is xbi in component Ca (uM_per_ms).
* ALGEBRAIC[85] is xbs in component Ca (uM_per_ms).
* STATES[24] is tropi in component Ca (uM).
* STATES[25] is trops in component Ca (uM).
* CONSTANTS[71] is taud in component Ca (ms).
* CONSTANTS[72] is taups in component Ca (ms).
* ALGEBRAIC[84] is jd in component Ca (uM_per_ms).
* CONSTANTS[73] is K_i in component reversal_potentials (mM).
* CONSTANTS[74] is prNaK in component reversal_potentials (dimensionless).
* RATES[0] is d/dt V in component cell (mV).
* RATES[2] is d/dt xh in component INa (dimensionless).
* RATES[3] is d/dt xj in component INa (dimensionless).
* RATES[1] is d/dt xm in component INa (dimensionless).
* RATES[5] is d/dt c1 in component ICaL (dimensionless).
* RATES[6] is d/dt c2 in component ICaL (dimensionless).
* RATES[7] is d/dt xi1ca in component ICaL (dimensionless).
* RATES[8] is d/dt xi1ba in component ICaL (dimensionless).
* RATES[9] is d/dt xi2ca in component ICaL (dimensionless).
* RATES[10] is d/dt xi2ba in component ICaL (dimensionless).
* RATES[11] is d/dt xr in component IKr (dimensionless).
* RATES[13] is d/dt xs1 in component IKs (dimensionless).
* RATES[14] is d/dt xs2 in component IKs (dimensionless).
* RATES[15] is d/dt xtos in component Ito (dimensionless).
* RATES[16] is d/dt ytos in component Ito (dimensionless).
* RATES[17] is d/dt xtof in component Ito (dimensionless).
* RATES[18] is d/dt ytof in component Ito (dimensionless).
* RATES[22] is d/dt Ca_JSR in component Irel (uM).
* RATES[23] is d/dt xir in component Irel (uM_per_ms).
* RATES[19] is d/dt Na_i in component Na (mM).
* RATES[4] is d/dt Ca_dyad in component Ca (uM).
* RATES[20] is d/dt Ca_submem in component Ca (uM).
* RATES[12] is d/dt Ca_i in component Ca (uM).
* RATES[21] is d/dt Ca_NSR in component Ca (uM).
* RATES[24] is d/dt tropi in component Ca (uM).
* RATES[25] is d/dt trops in component Ca (uM).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 8.314472;
CONSTANTS[1] = 308;
CONSTANTS[2] = 96.4853415;
CONSTANTS[3] = 5.4;
CONSTANTS[4] = 1.8;
CONSTANTS[5] = 136;
STATES[0] = -87.169816169406;
CONSTANTS[6] = 8;
CONSTANTS[7] = 0;
CONSTANTS[8] = 400;
CONSTANTS[9] = 3;
CONSTANTS[10] = -15;
STATES[1] = 0.001075453357;
STATES[2] = 0.990691306716;
STATES[3] = 0.993888937283;
CONSTANTS[11] = 12;
STATES[4] = 1.716573130685;
STATES[5] = 0.000018211252;
STATES[6] = 0.979322592773;
STATES[7] = 0.001208153482;
STATES[8] = 0.000033616596;
STATES[9] = 0.004173008466;
STATES[10] = 0.015242594688;
CONSTANTS[12] = 182;
CONSTANTS[13] = 0.00054;
CONSTANTS[14] = 0;
CONSTANTS[15] = 8;
CONSTANTS[16] = -40;
CONSTANTS[17] = 3;
CONSTANTS[18] = -40;
CONSTANTS[19] = 4;
CONSTANTS[20] = -40;
CONSTANTS[21] = 11.32;
CONSTANTS[22] = 3;
CONSTANTS[23] = 6.09365;
CONSTANTS[24] = 1.03615e-4;
CONSTANTS[25] = 0.00413;
CONSTANTS[26] = 0.00224;
CONSTANTS[27] = 0.3;
CONSTANTS[28] = 3;
CONSTANTS[29] = 0.00195;
CONSTANTS[30] = 78.0329;
CONSTANTS[31] = 1;
CONSTANTS[32] = 3;
CONSTANTS[33] = 0.3;
STATES[11] = 0.007074239331;
CONSTANTS[34] = 0.0125;
STATES[12] = 0.256752008084;
CONSTANTS[35] = 0.1386;
STATES[13] = 0.048267587131;
STATES[14] = 0.105468807033;
STATES[15] = 0.00364776906;
STATES[16] = 0.174403618112;
STATES[17] = 0.003643592594;
STATES[18] = 0.993331326442;
CONSTANTS[36] = 0.04;
CONSTANTS[37] = 0.11;
STATES[19] = 11.441712311614;
CONSTANTS[38] = 1.5;
CONSTANTS[39] = 1.5;
CONSTANTS[40] = 12;
STATES[20] = 0.226941113355;
CONSTANTS[41] = 0.84;
CONSTANTS[42] = 0.3;
CONSTANTS[43] = 1.3;
CONSTANTS[44] = 87.5;
CONSTANTS[45] = 12.3;
CONSTANTS[46] = 0.0036;
STATES[21] = 104.450004990523;
CONSTANTS[47] = 90;
STATES[22] = 97.505463697266;
CONSTANTS[48] = 2.58079;
CONSTANTS[49] = 26841.8;
CONSTANTS[50] = 9000;
CONSTANTS[51] = 0.3576;
CONSTANTS[52] = 0.05;
CONSTANTS[53] = 11.3;
CONSTANTS[54] = 100;
CONSTANTS[55] = 30;
STATES[23] = 0.006679257264;
CONSTANTS[56] = 0.5;
CONSTANTS[57] = 50;
CONSTANTS[58] = 0.4;
CONSTANTS[59] = 0.00002069;
CONSTANTS[60] = 24;
CONSTANTS[61] = 7;
CONSTANTS[62] = 47;
CONSTANTS[63] = 0.6;
CONSTANTS[64] = 15;
CONSTANTS[65] = 0.3;
CONSTANTS[66] = 42;
CONSTANTS[67] = 13;
CONSTANTS[68] = 0.0327;
CONSTANTS[69] = 0.0196;
CONSTANTS[70] = 70;
STATES[24] = 22.171689894953;
STATES[25] = 19.864701949854;
CONSTANTS[71] = 4;
CONSTANTS[72] = 0.5;
CONSTANTS[73] = 140;
CONSTANTS[74] = 0.01833;
CONSTANTS[75] = CONSTANTS[2]/( CONSTANTS[0]*CONSTANTS[1]);
CONSTANTS[76] = ( (( CONSTANTS[29]*CONSTANTS[27])/CONSTANTS[28])*CONSTANTS[26])/CONSTANTS[25];
CONSTANTS[77] = (exp(CONSTANTS[5]/67.3000) - 1.00000)/7.00000;
CONSTANTS[78] = (1.00000 - CONSTANTS[53])*CONSTANTS[47] - 50.0000;
CONSTANTS[79] = (1.00000/CONSTANTS[75])*log(CONSTANTS[3]/CONSTANTS[73]);
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[22] = (STATES[21] - STATES[22])/CONSTANTS[54];
ALGEBRAIC[1] = (STATES[0]<- 40.0000 ? 0.135000*exp((80.0000+STATES[0])/- 6.80000) : 0.00000);
ALGEBRAIC[8] = (STATES[0]<- 40.0000 ? 3.56000*exp( 0.0790000*STATES[0])+ 310000.*exp( 0.350000*STATES[0]) : 1.00000/( 0.130000*(1.00000+exp((STATES[0]+10.6600)/- 11.1000))));
RATES[2] = ALGEBRAIC[1]*(1.00000 - STATES[2]) - ALGEBRAIC[8]*STATES[2];
ALGEBRAIC[2] = (STATES[0]<- 40.0000 ? ( ( - 127140.*exp( 0.244400*STATES[0]) - 3.47400e-05*exp( - 0.0439100*STATES[0]))*1.00000*(STATES[0]+37.7800))/(1.00000+exp( 0.311000*(STATES[0]+79.2300))) : 0.00000);
ALGEBRAIC[9] = (STATES[0]<- 40.0000 ? ( 0.121200*exp( - 0.0105200*STATES[0]))/(1.00000+exp( - 0.137800*(STATES[0]+40.1400))) : ( 0.300000*exp( - 2.53500e-07*STATES[0]))/(1.00000+exp( - 0.100000*(STATES[0]+32.0000))));
RATES[3] = ALGEBRAIC[2]*(1.00000 - STATES[3]) - ALGEBRAIC[9]*STATES[3];
ALGEBRAIC[0] = (fabs(STATES[0]+47.1300)>0.00100000 ? ( 0.320000*1.00000*(STATES[0]+47.1300))/(1.00000 - exp( - 0.100000*(STATES[0]+47.1300))) : 3.20000);
ALGEBRAIC[7] = 0.0800000*exp(- STATES[0]/11.0000);
RATES[1] = ALGEBRAIC[0]*(1.00000 - STATES[1]) - ALGEBRAIC[7]*STATES[1];
ALGEBRAIC[5] = 1.00000/(1.00000+exp(- (STATES[0] - 1.50000)/16.7000));
ALGEBRAIC[16] = (fabs(STATES[0]+30.0000)<0.00100000/0.0687000 ? 1.00000/(7.19000e-05/0.148000+0.000131000/0.0687000) : 1.00000/(( 7.19000e-05*(STATES[0]+30.0000))/(1.00000 - exp( - 0.148000*(STATES[0]+30.0000)))+( 0.000131000*(STATES[0]+30.0000))/(exp( 0.0687000*(STATES[0]+30.0000)) - 1.00000)));
RATES[13] = (ALGEBRAIC[5] - STATES[13])/ALGEBRAIC[16];
ALGEBRAIC[4] = (fabs(STATES[0]+7.00000)>0.00100000 ? ( 0.00138000*1.00000*(STATES[0]+7.00000))/(1.00000 - exp( - 0.123000*(STATES[0]+7.00000))) : 0.00138000/0.123000);
ALGEBRAIC[11] = (fabs(STATES[0]+10.0000)>0.00100000 ? ( 0.000610000*1.00000*(STATES[0]+10.0000))/(exp( 0.145000*(STATES[0]+10.0000)) - 1.00000) : 0.000610000/0.145000);
ALGEBRAIC[15] = 1.00000/(ALGEBRAIC[4]+ALGEBRAIC[11]);
ALGEBRAIC[19] = 1.00000/(1.00000+exp(- (STATES[0]+50.0000)/7.50000));
RATES[11] = (ALGEBRAIC[19] - STATES[11])/ALGEBRAIC[15];
ALGEBRAIC[12] = ALGEBRAIC[5];
ALGEBRAIC[20] = 4.00000*ALGEBRAIC[16];
RATES[14] = (ALGEBRAIC[12] - STATES[14])/ALGEBRAIC[20];
ALGEBRAIC[6] = - (STATES[0]+3.00000)/15.0000;
ALGEBRAIC[17] = 1.00000/(1.00000+exp(ALGEBRAIC[6]));
ALGEBRAIC[22] = 9.00000/(1.00000+exp(- ALGEBRAIC[6]))+0.500000;
RATES[15] = (ALGEBRAIC[17] - STATES[15])/ALGEBRAIC[22];
ALGEBRAIC[21] = ALGEBRAIC[17];
ALGEBRAIC[13] = ( (- STATES[0]/30.0000)*STATES[0])/30.0000;
ALGEBRAIC[24] = 3.50000*exp(ALGEBRAIC[13])+1.50000;
RATES[17] = (ALGEBRAIC[21] - STATES[17])/ALGEBRAIC[24];
ALGEBRAIC[18] = 1.00000/(1.00000+exp(- (STATES[0] - CONSTANTS[14])/CONSTANTS[15]));
ALGEBRAIC[23] = ALGEBRAIC[18]/CONSTANTS[31];
ALGEBRAIC[25] = (1.00000 - ALGEBRAIC[18])/CONSTANTS[31];
ALGEBRAIC[26] = 1.00000/(1.00000+pow(CONSTANTS[22]/STATES[4], 3.00000));
ALGEBRAIC[34] = 10.0000+ 4954.00*exp(STATES[0]/15.6000);
ALGEBRAIC[35] = CONSTANTS[30]/(1.00000+pow(STATES[4]/CONSTANTS[23], 4.00000))+0.100000;
ALGEBRAIC[33] = 1.00000 - 1.00000/(1.00000+exp(- (STATES[0] - CONSTANTS[18])/CONSTANTS[19]));
ALGEBRAIC[36] = (ALGEBRAIC[34] - ALGEBRAIC[35])*ALGEBRAIC[33]+ALGEBRAIC[35];
ALGEBRAIC[38] = 1.00000/(1.00000+exp(- (STATES[0] - CONSTANTS[20])/CONSTANTS[21]));
ALGEBRAIC[39] = ( ALGEBRAIC[26]*ALGEBRAIC[38])/ALGEBRAIC[36];
ALGEBRAIC[40] = (1.00000 - ALGEBRAIC[38])/ALGEBRAIC[36];
ALGEBRAIC[37] = (ALGEBRAIC[34] - 450.000)*ALGEBRAIC[33]+450.000;
ALGEBRAIC[41] = ALGEBRAIC[38]/ALGEBRAIC[37];
ALGEBRAIC[42] = (1.00000 - ALGEBRAIC[38])/ALGEBRAIC[37];
RATES[6] = ( ALGEBRAIC[25]*STATES[5]+ ALGEBRAIC[40]*STATES[9]+ ALGEBRAIC[42]*STATES[10]) - (ALGEBRAIC[39]+ALGEBRAIC[41]+ALGEBRAIC[23])*STATES[6];
ALGEBRAIC[30] = 1.00000/(1.00000+exp(- (STATES[0] - CONSTANTS[16])/CONSTANTS[17]));
ALGEBRAIC[31] = (1.00000 - ALGEBRAIC[30])/CONSTANTS[32];
ALGEBRAIC[28] = 0.0241680*ALGEBRAIC[26];
ALGEBRAIC[43] = ( (( (( ALGEBRAIC[31]*ALGEBRAIC[23])/ALGEBRAIC[25])*ALGEBRAIC[28])/CONSTANTS[24])*ALGEBRAIC[40])/ALGEBRAIC[39];
RATES[9] = ( ALGEBRAIC[31]*STATES[7]+ ALGEBRAIC[39]*STATES[6]) - (ALGEBRAIC[40]+ALGEBRAIC[43])*STATES[9];
ALGEBRAIC[32] = ALGEBRAIC[31];
ALGEBRAIC[44] = ( (( (( ALGEBRAIC[32]*ALGEBRAIC[23])/ALGEBRAIC[25])*CONSTANTS[25])/CONSTANTS[26])*ALGEBRAIC[42])/ALGEBRAIC[41];
RATES[10] = ( ALGEBRAIC[32]*STATES[8]+ ALGEBRAIC[41]*STATES[6]) - (ALGEBRAIC[42]+ALGEBRAIC[44])*STATES[10];
ALGEBRAIC[45] = (((((1.00000 - STATES[7]) - STATES[9]) - STATES[8]) - STATES[10]) - STATES[5]) - STATES[6];
RATES[5] = ( ALGEBRAIC[23]*STATES[6]+ CONSTANTS[24]*STATES[7]+ CONSTANTS[26]*STATES[8]+ CONSTANTS[28]*ALGEBRAIC[45]) - (ALGEBRAIC[25]+CONSTANTS[27]+CONSTANTS[25]+ALGEBRAIC[28])*STATES[5];
ALGEBRAIC[27] = 0.0182688*ALGEBRAIC[26];
ALGEBRAIC[29] = ( (( ALGEBRAIC[27]*CONSTANTS[27])/CONSTANTS[28])*CONSTANTS[24])/ALGEBRAIC[28];
RATES[7] = ( ALGEBRAIC[28]*STATES[5]+ ALGEBRAIC[43]*STATES[9]+ ALGEBRAIC[27]*ALGEBRAIC[45]) - (ALGEBRAIC[31]+CONSTANTS[24]+ALGEBRAIC[29])*STATES[7];
RATES[8] = ( CONSTANTS[25]*STATES[5]+ ALGEBRAIC[44]*STATES[10]+ CONSTANTS[29]*ALGEBRAIC[45]) - (ALGEBRAIC[32]+CONSTANTS[26]+CONSTANTS[76])*STATES[8];
ALGEBRAIC[53] = (STATES[0]+33.5000)/10.0000;
ALGEBRAIC[59] = 1.00000/(1.00000+exp(ALGEBRAIC[53]));
ALGEBRAIC[55] = (STATES[0]+60.0000)/10.0000;
ALGEBRAIC[62] = 3000.00/(1.00000+exp(ALGEBRAIC[55]))+30.0000;
RATES[16] = (ALGEBRAIC[59] - STATES[16])/ALGEBRAIC[62];
ALGEBRAIC[61] = ALGEBRAIC[59];
ALGEBRAIC[57] = (STATES[0]+33.5000)/10.0000;
ALGEBRAIC[64] = 20.0000/(1.00000+exp(ALGEBRAIC[57]))+20.0000;
RATES[18] = (ALGEBRAIC[61] - STATES[18])/ALGEBRAIC[64];
ALGEBRAIC[72] = ( CONSTANTS[58]*STATES[12]*STATES[12])/( STATES[12]*STATES[12]+ CONSTANTS[56]*CONSTANTS[56]);
ALGEBRAIC[73] = (( CONSTANTS[59]*STATES[21]*STATES[21])/( STATES[21]*STATES[21]+ CONSTANTS[57]*CONSTANTS[57]))*( STATES[21]*16.6670 - STATES[12]);
ALGEBRAIC[79] = ( CONSTANTS[60]*CONSTANTS[61])/( (CONSTANTS[61]+STATES[12])*(CONSTANTS[61]+STATES[12]));
ALGEBRAIC[80] = ( CONSTANTS[62]*CONSTANTS[63])/( (CONSTANTS[63]+STATES[12])*(CONSTANTS[63]+STATES[12]));
ALGEBRAIC[81] = ( CONSTANTS[64]*CONSTANTS[65])/( (CONSTANTS[65]+STATES[12])*(CONSTANTS[65]+STATES[12]));
ALGEBRAIC[82] = ( CONSTANTS[66]*CONSTANTS[67])/( (CONSTANTS[67]+STATES[12])*(CONSTANTS[67]+STATES[12]));
ALGEBRAIC[83] = 1.00000/(1.00000+ALGEBRAIC[79]+ALGEBRAIC[80]+ALGEBRAIC[81]+ALGEBRAIC[82]);
ALGEBRAIC[86] = CONSTANTS[68]*STATES[12]*(CONSTANTS[70] - STATES[24]) - CONSTANTS[69]*STATES[24];
ALGEBRAIC[84] = (STATES[20] - STATES[12])/CONSTANTS[71];
RATES[12] = ALGEBRAIC[83]*(((ALGEBRAIC[84] - ALGEBRAIC[72])+ALGEBRAIC[73]) - ALGEBRAIC[86]);
RATES[24] = ALGEBRAIC[86];
ALGEBRAIC[85] = CONSTANTS[68]*STATES[20]*(CONSTANTS[70] - STATES[25]) - CONSTANTS[69]*STATES[25];
RATES[25] = ALGEBRAIC[85];
ALGEBRAIC[87] = (- STATES[23]+ALGEBRAIC[72]) - ALGEBRAIC[73];
RATES[21] = ALGEBRAIC[87];
ALGEBRAIC[69] = (STATES[22]>50.0000&&STATES[22]<CONSTANTS[47] ? (STATES[22] - 50.0000)/1.00000 : STATES[22]>=CONSTANTS[47] ? CONSTANTS[53]*STATES[22]+CONSTANTS[78] : 0.00000);
ALGEBRAIC[70] = ( STATES[21]*ALGEBRAIC[69])/CONSTANTS[47];
ALGEBRAIC[88] = STATES[20]/1000.00;
ALGEBRAIC[14] = STATES[0]*2.00000*CONSTANTS[75];
ALGEBRAIC[89] = (fabs(ALGEBRAIC[14])<0.00100000 ? ( 4.00000*CONSTANTS[13]*CONSTANTS[2]*CONSTANTS[75]*( ALGEBRAIC[88]*exp(ALGEBRAIC[14]) - 0.341000*CONSTANTS[4]))/( 2.00000*CONSTANTS[75]) : ( 4.00000*CONSTANTS[13]*STATES[0]*CONSTANTS[2]*CONSTANTS[75]*( ALGEBRAIC[88]*exp(ALGEBRAIC[14]) - 0.341000*CONSTANTS[4]))/(exp(ALGEBRAIC[14]) - 1.00000));
ALGEBRAIC[71] = exp( - CONSTANTS[52]*(STATES[0]+30.0000))/(1.00000+exp( - CONSTANTS[52]*(STATES[0]+30.0000)));
ALGEBRAIC[91] = (CONSTANTS[48]/1.00000)*ALGEBRAIC[45]*fabs(ALGEBRAIC[89])*ALGEBRAIC[71];
RATES[23] = ALGEBRAIC[91]*ALGEBRAIC[70] - ( STATES[23]*(1.00000 - ( CONSTANTS[55]*ALGEBRAIC[87])/STATES[21]))/CONSTANTS[55];
ALGEBRAIC[92] = ( (( ALGEBRAIC[45]*ALGEBRAIC[70]*fabs(ALGEBRAIC[89])*CONSTANTS[49])/1.00000)*exp( - CONSTANTS[51]*(STATES[0]+30.0000)))/(1.00000+exp( - CONSTANTS[51]*(STATES[0]+30.0000)));
ALGEBRAIC[94] = ALGEBRAIC[45]*CONSTANTS[50]*fabs(ALGEBRAIC[89]);
ALGEBRAIC[96] = ALGEBRAIC[92]+ALGEBRAIC[94];
RATES[4] = ALGEBRAIC[96] - (STATES[4] - STATES[20])/CONSTANTS[72];
ALGEBRAIC[90] = CONSTANTS[12]*ALGEBRAIC[45]*ALGEBRAIC[89];
ALGEBRAIC[67] = 1.00000/(1.00000+pow(CONSTANTS[42]/STATES[20], 3.00000));
ALGEBRAIC[95] = pow(STATES[19], 3.00000)*CONSTANTS[4]*exp( STATES[0]*0.350000*CONSTANTS[75]) - pow(CONSTANTS[5], 3.00000)*ALGEBRAIC[88]*exp( STATES[0]*(0.350000 - 1.00000)*CONSTANTS[75]);
ALGEBRAIC[66] = 1.00000+ 0.200000*exp( STATES[0]*(0.350000 - 1.00000)*CONSTANTS[75]);
ALGEBRAIC[97] = CONSTANTS[43]*pow(STATES[19], 3.00000)+ pow(CONSTANTS[44], 3.00000)*ALGEBRAIC[88];
ALGEBRAIC[98] = pow(CONSTANTS[45], 3.00000)*CONSTANTS[4]*(1.00000+ALGEBRAIC[88]/CONSTANTS[46]);
ALGEBRAIC[68] = CONSTANTS[46]*pow(CONSTANTS[5], 3.00000)*(1.00000+pow(STATES[19]/CONSTANTS[45], 3.00000));
ALGEBRAIC[99] = pow(STATES[19], 3.00000)*CONSTANTS[4]+ pow(CONSTANTS[5], 3.00000)*ALGEBRAIC[88];
ALGEBRAIC[100] = ALGEBRAIC[97]+ALGEBRAIC[98]+ALGEBRAIC[68]+ALGEBRAIC[99];
ALGEBRAIC[101] = ( CONSTANTS[41]*ALGEBRAIC[67]*ALGEBRAIC[95])/( ALGEBRAIC[66]*ALGEBRAIC[100]);
ALGEBRAIC[74] = ( CONSTANTS[60]*CONSTANTS[61])/( (CONSTANTS[61]+STATES[20])*(CONSTANTS[61]+STATES[20]));
ALGEBRAIC[75] = ( CONSTANTS[62]*CONSTANTS[63])/( (CONSTANTS[63]+STATES[20])*(CONSTANTS[63]+STATES[20]));
ALGEBRAIC[76] = ( CONSTANTS[64]*CONSTANTS[65])/( (CONSTANTS[65]+STATES[20])*(CONSTANTS[65]+STATES[20]));
ALGEBRAIC[77] = ( CONSTANTS[66]*CONSTANTS[67])/( (CONSTANTS[67]+STATES[20])*(CONSTANTS[67]+STATES[20]));
ALGEBRAIC[78] = 1.00000/(1.00000+ALGEBRAIC[74]+ALGEBRAIC[75]+ALGEBRAIC[76]+ALGEBRAIC[77]);
RATES[20] = ALGEBRAIC[78]*( 50.0000*(((STATES[23] - ALGEBRAIC[84]) - ALGEBRAIC[90])+ALGEBRAIC[101]) - ALGEBRAIC[85]);
ALGEBRAIC[63] = 1.00000/(1.00000+ 0.124500*exp( - 0.100000*STATES[0]*CONSTANTS[75])+ 0.0365000*CONSTANTS[77]*exp( - STATES[0]*CONSTANTS[75]));
ALGEBRAIC[65] = ( (( CONSTANTS[38]*ALGEBRAIC[63]*STATES[19])/(STATES[19]+CONSTANTS[40]))*CONSTANTS[3])/(CONSTANTS[3]+CONSTANTS[39]);
ALGEBRAIC[102] = CONSTANTS[6]*ALGEBRAIC[101];
ALGEBRAIC[105] = (1.00000/CONSTANTS[75])*log(CONSTANTS[5]/STATES[19]);
ALGEBRAIC[106] = CONSTANTS[11]*STATES[2]*STATES[3]*STATES[1]*STATES[1]*STATES[1]*(STATES[0] - ALGEBRAIC[105]);
RATES[19] = - (ALGEBRAIC[106]+ 3.00000*ALGEBRAIC[65]+ 3.00000*ALGEBRAIC[102])/( CONSTANTS[6]*1000.00);
ALGEBRAIC[46] = 1.02000/(1.00000+exp( 0.238500*((STATES[0] - CONSTANTS[79]) - 59.2150)));
ALGEBRAIC[47] = ( 0.491240*exp( 0.0803200*((STATES[0] - CONSTANTS[79])+5.47600))+ 1.00000*exp( 0.0617500*((STATES[0] - CONSTANTS[79]) - 594.310)))/(1.00000+exp( - 0.514300*((STATES[0] - CONSTANTS[79])+4.75300)));
ALGEBRAIC[48] = ALGEBRAIC[46]/(ALGEBRAIC[46]+ALGEBRAIC[47]);
ALGEBRAIC[49] = CONSTANTS[33]* pow((CONSTANTS[3]/5.40000), 1.0 / 2)*ALGEBRAIC[48]*(STATES[0] - CONSTANTS[79]);
ALGEBRAIC[54] = 1.00000/(1.00000+exp(ALGEBRAIC[53]));
ALGEBRAIC[56] = CONSTANTS[36]*STATES[15]*(STATES[16]+ 0.500000*ALGEBRAIC[54])*(STATES[0] - CONSTANTS[79]);
ALGEBRAIC[58] = CONSTANTS[37]*STATES[17]*STATES[18]*(STATES[0] - CONSTANTS[79]);
ALGEBRAIC[60] = ALGEBRAIC[56]+ALGEBRAIC[58];
ALGEBRAIC[93] = 2.00000*CONSTANTS[6]*ALGEBRAIC[90];
ALGEBRAIC[50] = 1.00000/(1.00000+exp((STATES[0]+33.0000)/22.4000));
ALGEBRAIC[51] = CONSTANTS[34]* pow((CONSTANTS[3]/5.40000), 1.0 / 2)*STATES[11]*ALGEBRAIC[50]*(STATES[0] - CONSTANTS[79]);
ALGEBRAIC[103] = (1.00000/CONSTANTS[75])*log((CONSTANTS[3]+ CONSTANTS[74]*CONSTANTS[5])/(CONSTANTS[73]+ CONSTANTS[74]*STATES[19]));
ALGEBRAIC[52] = 1.00000+0.800000/(1.00000+pow(0.500000/STATES[12], 3.00000));
ALGEBRAIC[104] = CONSTANTS[35]*ALGEBRAIC[52]*STATES[13]*STATES[14]*(STATES[0] - ALGEBRAIC[103]);
ALGEBRAIC[3] = floor(VOI/CONSTANTS[8])*CONSTANTS[8];
ALGEBRAIC[10] = (VOI - ALGEBRAIC[3]>=CONSTANTS[7]&&VOI - ALGEBRAIC[3]<=CONSTANTS[7]+CONSTANTS[9] ? CONSTANTS[10] : 0.00000);
ALGEBRAIC[107] = - (ALGEBRAIC[106]+ALGEBRAIC[49]+ALGEBRAIC[51]+ALGEBRAIC[104]+ALGEBRAIC[60]+ALGEBRAIC[102]+ALGEBRAIC[93]+ALGEBRAIC[65]+ALGEBRAIC[10]);
RATES[0] = ALGEBRAIC[107];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = (STATES[0]<- 40.0000 ? 0.135000*exp((80.0000+STATES[0])/- 6.80000) : 0.00000);
ALGEBRAIC[8] = (STATES[0]<- 40.0000 ? 3.56000*exp( 0.0790000*STATES[0])+ 310000.*exp( 0.350000*STATES[0]) : 1.00000/( 0.130000*(1.00000+exp((STATES[0]+10.6600)/- 11.1000))));
ALGEBRAIC[2] = (STATES[0]<- 40.0000 ? ( ( - 127140.*exp( 0.244400*STATES[0]) - 3.47400e-05*exp( - 0.0439100*STATES[0]))*1.00000*(STATES[0]+37.7800))/(1.00000+exp( 0.311000*(STATES[0]+79.2300))) : 0.00000);
ALGEBRAIC[9] = (STATES[0]<- 40.0000 ? ( 0.121200*exp( - 0.0105200*STATES[0]))/(1.00000+exp( - 0.137800*(STATES[0]+40.1400))) : ( 0.300000*exp( - 2.53500e-07*STATES[0]))/(1.00000+exp( - 0.100000*(STATES[0]+32.0000))));
ALGEBRAIC[0] = (fabs(STATES[0]+47.1300)>0.00100000 ? ( 0.320000*1.00000*(STATES[0]+47.1300))/(1.00000 - exp( - 0.100000*(STATES[0]+47.1300))) : 3.20000);
ALGEBRAIC[7] = 0.0800000*exp(- STATES[0]/11.0000);
ALGEBRAIC[5] = 1.00000/(1.00000+exp(- (STATES[0] - 1.50000)/16.7000));
ALGEBRAIC[16] = (fabs(STATES[0]+30.0000)<0.00100000/0.0687000 ? 1.00000/(7.19000e-05/0.148000+0.000131000/0.0687000) : 1.00000/(( 7.19000e-05*(STATES[0]+30.0000))/(1.00000 - exp( - 0.148000*(STATES[0]+30.0000)))+( 0.000131000*(STATES[0]+30.0000))/(exp( 0.0687000*(STATES[0]+30.0000)) - 1.00000)));
ALGEBRAIC[4] = (fabs(STATES[0]+7.00000)>0.00100000 ? ( 0.00138000*1.00000*(STATES[0]+7.00000))/(1.00000 - exp( - 0.123000*(STATES[0]+7.00000))) : 0.00138000/0.123000);
ALGEBRAIC[11] = (fabs(STATES[0]+10.0000)>0.00100000 ? ( 0.000610000*1.00000*(STATES[0]+10.0000))/(exp( 0.145000*(STATES[0]+10.0000)) - 1.00000) : 0.000610000/0.145000);
ALGEBRAIC[15] = 1.00000/(ALGEBRAIC[4]+ALGEBRAIC[11]);
ALGEBRAIC[19] = 1.00000/(1.00000+exp(- (STATES[0]+50.0000)/7.50000));
ALGEBRAIC[12] = ALGEBRAIC[5];
ALGEBRAIC[20] = 4.00000*ALGEBRAIC[16];
ALGEBRAIC[6] = - (STATES[0]+3.00000)/15.0000;
ALGEBRAIC[17] = 1.00000/(1.00000+exp(ALGEBRAIC[6]));
ALGEBRAIC[22] = 9.00000/(1.00000+exp(- ALGEBRAIC[6]))+0.500000;
ALGEBRAIC[21] = ALGEBRAIC[17];
ALGEBRAIC[13] = ( (- STATES[0]/30.0000)*STATES[0])/30.0000;
ALGEBRAIC[24] = 3.50000*exp(ALGEBRAIC[13])+1.50000;
ALGEBRAIC[18] = 1.00000/(1.00000+exp(- (STATES[0] - CONSTANTS[14])/CONSTANTS[15]));
ALGEBRAIC[23] = ALGEBRAIC[18]/CONSTANTS[31];
ALGEBRAIC[25] = (1.00000 - ALGEBRAIC[18])/CONSTANTS[31];
ALGEBRAIC[26] = 1.00000/(1.00000+pow(CONSTANTS[22]/STATES[4], 3.00000));
ALGEBRAIC[34] = 10.0000+ 4954.00*exp(STATES[0]/15.6000);
ALGEBRAIC[35] = CONSTANTS[30]/(1.00000+pow(STATES[4]/CONSTANTS[23], 4.00000))+0.100000;
ALGEBRAIC[33] = 1.00000 - 1.00000/(1.00000+exp(- (STATES[0] - CONSTANTS[18])/CONSTANTS[19]));
ALGEBRAIC[36] = (ALGEBRAIC[34] - ALGEBRAIC[35])*ALGEBRAIC[33]+ALGEBRAIC[35];
ALGEBRAIC[38] = 1.00000/(1.00000+exp(- (STATES[0] - CONSTANTS[20])/CONSTANTS[21]));
ALGEBRAIC[39] = ( ALGEBRAIC[26]*ALGEBRAIC[38])/ALGEBRAIC[36];
ALGEBRAIC[40] = (1.00000 - ALGEBRAIC[38])/ALGEBRAIC[36];
ALGEBRAIC[37] = (ALGEBRAIC[34] - 450.000)*ALGEBRAIC[33]+450.000;
ALGEBRAIC[41] = ALGEBRAIC[38]/ALGEBRAIC[37];
ALGEBRAIC[42] = (1.00000 - ALGEBRAIC[38])/ALGEBRAIC[37];
ALGEBRAIC[30] = 1.00000/(1.00000+exp(- (STATES[0] - CONSTANTS[16])/CONSTANTS[17]));
ALGEBRAIC[31] = (1.00000 - ALGEBRAIC[30])/CONSTANTS[32];
ALGEBRAIC[28] = 0.0241680*ALGEBRAIC[26];
ALGEBRAIC[43] = ( (( (( ALGEBRAIC[31]*ALGEBRAIC[23])/ALGEBRAIC[25])*ALGEBRAIC[28])/CONSTANTS[24])*ALGEBRAIC[40])/ALGEBRAIC[39];
ALGEBRAIC[32] = ALGEBRAIC[31];
ALGEBRAIC[44] = ( (( (( ALGEBRAIC[32]*ALGEBRAIC[23])/ALGEBRAIC[25])*CONSTANTS[25])/CONSTANTS[26])*ALGEBRAIC[42])/ALGEBRAIC[41];
ALGEBRAIC[45] = (((((1.00000 - STATES[7]) - STATES[9]) - STATES[8]) - STATES[10]) - STATES[5]) - STATES[6];
ALGEBRAIC[27] = 0.0182688*ALGEBRAIC[26];
ALGEBRAIC[29] = ( (( ALGEBRAIC[27]*CONSTANTS[27])/CONSTANTS[28])*CONSTANTS[24])/ALGEBRAIC[28];
ALGEBRAIC[53] = (STATES[0]+33.5000)/10.0000;
ALGEBRAIC[59] = 1.00000/(1.00000+exp(ALGEBRAIC[53]));
ALGEBRAIC[55] = (STATES[0]+60.0000)/10.0000;
ALGEBRAIC[62] = 3000.00/(1.00000+exp(ALGEBRAIC[55]))+30.0000;
ALGEBRAIC[61] = ALGEBRAIC[59];
ALGEBRAIC[57] = (STATES[0]+33.5000)/10.0000;
ALGEBRAIC[64] = 20.0000/(1.00000+exp(ALGEBRAIC[57]))+20.0000;
ALGEBRAIC[72] = ( CONSTANTS[58]*STATES[12]*STATES[12])/( STATES[12]*STATES[12]+ CONSTANTS[56]*CONSTANTS[56]);
ALGEBRAIC[73] = (( CONSTANTS[59]*STATES[21]*STATES[21])/( STATES[21]*STATES[21]+ CONSTANTS[57]*CONSTANTS[57]))*( STATES[21]*16.6670 - STATES[12]);
ALGEBRAIC[79] = ( CONSTANTS[60]*CONSTANTS[61])/( (CONSTANTS[61]+STATES[12])*(CONSTANTS[61]+STATES[12]));
ALGEBRAIC[80] = ( CONSTANTS[62]*CONSTANTS[63])/( (CONSTANTS[63]+STATES[12])*(CONSTANTS[63]+STATES[12]));
ALGEBRAIC[81] = ( CONSTANTS[64]*CONSTANTS[65])/( (CONSTANTS[65]+STATES[12])*(CONSTANTS[65]+STATES[12]));
ALGEBRAIC[82] = ( CONSTANTS[66]*CONSTANTS[67])/( (CONSTANTS[67]+STATES[12])*(CONSTANTS[67]+STATES[12]));
ALGEBRAIC[83] = 1.00000/(1.00000+ALGEBRAIC[79]+ALGEBRAIC[80]+ALGEBRAIC[81]+ALGEBRAIC[82]);
ALGEBRAIC[86] = CONSTANTS[68]*STATES[12]*(CONSTANTS[70] - STATES[24]) - CONSTANTS[69]*STATES[24];
ALGEBRAIC[84] = (STATES[20] - STATES[12])/CONSTANTS[71];
ALGEBRAIC[85] = CONSTANTS[68]*STATES[20]*(CONSTANTS[70] - STATES[25]) - CONSTANTS[69]*STATES[25];
ALGEBRAIC[87] = (- STATES[23]+ALGEBRAIC[72]) - ALGEBRAIC[73];
ALGEBRAIC[69] = (STATES[22]>50.0000&&STATES[22]<CONSTANTS[47] ? (STATES[22] - 50.0000)/1.00000 : STATES[22]>=CONSTANTS[47] ? CONSTANTS[53]*STATES[22]+CONSTANTS[78] : 0.00000);
ALGEBRAIC[70] = ( STATES[21]*ALGEBRAIC[69])/CONSTANTS[47];
ALGEBRAIC[88] = STATES[20]/1000.00;
ALGEBRAIC[14] = STATES[0]*2.00000*CONSTANTS[75];
ALGEBRAIC[89] = (fabs(ALGEBRAIC[14])<0.00100000 ? ( 4.00000*CONSTANTS[13]*CONSTANTS[2]*CONSTANTS[75]*( ALGEBRAIC[88]*exp(ALGEBRAIC[14]) - 0.341000*CONSTANTS[4]))/( 2.00000*CONSTANTS[75]) : ( 4.00000*CONSTANTS[13]*STATES[0]*CONSTANTS[2]*CONSTANTS[75]*( ALGEBRAIC[88]*exp(ALGEBRAIC[14]) - 0.341000*CONSTANTS[4]))/(exp(ALGEBRAIC[14]) - 1.00000));
ALGEBRAIC[71] = exp( - CONSTANTS[52]*(STATES[0]+30.0000))/(1.00000+exp( - CONSTANTS[52]*(STATES[0]+30.0000)));
ALGEBRAIC[91] = (CONSTANTS[48]/1.00000)*ALGEBRAIC[45]*fabs(ALGEBRAIC[89])*ALGEBRAIC[71];
ALGEBRAIC[92] = ( (( ALGEBRAIC[45]*ALGEBRAIC[70]*fabs(ALGEBRAIC[89])*CONSTANTS[49])/1.00000)*exp( - CONSTANTS[51]*(STATES[0]+30.0000)))/(1.00000+exp( - CONSTANTS[51]*(STATES[0]+30.0000)));
ALGEBRAIC[94] = ALGEBRAIC[45]*CONSTANTS[50]*fabs(ALGEBRAIC[89]);
ALGEBRAIC[96] = ALGEBRAIC[92]+ALGEBRAIC[94];
ALGEBRAIC[90] = CONSTANTS[12]*ALGEBRAIC[45]*ALGEBRAIC[89];
ALGEBRAIC[67] = 1.00000/(1.00000+pow(CONSTANTS[42]/STATES[20], 3.00000));
ALGEBRAIC[95] = pow(STATES[19], 3.00000)*CONSTANTS[4]*exp( STATES[0]*0.350000*CONSTANTS[75]) - pow(CONSTANTS[5], 3.00000)*ALGEBRAIC[88]*exp( STATES[0]*(0.350000 - 1.00000)*CONSTANTS[75]);
ALGEBRAIC[66] = 1.00000+ 0.200000*exp( STATES[0]*(0.350000 - 1.00000)*CONSTANTS[75]);
ALGEBRAIC[97] = CONSTANTS[43]*pow(STATES[19], 3.00000)+ pow(CONSTANTS[44], 3.00000)*ALGEBRAIC[88];
ALGEBRAIC[98] = pow(CONSTANTS[45], 3.00000)*CONSTANTS[4]*(1.00000+ALGEBRAIC[88]/CONSTANTS[46]);
ALGEBRAIC[68] = CONSTANTS[46]*pow(CONSTANTS[5], 3.00000)*(1.00000+pow(STATES[19]/CONSTANTS[45], 3.00000));
ALGEBRAIC[99] = pow(STATES[19], 3.00000)*CONSTANTS[4]+ pow(CONSTANTS[5], 3.00000)*ALGEBRAIC[88];
ALGEBRAIC[100] = ALGEBRAIC[97]+ALGEBRAIC[98]+ALGEBRAIC[68]+ALGEBRAIC[99];
ALGEBRAIC[101] = ( CONSTANTS[41]*ALGEBRAIC[67]*ALGEBRAIC[95])/( ALGEBRAIC[66]*ALGEBRAIC[100]);
ALGEBRAIC[74] = ( CONSTANTS[60]*CONSTANTS[61])/( (CONSTANTS[61]+STATES[20])*(CONSTANTS[61]+STATES[20]));
ALGEBRAIC[75] = ( CONSTANTS[62]*CONSTANTS[63])/( (CONSTANTS[63]+STATES[20])*(CONSTANTS[63]+STATES[20]));
ALGEBRAIC[76] = ( CONSTANTS[64]*CONSTANTS[65])/( (CONSTANTS[65]+STATES[20])*(CONSTANTS[65]+STATES[20]));
ALGEBRAIC[77] = ( CONSTANTS[66]*CONSTANTS[67])/( (CONSTANTS[67]+STATES[20])*(CONSTANTS[67]+STATES[20]));
ALGEBRAIC[78] = 1.00000/(1.00000+ALGEBRAIC[74]+ALGEBRAIC[75]+ALGEBRAIC[76]+ALGEBRAIC[77]);
ALGEBRAIC[63] = 1.00000/(1.00000+ 0.124500*exp( - 0.100000*STATES[0]*CONSTANTS[75])+ 0.0365000*CONSTANTS[77]*exp( - STATES[0]*CONSTANTS[75]));
ALGEBRAIC[65] = ( (( CONSTANTS[38]*ALGEBRAIC[63]*STATES[19])/(STATES[19]+CONSTANTS[40]))*CONSTANTS[3])/(CONSTANTS[3]+CONSTANTS[39]);
ALGEBRAIC[102] = CONSTANTS[6]*ALGEBRAIC[101];
ALGEBRAIC[105] = (1.00000/CONSTANTS[75])*log(CONSTANTS[5]/STATES[19]);
ALGEBRAIC[106] = CONSTANTS[11]*STATES[2]*STATES[3]*STATES[1]*STATES[1]*STATES[1]*(STATES[0] - ALGEBRAIC[105]);
ALGEBRAIC[46] = 1.02000/(1.00000+exp( 0.238500*((STATES[0] - CONSTANTS[79]) - 59.2150)));
ALGEBRAIC[47] = ( 0.491240*exp( 0.0803200*((STATES[0] - CONSTANTS[79])+5.47600))+ 1.00000*exp( 0.0617500*((STATES[0] - CONSTANTS[79]) - 594.310)))/(1.00000+exp( - 0.514300*((STATES[0] - CONSTANTS[79])+4.75300)));
ALGEBRAIC[48] = ALGEBRAIC[46]/(ALGEBRAIC[46]+ALGEBRAIC[47]);
ALGEBRAIC[49] = CONSTANTS[33]* pow((CONSTANTS[3]/5.40000), 1.0 / 2)*ALGEBRAIC[48]*(STATES[0] - CONSTANTS[79]);
ALGEBRAIC[54] = 1.00000/(1.00000+exp(ALGEBRAIC[53]));
ALGEBRAIC[56] = CONSTANTS[36]*STATES[15]*(STATES[16]+ 0.500000*ALGEBRAIC[54])*(STATES[0] - CONSTANTS[79]);
ALGEBRAIC[58] = CONSTANTS[37]*STATES[17]*STATES[18]*(STATES[0] - CONSTANTS[79]);
ALGEBRAIC[60] = ALGEBRAIC[56]+ALGEBRAIC[58];
ALGEBRAIC[93] = 2.00000*CONSTANTS[6]*ALGEBRAIC[90];
ALGEBRAIC[50] = 1.00000/(1.00000+exp((STATES[0]+33.0000)/22.4000));
ALGEBRAIC[51] = CONSTANTS[34]* pow((CONSTANTS[3]/5.40000), 1.0 / 2)*STATES[11]*ALGEBRAIC[50]*(STATES[0] - CONSTANTS[79]);
ALGEBRAIC[103] = (1.00000/CONSTANTS[75])*log((CONSTANTS[3]+ CONSTANTS[74]*CONSTANTS[5])/(CONSTANTS[73]+ CONSTANTS[74]*STATES[19]));
ALGEBRAIC[52] = 1.00000+0.800000/(1.00000+pow(0.500000/STATES[12], 3.00000));
ALGEBRAIC[104] = CONSTANTS[35]*ALGEBRAIC[52]*STATES[13]*STATES[14]*(STATES[0] - ALGEBRAIC[103]);
ALGEBRAIC[3] = floor(VOI/CONSTANTS[8])*CONSTANTS[8];
ALGEBRAIC[10] = (VOI - ALGEBRAIC[3]>=CONSTANTS[7]&&VOI - ALGEBRAIC[3]<=CONSTANTS[7]+CONSTANTS[9] ? CONSTANTS[10] : 0.00000);
ALGEBRAIC[107] = - (ALGEBRAIC[106]+ALGEBRAIC[49]+ALGEBRAIC[51]+ALGEBRAIC[104]+ALGEBRAIC[60]+ALGEBRAIC[102]+ALGEBRAIC[93]+ALGEBRAIC[65]+ALGEBRAIC[10]);
}