na-lab-gauss seidel

code for gauss seidel

a = [4.63 -1.21 3.22 ;-3.07 5.48 2.11 ; 1.26 3.11 4.57;];

b = [2.22 ; -3.17 ; 5.11;];

x = [0;0;0];

n= input("Enter no of iterations :");

for k=1:n

xold = x;

for i=1:size(a,1)

sum = 0;

for j=1:i-1

sum = sum+a(i,j)*x(j);

end

for j=i+1:size(a,1)

sum = sum + a(i,j)*xold(j);

end

x(i)= (1/a(i,i))*(b(i)-sum);

end

end

disp(x)