In mathematics, Hölder's theorem states that the gamma function does not satisfy any algebraic differential equation whose coefficients are rational functions. This result was first proved by Otto Hölder in 1887; several alternative proofs have subsequently been found.
The theorem also generalizes to the                     q                -gamma function.
For every                     n        ∈                              N                                0                                  , there is no non-zero polynomial                     P        ∈                  C                [        X        ;                  Y                      0                          ,                  Y                      1                          ,        …        ,                  Y                      n                          ]                 such that
                    ∀        z        ∈                  C                ∖                              Z                                ≤            0                          :                P                          (          z          ;          Γ          (          z          )          ,                      Γ            ′                    (          z          )          ,          …          ,                                    Γ                              (                n                )                                              (          z          )          )                =        0        ,                where                     Γ                 is the gamma function.                             ◼                
For example, define                     P        ∈                  C                [        X        ;                  Y                      0                          ,                  Y                      1                          ,                  Y                      2                          ]                 by                     P                                                                     =                                            df                                                                       X                      2                                    Y                      2                          +        X                  Y                      1                          +        (                  X                      2                          −                  ν                      2                          )                  Y                      0                                  . Then the equation
                    P        (        z        ;        f        (        z        )        ,                  f          ′                (        z        )        ,                  f          ″                (        z        )        )        =                  z                      2                                    f          ″                (        z        )        +        z                  f          ′                (        z        )        +        (                  z                      2                          −                  ν                      2                          )        f        (        z        )        ≡        0                is called an algebraic differential equation, which, in this case, has the solutions                     f        =                  J                      ν                                   and                     f        =                  Y                      ν                                   — the Bessel functions of the first and second kind respectively. Hence, we say that                               J                      ν                                   and                               Y                      ν                                   are differentially algebraic (also algebraically transcendental). Most of the familiar special functions of mathematical physics are differentially algebraic. All algebraic combinations of differentially algebraic functions are differentially algebraic. Furthermore, all compositions of differentially algebraic functions are differentially algebraic. Hölder’s Theorem simply states that the gamma function,                     Γ                , is not differentially algebraic and is therefore transcendentally transcendental.
Let                     n        ∈                              N                                0                                  , and assume that a non-zero polynomial                     P        ∈                  C                [        X        ;                  Y                      0                          ,                  Y                      1                          ,        …        ,                  Y                      n                          ]                 exists such that
                    ∀        z        ∈                  C                ∖                              Z                                ≤            0                          :                P                          (          z          ;          Γ          (          z          )          ,                      Γ            ′                    (          z          )          ,          …          ,                                    Γ                              (                n                )                                              (          z          )          )                =        0.                As a non-zero polynomial in                               C                [        X        ]                 can never give rise to the zero function on any non-empty open domain of                               C                         (by the Fundamental Theorem of Algebra), we may suppose, without loss of generality, that                     P                 contains a monomial term having a non-zero power of one of the indeterminates                               Y                      0                          ,                  Y                      1                          ,        …        ,                  Y                      n                                  .
Assume also that                     P                 has the lowest possible overall degree with respect to the lexicographic ordering                               Y                      0                          <                  Y                      1                          <        …        <                  Y                      n                          <        X                . For example,
                    deg                          (          −          3                      X                          10                                            Y                          0                                      2                                            Y                          1                                      4                                +          i                      X                          2                                            Y                          2                                )                <        deg                          (          2          X                      Y                          0                                      3                                −                      Y                          1                                      4                                )                        because the highest power of                               Y                      0                                   in any monomial term of the first polynomial is smaller than that of the second polynomial.
