Given a linear system i want to do Gauss-Seidel iteration.But getting error.Would appreciate if anyone fix me the code. Thanks.

Question

Abid Shahriar on 30 May 2021

0

Edited: KALYAN ACHARJYA on 30 May 2021

>> x(1)=1.5;

>> y(1)=-1;

>> z(1)=1;

>> iter(1)=0;

>> for n=1:4

iter(n+1)=n;

x(n+1)=(16+3*y(n)-2*z(n))/8;

y(n+1)=(-9+2*x(n+1)-z(n))/8;

z(n+1)=(1+2*x(n+1)-4*y(n+1))/10;

end

>> Solution = [iter',x',y',z’]

0 Comments
ShowHide -1 older comments

Sign in to comment.

Sign in to answer this question.

Answer 1

KALYAN ACHARJYA on 30 May 2021

0
Link

Direct link to this answer

https://www.mathworks.com/matlabcentral/answers/843705-given-a-linear-system-i-want-to-do-gauss-seidel-iteration-but-getting-error-would-appreciate-if-anyo#answer_712810

No coding syntax error

x(1)=1.5;
y(1)=-1;
z(1)=1;
iter(1)=0;
for n=1:4
 iter(n+1)=n;   
 x(n+1)=(16+3*y(n)-2*z(n))/8;
 y(n+1)=(-9+2*x(n+1)-z(n))/8;
 z(n+1)=(1+2*x(n+1)-4*y(n+1))/10;
end
Solution = [iter',x',y',z']

Result:

Solution =
            0          1.5           -1            1
            1        1.375     -0.90625       0.7375
            2       1.4758     -0.84824      0.73445
            3       1.4983     -0.84223      0.73655
            4          1.5     -0.84206      0.73683