In this paper, numerical solution of the initial-boundary value problem for Benjamin-Bona-Mahony equation under homogeneous boundary is considered. A linearized three-level difference scheme is introduced. The difference scheme simulates a conservative quantity of the problem. Furthermore, convergence and stability of the difference scheme are proved by discrete energy method. Numerical experiments verify the reliability of this method.
