Méthode de Jacobi

La méthode de Jacobi est une technique itérative pour résoudre un système linéaire.

function [x]=Jacobi(A,b,x0,k)
//
// On veut resoudre le systeme Ax=b
// de maniere recursive
// x0 condition initiale de l'algorithme de Jacobi
// k nombre d'itération
//
  n=size(A,'c')
  x=x0;
  for p=1:k
    for i=1:n
      x(i)=(b(i)-A(i,1:n)*x+A(i,i)*x(i))/A(i,i)
    end
  end
   
endfunction