This function solves the equation a %*% x = b
for x
,
where a
is tridiagonal and b
can be either a vector or a
matrix.
# S3 method for tridiag
solve(a, b, …)
a tridiag
object: a square tridiagonal (all zeroes
except for the main diagonal and the diagonals immediately above and
below it) matrix containing the coefficients of the linear system.
a vector or matrix giving the right-hand side(s) of the
linear system. If missing, b
is taken to be an identity matrix
and the function will return the inverse of a
.
ignored
Uses the LINPACK dgtsv
routine.