Re: WorkingPrecision and NDSolve
- To: mathgroup at smc.vnet.net
- Subject: [mg104206] Re: [mg104036] WorkingPrecision and NDSolve
- From: Dan Dubin <ddubin at ucsd.edu>
- 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.
--