Abstract: The precise integration method is an accurate numerical method for solving first-order differentialequations． The fundamental governing equations of layered elastic foundation were introduced into the duality system，and governing equations were represented as Hamilton canonical equations in the frequency-wave number domain．Based on which，a new method was proposed to solve dynamic flexibility of layered elastic foundation． The segmentmatrix of micro layer was obtained by using the precise integration method，and then，through merging segment matrixof micro layer，it was evaluated the segment matrix of the whole foundation． Finally，the dynamic flexibility of layeredelastic foundation was calculated by incorporating the radiation condition of boundary into the segment matrix of thewhole foundation． In the end，accuracy of the method was verified by numerical examples．