Abstract
We present a revised mathematical model of the immune response to Bacillus Calmette-Guérin (BCG) treatment of bladder cancer, optimized according to biological and clinical data accumulated during the last decade. The improved model accounts for cytotoxic T lymphocyte differentiation as an integral element of the delayed immune response, as well as the logistic growth terms for cancer cell proliferation. Three equilibria are demonstrated for the proposed model, which is assumed to be influenced by white noise stochastic perturbations that are directly proportional to the system state deviation from an equilibrium. Stability conditions for all equilibria are analyzed using the Kolmanovskii-Shaikhet general method of Lyapunov functionals construction.