Weighted Sobolev orthogonal polynomials on the unit ball
  For the weight function $W_\mu(x) = (1-|x|^2)^\mu$, $\mu > -1$, $\lambda > 0$ and $b_\mu$ a normalizing constant, a family of mutually orthogonal polynomials on the unit ball with respect to the inner product $$   \la f,g \ra = {b_\mu [\int_{\BB^d} f(x) g(x) W_\mu(x) dx +   \lambda \int_{\BB^d} \nabla f(x) \cdot \nabla g(x) W_\mu(x) dx]} $$ are constructed in terms of spherical harmonics and a sequence of Sobolev orthog onal polynomials of one variable. The latter ones, hence, the orthogonal polynomials with respect to $\la \cdot,\cdot\ra$, can be generated through a recursive formula. 
Teresa E. Perez, Miguel A. Pinar, Yuan Xu
math.CA