In this paper, a matrix formulated approach is presented for numerical solution of a generalized wave equation with a nonlocal boundary condition. Some theoretical results are developed to simplify and reduce computational costs. Also, convergence analysis and error estimation of the presented method are investigated. Finally, some illustrative examples are compared with the closed-form solutions and existing numerical methods to demonstrate the efficiency and accuracy of the proposed method.