Monte Carlo simulations of parapatric speciation