longfloat library HP PRIME

[robmio]

Thanks for the reply. Part of the problem is related to the high values to use with the Gamma function. One method is to use the LogGamma function (see the code below) and translate the result with Exp (LogGamm). It works. There remains the problem of the incomplete Beta function, which cannot be "transformed" into other functions (at least I think).
Best regards, Roberto.

Code:

#cas
LogGamma(ασ):=
BEGIN
RETURN int(Psi(x),x,1,ασ);
END;
#end