b=2.718281828;x=1:1:38000;
a=73;
for i=1:length(x)

y(i)=(1/(a*sqrt(pi)))*(a/x(i))^1.5*exp(-a/x(i));
end
plot(x,y)