Abstract
This paper presents an axisymmetric numerical integration by using the boundary integral equation and the polyharmonia function. In conventional numerical integration, a given region is divided into several standard regions, for which a rule of approximate integration is available. However, it is difficult to decomposite a boundary region into elementary standard regions. Presented method does not require decomposition. This method requires a boundary geometry of the region and arbitrary internal points. The integral value is calculated after solving the discretized boundary integal equation. In order to investigate the efficiency of this method, several examples are given.