  We are concerned with the tensor equations whose coefficient tensor is an
M-tensor. We first propose a Newton method for solving the equation with a
positive constant term and establish its global and quadratic convergence. Then
we extend the method to solve the equation with a nonnegative constant term and
establish its convergence. At last, we do numerical experiments to test the
proposed methods. The results show that the proposed method is quite efficient.
