Module SSYR2K in IMSLMPerforms one of the symmetric rank 2k operations C = alpha*A*trans(B) + alpha*B*trans(A) + beta*C, or C = alpha*trans(A)*B + alpha*trans(B)*A + beta*C, where alpha and beta are scalars, C is an n by n symmetric matrix and A and B are n by k matrices in the first case and k by n matrices in the second case. Trans represents the transpose of the matrix.