A fourth order in time and second order in space scheme using a finite-difference method is developed for the non-linear Boussinesq equation. For the solution of the resulting non-linear system a predictor-corrector pair is proposed. The method is analyzed for local truncation error and stability. The results of a number of numerical experiments for both the single and the double-soliton waves are given.