function B=matrixinverse(A)
[a b]=size(A);
K=[A,eye(a,b)];
J=rref(K);
B=J(:,(b+1):2*b);
return;