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