Improves the computed solution to a symmetric positive definite tridiagonal system of linear equations AX=B, and provides forward and backward error bounds for the solution.