Principal pivoting algorithm of Lemke