Sobolev Orthogonal Polynomials: Asymptotics and Symbolic Computation

. The Sobolev polynomials, which are orthogonal with respect to an inner product involving derivatives, are considered. The theory about these nonstandard polynomials has been developed along the last 40 years. The local asymptotics of these polynomials can be described by the Mehler-Heine formulae, which connect the polynomials with the Bessel functions of the first kind. In recent years, the formulae have been computed for discrete Sobolev orthogonal polynomials in several particular cases. We improve various known results by unifying them. Besides, an algorithm to compute these formulae effectively is presented. The algorithm allows to construct a computer program based on Mathematica ® language, where the corresponding Mehler-Heine formulae are automatically obtained. Applications and examples show the efficiency of the approach developed.