Incomplete Gamma Function of the Second Kind

[toml_12953]

I'm trying to compute the Gamma function by computing the incomplete gamma function of the first kind and the second kind and adding them together.

(no jokes about Close Encounters of the Third Kind, please! )

The first kind computes with no problem but the second kind gives me an error:

"This function works only in radian mode Error: Bad Argument Value"

It doesn't use any trig functions that I can see so why the error? I know it does work in radian mode but I'm trying to understand what's happening.

Code:

EXPORT γ(A,X)
BEGIN
  RETURN approx(int(T^(A-1)*e^(-T),T,0,X));
END;

EXPORT Γ(A,X)
BEGIN
  RETURN approx(int(T^(A-1)*e^(-T),T,X,∞));
END;

EXPORT Γ(A)
BEGIN
  RETURN approx(γ(A,30) + Γ(A,30));
END;

[Albert Chan]

(05-25-2021 10:26 PM)toml_12953 Wrote:  "This function works only in radian mode Error: Bad Argument Value"

It doesn't use any trig functions that I can see so why the error?

Depending on user input, CAS might not recognize it related to gamma function.
It may tried to remove infinity by substitution, x = tan(t), so that x→∞ map to t→pi/2

CAS> gamma1(a,x) := int(t^(a-1)*e^-t, t=x .. inf)

CAS> gamma1(1, i)       → e^-i
CAS> gamma1(1.0, i)    → integrate(tan(t)^0.*e^(-tan(t))*(1+tan(t)^2), t, atan(i), π/2)

Actually, x=tan(t) here just shift the infinity ... atan(i)

Does atan(i) even have a limit ? CAS seems to think it is just a value.

CAS> limit(atan(x),x=i)    → ATAN(i)

Ref: DLMF: Chapter 8 Incomplete Gamma and Related Functions