function x=gaussjordan(A,b)
[n,m]=size(A)
for j=1:n
  for i=1:n
    if i~=j
      t=A(i,j)/A(j,j)
      A(i,:)=A(i,:)-t*A(j,:)
      b(i)=b(i)-t*b(j)
    endif
  endfor
endfor
x=b./diag(A)
endfunction