
    'MaxFuncEvals',10000,'Display','final')
%%
%
% Verifying the solution, we see that we find a matrix decomposition which
% fits the matrix with very small relative error (given the stopping
% tolerance of |1e-5| used by the optimizer).
[U,V] = example2_extract(m,n,k,out.X);
norm(Data.A-U*V')/norm(Data.A)
%%
%
% <html><hr></html>
%% References
%
% [1] Fletcher, R. & Reeves, C.M. (1964). 
% Function minimization by conjugate gradients. 
% _The Computer Journal_, 7, 149-154.
%
% [2] Hestenes, M. R. and Sti