In this paper we propose a mathematical model to describe the evolution of leukemiain the bone marrow. The model is based on a reaction-diffusion system of equations in a porousmedium. We show the existence of two stationary solutions, one of them corresponds to the normalcase and another one to the pathological case. The leukemic state appears as a result of a bifurcationwhen the normal state loses its stability. The critical conditions of leukemia developmentare determined by the proliferation rate of leukemic cells and by their capacity to diffuse. Theanalytical results are confirmed and illustrated by numerical simulations.