fit_hbd_psr_on_best_grid_size function

Fit pulled speciation rates of birth-death models on a time grid with optimal size.