The Sobolev orthogonal polynomials, which are orthogonal with respect to an inner product involving derivatives, are nonstandard orthogonal polynomials that 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. Obtaining this type of asymptotics requires many calculations. Hence, we construct this algorithm that allows to compute the corresponding Mehler-Heine formulae automatically. Theoretical details, applications and examples are shown in "Sobolev Orthogonal Polynomials: Asymptotics and Symbolic Computation", East Asian J. Appl. Math. 12(3) (2022), pp. 535-563, doi: 10.4208/eajam.240221.130921.