fitdistr function

Maximum-likelihood Fitting of Univariate Distributions