Date: Oct 16, 2009 7:04 AM
Author: André Hautot
Subject: WorkingPrecision in ODE

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@ulg.ac.be)