Abstract: The method for analyzing a single pile and pile groups in non-homogeneous soil is presented. The pile is discretized into cylinder and circular elements, and the corresponding soil is also discretized into the elements of the same amount and type. The equations of stress and displacements are established for piles and soil, respectively. The beam theory is used to establish the pile equation of stress and displacements, and the fundamental solutions of layered solids are used to establish the pile equation. The global equations for obtaining loads and displacements of pile-foundation interaction are presented by using displacement compatibility. The slipping problem of a pile and foundation is also considered. Finally, three examples of a single pile or pile groups are analyzed to demonstrate this numerical method. Because of using the fundamental solutions of layered elastic solids, the proposed method can be used for analysis of the pile-foundation interaction in non-homogeneous soil.