A computer program based on Bishop's simplified method of slices (1954) and capable of analyzing the slope stability of a multilayered soil mass is described. The computer program was specifically developed for analyzing the slope stability of highway bridge approach embankments; however, it can be applied to a broad spectrum of practical slope configurations and bearing capacity problems. Details of the use, applications, and accuracy of the program are presented. Important features of the computer program include a grid type, search operation for locating the critical shear surface and a ledger printout of the forces acting on each individual slice. The latter feature was included so that results of the computer program could be compared to those obtained from manual computations. Pore pressures in the computer program are handled in a manner described by Bishop (1954). Additionally, for seepage cases, infinite slope conditions are assumed and used to simulate a flow net.