Abstract
Analyzing the dynamics of tumor-immune systems can play an important role in the fight against cancer, since it can foster the development of more effective medical treatments. This paper was aimed at making a contribution to the study of tumor-immune dynamics by presenting a new model of cancer growth based on fractional-order differential equations. By investigating the system dynamics, the manuscript highlights the chaotic behaviors of the proposed cancer model for both the commensurate and the incommensurate cases. Bifurcation diagrams, the Lyapunov exponents, and phase plots confirm the effectiveness of the conceived approach. Finally, some considerations regarding the biological meaning of the obtained results are reported through the manuscript.