Solver of linear systems with gauss method

Solver of linear systems with gauss method

Linear system resolution function
Gauss function is a function created for the arbitrary size of compatible
linear systems. This function receives in input a square matrix of
dimension nxn and the column vector of known terms [X, P, U, L]=gauss(A,b)
returns 4 results. X is the vector of the system solution P is the vector
of permutations performed by the algorithm U is the upper triangular
matrix and L is the lower triangular matrix.
The output parameters U L and p are optional