3s is maximum likelihood program that implements the multispecies coalescent model with and without gene flow.
Yang Z. 2002. Likelihood and Bayes estimation of ancestral population sizes in Hominoids using data from multiple loci. Genetics 162:1811-1823.
Yang Z. 2010. A likelihood ratio test of speciation with gene flow using genomic sequence data. Genom Biol Evol 2:200-211.
Zhu T, Yang Z. 2012. Maximum likelihood implementation of an isolation-with-migration model with three species for testing speciation with gene flow. Mol Biol Evol 29:3131-3142.
Dalquen D, Zhu T, Yang Z. 2017. Maximum likelihood implementation of an isolation-with-migration model for three species. Syst Biol 66:379-398.