I_{ii}=\int dV\rho(\vec{r})(r^2-x_i^2)