The functions on boundary are expanded into wavelet series，and the coefficient matrices are obtained from the values of kernal functions on the boundary using a circulant matrix method. The approximate solution of boundary integral equation is worked out by solving the unknown wavelet coefficients. An example of 2D Laplace equation is given.