Can you also explain what each line of the code does?

Write a user-defined MATLAB function that decomposes an n times n matrix [A] into a lower triangular matrix [L] and an upper triangular matrix [U] (such that [A]=[L] [U] using the Gauss elimination method (without pivoting). For the function name and arguments, use [L, U] = LUdecompGauss (A), where the input argument A is the matrix to be decomposed and the output arguments L and U are the corresponding upper and lower triangular matrices. Use LUdecompGauss to determine the LU decomposition of the following matrix: [4 -8 2 -8 -1 0 -3.5 -4 3 -3 10 1 2 -3.5 3.75 -0.5]