routine:
%
% * DBDSDC computes the singular value decomposition (SVD) of a real
% N-by-N (upper or lower) bidiagonal matrix B = U * S * VT,
% using a divide and conquer method, where S is a diagonal matrix
% with non-negative diagonal elements (the singular values of B), and
% U and VT are orthogonal matrices of left and right singular vectors,
% respectively.

%%
% Matlab's built-in function svd seems to use the lapack subroutine xGESVD.
% Meanwhile, Octave's built-in function svd support both a