“…They include the BBM-Burgers equation, a modification of the Benjamin-Bona-Mahony equation, [12], which includes a dissipative term. On the other hand, most of the work presented in the literature about the numerical approximation of pseudo-parabolic equations seems to be focused on the use of finite differences, [5,6,80,7,26,32,8], as well as finite elements, [9,56], and finite volumes of different type for the discretization in space, sometimes combined with a domain decomposition method for a more accurate approximation of convective and diffusive effects, [4,82]. This different numerical treatment of the advective and diffusive fluxes, mentioned above, also makes influence in some choices of time integrators for pseudo-parabolic problems.…”