Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2009

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Problem with an integral

  • To: mathgroup at smc.vnet.net
  • Subject: [mg95645] Re: [mg95571] Problem with an integral
  • From: "David Park" <djmpark at comcast.net>
  • Date: Fri, 23 Jan 2009 05:09:56 -0500 (EST)
  • References: <18607536.1232540493860.JavaMail.root@m02>

When I see an expression like:

fun = x^(n + 1)*E^(-x + (I*k)/x)

it always makes me cringe a little. That is because it has parameters and a
variable and it is always better to make these a part of the definition. So
I would write:

fun[n_, k_][x_] := x^(n + 1)*E^(-x + (I*k)/x)

Then if I wanted to see precisely what I was doing, I would use the
Student's Integral section of the Presentation package.

Needs["Presentations`Master`"]

integrate[fun[1, 1][x], {x, 0, \[Infinity]}]
% // UseIntegrate[]
% // N

giving...

\!\(
\*SubsuperscriptBox[\(\[Integral]\), \(0\), \(\[Infinity]\)]\(\(
\*SuperscriptBox[\(E\), \(
\*FractionBox[\(I\), \(x\)] - x\)]\ 
\*SuperscriptBox[\(x\), \(2\)]\) \[DifferentialD]x\)\)
4 (MeijerG[{{}, {}}, {{0, 3/2, 2}, {1/2}}, 1/16] + 
   I MeijerG[{{}, {}}, {{1/2, 3/2, 2}, {0}}, 1/16])
1.67562 + 0.832688 I

If one wanted to do simplifications on the integrand or use other
integration techniques after the first step, that could also be done.


David Park
djmpark at comcast.net
http://home.comcast.net/~djmpark/  




From: Jepessen [mailto:jepessen at gmail.com] 

Hi.

I'm working with a little problem in Mathematica 7.0.0.
I want to integrate this function

fun = x^(n + 1)*E^(-x + (I*k)/x)

Assuming that both n and k are greater than zero, I write

integral = FullSimplify[Assuming[{k > 0, n > 0}, Integrate[fun, {x, 0,
\[Infinity]}]]]

And I obtain a symbolic result.
But, when I want to put some specific value, like this

integral /. {n -> 1, k -> 1}

I obtain always ComplexInfinity, and/or other errors.
So, I've tried to evaluate numerically the integral for the same
specific values, in this way

NIntegrate[Evaluate[fun /. {n -> 1, k -> 1}], {x, 0, \[Infinity]}]

And I obtain a finite numeric result. So, there's some error in
symbolic computation, or I miss something when I try to integrate the
formula?

Thanks for answers




  • Prev by Date: Re: message-driven function: more explanation
  • Next by Date: Re: Re: Re: Which editor do you use for math articles
  • Previous by thread: Re: Problem with an integral
  • Next by thread: message-driven function: more explanation