Under the assumption of potential flow, free-surface flow around a hydrofoil is calculated by the high-order spectra1!boundary-integral method, This method is one of the most efficient numerical methods by which the nonlinear interactions between hydrofoil and free-surface can be simulated in time-domain. In this method. the wave potential which represents the nonlinear evolution of free-surface is solved by the high-order spectral method and the body potential which provides the effects of hydrofoil and shed vortex is solved by the boundary-integral method. The calculated free-surface profiles which are generated by a uniformly translating hydrofoil are compared with other experimental results. And they show relatively good agreements each other. As another example, free-surface flow generated by a heaving and translating hydrofoil is calculated and discussed.