elim.m uses recursive orthogonalization now