[Matlab] Gauss-Seidel et Cholesky

Résolu/Fermé
felecity - 15 mai 2009 à 23:37
 sinox - 17 mai 2011 à 18:29
bsr à tous je cherche un programme en matlab pour résoudre par la méthode de Gauss seidel et la méthode directe de cholesky un système d'équations AX=b merci

2 réponses

function B=cholesky(A)


[n,m]=size(A);


if m ~=n
Error('A n est pas une matrice carre')
end
B=zeros(size(A));
for j=1:n
S(j,j)=A(j,j)
for k=1:j-1
S(j,j)=S(j,j)-B(j,k)^2;
end


B(j,j)=sqrt(S(j,j));

for i=j+1:n
B(i,j)=A(i,j)
for k=1:j-1
B(i,j)=A(i,j)-B(i,k)*B(j,k);
end
B(i,j)=B(i,j)/B(j,j)
end
end
29
merci bcp
0
thank ya !
0