MathGroup Archive 2004

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

Search the Archive

Conflict between Mathematica 4.1 and 5.01?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg46051] Conflict between Mathematica 4.1 and 5.01?
  • From: Hugh Walker <hwalker at gvtc.com>
  • Date: Tue, 3 Feb 2004 18:20:57 -0500 (EST)
  • Sender: owner-wri-mathgroup at wolfram.com

I run under Mac OS X 10.2.8 on an iMac.

The code below was extracted from a package I wrote some time ago using 
Mathematica 4.1. Recently I upgraded to 5.01.

If I execute the last line in 4.1, everything works as intended. 
However, the identical code executed in 5.01 results in numerous fatal 
errors. I find nothing in the 5.01 documentation that indicates any 
such incompatibility. My frustration has boiled over. When pasted in a 
4.1 notebook, this code works; but the same code bombs out when pasted 
in a 5.01 notebook. Can anyone spot the cause of this very strange 
behavior? Can anyone reproduce these results?

The job of the code snippets is to solve an eigenvalue problem by the 
shooting method. The particular example here ought to produce the 
ground state energy of the simple oscillator. Any comments would be 
most welcome.

sol[V_, e_, {q_, qLeft_, qRight_}] :=
   First[psi /. NDSolve[{(-(1/2))*Derivative[2][psi][q] +
         V*psi[q] == e*psi[q], psi[qLeft] == 0,
       Derivative[1][psi][qLeft] == 1}, psi,
      {q, qLeft, qRight}]]

pistol[V_, e_, {q_, qLeft_, qRight_}] :=
   sol[V, e, {q, qLeft, qRight}][qRight]

eigenValue[V_, eStart_, {q_, qLeft_, qRight_}] :=
   e /. FindRoot[pistol[V, e, {q, qLeft, qRight}] == 0,
     {e, eStart, eStart + 0.01*Abs[eStart]}]

eigenValue[0.5*q^2, 0.4, {q, -4., 4.}]


==========
Hugh Walker


  • Prev by Date: RE: How to plot...
  • Next by Date: Re: Re: Nasty bug in Integrate (version 5.0)
  • Previous by thread: Re: arnoldi method with Mathematica
  • Next by thread: Re: Conflict between Mathematica 4.1 and 5.01?