suivant: Le programme de factorisation
monter: La décomposition de Cholesky
précédent: La méthode par identification
Table des matières
Index
//utilise decomplu et splitmat ci-dessus
//A:=[[1,0,-1],[0,2,4],[-1,4,11]]
//A:=[[1,1,1],[1,2,4],[1,4,11]]
//A:=[[1,0,-2],[0,2,6],[0,2,11]]
//A:=[[1,-2,4],[-2,13,-11],[4,-11,21]]
//A:=[[-1,-2,4],[-2,13,-11],[4,-11,21]] (pas def pos)
//A:=[[24,66,13],[66,230,-11],[13,-11,210]]
choles(A):={
local j,n,p,L,U,D,p0;
n:=size(A);
A:=1/2*(A+tran(A));
(p,L,U):=decomplu(A);
p0:=makelist(x->x,0,n-1);
if (p!=p0) {return "pas definie positive ";}
D:=makemat(0,n,n);
for (j:=0;j<n;j++) {
//if (U[j,j]<0) {return "pas def positive";}
D[j,j]:=sqrt(U[j,j]);
}
return normal(L*D);
}
Documentation de giac écrite par Renée De Graeve