Homopolymer adsorption from a dilute solution on an interacting (attractive) surface under static equilibrium conditions is studied in the framework of a Hamiltonian model. The model makes use of the density of chain ends n(1,e) and utilizes the concept of the propagator G describing conformational probabilities to locally define the polymer segment density or volume fraction phi; both n(1,e) and phi enter into the expression for the system free energy. The propagator G obeys the Edwards diffusion equation for walks in a self-consistent potential field. The equilibrium distribution of chain ends and, consequently, of chain conformational probabilities is found by minimizing the system free energy. This results in a set of model equations that constitute the exact continuum-space analog of the Scheutjens-Fleer (SF) lattice statistical theory for the adsorption of interacting chains. Since for distances too close to the surface the continuum formulation breaks down, the continuum model is here employed to describe the probability of chain configurations only for distances z greater than 2l, where l denotes the segment length, from the surface; instead, for distances z < or = 2l, the SF lattice model is utilized. Through this novel formulation, the lattice solution at z = 2l provides the boundary condition for the continuum model. The resulting hybrid (lattice for distances z < or = 2l, continuum for distances z > 2l) model is solved numerically through an efficient implementation of the pseudospectral collocation method. Representative results obtained with the new model and a direct application of the SF lattice model are extensively compared with each other and, in all cases studied, are found to be practically identical.