Summary: In this paper we develop new tensor methods for unconstrained convex optimization, which solve at each iteration an auxiliary problem of minimizing convex multivariate polynomial. We analyze the simplest scheme, based on minimization of a regularized local model of the objective function, and its accelerated version obtained in the framework of estimating sequences. Their rates of convergence are compared with the worst-case lower complexity bounds for corresponding problem classes. Finally, for the third-order methods, we suggest an efficient technique for solving the auxiliary problem, which is based on the recently developed relative smoothness condition [H. H. Bauschke et al., Math. Oper. Res. 42, No. 2, 330–348 (2017; Zbl 1364.90251); H. Lu et al., SIAM J. Optim. 28, No. 1, 333–354 (2018; Zbl 1392.90090)]. With this elaboration, the third-order methods become implementable and very fast. The rate of convergence in terms of the function value for the accelerated third-order scheme reaches the level \(O\left(\frac{1}{k^4}\right)\), where \(k\) is the number of iterations. This is very close to the lower bound of the order \(O\left(\frac{1}{k^5}\right)\), which is also justified in this paper. At the same time, in many important cases the computational cost of one iteration of this method remains on the level typical for the second-order methods.