x-A'*b)
tic ; x = pinv(A)*b      ; toc, norm (A'*A*x-A'*b)
tic ; x = inverse(A)*b   ; toc, norm (A'*A*x-A'*b)
tic ; x = factorize(A)\b ; toc, norm (A'*A*x-A'*b)
</pre><pre class="codeoutput">Elapsed time is 0.056897 seconds.
ans =
   2.7016e-12
Elapsed time is 0.338749 seconds.
ans =
   2.3034e-12
Elapsed time is 0.079205 seconds.
ans =
   3.0968e-12
Elapsed time is 0.079404 seconds.
ans =
   3.0968e-12
</pre><p>FACTORIZE is better than BACKSLASH because you can reuse the factorization for different right-h