code 41764;
real procedure LAPLACEDENS(X,LOC,SCALE);
    value X, LOC, SCALE; real X, LOC, SCALE;
begin
    if SCALE ≤ 0
    then STATAL3 ERROR(“LAPLACEDENS”,3,SCALE);
    X := (X - LOC) / SCALE;
    LAPLACEDENS := .5 / SCALE × EXP(-ABS(X));
end LAPLACEDENS;
eop