fit_hbd_pdr_on_grid function

Fit pulled diversification rates of birth-death models on a time grid.