%Gauss-seidel method
f=@(x) [10*x(1)+2*x(2)-x(3)-27;
    -3*x(1)-6*x(2)+2*x(3)+61.5;
    x(1)+x(2)+5*x(3)+21.5];
tol=10^(-8);
%초기값은 (0,0,0)으로 설정, x(1)=a, x(2)=b, x(3)=c로 설정
i=1;
a(1)=27/10;
b(1)=(-61.5+3*a(1))/-6;
c(1)=(-21.5-a(1)-5*b(1));
while max([abs((a(i+1)-a(i))/a(i+1)) abs((b(i+1)-b(i))/b(i+1)) abs((c(i+1)-c(i))/c(i+1))])>tol
    a(i+2)=(27-2*b(i+1)-c(i+1))/10;
    b(i+2)=(-61.5+3*(a(i+2))-2*c(i+1))/-6;
    c(i+2)=-21.5-a(i+2)-5*b(i+2);
    i=i+1;
end
fprintf('수치해는 %f %f %f 이다.\n',a(i+1),b(i+1),c(i+1))

코딩을 이렇게 했는데 while absmax([abs((a(i+1)-a(i))/a(i+1)) abs((b(i+1)-b(i))/b(i+1)) abs((c(i+1)-c(i))/c(i+1))])>tol 이 라인이 인덱스가 행렬의 차원을 초과한다는데 어떻게 고쳐야하는지 알려주시면 감사하겠습니다.