% Example 3.20 svd
A = [ 0.913 0.659
      0.780 0.563
      0.457 0.330 ];
  
[U, S, V] = svd(A);

% apprimate matrix A using only 
the first singular value
B = S(1,1)*U(:,1)*(V(:,1).'),