#title 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")