a = [2 3 1 5; 1 0 3 1; 0 2 -3 2; 0 2 3 1];
b = inv(a);
disp(a * b); % result is (approx.) identity matrix