C++ code to simulate a Wright-Fisher model according to Simon-Park-Krug algorithm with a decreasing mutation rate and a sub-linear increasing fitness.