function k=gaussjordan(A,B,m,n);
for i=1:n
      s=A(i,i);
       for j=1:n
            A(i,j)=A(i,j)/s
        end
       B(i)=B(i)/s;
        for k=1:m
            if k~=i
                 s=A(k,i);
                 for j=1:n
                     A(k,j)=A(k,j)-A(i,j)*s;
                 end
                 B(k)=B(k)-B(i)*s;
             end
         end
     end