Next, observe that
                                                                        ∀                z                ∈                                  C                                ∖                                                      Z                                                        ≤                    0                                                  :                                                                             P                                                  (                  z                  +                  1                  ;                  Γ                  (                  z                  +                  1                  )                  ,                                                            Γ                                              (                        1                        )                                                                              (                  z                  +                  1                  )                  ,                                                            Γ                                              (                        2                        )                                                                              (                  z                  +                  1                  )                  ,                  …                  ,                                                            Γ                                              (                        n                        )                                                                              (                  z                  +                  1                  )                  )                                                                                    =                                                             P                                                  (                  z                  +                  1                  ;                  z                  Γ                  (                  z                  )                  ,                  [                  z                  Γ                  (                  z                  )                                      ]                                          (                      1                      )                                                        ,                  [                  z                  Γ                  (                  z                  )                                      ]                                          (                      2                      )                                                        ,                  …                  ,                  [                  z                  Γ                  (                  z                  )                                      ]                                          (                      n                      )                                                        )                                                                                    =                                                             P                                                  (                  z                  +                  1                  ;                  z                  Γ                  (                  z                  )                  ,                  z                                                            Γ                                              (                        1                        )                                                                              (                  z                  )                  +                  Γ                  (                  z                  )                  ,                  z                                                            Γ                                              (                        2                        )                                                                              (                  z                  )                  +                  2                                                            Γ                                              (                        1                        )                                                                              (                  z                  )                  ,                  …                  ,                  z                                                            Γ                                              (                        n                        )                                                                              (                  z                  )                  +                  n                                                            Γ                                              (                        n                        −                        1                        )                                                                              (                  z                  )                  )                                .                                                            If we define a second polynomial                     Q        ∈                  C                [        X        ;                  Y                      0                          ,                  Y                      1                          ,        …        ,                  Y                      n                          ]                 by the transformation
                    Q                                                                     =                                            df                                                             P        (        X        +        1        ;        X                  Y                      0                          ,        X                  Y                      1                          +                  Y                      0                          ,        X                  Y                      2                          +        2                  Y                      1                          ,        …        ,        X                  Y                      n                          +        n                  Y                      n            −            1                          )        ,                then we obtain the following algebraic differential equation for                     Γ                :
                    ∀        z        ∈                  C                ∖                              Z                                ≤            0                          :                Q                          (          z          ;          Γ          (          z          )          ,                      Γ            ′                    (          z          )          ,          …          ,                                    Γ                              (                n                )                                              (          z          )          )                ≡        0.                Furthermore, if                               X                      h                                    Y                      0                                              h                              0                                                              Y                      1                                              h                              1                                                    ⋯                  Y                      n                                              h                              n                                                             is the highest-degree monomial term in                     P                , then the highest-degree monomial term in                     Q                 is                               X                      h            +                          h                              0                                      +                          h                              1                                      +            ⋯            +                          h                              n                                                              Y                      0                                              h                              0                                                              Y                      1                                              h                              1                                                    ⋯                  Y                      n                                              h                              n                                                            . Consequently, the polynomial
                    Q        −                  X                                    h                              0                                      +                          h                              1                                      +            ⋯            +                          h                              n                                                    P                has a smaller overall degree than                     P                , and as it clearly gives rise to an algebraic differential equation for                     Γ                , it must be the zero polynomial by the minimality assumption on                     P                . Hence, defining                     R        ∈                  C                [        X        ]                 by                     R                                                                     =                                            df                                                                       X                                    h                              0                                      +                          h                              1                                      +            ⋯            +                          h                              n                                                            , we get
                    Q        =        P        (        X        +        1        ;        X                  Y                      0                          ,        X                  Y                      1                          +                  Y                      0                          ,        X                  Y                      2                          +        2                  Y                      1                          ,        …        ,        X                  Y                      n                          +        n                  Y                      n            −            1                          )        =        R        (        X        )        ⋅        P        (        X        ;                  Y                      0                          ,                  Y                      1                          ,        …        ,                  Y                      n                          )        .                Now, let                     X        =        0                 in                     Q                 to obtain
                    Q        (        0        ;                  Y                      0                          ,                  Y                      1                          ,        …        ,                  Y                      n                          )        =        P        (        1        ;        0        ,                  Y                      0                          ,        2                  Y                      1                          ,        …        ,        n                  Y                      n            −            1                          )        =        R        (        0        )        ⋅        P        (        0        ;                  Y                      0                          ,                  Y                      1                          ,        …        ,                  Y                      n                          )        =                  0                                    C                        [                          Y                              0                                      ,                          Y                              1                                      ,            …            ,                          Y                              n                                      ]                          .                A change of variables then yields                     P        (        1        ;        0        ,                  Y                      1                          ,                  Y                      2                          ,        …        ,                  Y                      n                          )        =                  0                                    C                        [                          Y                              0                                      ,                          Y                              1                                      ,            …            ,                          Y                              n                                      ]                                  , and an application of mathematical induction (along with a change of variables at each induction step) to the earlier expression
                    P        (        X        +        1        ;        X                  Y                      0                          ,        X                  Y                      1                          +                  Y                      0                          ,        X                  Y                      2                          +        2                  Y                      1                          ,        …        ,        X                  Y                      n                          +        n                  Y                      n            −            1                          )        =        R        (        X        )        ⋅        P        (        X        ;                  Y                      0                          ,                  Y                      1                          ,        …        ,                  Y                      n                          )                reveals that
                    ∀        m        ∈                  N                :                P        (        m        ;        0        ,                  Y                      1                          ,                  Y                      2                          ,        …        ,                  Y                      n                          )        =                  0                                    C                        [                          Y                              0                                      ,                          Y                              1                                      ,            …            ,                          Y                              n                                      ]                          .                This is possible only if                     P                 is divisible by                               Y                      0                                  , which contradicts the minimality assumption on                     P                . Therefore, no such                     P                 exists, and so                     Γ                 is not differentially algebraic. Q.E.D.