-function [x, k] = conjugate_gradient_method(A, b, x0, tolerance, max_iterations)
+function [x, k] = vanilla_cgm(A, b, x0, tolerance, max_iterations)
%
% Solve,
%
%
% All vectors are assumed to be *column* vectors.
%
- zero_vector = zeros(length(x0), 1);
-
k = 0;
x = x0; % Eschew the 'k' suffix on 'x' for simplicity.
rk = A*x - b; % The first residual must be computed the hard way.