The variation of the constants is a method from the theory of linear ordinary differential equations for the determination of a special solution of an inhomogeneous linear differential equation system of the first order or an inhomogeneous linear differential equation of any order. A complete solution ( fundamental system ) of the corresponding homogeneous differential equation is required for this .
Leonhard Euler used a forerunner of this method as early as 1748 in connection with astronomical problems. In its current form, the method was developed by the mathematician Joseph-Louis Lagrange .
motivation
First order linear differential equation
Let and be continuous functions, then the linear differential equation of the first order reads
![{\ displaystyle a: \ mathbb {R} \ rightarrow \ mathbb {R}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b43f9510a452746b2a5eb2348a0eec181681ec5a)
![{\ displaystyle b: \ mathbb {R} \ rightarrow \ mathbb {R}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fddec4013249dd608899f9f0fb746ab799404000)
![{\ displaystyle y '(x) = a (x) y (x) + b (x) \.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/25efc3b5fb12d7ee0ec565406183c96125df57dd)
Let also be an antiderivative of , then
![A.](https://wikimedia.org/api/rest_v1/media/math/render/svg/7daff47fa58cdfd29dc333def748ff5fa4c923e3)
![a](https://wikimedia.org/api/rest_v1/media/math/render/svg/ffd2487510aa438433a2579450ab2b3d557e5edc)
![{\ displaystyle A (x): = \ int _ {x_ {0}} ^ {x} a (t) {\ rm {d}} t \,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/437c803d763a286724a28a1703bb2885d2aae7ea)
where suitable boundary conditions must be satisfied. Then
![x_ {0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/86f21d0e31751534cd6584264ecf864a6aa792cf)
![{\ displaystyle \ left \ {y (x) = ce ^ {A (x)} \ | \ c \ in \ mathbb {R} \ right \}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/80ec1ec53d15e88e24e4c106c5ae72b3fa95468d)
the set of all solutions to the homogeneous differential equation .
![y '(x) = a (x) y (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/7428296a80bb7ff904bee1083da5271a8df3368b)
To solve the inhomogeneous differential equation, the function is now introduced and the approach of varying the constants selected
![{\ displaystyle c (x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b20e1deb5bbfe90c658811e9635d865f02902791)
-
.
This results in a clear assignment between the functions and , because it is an always positive, continuously differentiable function. The derivation of this approach function is
![y](https://wikimedia.org/api/rest_v1/media/math/render/svg/b8a6208ec717213d4317e666f1ae872e00620a0d)
![c](https://wikimedia.org/api/rest_v1/media/math/render/svg/86a67b81c2de995bd608d5b2df50cd8cd7d92455)
![{\ displaystyle \ exp (A (x))}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6fc3b1e6b4e8213c190626c017a141917dd16085)
![{\ displaystyle y '(x) = c (x) a (x) e ^ {A (x)} + c' (x) e ^ {A (x)} = a (x) y (x) + c '(x) e ^ {A (x)} \.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/32581122701c5341591b251568c28414c184787d)
So solves the inhomogeneous differential equation
![y](https://wikimedia.org/api/rest_v1/media/math/render/svg/b8a6208ec717213d4317e666f1ae872e00620a0d)
![{\ displaystyle y '(x) = a (x) y (x) + b (x) \}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a552f507575c1b91ba0a0698513610973fcc2d0b)
exactly when
![{\ displaystyle \ c '(x) = b (x) e ^ {- A (x)}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/44a80cb09a774f6b5ca9cc0ae3a53e999cf4d982)
applies. For example is
![{\ displaystyle c (x): = \ int _ {x_ {0}} ^ {x} b (t) e ^ {- A (t)} {\ rm {d}} t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bc02306436d84c4072fe991d6fe8fe9c1ee05447)
such a function and thus
![{\ displaystyle y_ {sp} (x): = e ^ {A (x)} \ cdot \ int _ {x_ {0}} ^ {x} b (t) e ^ {- A (t)} {\ rm {d}} t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/94fad3d6d72581a083f92d340f60c055d1efcc75)
the special solution with . So is
![y _ {{sp}} (x_ {0}) = 0](https://wikimedia.org/api/rest_v1/media/math/render/svg/563b0c3bb0dd6e82e2dbfce006cfc7f2a3bc2684)
![{\ displaystyle \ left \ {y (x) = e ^ {A (x)} \ cdot \ left [\ int _ {x_ {0}} ^ {x} b (t) e ^ {- A (t) } {\ rm {d}} t + c \ right] \ | \ c \ in \ mathbb {R} \ right \}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1733619e01924dbba9a4beee91ea3f3e7e0e58ed)
the set of all solutions to the inhomogeneous differential equation .
![{\ displaystyle y '(x) = a (x) y (x) + b (x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ad73154ce77549f6bbf1cf5f24a3b8a7a0f38cca)
example
If a DC voltage is applied to a coil with the inductance and the ohmic resistor , then the following applies to the voltage at the resistor
![R.](https://wikimedia.org/api/rest_v1/media/math/render/svg/4b0bfb3769bf24d80e15374dc37b0441e2616e33)
![U_ {0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/52116e300c6199496d7b4a60d417ba34a4e569dc)
![U (t) = U_ {0} -L {\ dot I} (t). \! \,](https://wikimedia.org/api/rest_v1/media/math/render/svg/9340122cae0d38c2cadfae7cca91ee3d0f463d52)
According to Ohm's law , this also applies
-
.
It is therefore an inhomogeneous linear differential equation of the first order with constant coefficients, which is now to be solved using the method of varying the constants.
For the corresponding homogeneous differential equation
![{\ displaystyle {\ dot {I}} _ {h} (t) = - {\ frac {R} {L}} I_ {h} (t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f67ed79daeaf861f229b8ab57a9364fa5e65398a)
is the solution for each constant
-
.
As a starting point for the solution of the inhomogeneous differential equation, replace the constant with a variable expression . So you bet
![c](https://wikimedia.org/api/rest_v1/media/math/render/svg/86a67b81c2de995bd608d5b2df50cd8cd7d92455)
![c (t)](https://wikimedia.org/api/rest_v1/media/math/render/svg/6e77d04462d26c4fdbfe8f988b182babd15059df)
![I (t): = c (t) e ^ {{- {\ frac {R} {L}} t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/be8736e0ae299201e42f5e676b1bcfc335e581fa)
and tries to determine a differentiable function that satisfies the inhomogeneous differential equation. It follows
![c (t)](https://wikimedia.org/api/rest_v1/media/math/render/svg/6e77d04462d26c4fdbfe8f988b182babd15059df)
![I.](https://wikimedia.org/api/rest_v1/media/math/render/svg/535ea7fc4134a31cbe2251d9d3511374bc41be9f)
![{\ displaystyle {\ begin {aligned} {\ frac {U_ {0}} {R}} - {\ frac {L} {R}} {\ dot {I}} (t) & = {\ frac {U_ {0}} {R}} - {\ frac {L} {R}} {\ dot {c}} (t) e ^ {- {\ frac {R} {L}} t} + {\ frac { L} {R}} c (t) {\ frac {R} {L}} e ^ {- {\ frac {R} {L}} t} \\ & = {\ frac {U_ {0}} { R}} - {\ frac {L} {R}} {\ dot {c}} (t) e ^ {- {\ frac {R} {L}} t} + c (t) e ^ {- { \ frac {R} {L}} t} \\ & = {\ frac {U_ {0}} {R}} - {\ frac {L} {R}} {\ dot {c}} (t) e ^ {- {\ frac {R} {L}} t} + I (t). \ end {aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d7511fbed99ec4b48dde3d33329a7c3e96aea95a)
Accordingly, the inhomogeneous differential equation is solved if and only if applies
-
.
This boundary value condition is equivalent to or after integration with . So the solution to the inhomogeneous differential equation is
![{\ displaystyle \ textstyle {\ dot {c}} (t) = {\ frac {U_ {0}} {L}} e ^ {{\ frac {R} {L}} t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2bfde6471f8381cda8bb2281bf29fbd52ac08ce0)
![\ textstyle c (t) = {\ frac {U_ {0}} {R}} e ^ {{{\ frac {R} {L}} t}} + d](https://wikimedia.org/api/rest_v1/media/math/render/svg/35ea0464f11f76f4290602b9eb05862c47ae5ef0)
-
.
The constant can be determined from the initial condition and results in the solution
![d](https://wikimedia.org/api/rest_v1/media/math/render/svg/e85ff03cbe0c7341af6b982e47e9f90d235c66ab)
![I (0) = 0](https://wikimedia.org/api/rest_v1/media/math/render/svg/fb32036dba98d7bf35253b755db5efbd89be0e29)
-
.
Inhomogeneous systems of linear differential equations of the first order
The above procedure can be generalized in the following way:
formulation
Let and be continuous functions and a fundamental matrix of the homogeneous problem as well as the matrix that results from by replacing the -th column with . Then
![A: {\ mathbb {R}} \ rightarrow {\ mathbb {R}} ^ {{n \ times n}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ed21463d3d239641f60814fbf4d9be06ae9bf0f)
![b: {\ mathbb {R}} \ rightarrow {\ mathbb {R}} ^ {n}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2c79169e826dfbfab9d8ab194c5eaf1cfdf62ead)
![\ Phi (x) = (y_ {1} (x) \ | \ \ cdots \ | \ y_ {n} (x))](https://wikimedia.org/api/rest_v1/media/math/render/svg/ea494ce84b1543f8a90dc54628e9f8ddd6c125c8)
![y '(x) = A (x) y (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/bfd3381168112b842fdd36a5d03427cd67dcc55c)
![\ Phi _ {k} (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/b378cdb87f9dac4a916d1da87800f3cba8df876f)
![\ Phi (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/79e4f01c93494fbb5dcd75761f4468121b00b294)
![k](https://wikimedia.org/api/rest_v1/media/math/render/svg/c3c9a2c7b599b37105512c5d570edc034056dd40)
![b (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/97e149d1ad9ac13aa08fe1230a622e06d2dde3a0)
![y _ {{sp}} (x): = \ sum _ {{k = 1}} ^ {n} c_ {k} (x) y_ {k} (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/5201495d03a5a3d04e114c6c42d1fdcc6f0e44a1)
With
![c_ {k} (x): = \ int _ {{x_ {0}}} ^ {x} {\ frac {\ det \ Phi _ {k} (s)} {\ det \ Phi (s)}} {{\ rm {d}}} s](https://wikimedia.org/api/rest_v1/media/math/render/svg/b297dbf69802a11554dccf4f5cc2b3f467c004ba)
the solution of the inhomogeneous initial value problem and .
![y '(x) = A (x) y (x) + b (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/7325840c8533b903cee0fe0a9ad0c3762c438cf8)
![y (x_ {0}) = 0](https://wikimedia.org/api/rest_v1/media/math/render/svg/18976558c36a4fe8bab5d480c92393d73009b95e)
proof
Set
![y _ {{sp}} (x): = \ Phi (x) \ int _ {{x_ {0}}} ^ {x} \ Phi (s) ^ {{- 1}} b (s) {{\ rm {d}}} s \.](https://wikimedia.org/api/rest_v1/media/math/render/svg/90cb03e27ceb5bd3b31c37bea303773f2d818b79)
It is , and because of differentiation, one sees that the differential equation is satisfied. Well solves
![y _ {{sp}} (x_ {0}) = 0](https://wikimedia.org/api/rest_v1/media/math/render/svg/563b0c3bb0dd6e82e2dbfce006cfc7f2a3bc2684)
![\ Phi '(x) = A (x) \ Phi (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/618c8f9138421a87577addd6635ed10513422b3d)
![y _ {{sp}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8fba466b202f2b2b71fa5f67b1330233f9a3aef2)
![y _ {{sp}} '(x) = A (x) y _ {{sp}} (x) + b (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/42cfc91ce332bffe799449fc9e800913d3c11727)
![a (s): = \ Phi ^ {{- 1}} (s) b (s) \ in {\ mathbb {R}} ^ {n}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2799d5c9646ebeffadb84928b0cd42467d33d556)
for fixed the system of linear equations
![s](https://wikimedia.org/api/rest_v1/media/math/render/svg/01d131dfd7673938b947072a13a9744fe997e632)
![\ Phi (s) \ cdot a (s) = b (s) \.](https://wikimedia.org/api/rest_v1/media/math/render/svg/a0c352572e3bf8917a3cf14b13cea4bd0243d6ea)
According to Cramer's rule is
![a_ {k} (s) = {\ frac {\ det \ Phi _ {k} (s)} {\ det \ Phi (s)}} \, \ k = 1, \ ldots, n \.](https://wikimedia.org/api/rest_v1/media/math/render/svg/6ce2154305feef0bc042d5d4816cf02ec1f832fb)
So it applies
![y _ {{sp}} (x) = \ int _ {{x_ {0}}} ^ {x} \ Phi (x) a (s) {{\ rm {d}}} s = \ sum _ {{ k = 1}} ^ {n} \ left [\ int _ {{x_ {0}}} ^ {x} {\ frac {\ det \ Phi _ {k} (s)} {\ det \ Phi (s )}} {{\ rm {d}}} s \ right] y_ {k} (x) \.](https://wikimedia.org/api/rest_v1/media/math/render/svg/80f4cb24c1653159831404b8e7f637f987f46efd)
Special case: resonance case
If the inhomogeneity itself is the solution of the homogeneous problem, ie , this is called a resonance case . In this case it is
![b](https://wikimedia.org/api/rest_v1/media/math/render/svg/f11423fbb2e967f986e36804a8ae4271734917c3)
![b '(x) = A (x) b (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/650ae7c81c2b3c015c4698cfa72ae3d9f7a98b40)
![\ y _ {{sp}} (x): = (x-x_ {0}) b (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/9bfe671e7526248d0a406bc93437a2c5afe12e19)
the solution of the inhomogeneous initial value problem and .
![y '(x) = A (x) y (x) + b (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/7325840c8533b903cee0fe0a9ad0c3762c438cf8)
![y (x_ {0}) = 0](https://wikimedia.org/api/rest_v1/media/math/render/svg/18976558c36a4fe8bab5d480c92393d73009b95e)
Inhomogeneous linear differential equations of higher order
Solving a higher order differential equation is equivalent to solving an appropriate first order differential equation system. In this way the above procedure can be used to construct a special solution for a higher order differential equation.
formulation
Let be continuous functions and a fundamental matrix of the homogeneous problem whose first row reads, as well as the matrix that results from by replacing the -th column with . Then
![a_ {0}, \ ldots, a _ {{n-1}}, b: {\ mathbb {R}} \ rightarrow {\ mathbb {R}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6db3e73274777850eb71077077b142561e1b1090)
![\ Phi (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/79e4f01c93494fbb5dcd75761f4468121b00b294)
![\ textstyle y ^ {{(n)}} (x) = \ sum _ {{k = 0}} ^ {{n-1}} a_ {k} (x) y ^ {{(k)}} ( x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/b4ed74dbd1140d6db7ecd26661d54b2ec69c0c4d)
![(y_ {1} (x) \ | \ \ cdots \ | \ y_ {n} (x))](https://wikimedia.org/api/rest_v1/media/math/render/svg/d30b409f0bf372b8715f3729fd308f46f7ad71f4)
![\ Phi _ {k} (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/b378cdb87f9dac4a916d1da87800f3cba8df876f)
![\ Phi (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/79e4f01c93494fbb5dcd75761f4468121b00b294)
![k](https://wikimedia.org/api/rest_v1/media/math/render/svg/c3c9a2c7b599b37105512c5d570edc034056dd40)
![\ textstyle {\ begin {pmatrix} 0 \\\ vdots \\ 0 \\ b (x) \\\ end {pmatrix}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f728ba6e06cac203327a22ce5b0ceff462bd8bc3)
![y _ {{sp}} (x): = \ sum _ {{k = 1}} ^ {n} c_ {k} (x) y_ {k} (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/5201495d03a5a3d04e114c6c42d1fdcc6f0e44a1)
With
![c_ {k} (x): = \ int _ {{x_ {0}}} ^ {x} {\ frac {\ det \ Phi _ {k} (s)} {\ det \ Phi (s)}} {{\ rm {d}}} s](https://wikimedia.org/api/rest_v1/media/math/render/svg/b297dbf69802a11554dccf4f5cc2b3f467c004ba)
the solution of the inhomogeneous initial value problem and .
![\ textstyle y ^ {{(n)}} (x) = \ sum _ {{k = 0}} ^ {{n-1}} a_ {k} (x) y ^ {{(k)}} ( x) + b (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/64b9e7a7af62548c7bfee3697d822aef7b823c24)
![y (x_ {0}) = 0](https://wikimedia.org/api/rest_v1/media/math/render/svg/18976558c36a4fe8bab5d480c92393d73009b95e)
proof
First consider the corresponding system of differential equations of the first order, consisting of equations
![n](https://wikimedia.org/api/rest_v1/media/math/render/svg/a601995d55609f2d9f5e233e36fbe9ea26011b3b)
-
With
The following applies: solves the scalar equation -th order if and only if the
solution of the above system is first order. By definition,
a fundamental matrix for this system is first order. Finally, we apply the method of varying the constants, which has been proven above.
![y (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/e871993bfd131a8b0c3591c26084cf8171a74dcd)
![n](https://wikimedia.org/api/rest_v1/media/math/render/svg/a601995d55609f2d9f5e233e36fbe9ea26011b3b)
![\ textstyle Y (x): = {\ begin {pmatrix} y (x) \\ y '(x) \\\ vdots \\ y ^ {{{(n-1)}} (x) \\\ end { pmatrix}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1115e845ffea71fe1b5ef55cce14ad666af62927)
![\ Phi](https://wikimedia.org/api/rest_v1/media/math/render/svg/aed80a2011a3912b028ba32a52dfa57165455f24)
Alternative: basic solution procedure
In the case of constant coefficients it is sometimes advantageous to use the basic solution method to construct a special solution: Is that homogeneous solution of which
![y_ {h}](https://wikimedia.org/api/rest_v1/media/math/render/svg/05e1b7d6e3a1b3e3261f7c1446b1269574ab45b5)
![\ textstyle y ^ {{(n)}} (x) = \ sum _ {{k = 0}} ^ {{n-1}} a_ {k} y ^ {{(k)}} (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/a0cc86d3dfceccbcd207a79d29719cadcd68c849)
![y_ {h} ^ {{(k)}} (x_ {0}) = 0 \, \ k = 0, \ ldots, n-2 \, \ y_ {h} ^ {{(n-1)}} (x_ {0}) = 1](https://wikimedia.org/api/rest_v1/media/math/render/svg/10671e8c5175ee7eecbfadd45c2c57aa5894a4ab)
met then is
![y _ {{sp}} (x): = \ int _ {{x_ {0}}} ^ {x} y_ {h} (x_ {0} + xt) b (t) {{\ rm {d}} } t](https://wikimedia.org/api/rest_v1/media/math/render/svg/93f54aa6bce2c235b4a367c9901acc1ec2f66323)
that special solution from with .
![\ textstyle y ^ {{(n)}} (x) = \ sum _ {{k = 0}} ^ {{n-1}} a_ {k} y ^ {{(k)}} (x) + b (x)](https://wikimedia.org/api/rest_v1/media/math/render/svg/f621d45cfe6fbf5016582008cb324ae97d16c8c4)
![y _ {{sp}} (x_ {0}) = 0](https://wikimedia.org/api/rest_v1/media/math/render/svg/563b0c3bb0dd6e82e2dbfce006cfc7f2a3bc2684)
proof
You check by differentiating
![y _ {{sp}} ^ {{(k)}} (x) = \ int _ {{x_ {0}}} ^ {x} y_ {h} ^ {{(k)}} (x_ {0} + xt) b (t) {{\ rm {d}}} t \, \ k = 0, \ ldots, n-1](https://wikimedia.org/api/rest_v1/media/math/render/svg/5c9a8b7ec11a0b96b8dbf23401e499d2aa5f4fd9)
and
![y _ {{sp}} ^ {{(n)}} (x) = b (x) + \ int _ {{x_ {0}}} ^ {x} y_ {h} ^ {{(n)}} (x_ {0} + xt) b (t) {{\ rm {d}}} t \.](https://wikimedia.org/api/rest_v1/media/math/render/svg/bee909cd696ebe6371d83b2696d5a25afdbe7bd2)
It turns out
![y _ {{sp}} ^ {{(n)}} (x) - \ sum _ {{k = 0}} ^ {{n-1}} a_ {k} y _ {{sp}} ^ {{( k)}} (x) = b (x) + \ int _ {{x_ {0}}} ^ {x} \ left [y_ {h} ^ {{(n)}} - \ sum _ {{k = 0}} ^ {{n-1}} a_ {k} y_ {h} ^ {{(k)}} \ right] (x_ {0} + xt) b (t) {{\ rm {d} }} t = b (x) \.](https://wikimedia.org/api/rest_v1/media/math/render/svg/efda58292c5fc8cf79893c3c5fe0edb8bff9239d)
Individual evidence
-
^ Forest Ray Moulton : An Introduction to Celestial Mechanics , 2nd ed. (First published by the Macmillan Company in 1914; reprinted in 1970 by Dover Publications, Inc., Mineola, New York), page 431
-
^ Leonhard Euler: Recherches sur la question des inégalités du mouvement de Saturne et de Jupiter, sujet proposé pour le prix de l'année 1748 par l'Académie Royale des Sciences de Paris, France: G. Martin, JB Coignard, & HL Guerin , 1749, online at: Google.com
-
^ Joseph-Louis Lagrange: (1766) “Solution de différensproblemèmes du calcul integral,” Mélanges de philosophie et de mathématique de la Société royale de Turin , vol. 3, pages 179-380.
-
↑ Wolfgang Walter : Ordinary differential equations. 3. Edition. Springer Verlag, 1986, ISBN 3-540-16143-0 , §2, Section II
-
↑ Wolfgang Walter : Ordinary differential equations. 3. Edition. Springer Verlag, 1986, ISBN 3-540-16143-0 , §16
-
↑ Differential equations of the nth order. In: Otto Forster : Analysis II. Vieweg Verlag, 1977, ISBN 3-499-27031-5 , Chapter II, §12.