fit_hbd_model_parametric function

Fit a parametric homogenous birth-death model to a timetree.