Prof. Frederic Chyzak,INRIA

We develop and compare two algorithms for computing first-order right-hand factors in the ring of linear Mahler operators l_r M^r+?+l_1 M+l_0 where l_0,…,l_r are polynomials in x and Mx=x^b M for some integer b≥2. In other words, we give algorithms for finding all formal infinite product solutions of linear functional equations l_r (x)f(x^(b^r ))+?+l_1 (x)f(x^b)+l_0 (x)f(x)=0.

The first of our algorithms is adapted from Petkov?ek's classical algorithm for the analogous problem in the case of linear recurrences. The second one proceeds by computing a basis of generalized power series solutions of the functional equation and by using Hermite--Padé approximants to detect those linear combinations of the solutions that correspond to first-order factors.

We present implementations of both algorithms and discuss their use in combination with criteria from the literature to prove the differential transcendance of power series solutions of Mahler equations.