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