In this paper a new mathematical model is developed for the dynamics between tumor cells, normal cells, immune cells, chemotherapy drug concentration and drug toxicity. Then, the theorem of Lyapunov stability is applied to design treatment strategies for drug protocols that ensure a desired rate of tumor cell kill and push the system to the area with smaller tumor cells. Using of this theorem a...