Fourier image methods for evolution equations of deep-water waves

The velocity potential of the fluid satisfies the Laplace equation with nonlocal boundary conditions on a free surface. This differential problem is transformed to an evolution equation in Fourier variables. The Fourier transform images of boundary functions are approximated by Picard's iterations and the method of lines on meshes related to roots of Hermite polynomials. Due to convolutions of sine and cosine functions the integral terms of Picard's iterations reveal unexpected instabilities for wave numbers in a neighborhood of zero