In this paper we calculate kinetic, potential and full energy with three-and four-particle direct correlations taken into account at wide temperature region on the base of the density matrix of the interacting Bose-particles [I. O. Vakarchuk, O. I. Hryhorchak, Journ. Phys. Stud. 3, 3005 (2009)]. In the low temperature limit the obtained expression for the full energy is equal to the wellknown expression for ground state energy in the approximation of "two sums over the wave vector". The results of this work can be applied for the numeric calculation of the heat capacity of liquid 4 He in order to check the theoretical and experimental results quantatively, especially in the λ-transition region. 1 Вступ Розрахунок внутрiшньої енергiї такої багатобозонної системи, як рiдкий 4 He має доволi давню iсторiю. Великою мiрою це пов'язано iз прагненням теоретично описати λ-подiбний хiд кривої теплоємностi в околi точки фазового переходу. Першi кроки були зробленi в напрямi наближеного обчислення енергiї основного стану i енергетичного спектру цiєї системи ще Боголюбовим [1] майже сiмдесят рокiв тому. Пiзнiше вiн разом iз Зубарєвим у роботi [2] знайшли хвильову функцiю основного стану, а також хвильовi функцiї нижнiх збуджених рiвнiв для слабонеiдеального бозе-газу з допомогою методу колективних змiнних. В границi зникаюче малої взаємодiї енергiя основного стану бозе-системи була вперше коректно знайдена в роботi [3]. Розрахунок проводився з використанням псевдотенцiалу.Дослiдження термодинамiчних функцiй основного стану рiдкого 4 He, зокрема внутрiшньої енергiї, в наближеннi вищому, нiж наближення Боголюбова, було проведено в рядi робiт [4,5,6,7]. Завдяки цим роботам було показано, що врахування три-та чотиричастинкових кореляцiй покращує значення для енергiї основного стану.Трохи згодом в роботах [8,9] з допомогою двочасових температурних функцiй Грiна були отриманi вирази для термодинамiчних функцiй рiдкого 4 He в широкотемпературному дiапазонi. В пiдходi колективних змiнних термодинамiчнi функцiї багатобозонної системи в широкотемпературнiй областi були знайденi в роботi [10]. Розрахунок проводився за допомогою уcереднення з матрицею густини взаємодiючих бозе-частинок в наближеннi парних кореляцiй. Узгодження з експериментальними даними отриманих результатiв для кiнетичної, потенцiальної та повної енергiї є досить добрим [11,12,13,14], однак неповним. Це великою мiрою пов'язано з тим, що для матрицi густини було взято лише наближення парних кореляцiї. Для бiльш точних результатiв потрiбно врахувати три-та чотиричастинковi кореляцiї. Як вiдомо [15,16], їх внесок в значення термодинамiчних функцiй може виявитися досить значним. Врахування внескiв три-i чотиричастинкових кореляцiй у внутрiшню енергiю багатобозонної системи i є предметом цiєї роботи. 1 arXiv:1506.03987v1 [cond-mat.quant-gas]