A parametric scheme for the numerical solution of the Boussinesq equation