Mathematical modeling has been used to simulate the interaction of chemotherapy and immunotherapy drugs intervention with the dynamics of tumor cells growth. This work studies the interaction of cells in the immune system, such as the natural killer, dendritic, and cytotoxic CD8+ T cells, with chemotherapy. Four different cases were considered in the simulation: no drug intervention, independent interventions (either chemotherapy or immunotherapy), and combined interventions of chemotherapy and immunotherapy. The system of ordinary differential equations was initially solved using the Runge-Kutta method and compared with two additional methods: the Explicit Euler and Heun’s methods. Results showed that the combined intervention is more effective compared to the other cases. In addition, when compared with Runge-Kutta, the Heun’s method presented a better accuracy than the Explicit Euler technique. The proposed mathematical model can be used as a tool to improve cancer treatments and targeted therapy.