Numerical calculation of the adaptive-dynamics model.

<p>Each panel shows evolutionary outcomes of the antagonistic coevolution after 100,000 generations under various ratios of the mortality coefficients to the coefficient of female fertilisation (<i>β_f</i>/<i>α</i> and <i>β_m</i>/<i>α</i>) The two panels differ in the value of male PRR (<b>A</b>: <i>μ</i> = 2.00, <b>B</b>: <i>μ</i> = 0.50). Each coloured region indicates a different evolutionary outcome (white: prevalence of sexual reproduction, light grey: co-occurrence of two reproductive modes; dark grey: oscillation of the occurrence of parthenogenesis). Other parameters are as in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0058141#pone-0058141-g002" target="_blank">Figure 2</a>.