Computing the probability of gene trees concordant with the species tree in the multispecies coalescent