Mathematica notebook with simulation methods

This notebook includes the simulation model for the paper "Introgression makes waves in inferred histories of effectivepopulation size". In the present state, the key elements of the simulations are expressed in functions, with the simulation of population structure with gene flow or introgression in the function "funkycoalescence". These functions are insufficiently commented as they stand, and comments need to be added. The function incorporates the coalescence approach as described by Hudson (1990), with mutation (added as "rate") and migration (or instantaneous separation into subpopulations).