LU Factorization without pivoting from Gaussian Elimination in Matlab

Question

I have successfully implemented Gaussian elimination in the following Matlab code. The code gives Gaussian elimination form, from which lower and upper triangles (L and U) can b