Discrete numerical solution of the generalized Boussinesq problem