loglikelihood_hbd function

Galculate the log-likelihood of a homogenous birth-death model.