clear, clc, close all

syms x y z

A = [2 5 10; 0 5 10; 0 0 2]
vecx = [x,y,z]'
A*vecx
b = [10 20 2]';
 
linsolve(A,b)
A^-1*b