code 41761;
real procedure FISHERDENS(X, DF1, DF2);
value X, DF1, DF2;
    real X, DF1, DF2;
FISHERDENS:= if DF1 ≤ 0
             then STATAL3 ERROR(“FISHERDENS”, 2, DF1)
             else
             if DF2 ≤ 0
             then STATAL3 ERROR(“FISHERDENS”, 3, DF2)
             else
             if X ≤ 0 then 0 else
             EXP(LOGGAMMA((DF1 + DF2) / 2) -
              LOGGAMMA(DF1 / 2) - LOGGAMMA(DF2 / 2) +
               (DF1 × LN(DF1) + DF2 × LN(DF2)) / 2 +
                (DF1 / 2 - 1) × LN(X) - (DF1 + DF2) / 2 ×
                  LN(DF1 × X + DF2));
eop