Robert Dodier wrote:
>
> I tried the program shown in Michel's message of Feb 20 with different
> values for s (5, 10, or 20, I have a slow computer) and the bump is
In fact my aim was to run the program for a spin as high as perhaps 1000,
but it may prove to be impossible ... In the "classical" limit, hbar = 0
this problem has solution by elliptic functions, so it is interesting to
understand how algebraic things go to elliptic limits, but this requires
apparently very high spin to work out. By the way the physical problem is
discussed by Senitzky in
http://prola.aps.org/abstract/PRA/v3/i1/p421_1
> present for the values I tried. Also, I replaced eigens_by_jacobi by dgeev
> (from LAPACK) and the bump is the same.
Well maybe the glitch is really here and needs to be understood, and
Mathematica is wrong.
>
> I think the next thing to consider is the construction of the matrix M.
> Maybe Mathematica and Maxima are computing different values for
> some radicals or something; just guessing.
Yes, perhaps the matrix itself is not computed exactly the same.
--
Michel Talon