The case of zero determinant was not treated correctly when dealing
with N*N matrices with N>4

For these matrices, a Gauss-pivot is used but it was incorrect in the
singular case