plot command is not work in Gauss Seidel Method

A = [

1 0 -1 0 0 0 0

-1 1 1 0 0 0 0

0 -1 2 -1 0 0 0

0 -1 0 2 0 1 0

0 -1 -1 0 2 1 1

0 0 1 -1 0 2 -1

0 0 0 -1 1 -1 2

];

b = [0 0 0 0 0 20 0]';

x=[1 1 1 1 1 1 1];

n=size(x,1);

normVal=Inf;

tol=1e-5; itr=0;

%%Gauss Seidel Method

while normVal>tol

x_old=x;

for i=1:n

sigma=0;

for j=1:i-1

sigma=sigma+A(i,j)*x(j);

end

for j=i+1:n

sigma=sigma+A(i,j)*x_old(j);

end

x(i)=(1/A(i,i))*(b(i)-sigma);

end

itr=itr+1;

normVal=norm(x_old-x);

end

%%

fprintf('system solution : \n%f\n%f\n%f\n%f in %d iterations',x,itr);

hold on

plot(itr,x, '+')

xlim([0 10])

ylim([0 2])

hold off

end

Answer 1

[John D'Errico]

When you tell us that something does not work, you need ot explain why you THINK it does not work. In fact, plot worked perfectly well there. At least, it tiried to work as well as it could.

itr
itr =
     2
     
x
x =
     0     1     1     1     1     1     1
     
plot(itr,x, '+')

So itr is the number 2. x is a vector. What did you expect it to do?

[suketu vaidya]

hello ,

i is not shosws full vactor x . i means x1 = 0 x2 = 1 x3=1 and so on . i want a graph that shows x1,x2,x3 in x axist and on y axis value of x1,x2,x3 and so on

[suketu vaidya]

easy way i try to plot ans of gauss seidel method

Answer 2

[suketu vaidya]

right ans is

hold on

plot(x,'o')