A numerical method to calculate the load-settlement of an axially loaded pile is presented. The method is based on an iterative scheme that tests the global equilibrium conditions and the load transfer (t-z) equations. The scheme is both accurate and computationally very economical.