I_{ik}=\int \rho\left(x_l^2\delta_{ik}-x_ix_k\right)\dd{V}