Fact-checked by Grok 2 weeks ago

Frobenius method

The Frobenius method is a mathematical procedure for finding power series solutions to second-order linear homogeneous ordinary differential equations of the form y'' + p(x) y' + q(x) y = 0 that possess a regular singular point, typically at x = 0.^[1]^[2] It assumes a solution of the form y(x) = x^r \sum_{n=0}^\infty a_n x^n, where the exponent r is determined by substituting the series into the equation and equating the coefficient of the lowest power of x to zero, yielding the indicial equation r(r-1) + p(0) r + q(0) = 0.^[3]^[4] This approach guarantees at least one analytic solution near the singular point when the conditions for a regular singularity hold, namely that x p(x) and x^2 q(x) are analytic at the point.^[1]^[2] Named after the German mathematician Ferdinand Georg Frobenius (1849–1917), the method was introduced in his 1873 paper "Ueber die Integration der linearen Differentialgleichungen durch Reihen," published in the Journal für die reine und angewandte Mathematik.^[2] Frobenius developed it as an extension of earlier power series techniques, building on work by mathematicians like Fuchs to handle singularities where standard Taylor series expansions fail due to non-analytic coefficients.^[2] The procedure involves solving the indicial equation for roots r_1 and r_2, then deriving a recurrence relation for the coefficients a_n from higher-order terms; depending on the roots—whether distinct and non-integer separated, equal, or differing by an integer—the second linearly independent solution may take a similar series form, require a logarithmic term, or involve reduction of order.^[3]^[4] The Frobenius method is fundamental in the theory of special functions, providing series representations for solutions to equations like Bessel's differential equation, which arises in problems involving cylindrical symmetry in physics, such as wave propagation and heat conduction.^[2]^[1] It also applies to Legendre's equation and other hypergeometric-type equations, yielding functions essential in quantum mechanics, electrostatics, and approximation theory.^[5] By enabling explicit series solutions near singularities, the method bridges analytic and asymptotic analysis, influencing modern computational approaches to differential equations in applied mathematics.^[5]

Introduction

Definition and Scope

The Frobenius method is a mathematical technique used to find series solutions for second-order linear homogeneous ordinary differential equations (ODEs) that have a regular singular point at x = 0. It assumes a solution of the form

y(x) = x^r \sum_{n=0}^\infty a_n x^n,

where r is a real exponent to be determined (often non-integer), and the coefficients a_n (with a_0 \neq 0) are calculated via recurrence relations. This method applies to ODEs expressed in the standard form

x^2 y'' + x p(x) y' + q(x) y = 0,

where p(x) and q(x) are analytic functions at x = 0.^[6]^[7] The scope of the Frobenius method is limited to linear homogeneous second-order ODEs with coefficients that are analytic throughout the complex plane except at isolated regular singular points. It provides solutions valid in a punctured disk around the singularity, typically for x > 0 or x < 0, and can yield one or two linearly independent series solutions depending on the roots of the associated indicial equation. Unlike methods for irregular singular points, it does not guarantee convergence on both sides of the singularity but excels in handling equations arising in physics, such as those in quantum mechanics and wave propagation.^[6]^[7] Central assumptions for the method include the ODE being in normal form y'' + P(x) y' + Q(x) y = 0, where x = 0 is a regular singular point if x P(x) and x^2 Q(x) are both analytic at x = 0. Equivalently, in the multiplied form above, p(x) and q(x) must possess Taylor series expansions around x = 0. These conditions ensure that the singularity is mild enough for the series to capture the solution's asymptotic behavior without logarithmic terms in the leading cases.^[6]^[7]^[8] This method is motivated by the limitations of ordinary power series solutions, which converge only at ordinary points where coefficients are analytic and fail at singular points due to divergent or undefined terms. By incorporating the flexible exponent r, the Frobenius approach systematically resolves these issues, enabling the derivation of fundamental solutions for equations with variable coefficients that model real-world phenomena like heat conduction in irregular domains.^[6]^[7]

Historical Context

