I_{ij}=\int_S \sigma(\vec{r})(r^2\delta_{ij}-r_i r_j)dS