A . x = b



I have a coefficient matrix A = [a_11 ... a_1n] ... [a_m1 ... a_mn], a
column vector of variables x = [x_1 ... x_n], and a right-hand vector
b = [b_1 ... b_m].

Is there an elegant way to solve for A . x = b without doing something
like:

        matrixmap("=", x, invert(A) . b)

E.g., something akin to: solve(A . x = b, x).