Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
1996
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 1996

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

Search the Archive

Re: Integration Constants 2

  • To: mathgroup at smc.vnet.net
  • Subject: [mg5187] Re: Integration Constants 2
  • From: Daniel Lichtblau <danl>
  • Date: Sat, 9 Nov 1996 02:24:00 -0500
  • Organization: wolfram.com
  • Sender: owner-wri-mathgroup at wolfram.com

Dennis Michael Winters wrote:
> 
> I posted a message recently (Integration...Constants) that received
> several valuable responses requesting more information.  Well, here is
> the actual problem:
> 
> In[1]:=Integrate[(-2 x+ao)^-2 (-x+bo)^-1, x]
> 
> Mathematica gives the following output:
> 
> Out[1]=-(1/((ao-2 bo)(ao- 2 x))) - (Log(-bo+x)/(ao^2 - 4 ao bo + 4 bo^2)) +
> 
> Log(-ao + 2 x)/(ao^2 - 4 ao bo + 4 bo^2)
> 
> This answer is incorrect.  If performed by hand, the answer is shown to
> have positive values for ao and bo.  This is confirmed by the fact that
> this integration is a classical result of chemistry.  Years of chemical
> reaction rate analysis show the values to be positive.  Again, this
> calculation was performed using the unix version of 2.2.3.  Any help from
> those who have already responded would be most appreciated.
> 
> --
> ----------------------------------
> -Dennis M. Winters               -
> -P.O. BOX 21215                  -
> -EMORY UNIVERSITY                -
> -ATLANTA, GA 30322               -
> -dwint01 at emory.edu               -
> -or winters at euch6h.chem.emory.edu-
> ----------------------------------

Below I show the result from version 3.0, which is mathematically
equivalent to your result. I then show that the integrand and the
differentiated result differ by a constant (which is in fact zero).


Out[4]//InputForm=
  1/((ao - 2*bo)*(-ao + 2*x)) - Log[-bo + x]/(ao - 2*bo)^2 +
   Log[-ao + 2*x]/(ao - 2*bo)^2

In[5]:= ee = (-2 x+ao)^-2 (-x+bo)^-1;

In[6]:= Integrate[ee, x] // InputForm
Out[6]//InputForm=
  1/((ao - 2*bo)*(-ao + 2*x)) - Log[-bo + x]/(ao - 2*bo)^2 +
   Log[-ao + 2*x]/(ao - 2*bo)^2

In[7]:= D[%, x] - ee // Together
Out[7]= 0

I do not understand what you mean by the phrase, "answer ... ha[s]
positive values for ao and bo". The variables ao and bo are parameters.


Daniel Lichtblau
Wolfram Research
danl at wolfram.com


  • Prev by Date: WORLDWIDE MATHEMATICA CONFERENCE
  • Next by Date: Speeding Computations w/ NIntegrate[]
  • Previous by thread: Integration Constants 2
  • Next by thread: Re: LinearSolve[] leads to core dump; why?