This example solves the system of equations $AX=B$ using iterative refinement and to compute the forward and backward error bounds, where and

Here $A$ is symmetric, stored in packed form, and must first be factorized by zsptrf.