This work is concerned with the development and application of the p-adaptive boundary integral equation method (BIEM) in practical elastostatics engineering situations. Some basic concepts inherent to the p-adaptive technique are summarized and discussed. A pseudocode which illustrates the way for generating the p-adaptive system of equations in microcomputers is also provided. Two numerical examples, which show the accuracy of the method discussed herein are included.