MathGroup Archive 2008

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

Search the Archive

Re: A Problem with Simplify

This seems to be very strange:

Mathematica 6.0 for Microsoft Windows (32-bit)
Copyright 1988-2008 Wolfram Research, Inc.

In[1]:= Imn = Integrate[Sin[(m*Pi*x)/L]*Sin[(n*Pi*x)/L], {x, 0, L}];

In[2]:= Imn // InputForm

Out[2]//InputForm= (L*n*Cos[n*Pi]*Sin[m*Pi] - L*m*Cos[m*Pi]*Sin[n*Pi])/
(m^2*Pi - n^2*Pi)

In[3]:= Simplify[Imn, Assumptions -> n == m]

                           L n Cos[n Pi] Sin[m Pi] - L m Cos[m Pi]
Sin[n Pi]
Simplify::infd: Expression
------------------------------------------------- simplified to
                                              2       2
                                             m  Pi - n  Pi

Out[3]= Indeterminate


In[4]:= Integrate[Sin[(n*Pi*x)/L]*Sin[(n*Pi*x)/L], {x, 0, L}] //

Out[4]//InputForm= (L*(2 - Sin[2*n*Pi]/(n*Pi)))/4

This is a bug! :(

Recently I also have found another strange bug in Simplfy (not so
gives Sin[x]*Cos[x] instead Sin[2*x]/2. The workaround is

Kevin J. McCann 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.
> The following is just fine
> Imn = Simplify[Integrate[
>         Sin[(m*Pi*x)/L]*
>           Sin[(n*Pi*x)/L],
>         {x, 0, L}]]
> However, if I specify that m and n are integers, I only get the
> "general" solution of zero, i.e. when m and n are not equal.
> Imn = Simplify[Integrate[
>         Sin[(m*Pi*x)/L]*
>           Sin[(n*Pi*x)/L],
>         {x, 0, L}],
>       Element[m, Integers] &&
>         Element[n, Integers]]
> The workaround is obvious in this case, but shouldn't Mathematica give multiple
> answers? Perhaps something similar to what it already does with Integrate?
> Kevin
> --
> Kevin J. McCann
> Research Associate Professor
> JCET/Physics
> Physics Building
> University of Maryland, Baltimore County
> 1000 Hilltop Circle
> Baltimore, MD 21250

  • Prev by Date: Re: pdf and accents
  • Next by Date: Re: A Problem with Simplify
  • Previous by thread: Re: A Problem with Simplify
  • Next by thread: Re: Re: A Problem with Simplify