likelihood_coal_var function

Likelihood of a birth-death model using a coalescent approch