function x=gausselim(A,b)n=size(A,1);for i=1:n-1[foo,p]=max(abs(A(i:n,i)));p=p+i-1;A([i,p],:)=A([p,i],:);b([i,p])=b([p,i]);for j=i+1:nm=A(j,i)/A(i,i);for k=i:nA(j,k)=A(j,k)-m*A(i,k);endb(j)=b(j)-m*b(i);endendx=b;for i=n:-1:1for j=i+1:nx(i)=x(i)-A(i,j)*x(j);endx(i)=x(i)/A(i,i);end