MathGroup Archive 2009

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

Search the Archive

WorkingPrecision in ODE

  • To: mathgroup at smc.vnet.net
  • Subject: [mg104040] WorkingPrecision in ODE
  • From: Andre Hautot <ahautot at ulg.ac.be>
  • Date: Fri, 16 Oct 2009 07:21:25 -0400 (EDT)

Who has any idea why

NDSolve[{u''[t] + 1/2 u'[t] + Sin[u[t]] == 3/2 Sin[2/3 t], u[0] == 0, 
u'[0] == 0}, u, {t, 0, 200},
WorkingPrecision -> precision, MaxSteps -> Infinity]

works if precision=50 but not if precision=60 ?  (Error test failure at 
t=0.012758...)

(tested with Mathematica 7.0)

Thanks in advance

A. Hautot (ahautot at ulg.ac.be)



  • Prev by Date: How to find pattern with Mathematica?
  • Next by Date: Re: Streaming Real-Time Stock Quotes in Mathematica
  • Previous by thread: Re: How to find pattern with Mathematica?
  • Next by thread: Crash when first evaluation is ReadList[]