A non-local rare mutations model for quasispecies and Prisoner's dilemma: numerical assessment of qualitative behaviour

An integro-differential model for evolutionary dynamics with mutations is investigated by improving the understanding of its behavior using numerical simulations. The proposed numerical approach can handle also density dependent fitness, and gives new insights about the role of mutation in the preservation of cooperation.