< tolerance instead of == 0 for solsup
trisup
resolG
modified exec in resolG
LU factorisation with dolittle
LU factorisation with Gauss
Ax=b resolution with LU factorisation
solve Lx=b where L lower triangular
computes the inverse of a matrix with Gauss