“…Now a days, much attention has been given to the literature of the stable methods for the numerical solution of Benjamina-Bona-Mohany (BBM) equations. In addition to that, Considerable efforts have been made by many mathematicians to obtain exact and approximate solutions of partial differential equations such as Benjamina-Bona-Mohany equations and a number of efficient, accurate and powerful methods have been developed by those mathematicians such as, Backlund transformation method [3], Lie group method [4], Adomian's decomposition method [5], Integral method [6], Hirota's bilinear method [7], homotopy analysis method [8], He's Homotopy perturbation method [9], Exp-Function method [10], Haar wavelet method [11] and Cardinal B-Spline wavelets method [12]. The aim of the present work is to develop Laguerre wavelets collocation method, mutually for solving partial differential equations with initial and boundary conditions of the BBM equations, which is simple, fast and guarantees the necessary accuracy for a relative small number of grid points.…”