An integro-differential model for evolutionary dynamics with mutations is investigated by improving the understanding of its behaviour 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.