This ordinary differential equation model, simulating the tumor-immune interactions involved in BCG immunotherapy to treat superficial bladder cancer, is described by the publication:
Bunimovich-Mendrazitsky, S., Shochat, E., Stone, L. "Mathematical Model of BCG Immunotherapy in Superficial Bladder Cancer". Bull. Math. Biol. 69, 1847–1870 (2007). DOI: 10.1007/s11538-007-9195-z

Comment:
This model is based on the system of ODEs given in Equation 4 of the publication manuscript.
Reproduction of Figure 4 was achieved by setting p4 = 0.085.
Abstract:
Immunotherapy with Bacillus Calmette-Guérin (BCG)-an attenuated strain of Mycobacterium bovis (M. bovis) used for anti tuberculosis immunization-is a clinically established procedure for the treatment of superficial bladder cancer. However, the mode of action has not yet been fully elucidated, despite much extensive biological experience. The purpose of this paper is to develop a first mathematical model that describes tumor-immune interactions in the bladder as a result of BCG therapy. A mathematical analysis of the ODE model identifies multiple equilibrium points, their stability properties, and bifurcation points. Intriguing regimes of bistability are identified in which treatment has potential to result in a tumor-free equilibrium or a full-blown tumor depending only on initial conditions. Attention is given to estimating parameters and validating the model using published data taken from in vitro, mouse and human studies. The model makes clear that intensity of immunotherapy must be kept in limited bounds. While small treatment levels may fail to clear the tumor, a treatment that is too large can lead to an over-stimulated immune system having dangerous side effects for the patient.