I_{ij}=\int(|\vec{r}|^2\delta_{ij}-x_i x_j)dm