ABSTRACT: This paper introduces a cancer model for differentiation therapy and analyzes its mathematical properties. We explore non-negativity, boundedness, and equilibrium existence. Furthermore, we investigate bifurcations in the model, particularly transcritical and Hopf types. Using Sotomayor’s theorem, we demonstrate the presence of these bifurcations. Additionally, we validate our theoretical analysis by conducting numerical simulations and comparing them with the system’s phase diagram.