MathGroup Archive 2013

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

Search the Archive

Integrate bug in v 9.0.0

  • To: mathgroup at smc.vnet.net
  • Subject: [mg129647] Integrate bug in v 9.0.0
  • From: "Alexey Popkov" <lehin.p at gmail.com>
  • Date: Sat, 2 Feb 2013 01:17:27 -0500 (EST)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • Delivered-to: l-mathgroup@wolfram.com
  • Delivered-to: mathgroup-newout@smc.vnet.net
  • Delivered-to: mathgroup-newsend@smc.vnet.net

In version 9.0.0 the following integral is reported as divergent:

In[71]:= Integrate[x^2/(x^2 - x0^2), {x, x1, x2},
 Assumptions -> x0 > x2 > x1 >= 0]

During evaluation of In[71]:= Integrate::idiv: Integral of x^2/(x^2-x0^2) 
does not converge on {x1,x2}. >>

Out[71]= Integrate[x^2/(x^2 - x0^2), {x, x1, x2},
 Assumptions -> x0 > x2 > x1 >= 0]

Versions 8.0.4 and 5.2 give equivalent expressions:

In[2]:= Integrate[x^2/(x^2 - x0^2), {x, x1, x2},
 Assumptions -> x0 > x2 > x1 >= 0]

Out[2]= -x1 + x2 + x0 ArcTanh[(x0 (x1 - x2))/(x0^2 - x1 x2)]

In[9]:=
Integrate[x^2/(x^2-x0^2),{x,x1,x2},Assumptions->x0>x2>x1>=0]
Out[9]=
-x1+x2+1/2 x0 (Log[x0+x1]+Log[x0-x2]-Log[(x0-x1) (x0+x2)])

In[8]:= FullSimplify[-x1 + x2 +
   1/2 x0 (Log[x0 + x1] + Log[x0 - x2] -
      Log[(x0 - x1) (x0 + x2)]) == -x1 + x2 +
   x0 ArcTanh[(x0 (x1 - x2))/(x0^2 - x1 x2)],
 Assumptions -> x0 > x2 > x1 >= 0]

Out[8]= True

Alexey 




  • Prev by Date: Re: steady state concetration
  • Next by Date: Re: Using NDSolve solution as an initial condition for another
  • Previous by thread: Re: steady state concetration
  • Next by thread: Re: Integrate bug in v 9.0.0