[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
Prev by Date:
**Re: large integration result for simple problem: 1/x,, also BesselJ**
Next by Date:
**Re: large integration result for simple problem: 1/x,, also BesselJ**
Previous by thread:
**Re: large integration result for simple problem: 1/x,, also BesselJ**
Next by thread:
**Re: large integration result for simple problem: 1/x,, also BesselJ**
| |