fit_coal_var function

Fit birth-death model using a coalescent approch