pop_het_exp function

Compute the population expected heterozygosity