This example solves the system of equations $AX=B$ using iterative refinement and to compute the forward and backward error bounds, where Here $A$ is symmetric and positive definite, and is treated as a band matrix, which must first be factorized by dpbtrf.