This paper describes a method of solid modeling using the integral equation and the three-dimensional poly-harmonic function. This method requires the arbitrary points and lines on the free-form surface of solid model. The free-form surfaces of solid model are searched, after solving the discretized integral equation. These free-form surfaces approximately satisfy the Poisson equation. This method needs many hours to calculate the free-form surface, however it is easy to prepare the input data. The Complicated and smooth free-form surface of solid model can be expressed by Only one patch. In order to investigate the efficiency of this method, several examples are given.