Develop your own M-file to determine the LU factorization of a square matrix without partial pivoting. That is, develop a function that is passed the square matrix and returns the triangular matrices [L] and [U]. Test your function by using it to solve the system in Prob. 10.3. Confirm that your function is working properly by verifying that [L][U] = [A] and by using the built-in function lu.