Use C++ please

Develop your own program to determine the LU factorization of a square matrix. That is, develop a function that takes a square matrix (of arbitrary dimensions) and returns the triangular matrices [L] and [U]. Confirm that your function is working properly by verifying that [L] Test your function by using it to solve the following system of linear equations: 3x -2x2 +310 8x1-2x2 + 5x,--26