Let $S=XAX^T$, where $X$ is a random vector with independent coordinates and $A$ is a matrix with coefficients from a Banach space. We want to give a "simple" formula that accurately approximates $||S||_p$ ($p$-th moment of $S$). We will present known results, including the case when $X$ is the standard normal vector. We will also deal with the non-Gaussian case in which almost nothing is known.