Matlab

Gauss- Seidel
clc;
%
m=input('Numero maximo de iteraciones:');
n=input('Numero de elementos:');
tol=input('Tolerancia:');
fori=1:1:n
for j=1:1:n+1
a(i,j)=input('a=');
end
end
for i=1:1:n
x(i)=input('x=');
end
for k=1:1:m
co=true;for i=1:1:n
y=0;
for j=1:1:n
if i==j
else
y=y+a(i,j)*x(j);end
end
y=(a(i,n+1)-y)/a(i,i);
if abs((y-x(i))/y)n
disp('No hay solucion');
elsefor j=k:1:n+1
T=a(k,j);
a(k,j)=a(f,j);
a(f,j)=T;
endp=a(k,k);
end
end
for i=1:1:n
if i==k
else
c=a(i,k);
for j=k+1:1:n+1a(i,j)=(a(i,j)*(p-a(k,j))*a(k,j)*c)/PA;
end
return
end
end
PA=p;
end
fori=1:1:n
x(i)=a(i,j)/p;
x(i)
end

Numero de elementos: 5
a= -512
a= 500
a= 0
a= 0
a= 0
a= -9
a= 12
a= -512
a= 500
a=0
a= 0
a= 0
a= 0
a= 12
a= -512
a= 500
a= 0
a= 0
a= 0
a= 0
a= 12
a= -512
a= 500
a= 0
a= 0
a= 0
a= 0
a= 12
a= -512
a=0
ans =

0.0180


ans =

4.3200e-004


ans =

1.0368e-005


ans =

2.4869e-007


ans =

5.8286e-009