We present a two-variable integro-differential equation describing bound systems of unequal mass particles. The method is based on an extension of the two-variable integro-differential equations in the D = 3(A − 1)-dimensional space, known as IDEA, describing the bound states of A-body systems. This method has been successfully applied in the past to few-body systems with the results obtained being in good agreement to those of competing methods. In the present work we investigate whether the same is true for unequal mass particles. Therefore, we first employ the formalism to obtain binding energies for Λ- and ΛΛ-nuclear systems and compare the results with some other results in the field. Secondly, we apply it to the ϕ- and ϕϕ-nuclear systems. Using ϕ–nucleon and ϕ–ϕ interactions, available in the literature, we found that mesic–nuclear systems could exist.