TOPICS
Search

Leading Order Analysis


A procedure for determining the behavior of an nth order ordinary differential equation at a removable singularity without actually solving the equation. Consider

 (d^ny)/(dz^n)=F((d^(n-1)y)/(dz^(n-1)),...,(dy)/(dz),y,z),
(1)

where F is analytic in z and rational in its other arguments. Proceed by making the substitution

 y(z)=a(z-z_0)^alpha
(2)

with alpha<1. For example, in the equation

 (d^2y)/(dz^2)=6y^2+Ay,
(3)

making the substitution gives

 aalpha(alpha-1)(z-z_0)^(alpha-2)=6a^2(z-z_0)^(2alpha)+Aa(az-z_0)^alpha.
(4)

The most singular terms (those with the most negative exponents) are called the "dominant balance terms," and must balance exponents and coefficients at the singularity. Here, the first two terms are dominant, so

 alpha-2=2alpha=>alpha=-2
(5)
 6a=6a^2=>a=1,
(6)

and the solution behaves as y(z)=(z-z_0)^(-2). The behavior in the neighborhood of the singularity is given by expansion in a Laurent series, in this case,

 y(z)=sum_(j=0)^inftya_j(z-z_0)^(j-2).
(7)

Plugging this series in yields

 sum_(j=0)^inftya_j(j-2)(j-3)(z-z_0)^(j-4) 
 =6sum_(j=0)^inftysum_(k=0)^inftya_ja_k(z-z_0)^(j+k-4)+Asum_(j=0)^inftya_j(z-z_0)^(j-2).
(8)

This gives recurrence relations, in this case with a_6 arbitrary, so the (z-z_0)^6 term is called the resonance or Kovalevskaya exponent. At the resonances, the coefficient will always be arbitrary. If no resonance term is present, the pole present is not ordinary, and the solution must be investigated using a psi function.


See also

Psi Function

Explore with Wolfram|Alpha

References

Tabor, M. Chaos and Integrability in Nonlinear Dynamics: An Introduction. New York: Wiley, p. 330, 1989.

Referenced on Wolfram|Alpha

Leading Order Analysis

Cite this as:

Weisstein, Eric W. "Leading Order Analysis." From MathWorld--A Wolfram Web Resource. https://mathworld.wolfram.com/LeadingOrderAnalysis.html

Subject classifications