In this paper, we present a way of applying the so-called He's variational iteration method (VIM) to numerically solve the non linear autonomous third-order ordinary differential equation (ODE) y ′′′ = y −2 obtained by considering a traveling wave solution admitted by a lubrication equation modeling a twodimensional spreading of a thin viscous film on a inclined slope. Approximate analytical solution is derived and compared to the results obtained from the Adomian decomposition method (ADM) proposed in [20], to the exact analytical solution [7,8], to a fifth order Runge-Kutta method (DOPRI), a fourth order Runge-Kutta method (RK4), a three-stage fifth order Runge-Kutta method (RKD5) developed in [18]. A very good agreement and accuracy is observed. Comparisons are obtained using symbolic capabilities of Maple 18.0 package.