fisher_reaction_diffusion Paul Cochrane Example simulation of the reaction-diffusion equation of Fisher (1937). This equation can be interpreted as describing the frequency of a focal allele as a function of one-dimensional space in a biallelic diploid population of constant size. This simulation is adapted from equation 1 of "The diffusive spread of alleles in heterogeneous populations", G T Skalski, Evolution, 58 (3), 470-478 (2004). Thanks to John Novembre at UC Berkeley for pointing out this example. t yes yes yes yes main x 1024 (-1,1) 1 main complex P no RK4IP 1 1000 100 main yes L no 1024 POut