share/matrix/invert.usg

   1 BMT@MIT-MC 09/29/77 21:32:25
   2 The file SHARE;INVERT > finds the inverse of a matrix using the adjoint
   3 method.  This allows a user to compute the inverse of a matrix with bfloat
   4 entries or polynomials with floating pt. coefficients without converting to
   5 cre-form.  The DETERMINANT command is used to compute cofactors, so if
   6 RATMX is FALSE (the default) the inverse is computed without changing the
   7 representation of the elements.  The functions ADJOINT and INVERT are
   8 provided.  The current (temporary?) implementation is inefficient for
   9 matrices of high order.
  10
  11 The DETOUT flag if true keeps the determinant factored out of the inverse.
  12
  13 Note: the results are not automatically expanded.  If the matrix originally
  14 had polynomial entries, better appearing output can be generated by
  15 EXPAND(INVERT(mat)),DETOUT.  If it is desirable to then divide through by
  16 the determinant this can be accomplished by XTHRU(%) or alternatively
  17 from scratch by EXPAND(ADJOINT(mat))/EXPAND(DETERMINANT(mat)).
  18 INVERT(mat):=ADJOINT(mat)/DETERMINANT(mat).