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.