For this routine two examples are presented, both of which involve computing the singular value decomposition of a matrix
, where
in the first example and
in the second.
must first be reduced to tridiagonal form by
dgebrd. The program then calls
dorgbr twice to form
and
, and passes these matrices to
dbdsqr, which computes the singular value decomposition of
.