esidual error:
%          err = norm(A-U*S*V')/norm(A).
%

%%
% The lapack interface routines come from the work of Tim Toolan.
% See <http://www.mathworks.co.kr/matlabcentral/fileexchange/16777-lapack>


%% Examples

%%
% Generate a matrix with sigular values of [1 3 5 7]:
A = randn(7,4)+i*randn(7,4);
[U,S,V] = svd(A,0);
A = U*diag([1 3 5 7])*V';

%%
% Change format:
format short g

%%
% Compute singular values by selecting 'gesvd' implicitly:
svd_lapack(A)

%%
% Compute the econo