The function decomposes the matrix A into a lower triangular matrix L

and an upper triangular matrix U, using Crout s method such that A=LU.

Input variables:

A The matrix of coefficients.

Output variable:

L Lower triangular matrix.

U Upper triangular matrix.