OUTFLOWS:
AB_REMATE ¼ AB_GENOTYPE
A_ALLELES(t) ¼ A_ALLELES(t À dt) + (RELEASE_A + MUTATION À
A_MATE) * dt
INIT A_ALLELES ¼ 200
INFLOWS:
RELEASE_A ¼ 2 * AA_GENOTYPE +AB_GENOTYPE
MUTATION ¼ RANDOM(À0.05,0.05) * (A_ALLELES+B_ALLELES)
OUTFLOWS:
A_MATE ¼ A_ALLELES
BB_GENOTYPE(t) ¼ BB_GENOTYPE(t À dt) + (Δ_BB_GENOTYPE À
BB_REMATE) * dt
INIT BB_GENOTYPE ¼ 0
INFLOWS:
Δ_BB_GENOTYPE ¼ TOTAL_ALLELES/2 * (1 À A_FREQ)^2
OUTFLOWS:
BB_REMATE ¼ BB_GENOTYPE
B_ALLELES(t) ¼ B_ALLELES(t À dt) + (RELEASE_B À B_MATE À MUTATION) * dt
INIT B_ALLELES ¼ 300
INFLOWS:
RELEASE_B ¼ 2 * BB_GENOTYPE + AB_GENOTYPE
OUTFLOWS:
B_MATE ¼ B_ALLELES
MUTATION ¼ RANDOM(À0.05,0.05) * (A_ALLELES+B_ALLELES)
AA_FREQ ¼ IF TIME > 0 THEN AA_GENOTYPE/TOTAL_GENOTYPE ELSE 0
AB_FREQ ¼ IF TIME > 0 THEN AB_GENOTYPE/TOTAL_GENOTYPE ELSE 0
A_FREQ ¼ A_ALLELES/(A_ALLELES+ B_ALLELES)
BB_FREQ ¼ IF TIME > 0 THEN BB_GENOTYPE/TOTAL_GENOTYPE ELSE 0
TOTAL_ALLELES ¼ A_ALLELES + B_ALLELES
TOTAL_GENOTYPE ¼ AA_GENOTYPE+AB_GENOTYPE+BB_GENOTYPE
12.2 Mating and Mutation of Alleles Model Equations
103