Не сходится метод наискорейшего спуска для решения СЛАУ

function x = MNS(A,b,n,tolerance)
%MNS - метод наискорейшего спуска
%Вход:
 %   A - матрица системы
 %   b - вектор-столбец правых частей
 %   n - порядок матрицы
 %   tolerance - точность решения
%Выход:
 %   x - вектор решения
 
 if det(A) == 0 %Проверка равенства нулю определителя системы
    disp('Матрица вырождена');
     return
 end
 
  AA = (A')*A; %Приведение системы уравнений к эквивалентной
   bb = (A')*b; 
   
x = zeros(n,1); %Вектор-столбец начального приближения
 nrm = 1;
 
while nrm >= tolerance
      AAA = AA*x;
  for i = 1:n %Вычисление вектора невязки r
      r(i) = bb(i) - AAA(i);  
  end %for
      rr=r';
       lambda = ((rr')*rr)/((rr')*(AA*rr));
      x = x + lambda*rr; %Вычисление приближения вектора x
       nrm = norm(rr); %Вычисление нормы вектора невязки 
end %while

function tutuMNS
n=5;
L=zeros(n); x=zeros(n,1);
for i=1:n
    x(i) = i;
    for j=1:i
        L(i,j) = i+j;
    end
end
A = 0.001*L*L'
b = A*x
tolerance = 1e-3;
x = MNS(A,b,n,tolerance)
xMATLAB = A\b
end % tutuMNS