MathGroup Archive 2008

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

Search the Archive

Re: Re: A Problem with Simplify

  • To: mathgroup at smc.vnet.net
  • Subject: [mg87647] Re: [mg87603] Re: A Problem with Simplify
  • From: Daniel Lichtblau <danl at wolfram.com>
  • Date: Tue, 15 Apr 2008 05:48:31 -0400 (EDT)
  • References: <ftkb7f$a9m$1@smc.vnet.net> <200804140943.FAA08090@smc.vnet.net>

Alexey Popkov wrote:
> On Apr 10, 10:14 am, "Kevin J. McCann" <Kevin.McC... at umbc.edu> wrote:
> 
>>I have the following rather simple integral of two sines, which should
>>evaluate to zero if m is not equal to n and to L/2 if they are the same.
> 
> 
> Try the following:
> Integrate[Exp[(a - 1)*x], x] /. a -> 1
> Integrate[Cos[(a - 1)*x], x] /. a -> 1
> Integrate[(a - 1)^x, {x, -1, 0}] /. a -> 1
> Integrate[Cos[a x]/Sin[x], x] /. a -> 1
> 
> There is the ONE underlying BUG! In some complicated cases this bug
> may result in random partial answers.
> 
> http://forum.ru-board.com/topic.cgi?forum=5&topic=10291&start=80#9

Removing the replacements, here are the Integrate results.

In[17]:= InputForm[Integrate[Exp[(a - 1)*x], x]]
Out[17]//InputForm= E^((-1 + a)*x)/(-1 + a)

In[18]:= InputForm[Integrate[Cos[(a - 1)*x], x]]
Out[18]//InputForm= Sin[(-1 + a)*x]/(-1 + a)

In[19]:= InputForm[Integrate[(a - 1)^x, {x, -1, 0}]]
Out[19]//InputForm= (-2 + a)/((-1 + a)*Log[-1 + a])

In[20]:= InputForm[Integrate[Cos[a*x]/Sin[x], x]]
Out[20]//InputForm=
((1 + a)*Hypergeometric2F1[1/2 - a/2, 1, 3/2 - a/2, E^((2*I)*x)] -
   (-1 + a)*E^((2*I)*a*x)*Hypergeometric2F1[(1 + a)/2, 1, (3 + a)/2,
     E^((2*I)*x)])/((-1 + a^2)*E^(I*(-1 + a)*x))

As far as I am aware thse are correct. What is the bug?

Daniel Lichtblau
Wolfram Research





  • Prev by Date: Re: Mathematica Player Pro!
  • Next by Date: Re: Numerical integration and list of points
  • Previous by thread: Re: A Problem with Simplify
  • Next by thread: Re: A Problem with Simplify