fit_hbd_model_on_grid function

Fit a homogenous birth-death model on a discrete time grid.