fit_dist_start function

Find starting values for distribution parameters