Inverse: Inverse of a Matrix

In matlib: Matrix Functions for Teaching and Learning Linear Algebra and Multivariate Statistics

Inverse of a Matrix

Description

Uses gaussianElimination to find the inverse of a square, non-singular matrix, X.

Usage

Inverse(X, tol = sqrt(.Machine$double.eps), verbose = FALSE, ...)

Arguments

X	a square numeric matrix
tol	tolerance for checking for 0 pivot
verbose	logical; if TRUE, print intermediate steps
...	other arguments passed on

Details

The method is purely didactic: The identity matrix, I, is appended to X, giving X | I. Applying Gaussian elimination gives I | X^{-1}, and the portion corresponding to X^{-1} is returned.

Value

the inverse of X

Author(s)

John Fox

Examples

  A <- matrix(c(2, 1, -1,
               -3, -1, 2,
               -2,  1, 2), 3, 3, byrow=TRUE)
  Inverse(A)
  Inverse(A, verbose=TRUE, fractions=TRUE)

matlib documentation built on Oct. 3, 2024, 1:09 a.m.