As large waves enter a harbor, during their propagation, the motions a floating body are large and if may even be damaged by waves. This phenomenon may be caused by harbor resonance, resulting from large motion at low wave frequency, which is close to the natural frequency of a vessel. In order to calculate the motion of a floating body in a harbor, it is necessary to use the wave forces containing the body-harbor interference. The simulation program to predict the motions of a floating body by waves in a harbor is developed, and this program is based on the method of velocity potential contiuation method proposed by Ijima and Yoshida The calculated results are shown by the variation of wave frequency, wave angle, and the position of a floating body.