|
|
|
|
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
Tu stock tout dans M, ce qui est faut, il faut utiliser au moins une variable temporaire, si ce n'est plus.
|
Le mieux pour ne pas te tromper, c'est de faire tout pour éviter de modifier M, et de mettre le résultat dans une autre matrice, en plus cela permet de tester que le résultat est bon.
|