In this study, the bound state energy of a four-quark system was analytically calculated as a two heavy–heavy anti-quarks ˉbˉbˉbˉb and two light–light quarks udud. Tetraquark was assumed to be a bound state of two-body system consisting of two mesons, each containing a light quark and a heavy antiquark. Due to the presence of heavy mesons in the tetraquark, Born–Oppenheimer approximation was used to study its bound states. To assess the bounding energy, Schrödinger equation was solved using lattice QCD ˉbˉbˉbˉb potential, having expanded the tetraquark potential ˉbˉbudˉbˉbud up to 11th term. Binding energy state and wave function, however, were obtained in the scalar u/du/d channel. Graphical results for wave functions obtained versus antiquark–antiquark distance ˉbˉbˉbˉb confirmed the existence of the tetraquark ˉbˉbudˉbˉbud. Analytical bound state energy obtained here was in good agreement with several numerical ones published in the literature, confirming the accuracy of the approach taken here.