Approximate Taylor methods for ODEs