Jacobi

m=size(A,1) ;

n=size(A,2) ;

[ M, N ] = split( A , b, 1.0, 1 ) ;

for iter = 1:max_it,

x_1 = x;

y = N*x + b ;

x = M \ (y) ;

error = norm(x-x_1 )/ norm( x );

if ( error <= tol ), break, end

end

...

résolution de systèmes linéaires par la méthode de Jacobi
- Jacobi 1 : test diagonale dynamique
- Jacobi 2 : test statique

résolution de systèmes linéaires par la méthode de Jacobi