Помимо того, что написал Acl, вы можете использовать LinearSolve
:
Я заимствую данные Acl
Remove["Global`*"];
eqns = {x - y == 1, 2 x + 2 y == 3, 5*x - 3*y - z == 2}
vars = {x, y, z}
A = CoefficientArrays[eqns, vars];
sol = LinearSolve[A[[2]], -A[[1]]]
, что дает
{5/4, 1/4, 7/2}
In[135]:= Thread[vars->sol]
Out[135]= {x->5/4,y->1/4,z->7/2}