code 41502;
real procedure NORMAL(X, MU, SIGMA);
    value X, MU, SIGMA; real X, MU, SIGMA;
NORMAL:= if SIGMA ≤ 0
         then STATAL3 ERROR(“NORMAL”, 3, SIGMA)
         else PHI((X - MU) / SIGMA);
eop