code 41754;
real procedure LOGNORMALDENS(X, MU, SIGMA);
    value X, MU, SIGMA; real X, MU, SIGMA;
LOGNORMALDENS:= if SIGMA ≤ 0
                then
                STATAL3 ERROR(“LOGNORMALDENS”, 3, SIGMA)
                else if X ≤ 0 then 0
                else
                  EXP(-(((LN(X) - MU) / SIGMA) ⭡ 2) / 2)
                       × .39894228040143 / X / SIGMA;
eop