code 41541;
real procedure CAUCHY(X, LOC, SCALE);
value X, LOC, SCALE; real X, LOC, SCALE;
CAUCHY:= if SCALE ≤ 0
    then STATAL3 ERROR(“CAUCHY”, 3, SCALE)
    else ARCTAN((X - LOC) / SCALE) × .31830988618379 + .5;
eop