Compute the LU decomposition of A.
A will be a subset of L * U with lower triangular matrix L and upper triangular matrix U.
The result is returned in a permuted form, according to the optional return value P.
Accuracy: The result is a valid enclosure.
See also: @infsup/qr, @infsup/chol.
Package: interval