A numerical solution to the equilibria of the two-locus two-allele selection model

Examination of the equilibria of the standard two-locus two-allele selection model leads to the construction of a polynomial with coefficients derived from selective values in the genotypic fitness matrix. This polynomial can be partially factored algebraically and numerical techniques are available to extract the roots of the remainder. Each root provides a possible value of the disequilibrium coefficient and the gametic frequencies at equilibrium, and these can be readily checked for stability.