n = length(x0);
M = eye(n);
% The standard CGM is equivalent to the preconditioned CGM is you
% use the identity matrix as your preconditioner.
n = length(x0);
M = eye(n);
% The standard CGM is equivalent to the preconditioned CGM is you
% use the identity matrix as your preconditioner.