In this talk, we first introduce some theoretical results for the Z-spectrum, and then we propose a projected Newton iteration for computing the nonnegative Z-eigenpairs of nonnegative tensors. A local quadratic convergence is also given. More specially, the proposed algorithm can be employed to solve the tensor equation arising from the multilinear PageRank problem. Numerical experiments are provided to illustrate the effectiveness of the proposed approach.