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.