Metodo De Jacobi

Método de jacobi en matlab

clear all

clc

En este primer paso se imprime el menú en pantalla

fprintf('Resolución del sistema Ax = b por \n');

fprintf('Jácobi \n\n')

n=input('Ingrese el orden del sistema = ');

fprintf('\n');

fprintf('Matriz coeficientes del sistema \n\n')

Aquí nos pide la aproximación inicial y la tolerancia

for i=1:n

for j=1:n

fprintf('coeficiente A(%d,%d) = ', i,j ) % Definición de la matriz A

A(i,j)=input(' ');

end

end

fprintf('\nVector términos independientes \n\n')

for i=1:n

fprintf('término b(%d) = ', i ) % Definición del vector b

b(i)=input(' ');

end

Se aplica la primera interaction dada

b=b';

fprintf('\nVector aproximación inicial\n\n')

for i=1:n

fprintf('xo(%d) = ', i ) % Definición de la aproximación inicial

xo(i)=input(' ');

end

xo=xo';

e=input('\nIngrese el valor de la tolerancia = '); % Definición de la tolerancia

fprintf('\n');

Se cicla hasta obtener el valor del paro en la matriz

for i=1:n

c(i)=b(i)/A(i,i); % Cálculo del vector c

for j=1:n

if i==j T(i,j)=0; , else T(i,j)=-A(i,j)/A(i,i);, end % Cálculo de la matriz T

end

end

c=c';

Er=A*xo-b;

Aquí es utilizado como criterio de paro en las iteraciones

while norm(Er,inf)>=e

x=T*xo+c;

xo=x;

Er=A*xo-b;

end

disp(x)

...