Abstract: In this work, we derive two numerical schemes for
solving a class of nonlinear partial differential equations. The first
method is of second order accuracy in space and time directions, the
scheme is unconditionally stable using Von Neumann stability
analysis, the scheme produced a nonlinear block system where
Newton-s method is used to solve it. The second method is of fourth
order accuracy in space and second order in time. The method is
unconditionally stable and Newton's method is used to solve the
nonlinear block system obtained. The exact single soliton solution
and the conserved quantities are used to assess the accuracy and to
show the robustness of the schemes. The interaction of two solitary
waves for different parameters are also discussed.