Baixe o app para aproveitar ainda mais
Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original
function [xnew, it, dx, xhist] = GaussSeidel(A, b, x1, eps, maxit) xhist = x1'; n = length(b); x = x1; xnew = x; for it = 1:maxit for i = 1:n suml = b(i) - A(i,1:(i-1))*xnew(1:(i-1)) - A(i,(i+1):n)*x((i+1):n); xnew(i) = suml/A(i,i); end dx = abs(xnew-x); err = norm(dx); relerr = err/(norm(x)+eps) x = xnew; xhist = [xhist;x']; if (relerr < eps) break; end end end
Compartilhar