Leukaemia accounts for around 3% of all cancer types diagnosed in adults, and is the most common type of cancer in children of paediatric age (typically ranging from 0 to 14 years). There is increasing interest in the use of mathematical models in oncology to draw inferences and make predictions, providing a complementary picture to experimental biomedical models. In this paper we recapitulate the state of the art of mathematical modelling of leukaemia growth dynamics, in time and response to treatment. We intend to describe the mathematical methodologies, the biological aspects taken into account in the modelling, and the conclusions of each study. This review is intended to provide researchers in the field with solid background material, in order to achieve further breakthroughs in the promising field of mathematical biology.