Re: UH oh, integration of rational functions has a bug.
- To: mathgroup at smc.vnet.net
- Subject: [mg41680] Re: UH oh, integration of rational functions has a bug.
- From: carlos at colorado.edu (Carlos Felippa)
- Date: Fri, 30 May 2003 03:56:25 -0400 (EDT)
- References: <bb4vr4$491$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Richard Fateman <fateman at cs.berkeley.edu> wrote in message news:<bb4vr4$491$1 at smc.vnet.net>...
> Mathematica 4.1
>
> In[1]:= 1 + 5*x + a*x^2 + 10*x^3 + 5*x^4 + x^5
>
>
> Integrate[1/%, x]
>
> %1 /. a -> 10
>
> Integrate[1/%, x]
>
> %1 /. a -> c
>
> Integrate[1/%, x]
>
> InputForm[%] ===>
>
> RootSum[1 + 5*#1 + c*#1^2 + 10*#1^3 + 5*#1^4 + #1^5 & ,
> Log[x - #1]/(5 + 30*#1^2 + 20*#1^3 + 5*#1^4 +
> 2*#1*Integrate`V[1][6]) & ]
>
> WHAT is THIS ^^^^^^^^^^^^^^^ ?? and why did it not
> appear in the identical problem where there was an "a"
> instead of a "c" ??
>
> It is somehow dependent on the sequence of commands.
This might be related to a problem reported in the thread
"Integration Bug" of Dec 2002. My post in that thread was
=============================================================
The Integrate behavior seems to be connected to lexicographic
ordering of the variables. For example (from 4.2):
ClearAll[x,z,a,b];
Print[Integrate[Cos[x+z]/(x+z)^2,{x,-Infinity,Infinity}] //InputForm];
Print[Integrate[Cos[x+z]/(x+z)^2,{z,-Infinity,Infinity}] //InputForm];
Print[Integrate[Cos[x+z]/(x+z)^3,{x,-Infinity,Infinity}] //InputForm];
Print[Integrate[Cos[x+z]/(x+z)^3,{z,-Infinity,Infinity}] //InputForm];
Print[Integrate[Sin[x+a+b]/(x+a+b),{x,-Infinity,Infinity}] //InputForm];
Print[Integrate[Sin[x+a+b]/(x+a+b),{a,-Infinity,Infinity}] //InputForm];
Print[Integrate[Sin[x+a+b]/(x+a+b),{b,-Infinity,Infinity}] //InputForm];
-Pi
0
-CosIntegral[-z]/2 + CosIntegral[z]/2
0
0
Pi - SinIntegral[-b - x] - SinIntegral[b + x]
0