This paper is devoted to the phenomenon of superintegrability. This phenomenon is manifested in the existence of a formula for character averages, expressed through the same characters at special points and of its various generalization. In this paper we develop a method of proving such formulas from first principle from Virasoro constraints and W -representation. We apply it to prove the formula for the Jack functions averages -appropriate analogue of characters for the β-deformed Hermitian Gaussian matrix model. We also sketch the construction of W -operators from Calogero-Ruijsenaars Hamiltonians.