Problem $6$

In your

HW$2.$py file create a function called naive$\_$LU$()$ that takes the following arguments as input in this order:

numpy array $A$ : a matrix

And returns the following output:

numpy array $L$ : the lower triangular matrix of the $LU$ decomposition for A $($ignore pivoting$)$

numpy array $U$ : the upper triangular matrix of the LU decomposition for A

Hints: You will not receive credit if you just call any built$-$in $lu()$ function. Also have a check inside the function which verifies that $A$ is square before performing the decomposition.