I_{ij}=I_{ji}=-\int dV\rho(\vec{r})x_ix_j