Glioblastoma Multiforme (GBM) is a highly aggressive and malignant type of brain tumour. Besides the typical hallmarks of cancer, such as uncontrolled cellular proliferation and instability, GBM also exhibits dramatic invasive potential and resistance to common therapies: even with a complete treatment including neurosurgery, chemotherapy and radiotherapy, the median survival time is about 10-16 months. Hence, there is a critical need to understand and replicate the biological complexity of the brain, in order to predict tumour evolution and arrange therapeutic strategies accordingly. For that purpose, mathematical and computational models can provide powerful instruments for investigating GBM progression: in the last decades, several models that describe brain tumour growth have been proposed, using different frameworks and accounting for different characteristics.