Hello, Use > x=rnorm(100, mean=3, sd=1) > library(MASS) >fitdistr(x, "normal") mean sd 2.93666631 0.99673982 (0.09967398) (0.07048015) Hope this helps, Shrieb