Решение СЛАУ методом Гаусса-Зейделя

clear;
clc;
a=[-5 3 3;1 -4 3;-5 -4 2]   
b=[24;-1;15]
e=0.00001;
n=3;
miter=100;
for i=1:n
    x(i)=0;
end
for iter=1:miter
    big=0;
        for i=1:n
            sum=0;
                  if i==1
                      for j=i+1:n
                          sum=sum+a(i,j)*x(j);
                      end
                  elseif i==n
                          for j=1:n-1
                              sum=sum+a(i,j)*x(j);
                          end
                  else 
                          for j=1:i-1
                              sum=sum+a(i,j)*x(j);
                          end
                          for j=i+1:n
                              sum=sum+a(i,j)*x(j);
                          end
                      
                  end
               temp=(b(i)-sum)/a(i,i);
               if abs(temp-x(i))>big
                   big=abs(temp-x(i));
               end
               x(i)=temp;
        end
            if big < e 
                niter=iter;
                break;
            end
end
disp('Решение получено с помощью метода Гаусса-Зейделя');
x=x'
disp('Решение получено с помощью средств MatLab');
x1=(a^-1)*b

x = zeros(1, n);
cmp = false;
k = 0;
while cmp == false
    for i = 1 : n
        sum = 0;
        for j = 1 : n
            if (j ~= i)
                sum = sum + a(i,j) * x(j);
            end
        end    
        temp(i) = x(i);
        x(i) = (b(i) - sum) / a(i,i);
    end
    for i = 1 : n
        if abs(x(i) - temp(i)) >= e
            break;
        end
    end
    if i == n
        cmp = true
    end
    k = k + 1;
end