Hardy-Weinberg calculation