Abstract: In this paper, the initial functions of Biot consolidation of one layer subsoil and multilayer subsoils are established by using the Laplace and Hankel transformation and the recurrence law of matrices respectively. Then the equations for evaluate the vertical deformation due to the Biot consolidation are derived. These equations may be solved easily by using numerical integration. Only three simultaneous linear algebric equations are to be solved no matter how many layers are involved in the subsoils.