Use of quasi-normal form to examine stability of tumor-free equilibrium in a mathematical model of BCG treatment of bladder cancer

Understanding the dynamics of human hosts and tumors is of critical importance. A mathematical model was developed by Bunimovich-Mendrazitsky et al. ([10]), who explored the immune response in bladder cancer as an effect of BCG treatment. This treatment exploits the host's own immune system to boost a response that will enable the host to rid itself of the tumor. Although this model was extensively studied using numerical simulation, no analytical results on global tumor dynamics were originally presented. In this work, we analyze stability in a mathematical model for BCG treatment of bladder cancer based on the use of quasi-normal form and stability theory. These tools are employed in the critical cases, especially when analysis of the linearized system is insufficient. Our goal is to gain a deeper insight into the BCG treatment of bladder cancer, which is based on a mathematical model and biological considerations, and thereby to bring us one step closer to the design of a relevant clinical protocol.