A boundary Nevanlinna-Pick interpolation problem is posed and solved in the
quaternionic setting. Given nonnegative real numbers κ1,…,κN, quaternions p1,…,pN all of modulus 1, so that the
2-spheres determined by each point do not intersect and pu=1 for u=1,…,N, and quaternions s1,…,sN, we wish to find a slice
hyperholomorphic Schur function s so that r→1r∈(0,1)lims(rpu)=suforu=1,…,N, and
r→1r∈(0,1)lim1−r1−s(rpu)su≤κu,foru=1,…,N. Our arguments relies on the theory of slice hyperholomorphic
functions and reproducing kernel Hilbert spaces