    -0.43568 +    0.17504i
      0.35993 +     0.1078i     -0.42549 -    0.71005i


err =

   1.0157e-15

</pre><h2>Diagnostics<a name="11"></a></h2><p>Lapack, a fortran computational library, has two different subroutines for the Singular Value Decompostion (SVD): xGESVD and xGESDD. xGESVD is based on an implicit QR iteration and xGESDD uses a divide-and-conquer approach. See <a href="http://www.netlib.org/lapack/lug/node32.html">http://www.netlib.org/lapack/lug/node32.html</a> and <a href="http://www.netl