Résoudre une equation linéaire avec gauss

Bonjour,
je commence à resoudre une equation linéaire à l'aide de matlab alrs j'ai déja vu sur le net et j'ai trouver le pgrm de princess sarah qui ma donner les résultats voulu mais le seul problème c que ya des étapes que j'ai pa bien compris
alrs svp sk vous pouvez me faire des comentaires sur se pgm:
merci d'avance
svp faite le moi le plus rapidement possible



function x = solution_gauss(A,b)

n=size(A,1);

for p=1:n
vec=[(1:p-1) n (p:n-1)];
test=1;
while A(p,p)==0
if test==n
error('La matrice n''est pas inversible')
end
A=A(vec,:);
b=b(vec);
test=test+1;
end
b(p)=b(p)/A(p,p);
A(p,:)=A(p,:)/A(p,p);
for q=p+1:n
b(q)=b(q)-A(q,p)*b(p);
A(q,:)=A(q,:)-A(q,p)*A(p,:);
end
end
x=zeros(n,1);
x(n)=b(n);
for p=n-1:-1:1
x(p)=b(p);
for q=p+1:n
x(p)=x(p)-A(p,q)*x(q);
end
end