Résoudre un système linéaire par la méthode itérative de Jacobi : jacobi_linsolve

### 6.63.7 Résoudre un système linéaire par la méthode itérative de Jacobi : jacobi_linsolve

jacobi_linsolve a comme paramètres :
la matrice du système A,
le second membre b sous la forme d’une liste,
la tolérance eps et
le nombre maximum d’itération maxiter.
la matrice du système A,
le second membre b sous la forme d’une liste,
la tolérance eps et
le nombre maximum d’itération maxiter.
jacobi_linsolve permet de résoudre un système d’équations linéaires de matrice A et de second membre b, par la méthode itérative de Jacobi avec une tolérance de eps et un nombre maximal d’itérations de maxiter.
jacobi_linsolve renvoie une liste qui est solution du système d’équations.
On tape :

A:=[[100,2],[2,100]];
jacobi_linsolve(A,[0,1],1e-12);

On obtient :

[-0.000200080032,0.0100040016006]

Si on tape :

evalf(linsolve(A,[0.1]))

On obtient :

[-0.000200080032013,0.0100040016006]