=svd_lapack(A,0,'gesdd'); toc

%%
% Measure times when the complete SVD is calculated:
tic, [U,S,V]=svd(A,0); toc
tic, [U,S,V]=svd_lapack(A,0,'gesvd'); toc
tic, [U,S,V]=svd_lapack(A,0,'gesdd'); toc

%%
% In Octave, we can use svd_driver() to select subroutines:
%
%   svd_driver('gesvd'), tic, [U,S,V]=svd(A,0); toc
%   svd_driver('gesdd'), tic, [U,S,V]=svd(A,0); toc


%% Convergence test
%
% SVD has implemented in many numerical libraries such as LINPACK, EISPACK,
% and Lapack. See <http://