jacob求解实对称矩阵的特征值特征向量c代码