|
|
Line 1: |
Line 1: |
| In [[mathematics]], a '''Cauchy–Euler equation''' (also known as the '''Euler–Cauchy equation''', or simply '''Euler's equation''') is a [[linear differential equation|linear]] [[homogeneous differential equation|homogeneous]] [[ordinary differential equation]] with [[variable coefficient]]s. It is sometimes referred to as an equidimensional equation. Because of its simple structure the equation can be replaced with an equivalent equation with [[constant coefficient]]s which can then be solved explicitly.
| | Got nothing to write about myself really.<br>Hurrey Im here and a part of this community.<br>I really hope I'm useful in some way .<br><br>Here is my site; [http://www.exit44.com/2014/08/dog-training-train-pit-bull/ best dog training methods] |
| | |
| ==The equation==
| |
| Let ''y''<sup>(''n'')</sup>(''x'') be the ''n''th derivative of the unknown function ''y''(''x''). Then a '''Cauchy–Euler equation''' of order ''n'' has the form
| |
| | |
| :<math>a_{n} x^n y^{(n)}(x) + a_{n-1} x^{n-1} y^{(n-1)}(x) + \cdots + a_0 y(x) = 0.</math>
| |
| | |
| The substitution <math>x = e^u</math> reduces this equation to a linear differential equation with constant coefficients. Alternatively a trial solution <math>y = x^m</math> may be used to solve for the basis solutions.<ref name="kreyszig">{{cite book|last=Kreyszig|first=Erwin|title=Advanced Engineering Mathematics|publisher=Wiley|date=May 10, 2006|isbn=978-0-470-08484-7}}</ref>
| |
| | |
| ===Second order – solving through trial solution===
| |
| [[File:Euler-Cauchy equation solution curves real roots.png|thumb|right|Typical solution curves for a second-order Euler–Cauchy equation for the case of two real roots]]
| |
| [[File:Euler-Cauchy equation solution curves double root.png|thumb|right|Typical solution curves for a second-order Euler–Cauchy equation for the case of a double root]]
| |
| [[File:Euler-Cauchy equation solution curves complex roots.png|thumb|right|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:<ref name="kreyszig" />
| |
| | |
| :<math>x^2\frac{d^2y}{dx^2} + ax\frac{dy}{dx} + by = 0. \,</math>
| |
| | |
| We assume a trial solution given by<ref name="kreyszig" />
| |
| | |
| :<math>y = x^m. \,</math>
| |
| | |
| Differentiating, we have:
| |
| | |
| :<math>\frac{dy}{dx} = mx^{m-1} \,</math>
| |
| | |
| and
| |
| | |
| :<math>\frac{d^2y}{dx^2} = m(m-1)x^{m-2}. \,</math>
| |
| | |
| Substituting into the original equation, we have:
| |
| | |
| :<math>x^2( m(m-1)x^{m-2} ) + ax( mx^{m-1} ) + b( x^m ) = 0 \,</math>
| |
| | |
| Or rearranging gives:
| |
| | |
| :<math>m^2 + (a-1)m + b = 0. \,</math>
| |
| | |
| We then can solve for ''m''. There are three particular cases of interest:
| |
| | |
| * Case #1: Two distinct roots, ''m''<sub>1</sub> and ''m''<sub>2</sub>
| |
| * Case #2: One real repeated root, ''m''
| |
| * Case #3: Complex roots, α ± β''i''
| |
| | |
| In case #1, the solution is given by:
| |
| :<math>y = c_1 x^{m_{1}} + c_2 x^{m_2} \,</math>
| |
| | |
| In case #2, the solution is given by
| |
| | |
| :<math>y = c_1 x^m \ln(x) + c_2 x^m \,</math>
| |
| | |
| To get to this solution, the method of [[reduction of order]] must be applied after having found one solution ''y'' = ''x''<sup>''m''</sup>.
| |
| | |
| In case #3, the solution is given by:
| |
| | |
| :<math>y = c_1 x^\alpha \cos(\beta \ln(x)) + c_2 x^\alpha \sin(\beta \ln(x)) \,</math>
| |
| | |
| :<math>\alpha = \mathop{\rm Re}(m)\,</math>
| |
| :<math>\beta = \mathop{\rm Im}(m)\,</math>
| |
| | |
| For <math>c_1\,</math> and <math>c_2\,</math> in the real plane
| |
| | |
| This form of the solution is derived by setting ''x'' = ''e''<sup>''t''</sup> and using [[Euler's formula]]
| |
| | |
| ===Second order – solution through change of variables===
| |
| :<math>x^2\frac{d^2y}{dx^2} + ax\frac{dy}{dx} + by = 0 \,</math>
| |
| | |
| We operate the variable substitution defined by
| |
| | |
| :<math>t = \ln(x). \,</math>
| |
| :<math>y(x) = \phi(\ln(x)) = \phi(t). \,</math>
| |
| Differentiating:
| |
| :<math>\frac{dy}{dx}=\frac{1}{x}\frac{d\phi}{dt}</math>
| |
| | |
| :<math>\frac{d^2y}{dx^2}=\frac{1}{x^2}\bigg(\frac{d^2\phi}{dt^2}-\frac{d\phi}{dt}\bigg).</math>
| |
| | |
| Substituting <math>\phi(t)</math>, we have
| |
| :<math>\frac{d^2\phi}{dt^2} + (a-1)\frac{d\phi}{dt} + b\phi = 0. \,</math>
| |
| | |
| This equation in <math>\phi(t)</math> can be easily solved using its characteristic polynomial
| |
| :<math>\lambda^2 + (a-1)\lambda +b = 0.</math>
| |
| | |
| Now, if <math>\lambda_1</math> and <math>\lambda_2</math> 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
| |
| | |
| :<math>\phi(t)=c_1 e^{\lambda_1 t} + c_2 e^{\lambda_2 t}</math>, where the exponentials may be complex.
| |
| | |
| If the roots are equal, the general solution is given by
| |
| | |
| :<math>\phi(t)=c_1 e^{\lambda_1 t} + c_2 t e^{\lambda_1 t}.</math>
| |
| | |
| In both cases, the solution <math>y(x)</math> may be found by setting <math>t=\ln(x)</math>, hence <math>\phi(\ln(x)) = y(x)</math>.
| |
| | |
| Hence, in the first case,
| |
| :<math>y(x)=c_1 x^{\lambda_1} + c_2 x^{\lambda_2}</math>,
| |
| | |
| and in the second case,
| |
| :<math>y(x)=c_1 x^{\lambda_1} + c_2 \ln(x) x^{\lambda_1}.</math>
| |
| | |
| ===Example===
| |
| Given
| |
| | |
| : <math>x^2u''-3xu'+3u=0\,,</math>
| |
| | |
| we substitute the simple solution ''x''<sup>α</sup>:
| |
| | |
| : <math>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\,.</math>
| |
| | |
| For ''x''<sup>α</sup> to be a solution, either ''x'' = 0, which gives the [[Trivial (mathematics)|trivial]] solution, or the coefficient of ''x''<sup>α</sup> is zero. Solving the quadratic equation, we get ''α'' = 1, 3. The general solution is therefore
| |
| | |
| : <math>u=c_1 x+c_2 x^3\,.</math>
| |
| | |
| ==Difference equation analogue==
| |
| There is a [[difference equation]] analogue to the Cauchy–Euler equation. For a fixed ''m'' > 0, define the sequence ''ƒ''<sub>''m''</sub>(''n'') as
| |
| | |
| : <math>f_m(n) := n(n+1)\cdots (n+m-1).</math>
| |
| | |
| Applying the difference operator to <math>f_m</math>, we find that
| |
| | |
| : <math>
| |
| \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}
| |
| </math>
| |
| | |
| If we do this ''k'' times, we will find that
| |
| | |
| : <math>
| |
| \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}
| |
| </math>
| |
| | |
| where the superscript <sup>(''k'')</sup> denotes applying the difference operator ''k'' times. Comparing this to the fact that the ''k''-th derivative of ''x''<sup>''m''</sup> equals
| |
| | |
| : <math>m(m-1)\cdots(m-k+1)\frac{x^m}{x^k}</math>
| |
| | |
| suggests that we can solve the ''N''-th order difference equation
| |
| | |
| : <math>f_N(n) y^{(N)}(n) + a_{N-1} f_{N-1}(n) y^{(N-1)}(n) + \cdots + a_0 y(n) = 0,</math>
| |
| | |
| in a similar manner to the differential equation case. Indeed, substituting the trial solution
| |
| | |
| : <math>y(n) = f_m(n) \, </math>
| |
| | |
| brings us to the same situation as the differential equation case,
| |
| | |
| : <math>m(m-1)\cdots(m-N+1) + a_{N-1} m(m-1) \cdots (m-N+2) + \cdots + a_1 m + a_0 = 0.</math>
| |
| | |
| 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 ''m''<sub>1</sub> will yield expressions involving a discrete version of ln,
| |
| | |
| : <math>\varphi(n) = \sum_{k=1}^n \frac{1}{k - m_1}.</math>
| |
| | |
| (Compare with: <math>\ln (x - m_1) = \int_{1+m_1}^x \frac{1}{t - m_1} \, dt.</math>)
| |
| | |
| In cases where fractions become involved, one may use
| |
| | |
| : <math>f_m(n) := \frac{\Gamma(n+m)}{\Gamma(n)}</math>
| |
| | |
| instead (or simply use it in all cases), which coincides with the definition before for integer ''m''.
| |
| | |
| ==See also==
| |
| * [[Hypergeometric differential equation]]
| |
| * [[Cauchy–Euler operator]]
| |
| | |
| ==References==
| |
| <references />
| |
| | |
| ==Bibliography==
| |
| *{{Mathworld|EulerDifferentialEquation}}
| |
| | |
| {{DEFAULTSORT:Cauchy-Euler Equation}}
| |
| [[Category:Ordinary differential equations]]
| |