Write a Scilab function in the form function [L, U] = lufactor(A) that computes the LU factorization of a matrix A elementof R^n times n using Gaussian elimination (without partial pivoting). Input: -A: n times n matrix middot Output: (matrices L and U such that A = LU) -L: n times n lower triangular matrix -U: n times n upper triangular matrix Check if the matrix A is indeed upper triangular or lower triangular, respectively. If not, do nothing and set x = %nan