code 41759;
real procedure WEIBULLDENS(X,LOC,SCALE,ALPHA);
value X,LOC,SCALE,ALPHA; real X,LOC,SCALE,ALPHA;
begin
  if SCALE ≤ 0
  then STATAL3 ERROR(“WEIBULLDENS”,3,SCALE) else
  if ALPHA ≤ 0 
  then STATAL3 ERROR(“WEIBULLDENS”,4,ALPHA);
  
  WEIBULLDENS:= if X ≤ LOC then 0 else
    (ALPHA / SCALE) × EXP((ALPHA - 1) ×
      LN((X - LOC) / SCALE) - ((X - LOC) / SCALE) ⭡ ALPHA)
end WEIBULLDENS;
eop