function det = determ(A)

% determ  Matrix determinant from plu.
%
% det = determ(A) computes the determinant of the square matrix A
% as the sign of the permutation times the product of pivots.

[P, L, U, sign] = splu(A);
det = sign * prod(diag(U));