c矩阵计算