suivant: La décomposition LU d'une
monter: Résolution de MX =
précédent: L'algorithme
Table des matières
Index
Voici le programme de résolution de AX = b :
//M*res=b res renvoie le vecteur solution
//M:=[[1,1],[1,-1]]; b:=[4,2]
//M:=[[1,1,1],[1,1,1],[1,1,1]];b:=[0,0,0]
//M:=[[1,2,1],[1,2,5],[1,2,1]];b:=[0,1,0]
linsolvm(M,b):={
local A,B,N,n,k,d,res;
d:=ncols(M);
//A est la matrice du systeme +le 2nd membre
A:=border(M,-b);
//N contient une base du noyau
N:=-ker(A);
n:=size(N);
//res a d+1 composante (la derniere=1)
res:=makelist(0,0,d);
//C_(k) designe les constantes arbitraires
for (k:=0;k<n-1;k++){
res:=res+N[k]*C_(k);
}
res:=res+N[n-1];
res:=suppress(res,d);
return res;
}
Documentation de giac écrite par Renée De Graeve