Services & ResourcesWolfram ForumsMathGroup Archive

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

Re: large integration result for simple problem: 1/x,, also BesselJ

• To: mathgroup at smc.vnet.net
• Subject: [mg122843] Re: large integration result for simple problem: 1/x,, also BesselJ
• From: Andrzej Kozlowski <akoz at mimuw.edu.pl>
• Date: Sat, 12 Nov 2011 07:33:51 -0500 (EST)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com
• References: <201111110955.EAA08514@smc.vnet.net> <7999FCC4-7E77-4ED4-AEF6-AF3CCFAD0FA7@mimuw.edu.pl> <4EBD4CAB.1090505@eecs.berkeley.edu> <0CA89E70-9E70-4FC8-9E5D-4C16A4934A1F@mimuw.edu.pl> <4EBD5686.5030106@eecs.berkeley.edu> <0FB5E525-5B36-44A8-BB3C-D885A4902D1B@mimuw.edu.pl> <4EBD5BAF.1080300@eecs.berkeley.edu>

On 11 Nov 2011, at 18:30, Richard Fateman wrote:

> On 11/11/2011 9:16 AM, Andrzej Kozlowski wrote:
>> On 11 Nov 2011, at 18:08, Richard Fateman wrote:
>>
>>> On 11/11/2011 8:38 AM, Andrzej Kozlowski wrote:
>>>> Mathematica 8 returns:
>>>>
>>>>
>>>> Integrate[BesselJ[n, b*x], {x, 0, Infinity},
>>>>  Assumptions ->   {Re[n]>   -1}]
>>>>
>>>>  b^(n - 2)*(b^2)^(1/2 - n/2)
>>>>
>>>> Andrzej Kozlowski
>>> So by your previous note, this answer from version 8.0  is wrong 
since it does not exclude Im[b]==0.
>>> I note that the formula is also wrong unless it somehow excludes 
b==0, when the integral is infinite,
>>> but the formula is indeterminate.
>> I assume you must mean that it is wrong because it does not exclude 
Im[b]!=0 (in other words, a non-real number).
> yes, I mistyped that.  With your correction, the bug is now correctly 
reported.  I assume you agree that it is a bug, and that furthermore the 
response to Integrate[1/x,{x,a,b}]  is miserable.
>


In my opinion it is a bug, albeit a minor one. Mathematica 8 answer to 
the other integral is different form that in v. 7:

Integrate[1/x, {x, a, b}]

ConditionalExpression[Log[b] - Log[a], ((Im[a] >= Im[b] && Re[a]*Im[b] 
<= Im[a]*Re[b]) ||
    (Re[a]*Im[b] >= Im[a]*Re[b] && Im[a] <= Im[b])) && ((Re[a/(b - 
a)] >= 0 && a^2 != a*b) ||
    NotElement[a/(a - b), Reals] || Re[a/(a - b)] >= 1)]

I looks to me like a reasonable answer. ConditionalExpression has now 
replaced If, and it seems to be better integrated with Mathematica's 
Assumptions mechanisms than is the case with If. For example:

Refine[ConditionalExpression[Sqrt[x^2], x > 0]]

ConditionalExpression[x,x>0]

while

Refine[If[x > 0, Sqrt[x^2]]]

If[x>0,Sqrt[x^2]]


Andrzej

• Follow-Ups:
  • Re: large integration result for simple problem: 1/x,, also BesselJ
    • From: Andrzej Kozlowski <akoz@mimuw.edu.pl>

• References:
  • large integration result for simple problem: 1/x,, also BesselJ
    • From: Richard Fateman <fateman@cs.berkeley.edu>