The Frobenius method, a powerful technique for obtaining series solutions to linear ordinary differential equations near regular singular points, was developed by the German mathematician Ferdinand Georg Frobenius in 1873. In his foundational paper "Über die Integration der linearen Differentialgleichungen durch Reihen," published in the Journal für die reine und angewandte Mathematik, Frobenius outlined a streamlined procedure to construct convergent power series solutions, emphasizing the role of an indicial equation to determine the leading exponents. This work marked a key advancement in addressing equations with rational coefficients that exhibit singularities, providing explicit convergence proofs for the resulting series.^[9] Frobenius' approach was deeply influenced by prior developments in the theory of differential equations. Leonhard Euler's investigations into equidimensional (Cauchy-Euler) equations during the mid-18th century introduced the characteristic indicial equation, which assumes solutions of the form y = x^r and yields roots dictating the series behavior—ideas that Frobenius adapted and generalized for more variable coefficients. Complementing this, Lazarus Fuchs' 1866 papers established the classification of regular singular points, where solutions remain analytic up to algebraic growth, and demonstrated that linear equations with finitely many such points admit solutions expressible as series involving powers and logarithms. Frobenius explicitly built on Fuchs' framework, simplifying the derivation of these solutions while recovering and extending Fuchs' results through a more direct algebraic manipulation.^[10] Following its introduction, the method evolved rapidly within the broader landscape of 19th-century analysis, generalizing Fuchs' singularity theory by enabling uniform treatment of second-order equations and higher systems. By the late 1880s, extensions addressed cases with repeated indicial roots, incorporating logarithmic terms into the series solutions to handle resonant behaviors at singularities. This refinement proved instrumental for applications to special functions; for instance, Henri Poincaré and contemporaries employed the method in the 1880s and 1890s to derive series representations for Bessel functions—arising in wave propagation—and hypergeometric functions, which underpin many physical models and unify numerous classical solutions, thereby solidifying the method's role in theoretical physics and pure mathematics.^[11]

Prerequisites

Power Series Solutions

The power series method provides a systematic approach to finding solutions to linear ordinary differential equations (ODEs) with analytic coefficients around points where the equation behaves regularly. An ordinary point x_0 for the second-order linear ODE y'' + P(x) y' + Q(x) y = 0 is defined as a point where both coefficient functions P(x) and Q(x) are analytic, meaning they possess Taylor series expansions that converge in some neighborhood of x_0.^[12]^[13] At such points, solutions can be expressed as power series of the form y(x) = \sum_{n=0}^{\infty} a_n (x - x_0)^n, where the coefficients a_n are determined recursively, yielding two linearly independent solutions for the general second-order case.^[14]^[15] To apply the method, one assumes the series solution y(x) = \sum_{n=0}^{\infty} a_n (x - x_0)^n and its term-by-term derivatives y'(x) = \sum_{n=1}^{\infty} n a_n (x - x_0)^{n-1} and y''(x) = \sum_{n=2}^{\infty} n(n-1) a_n (x - x_0)^{n-2}. Substituting these into the ODE y'' + P(x) y' + Q(x) y = 0, where P(x) and Q(x) are expanded as their own power series \sum_{k=0}^{\infty} p_k (x - x_0)^k and \sum_{k=0}^{\infty} q_k (x - x_0)^k, results in a single power series equation in powers of (x - x_0). Equating the coefficient of each power (x - x_0)^m to zero for m \geq 0 produces a recurrence relation that expresses higher coefficients a_{n+r} in terms of previous ones, typically for second-order equations taking the form a_{n+2} = f(a_n, a_{n+1}, \dots) after shifting indices and collecting terms, with the initial coefficients a_0 and a_1 arbitrary to generate the two independent solutions.^[14]^[16]^[12] For instance, consider the general form y'' + P(x) y' + Q(x) y = 0 with P(x) and Q(x) analytic at x_0 = 0. The substitution yields the coefficient equation for the lowest power (often involving a_0 and a_1) and subsequent recurrences, such as in the simple case y'' - y = 0, where the recurrence simplifies to a_{n+2} = a_n for n \geq 0, leading to solutions involving hyperbolic functions when summed.^[15]^[14] The power series solutions are guaranteed to converge in a disk of radius at least equal to the distance from x_0 to the nearest singularity of P(x) or Q(x) in the complex plane, ensuring analyticity of the solution within that region.^[13]^[12] This method succeeds precisely because the analyticity of the coefficients allows the series manipulations to preserve the equation's structure, but it encounters difficulties at points where P(x) or Q(x) have singularities, necessitating extensions like the Frobenius method for such cases.^[16]^[15]

Singular Points in ODEs

