function [x, Aprime, bprime] = Gauss(A, b)
Aprime= A;
bprime= b;
Aprimebprime=[A,b];
[M,N]=size(Ab);
if (det (A))==0
disp ('Cannot solve the matrix');
end
for i=1:N-1
for j=i+1:N
if Aprime(i)==Aprime(j)==0
disp('The matrix cannot be solved');
Aprimebprime(i,j:M)=Aprimebprime(i,j:M)-Aprimebprime(i,i)/Aprimebprime(j,i)*Aprimebprime(j,j:M);
Aprimebprime(i)=Aprimebprime(i)+Aprimebrime(j);
x=Aprimebprime(i);
end
end
end