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

1 view (last 30 days)
>> 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’]

Answers (1)

KALYAN ACHARJYA
KALYAN ACHARJYA on 30 May 2021
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
  3 Comments

Sign in to comment.

Categories

Find more on Systems Of Linear Equations in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!