In the theory of linear ordinary differential equations (ODEs), particularly second-order equations of the form y'' + P(x) y' + Q(x) y = 0, a point x_0 is classified as a singular point if at least one of the coefficient functions P(x) or Q(x) fails to be analytic at x_0. Analyticity here means the function can be represented by a power series with a positive radius of convergence in a neighborhood of x_0. At ordinary points, where both P(x) and Q(x) are analytic, solutions can be expressed as convergent power series around x_0, allowing for straightforward local analysis. Singular points, however, disrupt this regularity and require specialized techniques to understand solution behavior nearby.^[17] Singular points are further subdivided into regular singular points and irregular singular points based on the nature of the singularities in the coefficients. A singular point x_0 is a regular singular point if the modified coefficients (x - x_0) P(x) and (x - x_0)^2 Q(x) are both analytic at x_0. This condition ensures that the singularities are "mild" or of pole type with limited order: P(x) has at most a simple pole, and Q(x) has at most a pole of order 2 at x_0. If these conditions are not satisfied—for instance, if P(x) has a pole of order greater than 1 or Q(x) has a pole of order greater than 2, or if there are essential singularities—the point is an irregular singular point. Irregular singularities lead to more complex solution behaviors, often involving essential singularities or rapid growth/decay that cannot be captured by simple series expansions.^[18] Fuchs' theorem provides a foundational result for handling regular singular points in second-order linear ODEs. Named after Lazarus Fuchs, the theorem states that if x_0 is a regular singular point, then the equation possesses at least one solution of the form y(x) = (x - x_0)^r \sum_{n=0}^\infty a_n (x - x_0)^n, where r is a constant (determined by an indicial equation) and the series converges in some punctured neighborhood of x_0; the second linearly independent solution may take a similar form or involve a logarithmic term if the indicial roots differ by an integer. This guarantees that solutions near regular singular points exhibit algebraic branching or mild logarithmic growth, facilitating explicit construction via series methods. For Fuchsian equations—those with only regular singular points (finite or infinite in number)—the theorem extends to global analytic continuation along paths avoiding singularities.^[17] To illustrate, consider the equation x y'' + y' + y = 0, which in standard form is y'' + \frac{1}{x} y' + \frac{1}{x} y = 0. Here, x = 0 is a singular point since P(x) = 1/x and Q(x) = 1/x are not analytic at 0, but x P(x) = 1 and x^2 Q(x) = x are both analytic, confirming it as a regular singular point. In contrast, the equation y'' + \frac{1}{x^3} y = 0 has P(x) = 0 (analytic) but Q(x) = 1/x^3, so x^2 Q(x) = 1/x, which is not analytic at 0 due to the pole of order 1; thus, x = 0 is an irregular singular point.^[18] The distinction between regular and irregular singular points is crucial for solution strategies and analytic properties. At regular singular points, the guaranteed form of solutions enables analytic continuation around the point (with possible branch cuts), preserving much of the structure seen at ordinary points. Irregular singular points, however, typically demand asymptotic analysis or more advanced tools like the method of asymptotic expansions, as series solutions may diverge or fail to capture the essential behavior, limiting predictability and continuation. This classification underpins the development of methods tailored to each case, ensuring rigorous treatment of ODEs arising in physics and engineering.^[17]

Core Method

Series Expansion and Indicial Equation

The Frobenius method seeks series solutions to second-order linear ordinary differential equations with a regular singular point at x = 0, assuming a solution of the form

y(x) = x^r \sum_{n=0}^{\infty} a_n x^n = \sum_{n=0}^{\infty} a_n x^{n+r},

where r is a constant exponent to be determined and the coefficients a_n satisfy a_0 \neq 0.^[2]^[19] This ansatz generalizes the ordinary power series method by incorporating the factor x^r to account for the singularity, as introduced by Ferdinand Georg Frobenius in his 1873 paper on integrating linear differential equations via series.^[20] The method applies to equations writable in the form

x^2 y'' + x p(x) y' + q(x) y = 0,

where p(x) and q(x) are analytic at x = 0, expanded as power series p(x) = \sum_{k=0}^{\infty} p_k x^k and q(x) = \sum_{k=0}^{\infty} q_k x^k.^[19] Substituting the ansatz yields

y'(x) = \sum_{n=0}^{\infty} (n + r) a_n x^{n + r - 1}, \quad y''(x) = \sum_{n=0}^{\infty} (n + r)(n + r - 1) a_n x^{n + r - 2}.

Plugging these into the equation and multiplying through by x^{-r} to align powers produces a series in non-negative powers of x, starting from the lowest term.^[2]^[19] The coefficient of the lowest power, x^0, must vanish for the series to satisfy the equation, giving the indicial equation

r(r - 1) + p_0 r + q_0 = 0,

a quadratic equation in r with roots r_1 and r_2 (assuming r_1 \geq r_2).^[21]^[19] This equation determines the possible leading exponents for the solutions. If r_1 \neq r_2 and r_1 - r_2 is not an integer, two linearly independent Frobenius series solutions exist: one is y_1(x) = x^{r_1} \sum_{n=0}^{\infty} a_n x^n with a_0 \neq 0, and the other is y_2(x) = x^{r_2} \sum_{n=0}^{\infty} b_n x^n with b_0 \neq 0.^[2]^[20] For convenience, the leading coefficient is often normalized to a_0 = 1 or b_0 = 1.^[19]

Recurrence Relations for Coefficients

Once the indicial equation has been solved to determine the possible values of the exponent r, the coefficients a_n for n \geq 1 in the Frobenius series solution y(x) = \sum_{n=0}^{\infty} a_n x^{r+n} are computed using a recurrence relation derived from substituting the series into the differential equation and equating the coefficients of x^{r+n} to zero.^[22]^[4] Consider the second-order linear ODE in standard form for a regular singular point at x = 0: x^2 y'' + x p(x) y' + q(x) y = 0, where p(x) = \sum_{k=0}^{\infty} p_k x^k and q(x) = \sum_{k=0}^{\infty} q_k x^k are analytic at x = 0. Substituting the assumed series yields the general recurrence relation for n \geq 1:

