
mu1 = 1; mu2=5; mu3 = 10;
x = 0:0.01:2;
y1 = 2*x.*exp(-mu1*x.^2);
y2 = 2*x.*exp(-mu2*x.^2);
y3 = 2*x.*exp(-mu3*x.^2);

plot(x,y1,'r-', x,y2,'b:',  x,y3,'m--')
%following 4 lins should be put after plot()
legend('mu = 1', 'mu = 5','mu = 10')         
title('Reighley')
xlabel('x axis')
ylabel('y axis')
grid
axis([0 2 0 1])