Abstract
The "good" Boussinesq equation is transformed into a first order differential system. A fourth order finite difference scheme is derived for this system. The resulting scheme is analyzed for accuracy and stability. Newton's method and linearization techniques are used to solve the resulting nonlinear system. The exact solution and the conserved quantity are used to assess the accuracy and the efficiency of the derived method. Head-on and overtaking interactions of two solitons are also considered. The numerical results reveal the good performance of the derived method.