(eng): We investigate sufficient conditions for the presence of coexistence states for different genotypes in a diploid diallelic population with dominance distributed on a heterogeneous habitat, considering also the interaction between genes at multiple loci. In mathematical terms, this corresponds to the study of the Neumann boundary value problem.