code 41763;
real procedure CAUCHYDENS(X,LOC,SCALE);
value X,LOC,SCALE; real X,LOC,SCALE;
begin real PI,Q;

  if SCALE ≤ 0
  then STATAL3 ERROR(“CAUCHYDENS”,3,SCALE);
  PI:= 3.1415926535898;
  Q:= (X - LOC) / SCALE;
  CAUCHYDENS:= 1 / (PI × SCALE × (1 + Q × Q))
end CAUCHYDENS;
eop