Group Abstract Group Abstract

Message Boards Message Boards

0
|
5.1K Views
|
9 Replies
|
3 Total Likes
View groups...
Share
Share this post:

Different results solving the same integral?

Enrique Rodriguez
Posted 7 years ago

I need to evaluate an integral as in the attached file, for general parameters a, b and c. I obtain a conditional 0 as result using the Integrate function.

However, when I evaluate the integral again (numerically and analitically) substituing the parameters a,b and c by numbers (which fulfill the conditions obtained before), the result is not zero.

Should I trust the result 0 for general parameters a, b and c, or not? What could be happening?

Thanks so much

Attachments:
POSTED BY: Enrique Rodriguez
9 Replies
Sort By: 
In[2]:= (I (MeijerG[{{1/2, 1, 1}, {}}, {{1/2, 1/2, 1}, {}}, -((
      I Sqrt[a])/b), 1/2] - 
    MeijerG[{{1/2, 1, 1}, {}}, {{1/2, 1/2, 1}, {}}, (I Sqrt[a])/b, 1/
     2]))/(2 Sqrt[a] \[Pi]^(3/2) (b + c + I ww)) /. {a -> 1.1, 
  b -> 0.5, c -> 0.5, ww -> 1.1}

Out[2]= 1.05409 - 1.1595 I
POSTED BY: Alexander Trounev
POSTED BY: Mariusz Iwaniuk

We need to make an assumption for the integral. For example,

In[41]:= Integrate[
 1/Sqrt[(w^2 + a)]/((b - I*w)*(c + I*(w + ww))), {w, -Infinity, 
  Infinity}, Assumptions -> {a > 0, c > 0, b > 0, ww > 0}]

Out[41]= (I (MeijerG[{{1/2, 1, 1}, {}}, {{1/2, 1/2, 1}, {}}, -((
     I Sqrt[a])/b), 1/2] - 
   MeijerG[{{1/2, 1, 1}, {}}, {{1/2, 1/2, 1}, {}}, (I Sqrt[a])/b, 1/
    2]))/(2 Sqrt[a] \[Pi]^(3/2) (b + c + I ww))
POSTED BY: Alexander Trounev

Thank you Mariusz. It's very good. I know that Mathematica contains the solution to any problem.
POSTED BY: Alexander Trounev

Challenge accomplished. See my attached file Support v2.nb

Regards, MI.
POSTED BY: Mariusz Iwaniuk

Hi

I found a simple workaround to your problem see attched file.

Regards,MI

Attachments:
POSTED BY: Mariusz Iwaniuk

This is a challenge, we need to get around this obstacle.
POSTED BY: Alexander Trounev

I send e-mail to WRI support with this issue and they gave me the answer:

  Hello,

  Thank you for contacting Wolfram Technical Support. I understand that the following is not returning the correct result when 
  compared with the NIntegrate result. 

  Integrate[1/(Sqrt[1+w^2] ((1+I)+w+w^2)),{w,-Infinity,Infinity}]

  The fundamental problem seems to be with the following term in the Integrate result:

  MeijerG[{{1/2,1,1},{}},{{1/2,1/2,1},{}},-I,1/2]

  For some reason, Mathematica is not able to evaluate this numerically. It just returns unevaluated. Also, if you try to evaluate 
  this using arbitrary-precision, you get a bunch of errors. I have filed a report with our developers regarding this behavior. Thank 
  you for bringing this to our attention. 

  Regards,

  Luke Titus
  Wolfram Technical Support
  Wolfram Research Inc.
  http://support.wolfram.com
POSTED BY: Mariusz Iwaniuk

You may have hit a singular point. With the values {ww -> 1., a -> 2., b -> 1., c -> 1.} we get a numerical value.
POSTED BY: Gianluca Gorni
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Reply Preview
Attachments
Remove
or Discard