clear; clc;

A = randi(10, 3, 3)
x = randi(10, 3, 1)
y = A * x

det(A)
x_ = A \ y

B = [1, 2, 3; 4, 5, 6; 7, 8, 9]
det(B)
inv(B)