We consider the problem of a ship advancing in waves. In this method, the zone of free surface in the vicinity of body is discretized. On the discretized surface, the first-order and second-order derivatives of ship waves are represented by the B-Spline formulae. Different ship waves are approximated by cubic B-spline and the first and second order derivates of incident waves are calculated and compared with analytical value. It approves that this numerical method has sufficient accuracy and can be also applied to approximate the velocity potential on the free surface.