In this thesis we construct a new model of the immune response to the growth of breast cancer cells and investigate the impact of certain drug therapies on the cancer. We use delay differential equations to model the interaction of breast cancer cells with the immune system. The new model is constructed by combining two previous models. The first model accounts for different cell cycles and includes terms to evaluate drug treatments, but ignores quiescent tumor cells. The second model includes quiescent cells, but ignores the immune response and drug treatments. The new model is obtained by combining and modifying these two models to account for quiescent cells, immune cells and includes drug intervention terms. This new model is used to evaluate the effects of pulsed applications of the drug Paclitaxel for models with and without quiescent cells. We use sensitivity equation methods to analyze the sensitivity of the model with respect to the initial number of immune cytotoxic T-cells. Numerical experiments are conducted to compare the model predictions to observed behavior.