function f=lognpdf(x,eta,sigma)
%lognormal p.d.f.  for mean eta and vaol parameter sigma
f=(log(x/eta)-0.5*sigma^2).^2;
f=(1./(sigma*sqrt(2*pi)*x)).*exp(-f/(2*sigma^2));