We describe an algorithm based on Gaussian elimination for solving an n×n system of linear equations. It uses a pivoting strategy which is particularly appropriate whenever the residuals can be made small by solving for fewer than n of the unknowns as may happen, for example, in interpolation problems.