A semilinear pseudoparabolic equation with nonlocal integral boundary conditions is studied in the present paper. Using Rothe method, which is based on backward Euler finitedifference schema, we designed a suitable semidiscretization in time to approximate the original problem by a sequence of standard elliptic problems. The questions of convergence of the approximation scheme as well as the existence and uniqueness of the solution are investigated. Moreover, the Legendre pseudospectral method is employed to discretize the time-discrete approximation scheme in the space direction. The main advantage of the proposed approach lies in the fact that the full-discretization schema leads to a symmetric linear algebraic system, which may be useful for theoretical and practical reasons. Finally, numerical experiments are included to illustrate the effectiveness and robustness of the presented algorithm.