Evolution of phenotypic polymorphism in symbiont-pairing in plant-fungal symbiosis

This is Mathematica code to model the evolution of polyphenism in plants
(i.e., the ability of a single genotype to randomly develop into
individuals with different phenotypes). We are interested in the trait of
roots allocated to two different fungal partners when the benefits of a
particular ratio of partners vary with the environment. We use both an
analytical model and a simulation to investigate the conditions for
polyphenism to evolve and evolutionary endpoints -- the root allocation
phenotypes that evolved to be expressed and the probability that each one
is expressed. Along with the analytical model and simulation code, we also
include saved simulation results as .csv files so that the results can be
analyzed without needing to re-run the simulations. The saved simulation
results include the evolved traits (root allocation and phenotype
frequency) and the time it took for populations to evolve away from
monophenism.