Octave - Ayuda con Metodo Jacobi Octave

METODO DE JACOBI

Ingrese la matriz de coeficiente a:12

Ingrese la matriz de coeficiente b:2

Ingrese el vector con las aproximacimaciones Iniciales:15

Ingrese el número máximo de iteraciones:3

Ingrese la tolerancia:2

condicional=

 1

La matriz de transicion de jacobi:

0

El vector constante es:   1.666666666666667e-01

error: x(2): out of bound 1

error: called from

    metodo_jacobi at line 47 column 2

>>

clear all

clc

format long

fprintf('METODO DE JACOBI\n\n');

%Ingresamos las variables a calcular:

a=input('Ingrese la matriz de coeficiente a:');

b=input('Ingrese la matriz de coeficiente b:');

x=input('Ingrese el vector con las aproximacimaciones Iniciales:');

iter=input('Ingrese el número máximo de iteraciones:');

total=input('Ingrese la tolerancia:');

%Se calcula el condicional de la matriz de coeficientes

cond=norm(a)*norm(a^-1);

disp('condicional=');

disp(cond);

%Se calcula el determinante de la matriz de coeficientes

determinante=det(a);

if determinante==0;

disp('El determinante es cero, el problema no tiene solución única')

return

end

n=length(b);%numero de elementos del vector b

d=diag(diag(a)); %obtencion de la matriz diagonal

l=d-tril(a); %obtencion de la matriz diagonal superior L

u=d-triu(a);%obtencion de la matriz diagonal inferior u

fprintf('La matriz de transicion de jacobi:\n');

T=d^-1*(l+u); % matriz de transicion de jacobi

disp(T);

re=max(abs(eig(T))); %calculo del radio espectral

if re>1

disp('Radio Espectral mayor que 1');

disp('el método no converge');

return

end

fprintf('El vector constante es:');

C=d^-1*b; %vector constante C, para el metodo

disp(C);

i=0;

err=total+1;

z=[i,x(1),x(2),x(3),err]; %vector que me permite graficar la tabla

while err>total & i<iter

xi=T*x+C;

%disp(xi)

err=norm(xi-x); %norma 2

%err=max(abs(xi-x)); %norma 1

%err=norm(xi-x)/norm(xi); %norma relativa

x=xi;

i=i+1;

z(i,1)=i;

z(i,2)=x(1);

z(i,3)=x(2);

z(i,4)=x(3);

z(i,5)=err;

end

fprintf('TABLA:\n\n n x1 x2 x3 Error\n\n ');

disp(z)%impresión de la tabla.