Function: invert (<M>) Returns the inverse of the matrix <M>. The inverse is computed by the adjoint method.
This allows a user to compute the inverse of a matrix with bfloat entries or polynomials with floating pt. coefficients without converting to cre-form.
Cofactors are computed by the
determinant function, so if
false the inverse is computed without changing the representation of the elements.
The current implementation is inefficient for matrices of high order.
true, the determinant is factored out of the inverse.
The elements of the inverse are not automatically expanded. If <M> has polynomial elements, better appearing output can be generated by
expand (invert (m)), detout. If it is desirable to then divide through by the determinant this can be accomplished by
xthru (%)or alternatively from scratch by
expand (adjoint (m)) / expand (determinant (m)) invert (m) := adjoint (m) / determinant (m)
^^ (noncommutative exponent) for another method of inverting a matrix.
There are also some inexact matches for
?? invert to see them.
(%o1) true (%i2)