A three-dimensional time-domain panel method is developed, which solves the Neumann-Kelvin linearised ship motions problem. The initial boundary value problem is expressed as an integral equation using the transient free surface Green’s function. This integral equation is discretised using plane elements, on which the potential assumed to be constant. The memory part of the Green’s function and its spatial derivatives are evaluated using fourth order ordinary differential equations, including forward speed effects. It is shown that this method is much more efficient, and as accurate, as using classical quadrature methods. The effect of the forward speed on the evaluation of the memory part of the Green’s function is also investigated. I...