HP Forums

Full Version: Incomplete Gamma Function of the Second Kind
You're currently viewing a stripped down version of our content. View the full version with proper formatting.
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! Big Grin)

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;
(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
Reference URL's