用C++求两个对称矩阵之和与乘积的算法

请高手指点

按照矩阵运算公式进行，只不过计算循环要考虑对称性，结果依然是对称阵，如

//sum

for(int i = 0; i < n; ++i)

{

for(int j = 0; j <= i; ++i)

{

c[i][j] = a[i][j] + b[i][j];

c[j][i] = c[i][j];

}

}

//

for(int i = 0; i < n; ++i)

{

for(int j = 0; j <= i; ++i)

{

for(int l = 0; l < n; ++l)

c[i][j] += a[i][k] * b[j][k];

c[j][i] = c[i][j];

}

}