B=[0 2/5 1/5;-1/2 0 1/3;-2/10 1/10 0] 

c=[1/5;-4/3;11/10] 

x=[2 2 2]' 

format compact 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=B*x+c 

x=[2 2 2]' 

x(1)=0*x(1)+2/5*x(2)+1/5*x(3) 

x(1)=0*x(1)+2/5*x(2)+1/5*x(3)+c(1) 

x(2)=-1/2*x(1)+0*x(2)+0/3*x(3)+c(2) 

x(1)=0*x(1)+2/5*x(2)+1/5*x(3)+c(1) 

x(2)=-1/2*x(1)+0*x(2)+1/3*x(3)+c(2) 

x(3)=-1/5*x(1)+1/10*x(2)+0*x(3)+c(3) 

x(1)=0*x(1)+2/5*x(2)+1/5*x(3)+c(1) 

x(2)=-1/2*x(1)+0*x(2)+1/3*x(3)+c(2) 

x(3)=-1/5*x(1)+1/10*x(2)+0*x(3)+c(3) 

x(1)=0*x(1)+2/5*x(2)+1/5*x(3)+c(1) 

x(2)=-1/2*x(1)+0*x(2)+1/3*x(3)+c(2) 

x(3)=-1/5*x(1)+1/10*x(2)+0*x(3)+c(3) 

A=[-5 2 1;3 6 -2;2 -1 10],b=[-1;-8;11],x0=[2 2 2]' 

seidel(A,b,x0,0) 

seidel(A,b,x0,0) 

seidel(A,b,x0,3)