fit_distr function

Fit distribution to data