function [X] = systemsolver(A,B) K=[A,B]; X=rref(K); return;