Matlab - metodo de gauss seidel

   
Vista:

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
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder