A generalized mathematical model of the breast and ovarian cancer is developed by considering the fractional differential equations with Caputo time-fractional derivatives. The use of the fractional model shows that the time-evolution of the proliferating cell mass, the quiescent cell mass, and the proliferative function are significantly influenced by their history. Even if the classical model, based on the derivative of integer order has been studied in many papers, its analytical solutions are presented in order to make the comparison between the classical model and the fractional model. Using the finite difference method, numerical schemes to the Caputo derivative operator and Riemann-Liouville fractional integral operator are obtained. Numerical solutions to the fractional differential equations of the generalized mathematical model are determined for the chemotherapy scheme based on the function of "on-off" type. Numerical results, obtained with the Mathcad software, are discussed and presented in graphical illustrations. The presence of the fractional order of the time-derivative as a parameter of solutions gives important information regarding the proliferative function, therefore, could give the possible rules for more efficient chemotherapy.
Numerical solutions of multi-order fractional differential equations by Boubaker polynomials
