  Multilinear systems play an important role in scientific calculations of
practical problems. In this paper, we consider a tensor splitting method with a
relaxed Anderson acceleration for solving multilinear systems. The new method
preserves nonnegativity for every iterative step and improves the existing
ones. Furthermore, the convergence analysis of the proposed method is given.
The new algorithm performs effectively for numerical experiments.
