Pivot de gauss Matlab

Bonjour,
je voudrais faire l'algorithme de pivot de qauss en matlab mais ça n'a pas marché, voici le code que j'ai fait:

function[]=gauss()
n=input('donner lla dimension de la matrice ')
for i=1:n
for j=1:n
sprintf('M[%i,%i]',i,j)
M(i,j)=input('')
end
end

for k=1:(n-1)
pivot=M(k,k)
if pivot ~= 0
for i=(k+1):n
for j=k+1:n
M(i,j)=M(i,j)-(M(i,k)/pivot)*M(k,j)
end
end
else printf('probleme')
end
end

quelqu'un peut me le corriger?? SVP.merci.
Configuration: Windows Vista
Opera 9.60

2 réponses

  1. Contributeur
    Tu stock tout dans M, ce qui est faut, il faut utiliser au moins une variable temporaire, si ce n'est plus.
    Regarde sur mon site, dans ma classe matrice, il y a l'inversion par pivot de Gauss de mit.
    Sinon, MATLAB inverse très bien les matrice, pas besoin de se prendre la tête ;)
    1
    1. SVP pouvez vous m'expliquer davantage? j'ai pas compris,comment je dois stocker tous dans M? et quelle variable temporaire dois-je utiliser? quel est votre site?
      0