mathematica中矩阵的总和

Question

我想对之前计算的矩阵求和。 对于r = 1，n = 3;下标[P，i]是3x3矩阵，如P1，P2，P3。 我的代码是这样的：

'Y = Inverse[S];
 Print["Y=", MatrixForm[Y]];
 For[i = 1, i <= n, i++,  
 Subscript[P, i] = MatrixForm[Outer[Times, S[[All, i]], Y[[i]]]];
 Print["CarpimS=", MatrixForm[S[[All, i]]]];
 Print["CarpimY=", MatrixForm[Y[[i]]]];
 Print["P=", MatrixForm[Outer[Times, S[[All, i]], Y[[i]]]]];  
 ];

 toplamP = MatrixForm[ConstantArray[0, {n, n}]];
 For[i = 2. r + 1, i <= n, i++,
 toplamP = toplamP + Subscript[P, i];
 ];
 Print["ToplamP=", toplamP];'

但是mathematica只给了我P3而且P3没有矩阵形式。