A parametric scheme is proposed for the numerical solution of the nonlinear Boussinesq equation. The numerical method is developed by approximating the time and the space partial derivatives by finite-difference re placements and the nonlinear term by an appropriate linearized scheme. The resulting finite-difference method is analyzed for local truncation error and stability. The results of a number of numerical experiments are given for both the single and the double-soliton wave.