functionx = gauss(A, b)n = length(b);for k = 1 : n-1
if A(k,k)==0
fprintf('Error: the %dth pivot element equal to zero!\n',k);
return;
end
index = [k+1:n];
m = -A(index,k)/A(k,k);
A(index,index) = A(index,index) + m*A(k,index);
b(index) = b(index) + m*... |
|