Extensions for more definite integrals with log and exp functions
Subject: Extensions for more definite integrals with log and exp functions
From: Dieter Kaiser
Date: Mon, 13 Jul 2009 23:16:44 +0200
I have already reported some work for more definite integrals. At first
I have committed the following integrals:
1. CASE: Integrals generated by the Gamma function:
integrate(z^w*log(z)^m*exp(-t^s),z,0,inf)
2. CASE: Integrals generated by the Incomplete Gamma function:
integrate(z^w*log(z)^m*exp(-t^s),z,x,inf)
3. CASE: Integrals generated by the Beta function:
integrate(z^r*log(z)^n*(1-z)^s*log(1-z)^m,z,0,1)
The paper I have used reports some more algorithm for even more general
integrals (thanks to R. F. Fateman for the paper). I will have a look at
these integrals too.
This is the result I like most, because on the German wikipedia someone
wrote, that Maxima is not even able to solve this integral:
(%i4) integrate(exp(-t)*log(t),t,0,inf);
(%o4) -%gamma
Because we have implemented the expansion of the Incomplete Gamma
function, we get nice results in terms of the Error and the Exponential
functions too:
(%i21) assume(x>0)$
(%i22) integrate(t^-2*exp(-t^-2),t,x,inf);
(%o22) (sqrt(%pi)-gamma_incomplete(1/2,1/x^2))/2
(%i23) integrate(t^-2*exp(-t^-2),t,x,inf),gamma_expand:true;
(%o23) -sqrt(%pi)*(erfc(1/x)-1)/2
(%i24) integrate(t^-3*exp(-t^-2),t,x,inf);
(%o24) -(gamma_incomplete(1,1/x^2)-1)/2
(%i25) integrate(t^-3*exp(-t^-2),t,x,inf),gamma_expand:true;
(%o25) %e^-(1/x^2)*(%e^(1/x^2)-1)/2
Dieter Kaiser