We describe a pseudo-spectral numerical method to solve the systems of one-dimensional evolution equations for free surface waves in a homogeneous layer of an ideal fluid. We use the method to solve a system of one-dimensional integro-differential equations, first proposed by Ovsjannikov and later derived by Dyachenko, Zakharov, and Kuznetsov, to simulate the exact evolution of nonlinear free surface waves governed by the two-dimensional Euler equations. These equations are written in the transformed plane where the free surface is mapped onto a flat surface and do not require the common assumption that the waves have small amplitude used in deriving the weakly nonlinear Korteweg-de Vries and Boussinesq long-wave equations. We compare the solution of the exact reduced equations with these weakly nonlinear long-wave models and with the nonlinear long-wave equations of Su and Gardner that do not assume the waves have small amplitude. The Su and Gardner solutions are in remarkably close agreement with the exact Euler solutions for large amplitude solitary wave interactions while the interactions of low-amplitude solitary waves of all four models agree. The simulations demonstrate that our method is an efficient and accurate approach to integrate all of these equations and conserves the mass, momentum, and energy of the Euler equations over very long simulations.