Write a MATLAB function that accepts matrix A of size n x n and computes the LU decomposition of A with partial pivoting. The function call must be [L, U, P] = function(A). Please remember that you are not allowed to use any built in MATLAB functions to compute the LU decomposition. you are allowed to use functions like max. HINT: Test your code by creating a random matrix, computing the LU decomposition and checking the norm of P*A-L*U.