This paper studies Gardner equation, which represents long nonlinear internal waves. The collocation method based on B-splines is applied to the equation. The stability of the proposed numerical scheme is analyzed by using von Neumann theory. To observe some physical properties of long nonlinear internal waves, three test problems which contain the propagation of solitary waves, the interaction of solitary waves and evolution of solitons are considered. Also, the effect of nonlinearity on physical problems is investigated. In order to see this effect clearly, the same parameters are used during the computation for different degrees of nonlinearity in each problem.