
clear

set(0,                           ...
   'defaultaxesfontsize', 20,   ...
   'defaultaxeslinewidth', 1.0, ...
   'defaultlinelinewidth', 1.2, ...
   'defaultpatchlinewidth', 0.7);

k = [0:.01:5];
f=sqrt(k).*exp(-k);
plot(k,f,'linewidth',2)
xlabel('\Lambda','fontsize',20)
ylabel('f(\Lambda)','fontsize',20)

x=[2.55:.001:2.65];
plot(x,erf(x/sqrt(2)))
