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.