Fuchsian theory

The Fuchsian theory of linear differential equations, which is named after Lazarus Immanuel Fuchs, provides a characterization of various types of singularities and the relations among them. At any ordinary point of a homogeneous linear differential equation of order [math]\displaystyle{ n }[/math] there exists a fundamental system of [math]\displaystyle{ n }[/math] linearly independent power series solutions. A non-ordinary point is called a singularity. At a singularity the maximal number of linearly independent power series solutions may be less than the order of the differential equation.

Generalized series solutions

The generalized series at [math]\displaystyle{ \xi\in\mathbb{C} }[/math] is defined by

[math]\displaystyle{ (z-\xi)^\alpha\sum_{k=0}^\infty c_k(z-\xi)^k, \text{ with } \alpha,c_k \in \mathbb{C} \text{ and } c_0\neq0, }[/math]

which is known as Frobenius series, due to the connection with the Frobenius series method. Frobenius series solutions are formal solutions of differential equations. The formal derivative of [math]\displaystyle{ z^\alpha }[/math], with [math]\displaystyle{ \alpha\in\mathbb{C} }[/math], is defined such that [math]\displaystyle{ (z^\alpha)'=\alpha z^{\alpha-1} }[/math]. Let [math]\displaystyle{ f }[/math] denote a Frobenius series relative to [math]\displaystyle{ \xi }[/math], then

[math]\displaystyle{ {d^nf \over d z^n} = (z-\xi)^{\alpha-n}\sum_{k=0}^\infty (\alpha+k)^{\underline{n}} c_k(z-\xi)^k, }[/math]

where [math]\displaystyle{ \alpha^{\underline{n}}:=\prod_{i=0}^{n-1}(\alpha-i) = \alpha(\alpha-1)\cdots(\alpha-n+1) }[/math] denotes the falling factorial notation.

Indicial equation

Let [math]\displaystyle{ f:=(z-\xi)^{\alpha}\sum_{k=0}^{\infty}c_k(z-\xi)^k }[/math] be a Frobenius series relative to [math]\displaystyle{ \xi \in \mathbb{C} }[/math]. Let [math]\displaystyle{ Lf=f^{(n)} + q_1f^{(n-1)} + \cdots + q_nf }[/math] be a linear differential operator of order [math]\displaystyle{ n }[/math] with one valued coefficient functions [math]\displaystyle{ q_1, \dots, q_n }[/math]. Let all coefficients [math]\displaystyle{ q_1,\dots,q_n }[/math] be expandable as Laurent series with finite principle part at [math]\displaystyle{ \xi }[/math]. Then there exists a smallest [math]\displaystyle{ N\in\mathbb{N} }[/math] such that [math]\displaystyle{ (z-\xi)^Nq_i }[/math] is a power series for all [math]\displaystyle{ i\in\{1,\dots, n\} }[/math]. Hence, [math]\displaystyle{ Lf }[/math] is a Frobenius series of the form [math]\displaystyle{ Lf=(z-\xi)^{\alpha-n-N}\psi(z) }[/math], with a certain power series [math]\displaystyle{ \psi(z) }[/math] in [math]\displaystyle{ (z-\xi) }[/math]. The indicial polynomial is defined by [math]\displaystyle{ P_{\xi}:=\psi(0) }[/math] which is a polynomial in [math]\displaystyle{ \alpha }[/math], i.e., [math]\displaystyle{ P_{\xi} }[/math] equals the coefficient of [math]\displaystyle{ Lf }[/math] with lowest degree in [math]\displaystyle{ (z-\xi) }[/math]. For each formal Frobenius series solution [math]\displaystyle{ f }[/math] of [math]\displaystyle{ Lf=0 }[/math], [math]\displaystyle{ \alpha }[/math] must be a root of the indicial polynomial at [math]\displaystyle{ \xi }[/math], i. e., [math]\displaystyle{ \alpha }[/math] needs to solve the indicial equation [math]\displaystyle{ P_{\xi}(\alpha) = 0 }[/math].

If [math]\displaystyle{ \xi }[/math] is an ordinary point, the resulting indicial equation is given by [math]\displaystyle{ \alpha^{\underline{n}}=0 }[/math]. If [math]\displaystyle{ \xi }[/math] is a regular singularity, then [math]\displaystyle{ \deg(P_{\xi}(\alpha))=n }[/math] and if [math]\displaystyle{ \xi }[/math] is an irregular singularity, [math]\displaystyle{ \deg(P_{\xi}(\alpha))\lt n }[/math] holds. This is illustrated by the later examples. The indicial equation relative to [math]\displaystyle{ \xi=\infty }[/math] is defined by the indicial equation of [math]\displaystyle{ \widetilde{L}f }[/math], where [math]\displaystyle{ \widetilde{L} }[/math] denotes the differential operator [math]\displaystyle{ L }[/math] transformed by [math]\displaystyle{ z=x^{-1} }[/math]which is a linear differential operator in [math]\displaystyle{ x }[/math], at [math]\displaystyle{ x=0 }[/math].

Example: Regular singularity

The differential operator of order [math]\displaystyle{ 2 }[/math], [math]\displaystyle{ Lf := f''+\frac{1}{z}f'+\frac{1}{z^2}f }[/math], has a regular singularity at [math]\displaystyle{ z=0 }[/math]. Consider a Frobenius series solution relative to [math]\displaystyle{ 0 }[/math], [math]\displaystyle{ f := z^\alpha(c_0 + c_1z + c_2 z^2 + \cdots) }[/math] with [math]\displaystyle{ c_0\neq0 }[/math].

[math]\displaystyle{ \begin{align} Lf & = z^{\alpha-2}(\alpha(\alpha-1)c_0 + \cdots) + \frac{1}{z}z^{\alpha-1}(\alpha c_0 + \cdots) + \frac{1}{z^2}z^{\alpha}(c_0 + \cdots) \\[5pt] & = z^{\alpha-2}c_0(\alpha(\alpha-1) + \alpha + 1) + \cdots. \end{align} }[/math]

This implies that the degree of the indicial polynomial relative to [math]\displaystyle{ 0 }[/math] is equal to the order of the differential equation, [math]\displaystyle{ \deg(P_0(\alpha)) = \deg(\alpha^2 + 1) = 2 }[/math].

Example: Irregular singularity

The differential operator of order [math]\displaystyle{ 2 }[/math], [math]\displaystyle{ Lf:=f''+\frac{1}{z^2}f' + f }[/math], has an irregular singularity at [math]\displaystyle{ z=0 }[/math]. Let [math]\displaystyle{ f }[/math] be a Frobenius series solution relative to [math]\displaystyle{ 0 }[/math].

[math]\displaystyle{ \begin{align} Lf & = z^{\alpha-2}(\alpha(\alpha-1)c_0 + \cdots) + \frac{1}{z^2}z^{\alpha-1}(\alpha c_0 + \cdots) + z^{\alpha}(c_0 + \cdots) \\[5pt] & = z^{\alpha-3} c_0 \alpha + z^{\alpha-2}(c_0\alpha(\alpha-1) + c_1) + \cdots. \end{align} }[/math]

Certainly, at least one coefficient of the lower derivatives pushes the exponent of [math]\displaystyle{ z }[/math] down. Inevitably, the coefficient of a lower derivative is of smallest exponent. The degree of the indicial polynomial relative to [math]\displaystyle{ 0 }[/math] is less than the order of the differential equation, [math]\displaystyle{ \deg(P_0(\alpha)) = \deg(\alpha) = 1 \lt 2 }[/math].

Formal fundamental systems

We have given a homogeneous linear differential equation [math]\displaystyle{ Lf=0 }[/math] of order [math]\displaystyle{ n }[/math] with coefficients that are expandable as Laurent series with finite principle part. The goal is to obtain a fundamental set of formal Frobenius series solutions relative to any point [math]\displaystyle{ \xi\in\mathbb{C} }[/math]. This can be done by the Frobenius series method, which says: The starting exponents are given by the solutions of the indicial equation and the coefficients describe a polynomial recursion. W.l.o.g., assume [math]\displaystyle{ \xi=0 }[/math].

Fundamental system at ordinary point

If [math]\displaystyle{ 0 }[/math] is an ordinary point, a fundamental system is formed by the [math]\displaystyle{ n }[/math] linearly independent formal Frobenius series solutions [math]\displaystyle{ \psi_1, z\psi_2, \dots, z^{n-1}\psi_{n} }[/math], where [math]\displaystyle{ \psi_i\in\mathbb{C}z }[/math] denotes a formal power series in [math]\displaystyle{ z }[/math] with [math]\displaystyle{ \psi(0)\neq0 }[/math], for [math]\displaystyle{ i\in\{1,\dots,n\} }[/math]. Due to the reason that the starting exponents are integers, the Frobenius series are power series.

Fundamental system at regular singularity

If [math]\displaystyle{ 0 }[/math] is a regular singularity, one has to pay attention to roots of the indicial polynomial that differ by integers. In this case the recursive calculation of the Frobenius series' coefficients stops for some roots and the Frobenius series method does not give an [math]\displaystyle{ n }[/math]-dimensional solution space. The following can be shown independent of the distance between roots of the indicial polynomial: Let [math]\displaystyle{ \alpha\in\mathbb{C} }[/math] be a [math]\displaystyle{ \mu }[/math]-fold root of the indicial polynomial relative to [math]\displaystyle{ 0 }[/math]. Then the part of the fundamental system corresponding to [math]\displaystyle{ \alpha }[/math] is given by the [math]\displaystyle{ \mu }[/math] linearly independent formal solutions

[math]\displaystyle{ \begin{align} & z^\alpha \psi_0 \\ & z^\alpha \psi_1 + z^\alpha\log(z)\psi_0\\ & z^\alpha \psi_2 + 2z^\alpha\log(z)\psi_1 + z^\alpha\log^2(z)\psi_0\\ & \qquad \vdots\\ & z^\alpha \psi_{\mu-1} + \cdots + \binom{\mu-1}{k} z^{\alpha}\log^k(z)\psi_{\mu-k} + \cdots + z^\alpha \log^{\mu-1}(z)\psi_0 \end{align} }[/math]

where [math]\displaystyle{ \psi_i\in\mathbb{C}z }[/math] denotes a formal power series in [math]\displaystyle{ z }[/math] with [math]\displaystyle{ \psi(0)\neq0 }[/math], for [math]\displaystyle{ i\in\{0,\dots,\mu-1\} }[/math]. One obtains a fundamental set of [math]\displaystyle{ n }[/math] linearly independent formal solutions, because the indicial polynomial relative to a regular singularity is of degree [math]\displaystyle{ n }[/math].

General result

One can show that a linear differential equation of order [math]\displaystyle{ n }[/math] always has [math]\displaystyle{ n }[/math] linearly independent solutions of the form

[math]\displaystyle{ \exp(u(z^{-1/s}))\cdot z^\alpha(\psi_0(z^{1/s}) + \cdots + \log^k(z) \psi_k(z^{1/s}) + \cdots + \log^{w}(z) \psi_w(z^{1/s})) }[/math]

where [math]\displaystyle{ s\in\mathbb{N}\setminus\{0\}, u(z)\in\mathbb{C}[z] }[/math] and [math]\displaystyle{ u(0)=0, \alpha\in\mathbb{C}, w\in\mathbb{N} }[/math], and the formal power series [math]\displaystyle{ \psi_0(z),\dots,\psi_w\in\mathbb{C}z }[/math].

[math]\displaystyle{ 0 }[/math] is an irregular singularity if and only if there is a solution with [math]\displaystyle{ u\neq 0 }[/math]. Hence, a differential equation is of Fuchsian type if and only if for all [math]\displaystyle{ \xi\in\mathbb{C}\cup\{\infty\} }[/math] there exists a fundamental system of Frobenius series solutions with [math]\displaystyle{ u=0 }[/math] at [math]\displaystyle{ \xi }[/math].

References

• Ince, Edward Lindsay (1956). Ordinary Differential Equations. New York, USA: Dover Publications. ISBN 9780486158211. 
• Poole, Edgar Girard Croker (1936). Introduction to the theory of linear differential equations. New York: Clarendon Press. 
• Tenenbaum, Morris; Pollard, Harry (1963). Ordinary Differential Equations. New York, USA: Dover Publications. pp. Lecture 40. ISBN 9780486649405. https://archive.org/details/ordinarydifferen00tene_0. 
• Horn, Jakob (1905). Gewöhnliche Differentialgleichungen beliebiger Ordnung. Leipzig, Germany: G. J. Göschensche Verlagshandlung. 
• Schlesinger, Ludwig Lindsay (1897). Handbuch der Theorie der linearen Differentialgleichungen (2. Band, 1. Teil). Leipzig, Germany: B. G.Teubner. pp. 241 ff.