+- HP Forums (https://www.hpmuseum.org/forum)
+-- Forum: HP Calculators (and very old HP Computers) (/forum-3.html)
+--- Forum: HP Prime (/forum-5.html)
+--- Thread: Incomplete Gamma Function of the Second Kind (/thread-17004.html)
Incomplete Gamma Function of the Second Kind - toml_12953 - 05-25-2021 10:26 PM
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;RE: Incomplete Gamma Function of the Second Kind - Albert Chan - 05-26-2021 12:32 AM
(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