Maxwell-Boltzmann distribution

http://ko.wikipedia.org/wiki/%EB%A7%A5%EC%8A%A4%EC%9B%B0-%EB%B3%BC%EC%B8%A0%EB%A7%8C_%EB%B6%84%ED%8F%AC

mb <- function(x,a) { sqrt(2/pi) * x^2 * exp(-x^2/(2*a^2))/a^3 }

x <- seq(0,20,len=1000)
plot(x,mb(x,1),type="l",lwd=3)
lines(x,mb(x,2),col="red",lwd=3)
lines(x,mb(x,5),col="blue",lwd=3)
legend(12,0.5,legend=c("a=1","a=2","a=5"),col=c("black","red","blue"),lwd=3,bty="n")

install.packages("shotGroups")