Mathematical analysis and simulations of a cancer model with interleukins and delayed immunotherapy