Prev: Region borders in WorldPlot
Next: WorkingPrecision in ODE
From: Andre Hautot on 16 Oct 2009 07:03
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


From: Dan Dubin on 22 Oct 2009 22:41
>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.
--

 | 
Pages: 1
Prev: Region borders in WorldPlot
Next: WorkingPrecision in ODE