Cauchy–Euler equation

In mathematics, a Cauchy–Euler equation (also known as the Euler–Cauchy equation, or simply Euler's equation) is a linear homogeneous ordinary differential equation with variable coefficients. It is sometimes referred to as an equidimensional equation. Because of the particularly simple equidimensional structure the equation can be replaced with an equivalent equation with constant coefficients which can then be solved explicitly.

The equation

Let y(n)(x) be the nth derivative of the unknown function y(x). Then a Cauchy–Euler equation of order n has the form

a_{n} x^n y^{(n)}(x) + a_{n-1} x^{n-1} y^{(n-1)}(x) + \cdots + a_0 y(x) = 0.

The substitution x = e^u reduces this equation to a linear differential equation with constant coefficients. Alternatively a trial solution y = x^m may be used to solve for the basis solutions.[1]

Second order – solving through trial solution

Typical solution curves for a second-order Euler–Cauchy equation for the case of two real roots
Typical solution curves for a second-order Euler–Cauchy equation for the case of a double root
Typical solution curves for a second-order Euler–Cauchy equation for the case of complex roots

The most common Cauchy–Euler equation is the second-order equation, appearing in a number of physics and engineering applications, such as when solving Laplace's equation in polar coordinates. It is given by the equation:[1]

x^2\frac{d^2y}{dx^2} + ax\frac{dy}{dx} + by = 0. \,

We assume a trial solution given by[1]

y = x^m. \,

Differentiating, we have:

\frac{dy}{dx} = mx^{m-1} \,

and

\frac{d^2y}{dx^2} = m(m-1)x^{m-2}. \,

Substituting into the original equation, we have:

x^2( m(m-1)x^{m-2} ) + ax( mx^{m-1} ) + b( x^m ) = 0 \,

Or rearranging gives:

m^2 + (a-1)m + b = 0. \,

We then can solve for m. There are three particular cases of interest:

In case #1, the solution is given by:

y = c_1 x^{m_{1}} + c_2 x^{m_2} \,

In case #2, the solution is given by

y = c_1 x^m \ln(x) + c_2 x^m \,

To get to this solution, the method of reduction of order must be applied after having found one solution y = xm.

In case #3, the solution is given by:

y = c_1 x^\alpha \cos(\beta \ln(x)) + c_2 x^\alpha \sin(\beta \ln(x)) \,
\alpha = \mathop{\rm Re}(m)\,
\beta = \mathop{\rm Im}(m)\,

For c_1\, and c_2\, in the real plane

This form of the solution is derived by setting x = et and using Euler's formula

Second order – solution through change of variables

x^2\frac{d^2y}{dx^2} + ax\frac{dy}{dx} + by = 0 \,

We operate the variable substitution defined by

t = \ln(x). \,
y(x) = \phi(\ln(x)) = \phi(t). \,

Differentiating:

\frac{dy}{dx}=\frac{1}{x}\frac{d\phi}{dt}
\frac{d^2y}{dx^2}=\frac{1}{x^2}\bigg(\frac{d^2\phi}{dt^2}-\frac{d\phi}{dt}\bigg).

Substituting \phi(t), we have

\frac{d^2\phi}{dt^2} + (a-1)\frac{d\phi}{dt} + b\phi = 0. \,

This equation in \phi(t) can be easily solved using its characteristic polynomial

\lambda^2 + (a-1)\lambda +b = 0.

Now, if \lambda_1 and \lambda_2 are the roots of this polynomial, we analyze the two main cases: distinct roots and double roots:

If the roots are distinct, the general solution is given by

\phi(t)=c_1 e^{\lambda_1 t} + c_2 e^{\lambda_2 t}, where the exponentials may be complex.

If the roots are equal, the general solution is given by

\phi(t)=c_1 e^{\lambda_1 t} + c_2 t e^{\lambda_1 t}.

In both cases, the solution y(x) may be found by setting t=\ln(x), hence \phi(\ln(x)) = y(x).

Hence, in the first case,

y(x)=c_1 x^{\lambda_1} + c_2 x^{\lambda_2},

and in the second case,

y(x)=c_1 x^{\lambda_1} + c_2 \ln(x) x^{\lambda_1}.

Example

Given

x^2u''-3xu'+3u=0\,,

we substitute the simple solution xα:

x^2(\alpha(\alpha-1)x^{\alpha-2})-3x(\alpha x^{\alpha-1})+3x^\alpha=\alpha(\alpha-1)x^\alpha-3\alpha x^\alpha+3x^\alpha = (\alpha^2-4\alpha+3)x^\alpha = 0\,.

For xα to be a solution, either x = 0, which gives the trivial solution, or the coefficient of xα is zero. Solving the quadratic equation, we get α = 1, 3. The general solution is therefore

u=c_1 x+c_2 x^3\,.

Difference equation analogue

There is a difference equation analogue to the CauchyEuler equation. For a fixed m > 0, define the sequence ƒm(n) as

f_m(n) := n(n+1)\cdots (n+m-1).

Applying the difference operator to f_m, we find that


\begin{align}
Df_m(n) & = f_{m}(n+1) - f_m(n) \\
& = m(n+1)(n+2) \cdots (n+m-1) = \frac{m}{n} f_m(n).
\end{align}

If we do this k times, we will find that


\begin{align}
f_m^{(k)}(n) & = \frac{m(m-1)\cdots(m-k+1)}{n(n+1)\cdots(n+k-1)} f_m(n) \\
& = m(m-1)\cdots(m-k+1) \frac{f_m(n)}{f_k(n)},
\end{align}

where the superscript (k) denotes applying the difference operator k times. Comparing this to the fact that the k-th derivative of xm equals

m(m-1)\cdots(m-k+1)\frac{x^m}{x^k}

suggests that we can solve the N-th order difference equation

f_N(n) y^{(N)}(n) + a_{N-1} f_{N-1}(n) y^{(N-1)}(n) + \cdots + a_0 y(n) = 0,

in a similar manner to the differential equation case. Indeed, substituting the trial solution

y(n) = f_m(n) \,

brings us to the same situation as the differential equation case,

m(m-1)\cdots(m-N+1) + a_{N-1} m(m-1) \cdots (m-N+2) + \cdots + a_1 m + a_0 = 0.

One may now proceed as in the differential equation case, since the general solution of an N-th order linear difference equation is also the linear combination of N linearly independent solutions. Applying reduction of order in case of a multiple root m1 will yield expressions involving a discrete version of ln,

\varphi(n) = \sum_{k=1}^n \frac{1}{k - m_1}.

(Compare with: \ln (x - m_1) = \int_{1+m_1}^x \frac{1}{t - m_1} \, dt.)

In cases where fractions become involved, one may use

f_m(n) := \frac{\Gamma(n+m)}{\Gamma(n)}

instead (or simply use it in all cases), which coincides with the definition before for integer m.

See also

References

  1. 1.0 1.1 1.2 Kreyszig, Erwin (May 10, 2006). Advanced Engineering Mathematics. Wiley. ISBN 978-0-470-08484-7.

Bibliography