7.9 Solving Systems of Equations
Returns the matrix X
for which (matrix* M X)
must have the same number of rows.
It is typical for B (and thus X) to be a column matrix, but not required.
If B is not a column matrix, matrix-solve solves for all the columns in B
matrix-solve does not solve overconstrained or underconstrained systems, meaning that
M must be invertible.
If M is not invertible, the result of applying the failure thunk fail is
matrix-solve is implemented using matrix-gauss-elim to preserve exactness in its
output, with partial pivoting for greater numerical stability when M is not exact.
See vandermonde-matrix for an example that uses matrix-solve to compute Legendre