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).