[(r + n)(r + n - 1) + p_0 (r + n) + q_0] a_n = -\sum_{k=1}^{n} \left[ p_k (r + n - k) + q_k \right] a_{n-k},

or equivalently,

a_n = -\frac{1}{F(r + n)} \sum_{k=1}^{n} \left[ p_k (r + n - k) + q_k \right] a_{n-k},

where F(s) = s(s-1) + p_0 s + q_0 is the indicial polynomial, and the sum runs over previous coefficients with the understanding that a_m = 0 for m < 0.^[22]^[4] This form arises from collecting terms contributing to the power x^{r+n}, where the leading term involves a_n multiplied by the Euler operator applied to x^{r+n}, and the sum accounts for contributions from higher-order terms in the expansions of p(x) and q(x).^[3] The recurrence is typically a multi-term relation, but for ODEs where p(x) and q(x) are polynomials of low degree, it simplifies to a two- or three-term recurrence, facilitating explicit computation. The coefficients are determined sequentially: a_1 is found in terms of a_0 (which is arbitrary and often set to 1 for normalization), then a_2 in terms of a_0 and a_1, and so on, assuming F(r + n) \neq 0 for all n \geq 1 to ensure the denominator is nonzero.^[22] This sequential process generates the full power series solution around the singular point. When the indicial roots r_1 and r_2 are distinct with a non-integer difference, the method yields two linearly independent solutions: one series with root r_1 and coefficients \{a_n^{(1)}\}, and another with root r_2 and coefficients \{a_n^{(2)}\}, each satisfying their respective recurrence relations starting from their own a_0.^[4]^[22] The general solution is then a linear combination of these two series.^[23]

Standard Applications

Solution for Non-Integer Root Differences

When the roots r_1 and r_2 (with r_1 > r_2) of the indicial equation differ by a non-integer value, the Frobenius method produces two distinct Frobenius series solutions that are valid and linearly independent. These solutions take the form

y_1(x) = x^{r_1} \sum_{n=0}^\infty a_n x^n, \quad y_2(x) = x^{r_2} \sum_{n=0}^\infty b_n x^n,

where a_0 and b_0 are arbitrary nonzero constants, and the higher coefficients a_n (n \geq 1) and b_n (n \geq 1) are determined recursively by substituting each series into the original second-order linear ODE with a regular singular point.^[24]^[3] The recurrence relations for the coefficients of each series, as derived in the core Frobenius procedure, operate without complications in this case. Specifically, the non-integer difference r_1 - r_2 ensures that the denominators in the recurrence formulas—typically involving terms like (n + r_1 - r_2) or similar—do not vanish for any nonnegative integer n, avoiding indeterminate forms or overlaps between the series terms.^[25]^[4] This allows the coefficients to be computed straightforwardly for both series, yielding well-defined power series expansions around the singular point.^[26] Linear independence of y_1 and y_2 follows from the distinct leading exponents r_1 and r_2, which prevent one series from being a scalar multiple of the other. The Wronskian W(y_1, y_2), computed as y_1 y_2' - y_2 y_1', is nonzero in a neighborhood of the singular point due to the non-integer separation, confirming that y_1 and y_2 form a fundamental set of solutions.^[24]^[3] Consequently, the general solution to the ODE is y(x) = c_1 y_1(x) + c_2 y_2(x), where c_1 and c_2 are arbitrary constants.^[25]^[4] To verify the solutions, substituting y_1 (or similarly y_2) and its derivatives into the ODE results in an identity that holds term by term after equating coefficients of like powers of x, with the indicial equation ensuring balance at the lowest order and the recurrences handling all higher orders.^[24]^[26] This term-by-term satisfaction confirms that both series are exact solutions within their radius of convergence.^[3]

Example: Bessel's Differential Equation

Bessel's differential equation is a second-order linear ordinary differential equation of the form

x^2 y'' + x y' + (x^2 - \nu^2) y = 0,

