function [x,k,erro]=seidel(A,b,x0,tol)
%n=lengh(A);
%%%%% Matriz Iterativa %%%%%
n=length(A);
for i=1:n
for j=1:n
if i==j
F(i,j)=0;
else
F(i,j)=-A(i,j)/A(i,i);
d(i,1)=b(i,1)/A(i,i);
end
end
end
k=0; erro=tol+1;
x=x0;
while erro>tol
for i=1:n
x(i,1)=F(i,:)*x+d(i,1);
end
erro=max(abs(x-x0))
x0=x;
k=k+1;
end
end
Resolução
clear, clc
%================================================
%Dados de entrada
A=[10 2 1; 1 -15 1; 2 3 10];
b=[7 32 6]';
x0=[0 0 0]';
tol=0.001;
x,k,erro]=seidel(A,b,x0,tol)