  Metastasis, the spread of cancer cells from a primary tumor to secondary location(s) in the human organism, is the ultimate cause of death for the majority of cancer patients. That is why, it is crucial to understand metastases evolution in order to successfully combat the disease. We consider a metastasized cancer cell population after medical treatment (e.g. chemotherapy). Arriving in a different environment the cancer cells may change their lifespan and reproduction, thus they may proliferate into different types. If the treatment is effective, in the context of branching processes it means, the reproduction of cancer cells is such that the mean offspring of each cell is less than one. However, it is possible mutations to occur during cell division cycle. These mutations can produce a new cancer cell type, which is resistant to the treatment. Cancer cells from this new type may lead to the rise of a non-extinction branching process. The above scenario leads us to the choice of a reducible multi-type age-dependent branching process as a relevant framework for studying the asymptotic behavior of such complex structures. Our previous theoretical results are related to the asymptotic behavior of the waiting time until the first occurrence of a mutant starting a non-extinction process and the modified hazard function as a measure of immediate recurrence of cancer disease. In the present paper these asymptotic results are used for developing numerical schemes and algorithms implemented in Python via the NumPy package for approximate calculation of the corresponding quantities. In conclusion, our conjecture is that this methodology can be advantageous in revealing the role of the lifespan distribution of the cancer cells in the context of cancer disease evolution and other complex cell population systems, in general. 
Computational modelling of cancer evolution by multi-type branching   processes