Important: The linalg package has been deprecated. Use the superseding packages, LinearAlgebra and VectorCalculus, instead.

- For information on migrating linalg code to the new packages, see examples/LinearAlgebraMigration.

The function inverse computes the matrix inverse of A. An error occurs if the matrix is singular.

This function uses Cramer's rule for matrices of dimension less than or equal to 4 by 4, and for matrices with the 'sparse' indexing function.

For other matrices, the inverse is computed by applying the operations for the Gauss-Jordan reduction of A to an identity matrix of the same shape.

The command with(linalg,inverse) allows the use of the abbreviated form of this command.