how to write a code for a linear system of equations using the gauss-seidel iteration solver in matlab

main.m
% define the matrix A and the vector b
A = [10 -1 2; -1 11 -1; 2 -1 10];
b = [6; 25; -11];

% initial guess for the solution
x0 = [0; 0; 0];

% maximum number of iterations and tolerance
maxIter = 100;
tol = 1e-6;

% initialize x and the iteration counter
x = x0;
iter = 0;

while iter < maxIter
    % Gauss-Seidel iteration
    for i = 1:size(A, 1)
        x(i) = (b(i) - A(i,1:i-1)*x(1:i-1) - A(i,i+1:end)*x0(i+1:end)) / A(i,i);
    end
    
    % check for convergence
    if norm(A*x - b) < tol
        break;
    end
    
    % update iteration counter
    iter = iter + 1;
    
    % update the previous solution
    x0 = x;
end

% print the solution and the number of iterations
disp(['Solution: ', num2str(x')])
disp(['Number of iterations: ', num2str(iter)])
768 chars
37 lines