In this paper, a mathematical model for chemotherapy that takes tumor immune-system interactions into account is considered for a strongly targeted agent. We use a classical model originally formulated by Stepanova, but replace exponential tumor growth with a generalised logistic growth model function depending on a parameter v. This growth function interpolates between a Gompertzian model (in the limit v → 0) and an exponential model (in the limit v → ∞). The dynamics is multi-stable and equilibria and their stability will be investigated depending on the parameter v. Except for small values of v, the system has both an asymptotically stable microscopic (benign) equilibrium point and an asymptotically stable macroscopic (malignant) equilibrium point. The corresponding regions of attraction are separated by the stable manifold of a saddle. The optimal control problem of moving an initial condition that lies in the malignant region into the benign region is formulated and the structure of optimal singular controls is determined.