16440783e-32
     -1.11022302462516e-16
      1.11022302462516e-16


err =

      2.25593997819691e-16

</pre><p>Calculate SVD using ZGESDD:</p><pre class="codeinput">[U,S,V,err] = svd_lapack(A,[],<span class="string">'gesdd'</span>)
</pre><pre class="codeoutput">
U =

  Columns 1 through 3

         0.707106781186547     -2.77555756156289e-17         0.707106781186547
         0.707106781186547      2.77555756156289e-17        -0.707106781186547
                         0                         1      3.9