
%      [U, S, V] = svd_lapack (A, ECON)
%      [U, S, V] = svd_lapack (A, ECON, OPT)
%      S = svd_lapack (A, ...)
%      [U, S, V, err] = svd_lapack (A, ...)

%% Description
% Compute the singular value decomposition of a matrix 
%   A = U*S*V'
% where A is complex or real (not necessarily square), 
% U and V are unitary or orthogonal, and
% S is always non-negative real diagonal.
% 
% The singular values reveal the rank of the matrix. Moreover,
% the singular value decomposition may be used