In this paper, we present a reliable algorithm to study the well-known model of nonlinear dispersive waves proposed by Boussinesq. We solve the Cauchy problem of Boussinesq equation using variational iteration method (VIM). The numerical results of this method are compared with the exact solution of an artificial model to show the efficiency of the method. The approximate solutions show that VIM is a powerful mathematical tool for solving nonlinear problems.