In this paper, we study the Martin boundary associated with a harmonic structure given by a coupled partial differential equations system. We give an integral representation for non negative harmonic functions of this structure. In particular, we obtain such results for biharmonic functions (i.e. ${▵}^2\varphi =0$) and for non negative solutions of the equation ${▵}^2\varphi =\varphi$.

We study the problem of construction of the smooth interpolation formula presented as the minimizer of suitable functionals subject to interpolation constraints. We present a procedure for determining the interpolation formula that in a natural way leads to a linear combination of polyharmonic splines complemented with lower order polynomial terms. In general, such formulae can be very useful e.g. in geographic information systems or computer aided geometric design. A simple computational example...