where \nu is a real parameter known as the order. This equation has a regular singular point at x = 0.^[27] To apply the Frobenius method, assume a solution of the form y(x) = x^r \sum_{n=0}^\infty a_n x^n, with a_0 \neq 0. Substituting this series into the differential equation yields the indicial equation r^2 - \nu^2 = 0.^[28] The roots of the indicial equation are r = \nu and r = -\nu. When \nu is such that the roots differ by a non-integer (i.e., $2\nu not integer), the method produces two linearly independent Frobenius series solutions. For half-integer orders (where the difference is integer), independent series solutions still exist.^[27] Consider the larger root r = \nu. The corresponding solution is y_1(x) = x^\nu \sum_{n=0}^\infty a_n x^n, where the coefficients satisfy the recurrence relation a_n = -\frac{a_{n-2}}{n(n + 2\nu)} for n \geq 2, with a_0 arbitrary and a_1 = 0 (hence all odd-indexed coefficients vanish).^[28] For the smaller root r = -\nu, the solution is y_2(x) = x^{-\nu} \sum_{n=0}^\infty b_n x^n, with the recurrence b_n = -\frac{b_{n-2}}{n(n - 2\nu)} for n \geq 2, b_0 arbitrary, and b_1 = 0 (odd coefficients vanish).^[27] These solutions are the Bessel functions of the first kind, denoted J_\nu(x) = y_1(x) (normalized with a_0 = \frac{1}{2^\nu \Gamma(\nu+1)}) and J_{-\nu}(x) = y_2(x) (with appropriate normalization). For fixed non-integer \nu, both J_\nu(x) and J_{-\nu}(x) are defined via power series in fractional powers that converge for all finite x, though they have a branch point at x=0.^[29]

Exceptional Cases

Integer Root Differences

When the roots r_1 and r_2 of the indicial equation satisfy r_1 - r_2 = N, where N is a nonnegative integer and r_1 > r_2, the Frobenius method encounters a special case that requires careful handling of the second solution.^[30]^[31] The series solution corresponding to the larger root r_1 always succeeds, yielding a valid Frobenius solution of the form

y_1(x) = x^{r_1} \sum_{n=0}^{\infty} a_n x^n,

with coefficients a_n determined by the recurrence relation without issue.^[4]^[30] For the smaller root r_2, an attempt is made to construct a second solution y_2(x) = x^{r_2} \sum_{n=0}^{\infty} b_n x^n using the same recurrence relation for the coefficients b_n.^[31]^[4] However, at the step n = N in the recurrence, the denominator vanishes because it includes a factor of (r_2 + n - r_1) = 0, potentially causing the computation of b_N to fail.^[30]^[31] If the corresponding numerator in the recurrence is also zero at this step, the indeterminacy allows b_N to be chosen freely, and the series can be continued indefinitely, resulting in two linearly independent power series solutions without logarithmic terms.^[4]^[31] If the numerator does not vanish or the recurrence leads to an inconsistency (such as division by zero with a nonzero numerator), the attempted series for y_2 cannot be completed independently and includes terms proportional to y_1, making it linearly dependent.^[30]^[31] In this failure scenario, the second linearly independent solution takes the form

y_2(x) = c \, y_1(x) \ln x + x^{r_2} \sum_{n=0}^{\infty} c_n x^n,

where c is a nonzero constant and the coefficients c_n are determined by substituting into the original equation.^[4]^[30] To detect whether the recurrence for y_2 succeeds or fails, one computes the coefficients up to n = N - 1 and examines the term for b_N; specifically, if this coefficient cannot be defined due to a zero denominator without a matching zero numerator, the case requires the logarithmic adjustment.^[31]^[4] This check ensures the appropriate form for the general solution is used while maintaining linear independence.^[30]

Logarithmic Solutions and Tandem Recurrences

When the roots of the indicial equation differ by a positive integer and the recurrence relation for the series solution corresponding to the smaller root r_2 encounters a division by zero (recurrence failure), the second linearly independent solution incorporates a logarithmic term.^[25]^[32] This arises because the power series for r_2 terminates prematurely or becomes inconsistent without the additional logarithmic factor to balance the equation.^[4] The general form of this second solution is

y_2(x) = A y_1(x) \ln x + x^{r_2} \sum_{n=0}^\infty c_n x^n,

