I_{ij}=\int_V \rho(\vec{r})(r^2\delta_{ij}-r_i r_j)dV