MathGroup Archive 2009

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

Search the Archive

Re: WorkingPrecision and NDSolve

  • To: mathgroup at
  • Subject: [mg104206] Re: [mg104036] WorkingPrecision and NDSolve
  • From: Dan Dubin <ddubin at>
  • Date: Thu, 22 Oct 2009 22:36:32 -0400 (EDT)

>Why does this instruction properly work if precision =60 but not if
>precision=70  (Message : error test failure at t=0.01316... ) :
>NDSolve[{x''[t] + x'[t] + Sin[x[t]] == Sin[t/3], x[0] == 0, x'[0] == 0},
>x, {t, 0, 10}, WorkingPrecision -> 70]         ?
>Thanks in advance,
>A. Hautot

Andre, I have been wondering the same thing, which I too observed 
when solving other ODEs. Evidently if WorkingPrecision is set too 
high in NDSolve, the method fails. Do the good people at WRI have any 
comments on this issue? It looks like a bug to me.

  • Prev by Date: Timer event
  • Next by Date: Calendar questions
  • Previous by thread: WorkingPrecision and NDSolve
  • Next by thread: How to find pattern with Mathematica?