Write a matlab function for gaussian elimination based off of the following psuedo code

function Gauss(n, (a_ij), (l_i))

integer i, j, k, n

real r, rmax, smax, xmult

real array (a_ij)_{1:n x 1:n}, (l_i)_1:n; real array allocate (s_i)_1:n

for i = 1 to n

l_ii

smax 0

for j =1 to n

smax max(smax, abs.value(a_ij)

end for

s_i smax

end for

for k = 1 to n-1

rmax 0

for i = k to n

r abs. value (a_ij/s_{l_i}) %l_i = l_i

if (r > rmax) then

rmax r

j i

end if

end for

l_j l_k

for i = k +1 to n

xmult a_{l_i,k}/a_{l_k,k}

a_{l_i,k} xmult

for j = k+1 to n

a_{l_i,j} a_{l_i,j} -(xmult)*a_{l_k,j}

end for

end for

end for

deallocate array (s_i)

end function Gauss