metodo de gauss seidel
Publicado por Jhon (2 intervenciones) el 19/05/2013 20:43:28
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
function GS(A,b,x,N)
clc;
n=length(x);
for k=1:N
for i=1:n
s=b(i);
for j=1:i-1
s=s-A(i,j)*x(j);
end
for j=1+i:n
s=s-A(i,j)*x(j);
end
x(i)=s/A(i,i);
end
fprintf('%3.0f %10.4f %10.4f %10.4f\n',k,x);
pause(0.2)
end
esto es el codigo de gauss seidel, se que me falta una cosa q se me paso, pero no logro saber que puede ser el error, por favor ayudenme :(
Valora esta pregunta
0