Compute the probability distribution of gene tree histories and gene tree topolgies for a fixed three-taxon species tree under the coalescent model with gene flow between both pairs of sister population.