A mathematical study has been conducted in this article combining
immunotherapy, radiation therapy, and estrogen inhibition  to treat 
breast cancer.  In this regard, a computational model has been 
developed consisting of eight nonlinear reaction-diffusion equations 
which have been solved numerically by perturbed functional iterations. 
The interactions between the immune system, estrogen-inhibitor, 
antioxidant, localized radiation therapy and progress of cancer have
been considered in the reaction part of the equations.  The diffusion
part shows the flux of cancer cells toward distant organs followed by
various leukocytes and medications which try to destroy cancer cells
and inhibit any metastasis.  Computational results have  strong
agreement with some of the findings recorded in the literature on this
topic.