A system of coupled nonlinear partial differential equations with convective and dispersive terms was modified from Boussinesq-type equations. Through a special formulation, a system of nonlinear partial differential equations was solved alternately and explicitly in time without linearizing the nonlinearity. Coupled compact schemes of sixth order accuracy in space were developed to obtain numerical solutions. Within couple compact schemes, variables and their first and second derivatives were solved altogether. The sixth order accuracy in space is achieved with a memory-saving arrangement of state variables so that the linear system is banded instead of blocked. This facilitates solving very large systems. The efficiency, simplicity, and accuracy make this coupled compact method viable as variational and weighted residual methods. Results were compared with exact solutions which were obtained via devised forcing terms. Error analyses were carried out, and the sixth order convergence in space and second order convergence in time were demonstrated. Long time integration was also studied to show stability and error convergence rates.