Computing the probability of gene trees concordant with the species tree in the multispecies coalescent. (arXiv:2001.06741v2 [q-bio.PE] UPDATED)