Abstract: For a distributive lattice L, we consider the problem of interpolating functions f : D → L defined on a finite set D ⊆ Ln , by means of lattice polynomial functions of L. Two instances of this problem have already been solved. In the case when L is a distributive lattice with least and greatest elements 0 and 1, Goodstein proved that a function f : {0, 1}n → L can be interpolated by a lattice polynomial function p : Ln → L if an only if f is monotone; in this case, the interpolating polynomial p was shown to be unique. The interpolation problem was also considered in the more general set- ting where L is a distributive lattice, not necessarily bounded, and where D ⊆ Ln is allowed to range over cuboids D = {a1 , b1 } × · · · × {an , bn } with ai , bi ∈ L and ai < bi . In this case, the class of such partial functions that can be interpolated by lattice polynomial functions was completely described. In this paper, we extend these results by completely characterizing the class of lattice functions that can be interpolated by polynomial functions on arbitrary finite subsets D ⊆ Ln . As in the latter setting, interpolating polynomials are not necessarily unique. We provide explicit descriptions of all possible lattice polynomial functions that interpolate these lattice functions, when such an interpolation is available.