Abstract: According to the Winkler foundation beam theory, the foundation function is constantized by the integral median theorem, and a linear equation system is constructed for the entire pile and all nodes through the condition that the internal forces and displacements of the nodes are continuous. This process is implemented by matlab programming, and two examples are used to verify the correctness of the calculation of the response of the horizontally loaded pile body in the case of single-layer foundation and multi-layer foundation. In the calculation example, multiple sets of parameters are changed, which shows that the location of the maximum value of the bending moment of the pile body, the zero point value of the bending moment and the zero point value of displacement are affected by the resistance of the soil layer near the pile top and the soil layer in the middle of the pile under the horizontal load of the pile top. The resistance is controlled by the resistance of the entire soil layer around the pile, and for the horizontal response of the pile body under the same situation, the position where the bending moment zero appears will lag behind the position where the displacement zero appears by about (2~3)d.