He’s variational iteration method [1, 2], which is a modified general Lagrange multiplier method [3], has been shown to solve effectively, easily and accurately a large class of nonlinear problems with approximations which converge quickly to accurate solutions. It was successfully applied to autonomous ordinary differential equations [4], nonlinear partial differential equations with variable ...