/*
   There are a total of 13 entries in the algebraic variable array.
   There are a total of 4 entries in each of the rate and state variable arrays.
   There are a total of 43 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (day).
 * STATES[0] is OB_p in component OB_p (dimensionless).
 * CONSTANTS[0] is D_OB_u in component OB_p (first_order_rate_constant).
 * CONSTANTS[1] is K_D1_TGF_beta in component OB_p (pM).
 * CONSTANTS[2] is K_D2_TGF_beta in component model_parameters (pM).
 * CONSTANTS[3] is D_OB_p in component model_parameters (first_order_rate_constant).
 * ALGEBRAIC[3] is TGF_beta in component TGF_beta (pM).
 * STATES[1] is OB_a in component OB_a (dimensionless).
 * CONSTANTS[4] is A_OB_a in component OB_a (first_order_rate_constant).
 * STATES[2] is OC_a in component OC_a (dimensionless).
 * ALGEBRAIC[1] is D_OC_p in component OC_a (first_order_rate_constant).
 * CONSTANTS[5] is A_OC_a in component OC_a (first_order_rate_constant).
 * CONSTANTS[6] is K_D3_TGF_beta in component OC_a (pM).
 * CONSTANTS[7] is OC_p in component OC_a (dimensionless).
 * ALGEBRAIC[7] is RANKL in component RANKL (pM).
 * ALGEBRAIC[8] is pi_RANKL_act_OC_p in component RANKL (dimensionless).
 * STATES[3] is BV in component BV (dimensionless).
 * ALGEBRAIC[0] is OC_a_tilde in component BV (dimensionless).
 * ALGEBRAIC[2] is OB_a_tilde in component BV (dimensionless).
 * CONSTANTS[8] is OC_a_initial in component BV (dimensionless).
 * CONSTANTS[9] is OB_a_initial in component BV (dimensionless).
 * CONSTANTS[10] is k_form in component BV (first_order_rate_constant).
 * CONSTANTS[11] is k_res in component model_parameters (first_order_rate_constant).
 * CONSTANTS[12] is alpha in component TGF_beta (pM).
 * CONSTANTS[13] is D_TGF_beta_tilde in component TGF_beta (first_order_rate_constant).
 * CONSTANTS[14] is S_TGF_beta in component TGF_beta (flux).
 * CONSTANTS[36] is PTH in component PTH (pM).
 * CONSTANTS[37] is pi_PTH_rep_OB_p in component PTH (dimensionless).
 * CONSTANTS[38] is pi_PTH_act_OB_p in component PTH (dimensionless).
 * CONSTANTS[15] is beta_PTH in component PTH (flux).
 * CONSTANTS[16] is D_PTH_tilde in component PTH (first_order_rate_constant).
 * CONSTANTS[17] is P_PTH_d in component PTH (flux).
 * CONSTANTS[18] is K_D4_PTH in component PTH (pM).
 * CONSTANTS[19] is K_D5_PTH in component PTH (pM).
 * CONSTANTS[20] is K_D6_PTH in component PTH (pM).
 * CONSTANTS[21] is K_D7_PTH in component PTH (pM).
 * CONSTANTS[39] is pi_PTH_act_OB_a in component PTH (dimensionless).
 * CONSTANTS[40] is pi_PTH_rep_OB_a in component PTH (dimensionless).
 * ALGEBRAIC[4] is OPG in component OPG (pM).
 * CONSTANTS[22] is beta1_OPG in component OPG (flux).
 * CONSTANTS[23] is beta2_OPG in component OPG (flux).
 * CONSTANTS[24] is OPG_max in component OPG (pM).
 * ALGEBRAIC[5] is P_OPG_e in component OPG (flux).
 * CONSTANTS[25] is P_OPG_d in component OPG (flux).
 * CONSTANTS[26] is D_OPG_tilde in component OPG (first_order_rate_constant).
 * CONSTANTS[41] is pi_PTH_rep_OB in component OPG (dimensionless).
 * CONSTANTS[27] is K_A2_RANKL in component RANKL (per_pM).
 * CONSTANTS[28] is K_A1_RANKL in component RANKL (per_pM).
 * CONSTANTS[29] is RANK in component RANKL (pM).
 * ALGEBRAIC[6] is RANKL_eff in component RANKL (pM).
 * ALGEBRAIC[9] is RANKL_tot in component RANKL (pM).
 * ALGEBRAIC[12] is P_RANKL in component RANKL (flux).
 * ALGEBRAIC[10] is P_RANKL_e in component RANKL (flux).
 * CONSTANTS[30] is P_RANKL_d in component RANKL (flux).
 * CONSTANTS[31] is R1_RANKL in component RANKL (pM).
 * CONSTANTS[32] is R2_RANKL in component RANKL (pM).
 * CONSTANTS[33] is beta_RANKL in component RANKL (flux).
 * ALGEBRAIC[11] is D_RANKL in component RANKL (flux).
 * CONSTANTS[34] is D_RANKL_tilde in component RANKL (first_order_rate_constant).
 * CONSTANTS[42] is pi_PTH_act_OB in component RANKL (dimensionless).
 * CONSTANTS[35] is K_D8_RANKL in component RANKL (pM).
 * RATES[0] is d/dt OB_p in component OB_p (dimensionless).
 * RATES[1] is d/dt OB_a in component OB_a (dimensionless).
 * RATES[2] is d/dt OC_a in component OC_a (dimensionless).
 * RATES[3] is d/dt BV in component BV (dimensionless).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 1.0;
CONSTANTS[0] = 7.000e-4;
CONSTANTS[1] = 4.545e-3;
CONSTANTS[2] = 1.416e-3;
CONSTANTS[3] = 5.348e-0;
STATES[1] = 1.0;
CONSTANTS[4] = 1.890e-1;
STATES[2] = 1.0;
CONSTANTS[5] = 7.000e-1;
CONSTANTS[6] = 1.416e-3;
CONSTANTS[7] = 1.0;
STATES[3] = 100.0;
CONSTANTS[8] = 1.0;
CONSTANTS[9] = 1.0;
CONSTANTS[10] = 1.571e0;
CONSTANTS[11] = 1.000e0;
CONSTANTS[12] = 1.000e0;
CONSTANTS[13] = 1.000e0;
CONSTANTS[14] = 0.0;
CONSTANTS[15] = 2.500e2;
CONSTANTS[16] = 8.600e1;
CONSTANTS[17] = 0.0;
CONSTANTS[18] = 1.5e2;
CONSTANTS[19] = 1.5e2;
CONSTANTS[20] = 2.226e-1;
CONSTANTS[21] = 2.226e-1;
CONSTANTS[22] = 1.0;
CONSTANTS[23] = 0.0;
CONSTANTS[24] = 2.000e8;
CONSTANTS[25] = 0.0;
CONSTANTS[26] = 3.500e-1;
CONSTANTS[27] = 3.412e-2;
CONSTANTS[28] = 1.000e-3;
CONSTANTS[29] = 1.000e1;
CONSTANTS[30] = 0.0;
CONSTANTS[31] = 0.0;
CONSTANTS[32] = 1.0;
CONSTANTS[33] = 1.684e4;
CONSTANTS[34] = 1.013e1;
CONSTANTS[35] = 1.306e1;
CONSTANTS[36] = (CONSTANTS[15]+CONSTANTS[17])/CONSTANTS[16];
CONSTANTS[37] = 1.00000/(1.00000+CONSTANTS[36]/CONSTANTS[20]);
CONSTANTS[38] = CONSTANTS[36]/(CONSTANTS[18]+CONSTANTS[36]);
CONSTANTS[39] = CONSTANTS[36]/(CONSTANTS[19]+CONSTANTS[36]);
CONSTANTS[40] = 1.00000/(1.00000+CONSTANTS[36]/CONSTANTS[21]);
CONSTANTS[41] = CONSTANTS[37];
CONSTANTS[42] = CONSTANTS[38];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[3] = ( CONSTANTS[12]*CONSTANTS[11]*STATES[2]+CONSTANTS[14])/CONSTANTS[13];
RATES[0] =  CONSTANTS[0]*(ALGEBRAIC[3]/(CONSTANTS[1]+ALGEBRAIC[3])) -  CONSTANTS[3]*STATES[0]*(1.00000/(1.00000+ALGEBRAIC[3]/CONSTANTS[2]));
RATES[1] =  CONSTANTS[3]*STATES[0]*(1.00000/(1.00000+ALGEBRAIC[3]/CONSTANTS[2])) -  CONSTANTS[4]*STATES[1];
ALGEBRAIC[0] = STATES[2] - CONSTANTS[8];
ALGEBRAIC[2] = STATES[1] - CONSTANTS[9];
RATES[3] =  CONSTANTS[10]*ALGEBRAIC[2] -  CONSTANTS[11]*ALGEBRAIC[0];
ALGEBRAIC[1] = (VOI>0.00000&&VOI<=100.000 ? 0.0210000 : 0.00210000);
ALGEBRAIC[4] = ( ( CONSTANTS[22]*STATES[0]+ CONSTANTS[23]*STATES[1])*CONSTANTS[41]+CONSTANTS[25])/(( ( CONSTANTS[22]*STATES[0]+ CONSTANTS[23]*STATES[1])*CONSTANTS[41])/CONSTANTS[24]+CONSTANTS[26]);
ALGEBRAIC[6] =  ( CONSTANTS[31]*STATES[0]+ CONSTANTS[32]*STATES[1])*CONSTANTS[42];
ALGEBRAIC[7] =  (ALGEBRAIC[6]/(1.00000+ CONSTANTS[28]*ALGEBRAIC[4]+ CONSTANTS[27]*CONSTANTS[29]))*((CONSTANTS[33]+CONSTANTS[30])/(CONSTANTS[33]+ CONSTANTS[34]*ALGEBRAIC[6]));
ALGEBRAIC[8] = ALGEBRAIC[7]/(CONSTANTS[35]+ALGEBRAIC[7]);
RATES[2] =  ALGEBRAIC[1]*CONSTANTS[7]*ALGEBRAIC[8] -  CONSTANTS[5]*STATES[2]*(ALGEBRAIC[3]/(CONSTANTS[6]+ALGEBRAIC[3]));
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[3] = ( CONSTANTS[12]*CONSTANTS[11]*STATES[2]+CONSTANTS[14])/CONSTANTS[13];
ALGEBRAIC[0] = STATES[2] - CONSTANTS[8];
ALGEBRAIC[2] = STATES[1] - CONSTANTS[9];
ALGEBRAIC[1] = (VOI>0.00000&&VOI<=100.000 ? 0.0210000 : 0.00210000);
ALGEBRAIC[4] = ( ( CONSTANTS[22]*STATES[0]+ CONSTANTS[23]*STATES[1])*CONSTANTS[41]+CONSTANTS[25])/(( ( CONSTANTS[22]*STATES[0]+ CONSTANTS[23]*STATES[1])*CONSTANTS[41])/CONSTANTS[24]+CONSTANTS[26]);
ALGEBRAIC[6] =  ( CONSTANTS[31]*STATES[0]+ CONSTANTS[32]*STATES[1])*CONSTANTS[42];
ALGEBRAIC[7] =  (ALGEBRAIC[6]/(1.00000+ CONSTANTS[28]*ALGEBRAIC[4]+ CONSTANTS[27]*CONSTANTS[29]))*((CONSTANTS[33]+CONSTANTS[30])/(CONSTANTS[33]+ CONSTANTS[34]*ALGEBRAIC[6]));
ALGEBRAIC[8] = ALGEBRAIC[7]/(CONSTANTS[35]+ALGEBRAIC[7]);
ALGEBRAIC[5] =  ( CONSTANTS[22]*STATES[0]+ CONSTANTS[23]*STATES[1])*CONSTANTS[41]*(1.00000 - ALGEBRAIC[4]/CONSTANTS[24]);
ALGEBRAIC[9] =  ALGEBRAIC[7]*(1.00000+ CONSTANTS[28]*ALGEBRAIC[4]+ CONSTANTS[27]*CONSTANTS[29]);
ALGEBRAIC[10] =  CONSTANTS[33]*(1.00000 - ALGEBRAIC[9]/ALGEBRAIC[6]);
ALGEBRAIC[11] =  - CONSTANTS[34]*ALGEBRAIC[9];
ALGEBRAIC[12] = ALGEBRAIC[10]+CONSTANTS[30];
}