create a decomp.m function in MatLab, that outputs LU decomposition of matrix A, it's condition number and its pivot vector, following the pseocode below:

Use the result of the decomp.m function to create a solve.m function that solves x from Ax=b for Gauss elimintation using the following pseucode:

procedure Gauss(n, (aij), (L)) integer i, j, k, n; real r, rmax, smax, xmult real array (aij)nxn n real array alloc aimx( real array allocate (S)1m for i= l to n do smax 0 for j = 1 to n do Smax max(smax, laijl) end for end for for k1 ton - 1 do rmax 0 for i = k to n do if (r > rmax) then rmax r end if end for for i k + 1 to n do for j = k +1 to n do end for end for end for deallocate array (si) end procedure Gauss procedure Gauss(n, (aij), (L)) integer i, j, k, n; real r, rmax, smax, xmult real array (aij)nxn n real array alloc aimx( real array allocate (S)1m for i= l to n do smax 0 for j = 1 to n do Smax max(smax, laijl) end for end for for k1 ton - 1 do rmax 0 for i = k to n do if (r > rmax) then rmax r end if end for for i k + 1 to n do for j = k +1 to n do end for end for end for deallocate array (si) end procedure Gauss