Abstract: A general and considerably simple analysis, i.e. presented the propagator matrix method in two systems of vector functions is to solve the layered isotropic elastostatic problem. This approach avoids the complicated nature of the problem. The analytical expressions of displacement and stress components at any point of the medium by general surface loads are obtained in these two systems. Numerical examples given verify that the present formulation is correct, and also show that when the modified Padé approximant technique is used in the numerical integration the present method is particularly efficient and easy to apply in practice.