T=\frac{1}{2}\{\omega_x^2\int dV\rho(\vec{r})(y^2+z^2)+\omega_y^2\int dV\rho(\vec{r})(x^2+z^2)+\omega_z^2\int dV\rho(\vec{r})(x^2+y^2)+\omega_x\omega_y(-\int dV\rho(\vec{r})xy)+\omega_x\omega_z(-\int dV\rho(\vec{r})xz)+\omega_y\omega_y(-\int dV\rho(\vec{r})yz)+\omega_y\omega_x(-\int dV\rho(\vec{r})yx)+\omega_z\omega_x(-\int dV\rho(\vec{r})zx)+\omega_z\omega_y(-\int dV\rho(\vec{r})zy)\}