//clc() A = [1,0.5,1/3;1,2/3,1/2;1,3/4,3/5]; B = [1.833333;2.166667;2.35]; U = inv(A); X = U*B; x = det(X(1,1)); y = det(X(2,1)); z = det(X(3,1)); disp(x,"x = ") disp(y,"y = ") disp(z,"z = ")