In this paper, the boundary element method is applied to solve the diffraction of waves by multiple vertical cylinders under the assumption of linear wave theory. A numerical analysis by boundary element method is based on Green's theorem and introduced to an integral equation for the fluid velocity potential around the cylinders. The numerical results obtained in this study are compared with the experimental data and the results of the theory using multiple scattering techniques. The comparisons show strong agreement. This numerical analysis method developed by using boundary element method could be used broadly for the design of various offshore structures to be constructed in coastal zones in the future.