function f=logn(x,eta,sigma) %lognormal pdf with mean eta and standard deviation sigma-see (2.26) of fincnacne notes. f=exp(-((log(x/eta)-sigma^2/2).^2/(2*sigma^2))); f=f./(x*sigma*sqrt(2*pi));