Create R script to play with the neutral Wright-Fisher model