In this paper, we will use a modified  variational iteration method (MVIM) for solving an inverse heat conduction problem (IHCP). The approximation of the temperature and the heat flux at  are considered. This method is based on the use of Lagrange multipliers for the identification of optimal values of parameters in a functional in Euclidian space. Applying this technique, a rapid convergent sequence to the exact solution is produced. Moreover, this method does not require any discretization, linearization or small perturbation, thus it can be considered as an efficient method to solve this problem. To show the strength  and capability of this method, some examples are given