fit_hbds_model_parametric function

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