function showLogNormal(mu,sigma)

int = 0.001;
x = int:int:20*mu;
f =    (  1./( x.*sigma*sqrt(2*pi)     )  ).*exp(-(((log(x) - mu).^2)/(2*sigma*sigma)));
plot(x,f)
title(['mu = ',num2str(mu),' sigma = ',num2str(sigma)]);

disp('Select left and right limits of the distribution to calculate probability...');
mousein = ginput(2);

x1 = mousein(1,1);
x2 = mousein(2,1);

ind1 = floor((x1 - x(1))/int); 
x(ind1);

ind2 = floor((x2 - x(1))/int); 
x(ind2);

sum(f(ind1:ind2)*int)