where y_1(x) = x^{r_1} \sum_{n=0}^\infty a_n x^n is the first solution (with r_1 > r_2 and r_1 - r_2 = m, a positive integer), and A is a constant determined by substituting into the original differential equation.^[25]^[4] The constant A is determined by substituting the assumed form of y_2 into the original differential equation, which introduces terms that resolve the recurrence inconsistency and ensure the equation holds order by order.^[32]^[33] To find the coefficients c_n, the form leads to tandem recurrence relations, which couple the unknown c_n with the known coefficients a_n of y_1. Substituting y_2 into the ODE produces terms involving derivatives of the logarithmic part: specifically, y_2' includes A y_1' \ln x + A \frac{y_1}{x} + series terms, and y_2'' includes additional A y_1'' \ln x + 2A \frac{y_1'}{x} - A \frac{y_1}{x^2} + series terms.^[25] Equating coefficients of like powers of x for n \geq 0 results in two coupled equations per order: one from the indicial-like balance and another from the recurrence, typically of the form

(r_2 + n)(r_2 + n - 1) c_n + \sum_{k=0}^{n-1} (p_{n-k} (r_2 + k) + q_{n-k}) c_k + A \left[ \text{terms from } y_1'' \ln x \text{ and derivatives} \right] = 0,

where p_j and q_j are the series coefficients of the ODE's normalized form x^2 y'' + x p(x) y' + q(x) y = 0.^[4]^[32] These relations are solved sequentially, with the logarithmic contributions acting as inhomogeneous terms derived from differentiating y_1. In the special case of repeated roots (r_1 = r_2 = r), the form simplifies to y_2(x) = y_1(x) \ln x + x^r \sum_{n=0}^\infty c_n x^n, where the constant A = 1 by convention, and the coefficients satisfy c_n = \frac{da_n}{dr} \big|_{r=r_1}, obtained by differentiating the original recurrence for a_n(r) with respect to r and evaluating at the repeated root.^[4]^[32] For instance, if the recurrence for y_1 is (r + n) a_n = - \sum_{k=0}^{n-1} b_{n-k} a_k, then the tandem relation for c_n becomes (r + n) c_n + a_n = - \sum_{k=0}^{n-1} b_{n-k} c_k, decoupling to yield the derivative form.^[32] Linear independence of y_1 and y_2 is verified using the Wronskian determinant, which for y_2 = y_1 \ln x + series evaluates to W(y_1, y_2) = \frac{y_1^2}{x} (up to constants), nonzero in the domain where y_1 \neq 0.^[25] This confirms the solutions span the two-dimensional solution space near the regular singular point.^[4]

Convergence and Limitations

Radius of Convergence

The Frobenius series solutions to a second-order linear differential equation with a regular singular point at x = 0, expressed as y(x) = x^r \sum_{n=0}^\infty a_n x^n, converge at least in the disk |x| < \rho, where \rho is the minimum of the radii of convergence of the power series expansions of x p(x) and x^2 q(x) for the equation y'' + p(x) y' + q(x) y = 0. This \rho > 0 ensures analyticity of these coefficient functions in that disk, guaranteeing the validity of the series substitution and coefficient recurrence relations. For a regular singular point at the origin, the radius extends precisely to the distance from 0 to the nearest other singularity of p(x) or q(x) in the complex plane, beyond which the solution may exhibit singularities that limit further convergence.^[34]^[4] The radius of convergence can be determined explicitly using the ratio test on the series coefficients. For large n, the recurrence relation for a_n typically yields \left| \frac{a_{n+1}}{a_n} \right| \sim \frac{1}{\rho} as n \to \infty, implying that the series \sum a_n x^n converges for |x| < \rho by the root test criterion \lim_{n \to \infty} \left| \frac{a_{n+1} x^{n+1}}{a_n x^n} \right| = \frac{|x|}{\rho} < 1. This behavior mirrors that of ordinary power series solutions around regular points, with the Frobenius factor x^r not affecting the radius since r is finite. In practice, for equations like the Legendre or Bessel equations, application of the ratio test confirms \rho = 1, corresponding to singularities at x = \pm 1.^[25]^[4] When the indicial roots differ by an integer or are equal, one solution may involve a logarithmic term, such as y_2(x) = y_1(x) \ln |x| + x^r \sum_{n=0}^\infty b_n x^n. The power series component converges in the same disk |x| < \rho as y_1(x), but the \ln |x| factor introduces a branch point at x = 0, rendering the solution multi-valued with a branch cut typically emanating from the origin along the negative real axis. This does not alter the radius of convergence but restricts the domain to slit planes excluding the cut, ensuring analytic continuation around other points within \rho.^[33] Fuchs' relation provides a global constraint on the convergence and analytic structure across all singular points: for a second-order Fuchsian equation on the Riemann sphere with s regular singularities (including infinity), the sum over all singularities of the sums of the indicial exponents r_1 + r_2 equals s - 2. This relation ensures compatibility of local Frobenius solutions for global meromorphic continuation, linking local radii at each singularity to the overall distribution of poles and zeros in the coefficient functions. Violations indicate irregular singularities, where Frobenius series fail to converge adequately.^[35]

Irregular Singular Points

In the theory of linear ordinary differential equations, an irregular singular point at x = x_0 for the equation P(x) y'' + Q(x) y' + R(x) y = 0 is a singular point where at least one of (x - x_0) Q(x)/P(x) or (x - x_0)^2 R(x)/P(x) fails to be analytic at x_0.^[4] This contrasts with regular singular points, where both expressions are analytic. A classic example is the equation y'' + \frac{1}{x^3} y' = 0 at x = 0, where the coefficient of y' introduces a pole of order higher than expected for regularity.^[36] The Frobenius method, which assumes solutions of the form y(x) = (x - x_0)^r \sum_{n=0}^\infty a_n (x - x_0)^n, fails at irregular singular points because the resulting recurrence relations for the coefficients a_n become infinite in extent or divergent, preventing the determination of a finite power series solution.^[37] Instead, solutions typically exhibit essential singularities, often involving exponential factors that lead to rapid growth, decay, or oscillations not captured by power series.^[37] The degree of this irregularity is quantified by the Poincaré rank p, defined for a system Y' = A(x) Y with singularity at x = 0 as the integer p where A(x) = x^{-p} (A_0 + A_1 x + \cdots) and A_0 \neq 0; a rank of p = 0 corresponds to a regular singular point.^[38] To address irregular singular points, alternative approaches are employed, such as asymptotic expansions of the form y(x) = \exp(F(x)) Y(x), where F(x) is a polynomial capturing the leading exponential behavior and Y(x) is a formal series.^[37] The WKB (Wentzel-Kramers-Brillouin) method provides approximate solutions near such points by constructing asymptotic series tailored to the singularity's rank, particularly useful for high-order poles.^[39] Formal power series with negative powers may also be used, though they generally diverge. These methods imply that no closed-form convergent series solutions exist via Frobenius, necessitating numerical integration or perturbative techniques for practical computation and analysis.^[37]

References

[1]
https://math.libretexts.org/Bookshelves/Differential_Equations/Differential_Equations_for_Engineers_(Lebl)/7%3A_Power_series_methods/7.3%3A_Singular_Points_and_the_Method_of_Frobenius
[2]
Frobenius Method -- from Wolfram MathWorld
If x_0 is an ordinary point of the ordinary differential equation, expand y in a Taylor series about x_0. Commonly, the expansion point can be taken as ...
[3]
[PDF] THE METHOD OF FROBENIUS
The first step in using the method of Frobenius is to determine the values of r that allow us to solve the differential equation. We do this by exploiting the ...Missing: ordinary | Show results with:ordinary
[4]
[PDF] Method of Frobenius | Math 537 - Ordinary Differential Equations
Definition (Ordinary and Singular Points) x0 is an ordinary point of Eqn. (1) if P(x0) 6= 0 and Q(x)/P(x), R(x)/P(x), and. F(x)/P(x) are analytic at x0.
[5]
[PDF] Power Series and Differential Equations: The Method of Frobenius∗
Jan 29, 2014 · It is a procedure for finding power series formulas for solutions to a great many differential equations of interest in mathematical physics.
[6]
7.3: Singular Points and the Method of Frobenius - Math LibreTexts
Feb 23, 2025 · While behavior of ODEs at singular points is more complicated, certain singular points are not especially difficult to solve.Examples · Method of Frobenius · Theorem $\PageIndex{1}$ · Bessel Functions
[7]
Section - 4.3 Singular points and the method of Frobenius
If the indicial equation has two real roots r 1 and r 2 such that r 1 − r 2 is not an integer, then we have two linearly independent Frobenius-type solutions.
[8]
[PDF] Overview of the Frobenius method
(ii) A point x。 is called a singular point of the ode if a, (z) or a₂(z) is. NOT analytic at xo. In this case,. Xo is either. (a.) a regular singular point.
[9]
A geometric-analytic study of linear differential equations of order two
Frobenius, Ueber die Integration der linearen Differentialgleichungen durch Reihen, J. Reine Angew. Math., 76 (1873), 214-235. doi: 10.1515/crll.1873.76.214 ...
[10]
What′s in a Name: Why Cauchy and Euler Share the Cauchy ...
Nov 27, 2017 · The fact that that Cauchy was born after Euler's death should give pause as to why they are both identified with a well-known differential ...
[11]
Tandem Recurrence Relations for Coefficients of Logarithmic ... - MDPI
In his 1873 paper Frobenius merely set himself the task to recover, in a more direct and simpler way, results about the solutions of differential equations ...
[12]
[PDF] Power Series Solutions and Equations with Analytic Coefficients
Jun 28, 2016 · A point x0 is called an ordinary point of equation (2) if both p and q are analytic at x0. If x0 is not an ordinary point, it is called a ...
[13]
[PDF] 640:527 POWER SERIES FALL 2014 Section 10. Ordinary and ...
We say that a point x0 is an ordinary point for the differential equation (14) (or equivalently for (15)) if both p(x) and q(x) are analytic at x0.
[14]
Differential Equations - Series Solutions - Pauls Online Math Notes
Nov 16, 2022 · The basic idea to finding a series solution to a differential equation is to assume that we can write the solution as a power series.
[15]
[PDF] Chapter 7 Power series methods
1: Let us start with a very simple example y00 − y = 0. Let us try a power series solution near x0 = 0, which is an ordinary point. Every point is an ordinary.
[16]
Power Series, Part 2: Analytic Solutions at Ordinary Points of ODEs
Consequently we cannot find the solutions in the usual way (via the characteristic polynomial). Instead, we use the Power Series Method.
[17]
[PDF] Chapter 2 Local analysis I: Linear di↵erential equations
Thus, every point x , 0 is an ordinary point of the first equation, the point x = 0 is a regular singular point of the second equation and an irregular.Missing: theorem | Show results with:theorem
[18]
[PDF] Regular points and singular points of second-order linear differential ...
In the case of a regular singular point at the origin, one can work out a systematic method for obtaining the corresponding Frobenius series solutions for y1(x) ...<|control11|><|separator|>
[19]
https://math.libretexts.org/Bookshelves/Differential_Equations/A_First_Course_in_Differential_Equations_for_Scientists_and_Engineers_(Herman)/04:_Series_Solutions/4.04:_The_Frobenius_Method/4.4.02:_Roots_of_Indicial_Equation
[20]
Ueber die Integration der linearen Differentialgleichungen durch ...
Frobenius, G.. "Ueber die Integration der linearen Differentialgleichungen durch Reihen.." Journal für die reine und angewandte Mathematik 76 (1873): 214-235.Missing: nach | Show results with:nach
[21]
Indicial Equation -- from Wolfram MathWorld
An indicial equation, also called a characteristic equation, is a recurrence equation obtained during application of the Frobenius method of solving a second- ...
[22]
https://math.libretexts.org/Courses/Mount_Royal_University/Mathematical_Methods/3%3A_Series_Solutions_of_Linear_Second_order_Equations/3.5%3A_The_Method_of_Frobenius_I
[23]
Frobenius Method for Solving Second-Order Ordinary Differential ...
In this paper we have studied how to solve second-order ordinary differential equation with variable coefficient at a singular point t = 0
[24]
[PDF] Notes on Diffy Qs: Differential Equations for Engineers
Page 1. Notes on Diffy Qs. Differential Equations for Engineers by Jiří Lebl. August 6, 2025. (version 6.9). Page 2. 2. Typeset in LATEX.
[25]
[PDF] Power Series and Differential Equations: The Method of Frobenius∗
Jan 29, 2014 · It is a procedure for finding power series formulas for solutions to a great many differential equations of interest in mathematical physics.
[26]
[PDF] The Method of Frobenius - Trinity University
Suppose the roots of the indicial equation are r1 and r2. If r1 − r2 6∈ Z, then both r = r1 and r = r2 yield (linearly independent) solutions. If r1 − r2 ∈ Z, ...
[27]
[PDF] Lecture 5: Examples of Frobenius Series: Bessel's Equation and ...
In this lecture we will consider the Frobenius series solution of the Bessel equation, which arises during the process.
[28]
[PDF] The method of Frobenius, MATH 3410 - Differential equations for ...
Mar 19, 2018 · Bessel functions are the canonical solutions of Bessel's differential equation: x. 2 y. 00. +xy. 0. + x. 2. −ν. 2 y = 0. (1) for an arbitrary ...
[29]
[PDF] Math 2280 - Lecture 35
The difference r1 −r2, with r1 ≥ r2, is an integer if r1 = r2, or if r1 = r2 +N. In the latter case there might, or might not, be two Frobenius solutions.
[30]
[PDF] 1 Introduction - Trinity University
The method is named after its inventor, Ferdinand George Frobenius (1849–1917), ... r2 is not an integer then the method of Frobenius yields solutions of the form.Missing: Georg original
[31]
[PDF] 642:527 fall 2010 summary of the method of frobenius
Radius of convergence. All the power series in (2)–(4) are guaranteed to have radius of convergence at least as big as the smaller of R1 and R2. 1. Page 2 ...
[32]
[PDF] The Big Theorem on the Frobenius Method, With Applications
The indicial equation arising in the method of Frobenius to solve this differential equation is α0r(r − 1) + β0r + γ0 = 0. Page 2. Chapter & Page: 34–2.
[33]
[PDF] Algorithms for Simplifying Differential Equations - FSU Math
Finding solutions of differential equations in terms of solutions of lower order differential equations is ... Lemma 3.26 (Fuchs's Relation). Let L ∈ C(x)[∂] of ...
[34]
[PDF] Frobenius Series Solutions - Math
1 1 .3 Frobenius Series Solutions. 663 and. = + r)(n + r —. I )C,Xhl+r 2 fl in Eq. (20)—the original differential equation. rather than Eq. (21 )—we get. 2 (n + ...
[35]
[PDF] Lecture Eight Irregular Singular Points of Ordinary Differential ...
Due to the factor of the exponential function, a solution near an irregular singular point behaves very differently from that near a regular singular point. It ...
[36]
[PDF] Symbolic Methods for Solving Systems of Linear Ordinary ...
Jul 25, 2010 · The integer p(A) := p is called the Poincaré rank of the system [A]. Def: - If p(A) = 0, the point x0 = 0 is called a singularity of first ...
[37]
[PDF] Lectures Nine and Ten The WKB Approximation
We will now apply the WKB method to obtain the asymptotic solutions near an irregular singular point of a second-order linear homogeneous equation. While we ...