matlab-高斯賽德爾迭代法

% 高斯-賽爾德迭代法
clc;
clear;
% 第一小題
% A=[1 0.4 0.4
% 0.4 1 0.8
% 0.4 0.8 1];
% b=[ 1
% 2
% 3];

% 第二小題
A=[ 1 2 -2
1 1 1
2 2 1];
b=[ 1
1
1];
% A=[ 4 1 -1
% 2 5 2
% 1 1 3];
% b=[ 5
% -4
% 3];
n=1; %迭代次數
x0=[ 1
1
1];%初始向量
m=length(b);
x=zeros(n+1,m);
x(1,1:m)=x0(1:m);
for k=2:n+1
for i=1:m
x(k,i)=(b(i)-sum(A(i,1:i-1).*x(k,1:i-1))-sum(A(i,i+1:m).*x(k-1,i+1:m)))/A(i,i);
end
end
x1=A\b %較准確的值
x %結果