gl.sim.Neconst function

Simulate a population with constant mutation rate