R Distribution

x<-scan("data")

hist(x,breaks=100)

plot(density(x))

plot(ecdf(x))

library(MASS)

fitdistr(x,"gamma")