Matlab - Metodo iterativo de Jacobi

% Metodo Iterativo de Jacobi)

function X=jacobi(A,B,P,delta, max1)

% Datos

% - A es una amtriz invertible de orden N x N

% - B es una matriz de orden N x 1

% - P es una matriz de orden N x 1: el punto inicial

% - delta es la tolerancia para P

% - max1 es el número máximo de iteraciones

% Resultados

% - X es una matriz de orden N x 1:

% la aproximación a la solución de AX=B

% generada por el método iterativo de Jacobi

N = length(B);

  X = 0;

for k=1: max1

    for j=1:N

        X(j)=(B(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j);

    end

    X = zeros(1, 100);

    err=abs(norm(X'-P));

    relerr=err/(norm(X)+eps);

    P=X';

    if(err<delta)||(relerr<delta)

        break

    end

end

X=X';

??? Error using ==> mtimes

Inner matrix dimensions must agree.

Error in ==> jacobi at 16

        X(j)=(B(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j);

% Metodo Iterativo de Jacobi)

function X=jacobi(A,B,P,delta, max1)

% Datos

% - A es una amtriz invertible de orden N x N

% - B es una matriz de orden N x 1

% - P es una matriz de orden N x 1: el punto inicial

% - delta es la tolerancia para P

% - max1 es el número máximo de iteraciones

% Resultados

% - X es una matriz de orden N x 1:

% la aproximación a la solución de AX=B

% generada por el método iterativo de Jacobi

N = length(B);

  X = 0;

for k=1: max1

    for j=1:N

        X(j)=(B(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j);

    end

    err=abs(norm(X'-P));

    relerr=err/(norm(X)+eps);

    P=X';

    if(err<delta)||(relerr<delta)

        break

    end

end

>> A=rand(4); B=rand(4,1); P=zeros(4,1); delta=0.0001; max1=100;X=jacobi_vectorizada(A,B,P,delta, max1)

X =

  1.0e+058 *

   -0.3254

   -0.2953

   -0.1737

   -1.2749