How do I Show Equal

[Thomas D. Dean]

I am trying to duplicate a Maxima proof.

#######################################
# Theorum
# for any real k,a,b
# integrate(exp(x)*x^k,x=a..b)
# =exp(1)^(-I*Pi*k)*(gamma_incomplete(k+1,-b)
# - gamma_incomplete(k+1,-a))
# maxima proof
# S:'integrate(exp(x)*x^k,x)';
# S1:ev(S, nouns);
# subst(-x=exp(%i*pi)*x,S1);
# /* From Newton-Leibnitz formula S is equal */
# subst(x=b,%)-subst(x=a,%);
# radcan(%);
# subst(pi=%pi,%);
restart;
S:=integrate(exp(x)*x^k,x); # this seems to include the ev(S, nouns)
soln:=exp(1)^(-I*Pi*k)*(GAMMA(k+1,-b)-GAMMA(k+1,-a));
# from I should be able to reduce to soln
subs(x=b,S) - subs(x=a,S);

But, my lack of Maple expertise prevents reduction of this expression to
soln.

How do I do this?

Tom Dean

[Axel Vogt]

Int(exp(x)*x^k,x); value(%);
eval(%, x=b) - eval(%,x=a); # need not be correct to do that!
simplify(%, symbolic);
simplify(%, size);

gives

-k*(GAMMA(k,-a)-GAMMA(k,-b))*(-1)^(-k)+b^k*exp(b)-a^k*exp(a);

The 2-argument GAMMA is Maple's incomplete Gamma function
and you have to compare the definition in both systems to
check why there