MathGroup Archive 2007

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

Search the Archive

Re: Mathematica 6.0 bug in computing MathieuC

  • To: mathgroup at smc.vnet.net
  • Subject: [mg80403] Re: Mathematica 6.0 bug in computing MathieuC
  • From: T.bakri at gmail.com
  • Date: Tue, 21 Aug 2007 05:07:20 -0400 (EDT)
  • References: <fabgsc$ien$1@smc.vnet.net>

On Aug 20, 9:46 am, chuck009 <dmili... at comcast.com> wrote:
> > Dear all,
> > Something is still worng with the algorithm behind
> > computing MathieuC. It
> > reprts, for instance that MathieuC is pi-periodic
> > when its not.
>
> Hey T, I was unable to cause the following code to fail to give a pi-periodic function "poly" for any natural number r and parameter q:
>
> r = 4;
> q = -3 - 3*I
> a = N[MathieuCharacteristicA[r, q]]
> poly = MathieuC[a, q, x]/Exp[I*r*x];
> p1 = Plot[Im[poly], {x, 0, 10}];
> linelist = (Line[{{#1*Pi, -3}, {#1*Pi, 3}}] & ) /@ Range[1, 3]
> Show[{p1, Graphics[linelist]}]
>
> where the division by Exp[irx] is a consequence of Floquet's Theorem.

Take the parameter c such that MathieuCharacteristicB[2, c]=1.  One
can easilty numerically compute c. We find
c=-6.433074585014675206670244658478726486156264965498297244484.
Clearly, we have MathieuS(1,c,t) is Pi-periodic. However
MathieuC(1,c,t) should oscillate but is NOT Pi-periodic. Moreover the
Wronski-determinant ( i.e. MathieuS(1,c,t)*MathieuCPrime(1,c,t) -
MathieuC(1,c,t)*MathieuSPrime(1,c,t)) is in this case time-independent
and must be equal to -1. Mathematica 6.0 does this all worng. You can
check it yourself.

Regards,
T. Bakri



  • Prev by Date: Re: Manipulate with Wrapping SetterBars
  • Next by Date: Re: listplot and legend (again)
  • Previous by thread: Re: Mathematica 6.0 bug in computing MathieuC
  • Next by thread: keyboard help navigation