Let $(X_i,i=1,2,...)$ be a Gaussian sequence with $X_i\in N(0,1)$ for each i and suppose its correlation matrix $R={\left({\rho }_{ij}\right)}_{i,j\ge 1}$ is the matrix of some linear operator R:l₂→ l₂. Then for $f_i\in L²\left(\mu \right)$, i=1,2,..., where μ is the standard normal distribution, we estimate the variation of the sum of the Gaussian functionals $f_i\left(X_i\right)$, i=1,2,... .