An approach allowing us to solve exactly a Schrodinger eigenvalue problem for a molecular potential of arbitrary shape is developed. The method is based on an integral equation. The conditions, under which the original eigenvalue problem may be transformed into an algebraic one, are discussed. Owing to a special partitioning of space into cells, only some absolutely and uniformly convergent Green's function expansions are involved. As a consequence, all the algebraic manipulations are rigorously justified, in contrast to some direct generalizations of multiple scattering technique. Nevertheless, the method offered results in much the same secular equation, as one previously obtained [R. G. Brown, J. Phys. B 21, L309 (1988)l. The relationship between the present method and other approaches is discussed. The matching conditions between any two local representations of a true wave function, which are believed to be some extra and independent ones [E. Badralexe and A. J. Freeman, Phys. Rev. B 37, 10469 (1988)], are shown to result from a corresponding set of algebraic equations. Some computational aspects are also discussed. A preliminary result of the H: energy spectrum calculation is reported.