Abstract: | In this study, a mathematical, fractional‐order model was developed for B cell chronic lymphocytic leukemia, with immune system, and then analyzed. Interactions between B leukemia cells, natural killer cells, cytotoxic T cells, and T‐helper cells are considered to be incorporated into a system consisting of four fractional differential equations. For estimation of the parameters, clinical data of six patients were used. By numerical solution of the system, the interactions between the leukemia cell population and the immune system cell populations for values of α ∈ (0,1) at different times were explained. By determining points of equilibrium and stability of the system were met. Bifurcation analysis showed that use of the fractional‐order model, figure out unpredictable behaviors of the system such as saddle‐node, bistability and hysteresis phenomenon occurred in the system by changing the values of some of the parameters, it was predictable. |