b=2.718281828;x=.08:.08:1000;n1=length(x);
a=100;
for i=1:n1

y(i)=(1/(a*sqrt(pi)))*(a/x(i))^1.5*exp(-a/x(i));
end
%semilogy(x,y)
plot(x,y)
%sum(y)=.9252 pentru n=17000 si.8836 pentru n=7000 si.95o6 pentru n=40000;
%n=1000 da sum(y)=.6986 si sum(y)=.8465;