I_{ik}=\sum m\left(x_l^2\delta_{ik}-x_ix_k\right)