We use a system biology approach to translate the interaction of Bacillus Calmette-Gurin (BCG) + interleukin 2 (IL-2) for the treatment of bladder cancer into a mathematical model. The model is presented as a system of differential equations with the following variables: number of tumor cells, bacterial cells, immune cells, and cytokines involved in the tumor-immune response. This work investigates the delay effect induced by the proliferation of tumor antigen-specific effector cells after the immune system destroys BCG-infected urothelium cells following BCG and IL-2 immunotherapy in the treatment of bladder cancer. For the proposed model, three equilibrium states are found analytically. The stability of all equilibria is analyzed using the method of Lyapunov functionals construction and the method of linear matrix inequalities (LMIs).