Millions of people are affected by leukemia. It is important to understand the progression dynamics of this disease to be able to minimize the damage that is caused by it. This article provides a mathematical framework to develop strategies to control leukemia. Bifurcation analysis is a powerful mathematical tool used to deal with the nonlinear dynamics of any process. Several factors must be considered, and multiple objectives must be met simultaneously. Bifurcation analysis and Multiobjective Nonlinear Model Predictive Control (MNLMPC) calculations are performed on three leukemia models. The MATLAB program MATCONT was used to perform the bifurcation analysis. The MNLMPC calculations were performed using the optimization language PYOMO in conjunction with the state-of-the-art global optimization solvers IPOPT and BARON. The bifurcation analysis revealed the existence of limit points and branch in the models. The limit and branch points were beneficial because they enabled the multiobjective nonlinear model predictive control calculations to converge to the Utopia point in both problems, which is the most beneficial solution. A combination of bifurcation analysis and multiobjective nonlinear model predictive control for leukemia models is the main contribution of this paper.