This paper studies the smoothness of solutions of the higher dimensional polynomial-like iterative equation. The methods given by Zhang Weinian[7] and by Kulczycki M, Tabor J.[3] are improved by constructing a new operator for the structure of the equation in order to apply fixed point theorems. Existence, uniqueness and stability of continuously differentiable solutions are given.