code 41503;
real procedure NORMALINV(PROB, MU, SIGMA);
    value PROB, MU, SIGMA; real PROB, MU, SIGMA;
NORMALINV:= if SIGMA ≤ 0
            then STATAL3 ERROR(“NORMALINV”, 3, SIGMA)
            else if PROB < 10-14 ∨ PROB > 1 - 10-14
            then STATAL3 ERROR(“NORMALINV”, 1, PROB)
            else MU + PHINV(PROB) × SIGMA;
eop