thank you for this response, I think that this is what I want to use, I will check it.
may be what you want is LUDecomposition
LUDecomposition
{lu, p, c} = LUDecomposition[H]; u = lu SparseArray[{i_, j_} /; j >= i -> 1, {4, 4}] // MatrixForm
