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)