A three-dimensional analysis method for a flow with liquid surfaces using the boundary element method (BEM) is developed. When the flow is treated as a potential flow and Bernoulli's equation is applicable the BEM is thought to be suitable for the analysis of flow with liquid surfaces. The method can be applied to sloshing problems in containers with liquid surfaces, and the collision of the liquid surfaces with the container ceilings can also be analyzed.