fit_hbd_pdr_parametric function

Fit parameterized pulled diversification rates of birth-death models.