De Pillis chemotherapy model (Muscod)

This is source code written in Muscod for the De Pillis chemotherapy model.

/* volume of tumor in absolute cell count*/
rhs[0] = - a*x[0] (1-b*x[0]) -c*x[1]*x[0] - D*x[0] - K_T*(1- exp(- x[4]))*x[0];

/* volume of unspecific immune cells */
rhs[1] = e*x[3] - f*x[1] + g*x[0]*x[0]/(h+x[0]*x[0])-p*x[1]*x[0] - K_N*(1- exp(- x[4]))*x[1];

/* volume of tumor-specific cytotoxic T-cells */
rhs[2] = -m*x[2] + j*D*D*x[0]*x[0]/(k+ D*D*x[0]*x[0])*x[2] - q*x[1]*x[2] + (r_1*x[1] + r_2*x[3])*x[0] ...
- v*x[1]* x[2]*x[2] - K_L* (1- exp(- x[4]))*x[2] + p_I*x[2]*x[5]/(g_I + x[5]) + u[2];

/* volume of circulating lymphocytes */
rhs[3] = alpha - beta *x[3] - K_C *(1- exp(- x[4])) *x[3];

/* amount of chemotherapeutic drug */
rhs[4] = u[0] - gamma*x[4];

/* amount of immunotherapeutic drug */
rhs[5] = u[1] - mu_I*x[5];