function x = qlognorm(p,lambda,zeta)%QLOGNORM  The inverse log-normal distribution function% %          x = plognorm(p, lambda, zeta)%        Copyright (c) Anders Holtsbergif any(any(abs(2*p-1)>1))   error('A probability should be 0<=p<=1, please!')endif any(any(zeta<=0))   error('Parameter zeta is wrong')endx = exp(erfinv(2*p-1).*sqrt(2).*zeta + lambda);