An algorithm for solving some hypergeometric ODE classes.

Lemmus Chan



A simple algorithm for solving a given linear ODE, when it can be
mapped into the (Gauss) hypergeometric equation

          x (1 - x) y'' + ((a + b + 1) x - c) y' + a b y = 0

by means of
                                 n
                              a x  + b
                 TR := { x -> -------- ,   y -> P(x) y }
                                 n
                              c x  + d

is presented. In above, P(x) is an arbitrary function, {a,b,c,d} are
arbitrary constants and n is a rational number. The algorithm is based
on the observation that, if such a transformation TR exists, then a
simpler transformation of the form


                                 n
            power_TR := { x -> x ,    y -> P(x) y },

mapping the given ODE into another one with 3 regular singular points,
also exists. At the core there is the fact that the transformation TR
can be decomposed into power_TR and a Mobius transformation, and the
latter preserves the structure of the singularities of the
hypergeometric equation. Although the idea is simple, it allows to
systematically solve - when such a TR exists - varied ODE problems not
solved by the present Maple or Mathematica ODE solvers.

This is joint work with Edgardo Cheb-Terrab.