code 41530;
real procedure STUDENT(X, DF); value X, DF;
    real X, DF;
begin realIB;

     if DF ≤ 0 then
         STUDENT:= STATAL3 ERROR(“STUDENT”,2,DF)
     else
     begin IB:=
         INCOMPLETE BETA(DF/(DF + X × X),DF/2,0.5,10-12);
         if IB < 0 then IB:= 0
         else if IB > 2 then IB:= 2;
         STUDENT:= if X < 0 then IB / 2 else 1 - IB / 2
     end
end STUDENT;
eop