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