Algorithm 432 C2: Solution of the matrix equation AX + XB = C F4 | Synapse