# 1911 Encyclopædia Britannica/Differential Equation

See also Differential equation on Wikipedia; and our 1911 Encyclopædia Britannica disclaimer.
​DIFFERENTIAL EQUATION, in mathematics, a relation between one or more functions and their differential coefficients. The subject is treated here in two parts: (1) an elementary introduction dealing with the more commonly recognized types of differential equations which can be solved by rule; and (2) the general theory.
Part I.—Elementary Introduction.
Of equations involving only one independent variable, x (known as ordinary differential equations), and one dependent variable, y, and containing only the first differential coefficient dy/dx (and therefore said to be of the first order), the simplest form is that reducible to the type
dy/dx＝ƒ(x)/F(y),
leading to the result ƒF(y)dy − ƒf(x)dx＝A, where A is an arbitrary constant; this result is said to solve the differential equation, the problem of evaluating the integrals belonging to the integral calculus.
​Another simple form is
dy/dx + y P＝Q,
where P, Q are functions of x only; this is known as the linear equation, since it contains y and dy/dx only to the first degree. If ƒPdx = u, we clearly have
${\frac {d}{dx}}(ye^{u})=e^{u}\left({\frac {dy}{dx}}+\mathrm {P} y\right)=e^{u}\mathrm {Q} ,$
so that y = e$^{−u}$(ƒe$^u$Qdx + A) solves the equation, and is the only possible solution, A being an arbitrary constant. The rule for the solution of the linear equation is thus to multiply the equation by e$^u$, where u = ƒPdx.
A third simple and important form is that denoted by
y＝px + ƒ(p),
where p is an abbreviation for dy/dx; this is known as Clairaut’s form. By differentiation in regard to x it gives
$p=p+x{\frac {dp}{dx}}+f'(p){\frac {dp}{dx}},$
where
$f'(p)={\frac {d}{dp}}f(p);$
thus, either (i.) dp/dx＝ 0, that is, p is constant on the curve satisfying the differential equation, which curve is thus any one of the straight lines y = cx = ƒ(c), where c is an arbitrary constant, or else, (ii.) x + ƒ′(p) = 0; if this latter hypothesis be taken, and p be eliminated between x + ƒ′(p) = 0 and y = px + ƒ(p), a relation connecting x and y, not containing an arbitrary constant, will be found, which obviously represents the envelope of the straight lines y = cx + ƒ(c).
In general if a differential equation φ(x, y, dy/dx) = 0 be satisfied by any one of the curves F(x, y, c) = 0, where c is an arbitrary constant, it is clear that the envelope of these curves, when existent, must also satisfy the differential equation; for this equation prescribes a relation connecting only the co-ordinates x, y and the differential coefficient dy/dx, and these three quantities are the same at any point of the envelope for the envelope and for the particular curve of the family which there touches the envelope. The relation expressing the equation of the envelope is called a singular solution of the differential equation, meaning an isolated solution, as not being one of a family of curves depending upon an arbitrary parameter.
An extended form of Clairaut’s equation expressed by
y＝xF(p) + ƒ(p)
may be similarly solved by first differentiating in regard to p, when it reduces to a linear equation of which x is the dependent and p the independent variable; from the integral of this linear equation, and the original differential equation, the quantity p is then to be eliminated.
Other types of solvable differential equations of the first order are (1)
M dy/dx＝N,
where M, N are homogeneous polynomials in x and y, of the same order; by putting v = y/x and eliminating y, the equation becomes of the first type considered above, in v and x. An equation (aB≷bA)
(ax+by+c)dy/dx＝Ax+By+C
may be reduced to this rule by first putting x+h, y+k for x and y, and determining h, k so that ah+bk+c = 0, Ah+Bk+C = 0.
(2) An equation in which y does not explicitly occur,
ƒ(x, dy/dx)＝0,
may, theoretically, be reduced to the type dy/dx = F(x); similarly an equation F(y, dy/dx) = 0.
(3) An equation
ƒ(dy/dx, x, y)＝0,
which is an integral polynomial in dy/dx, may, theoretically, be solved for dy/dx, as an algebraic equation; to any root dy/dx = F$_1$(x, y) corresponds, suppose, a solution φ$_1$(x, y, c) = 0, where c is an arbitrary constant; the product equation φ$_1$(x, y, c)φ$_2$(x, y, c) . . . = 0, consisting of as many factors as there were values of dy/dx, is effectively as general as if we wrote φ$_1$(x, y, c$_1$)φ$_2$(x, y, c$_2$) . . . = 0; for, to evaluate the first form, we must necessarily consider the factors separately, and nothing is then gained by the multiple notation for the various arbitrary constants. The equation φ$_1$(x, y, c)φ$_2$(x, y, c) . . . = 0 is thus the solution of the given differential equation.
In all these cases there is, except for cases of singular solutions, one and only one arbitrary constant in the most general solution of the differential equation; that this must necessarily be so we may take as obvious, the differential equation being supposed to arise by elimination of this constant from the equation expressing its solution and the equation obtainable from this by differentiation in regard to x.
A further type of differential equation of the first order, of the form
dy/dx＝A + By + Cy $^2$
in which A, B, C are functions of x, will be briefly considered below under differential equations of the second order.
When we pass to ordinary differential equations of the second order, that is, those expressing a relation between x, y, dy/dx and d $^2$y/dx$^2$, the number of types for which the solution can be found by a known procedure is very considerably reduced. Consider the general linear equation
${\frac {d^{2}y}{dx^{2}}}+\mathrm {P} {\frac {dy}{dx}}+\mathrm {Q} =\mathrm {R} ,$
where P, Q, R are functions of x only. There is no method always effective; the main general result for such a linear equation is that if any particular function of x, say y $_1$, can be discovered, for which
${\frac {d^{2}y_{1}}{dx^{2}}}+\mathrm {P} {\frac {dy_{1}}{dx}}+\mathrm {Q} y_{1}=0,$
then the substitution y = y $_1$η in the original equation, with R on the right side, reduces this to a linear equation of the first order with the dependent variable dη/dx. In fact, if y = y $_1$η we have
${\frac {dy}{dx}}=y_{1}{\frac {d\eta }{dx}}+\eta {\frac {dy_{1}}{dx}}$ and ${\frac {d^{2}y}{dx^{2}}}=y_{1}{\frac {d^{2}\eta }{dx^{2}}}+2{\frac {dy_{1}}{dx}}{\frac {d\eta }{dx}}+\eta {\frac {d^{2}y_{1}}{dx^{2}}},$
and thus
${\frac {d^{2}y}{dx^{2}}}+\mathrm {P} {\frac {dy}{dx}}+\mathrm {Q} =y_{1}{\frac {d^{2}\eta }{dx^{2}}}+\left(2{\frac {dy_{1}}{dx}}+\mathrm {P} y_{1}\right){\frac {d\eta }{dx}}+\left({\frac {d^{2}y_{1}}{dx^{2}}}+\mathrm {P} {\frac {dy_{1}}{dx}}+\mathrm {Q} y_{1}\right)\eta ;$
if then
${\frac {d^{2}y_{1}}{dx^{2}}}+\mathrm {P} {\frac {dy_{1}}{dx}}+\mathrm {Q} y_{1}=0,$
and z denote dη/dx, the original differential equation becomes
$y_{1}{\frac {dz}{dx}}+\left(2{\frac {dy_{1}}{dx}}\mathrm {P} y_{1}\right)=\mathrm {R} .$
From this equation z can be found by the rule given above for the linear equation of the first order, and will involve one arbitrary constant; thence y = y $_1$ η = y $_1$ ∫ zdx + Ay $_1$, where A is another arbitrary constant, will be the general solution of the original equation, and, as was to be expected, involves two arbitrary constants.
The case of most frequent occurrence is that in which the coefficients P, Q are constants; we consider this case in some detail. If θ be a root of the quadratic equation θ$^2$ + θP + Q = 0, it can be at once seen that a particular integral of the differential equation with zero on the right side is y $_1$ = e$^{θx}$. Supposing first the roots of the quadratic equation to be different, and φ to be the other root, so that φ + θ = −P, the auxiliary differential equation for z, referred to above, becomes dz/dx + (θ − φ)z = Re$^{−θx}$ which leads to ze$^{(θ−φ)}$ = B + ∫ Re$^{−θx}$dx, where B is an arbitrary constant, and hence to
$y=Ae^{\theta x}+e^{\theta x}\int \mathrm {B} e^{(\phi -\theta )x}dx+e^{\theta x}\int e^{(\phi -\theta )x}\int Re^{-\phi x}dxdx,$
or say to $y=\mathrm {A} e^{\theta x}+\mathrm {C} e^{\phi x}+\mathrm {U}$, where A, C are arbitrary constants and U is a function of x, not present at all when R = 0. If the quadratic equation θ$^2$ + Pθ + Q = 0 has equal roots, so that 2θ = −P, the auxiliary equation in z becomes ${\frac {dz}{dx}}=Re^{\theta x}$ giving $z=B+\int \mathrm {R} e^{\theta x}dx$, where B is an arbitrary constant, and hence
$y=(\mathrm {A} +\mathrm {B} x)e^{\theta x}+e^{\theta x}\iint \mathrm {R} e^{-\theta x}dxdx,$
or, say, $y=(\mathrm {A} +\mathrm {B} x)e^{\theta x}+U$, where A, B are arbitrary constants, and U is a function of x not present at all when R = 0. The portion $\mathrm {A} e^{\theta x}+\mathrm {B} e^{\phi x}$ or $(\mathrm {A} +\mathrm {B} x)e^{\theta x}$ of the solution, which is known as the complementary function, can clearly be written down at once by inspection of the given differential equation. The remaining portion U may, by taking the constants in the complementary function properly, be replaced by any particular solution whatever of the differential equation
${\frac {d^{2}y}{dx^{2}}}+\mathrm {P} {\frac {dy}{dx}}+\mathrm {Q} =\mathrm {R} ;$
for if u be any particular solution, this has a form
$\displaystyle u=A_{0}e^{\theta x}+B_{0}e^{\phi x}+U,$
or a form
$\displaystyle u=(A_{0}+B_{0}x)e^{\theta x}+U;$
thus the general solution can be written
$\displaystyle (A-A_{0})e^{\theta x}+(B-B_{0})e^{\theta x}+u,$ or $\displaystyle \{A-A_{0}+(B-B_{0})x\}e^{\theta x}+u,$
where A − A$_0$, B − B$_0$, like A, B, are arbitrary constants.
A similar result holds for a linear differential equation of any order, say
${\frac {d^{n}y}{dx^{n}}}+\mathrm {P} _{1}{\frac {d^{n-1}y}{dx^{n-1}}}+\ \ldots \ +\mathrm {P} _{n}=\mathrm {R} ,$
where P$_1$, P$_2$, . . . P$_n$ are constants, and R is a function of x. If we form the algebraic equation θ$^n$ + P$_1$θ$^{n−1}$ + . . . + P$_n$ = 0, and all the roots of this equation be different, say they are θ$_1$, θ$_2$, . . . θ$_n$, the general solution of the differential equation is
$\displaystyle y=A_{1}e^{\theta _{1}x}+A_{2}e^{\theta _{2}x}+\ ...\ +A_{n}e^{\theta _{n}x}+u,$
where A$_1$, A$_2$, . . . A$_n$ are arbitrary constants, and u is any ​particular solution whatever; but if there be one root θ$_1$ repeated r times, the terms A$_1$ e$^{θ1x}$ + ... + A$_r$ e$^{θrx}$ must be replaced by (A$_1$ + A$_2$x + ... + A$_r$x$^{r−1}$)e$^{θ1x}$ where A$_1$, ... A$_n$ are arbitrary constants; the remaining terms in the complementary function will similarly need alteration of form if there be other repeated roots.
To complete the solution of the differential equation we need some method of determining a particular integral u; we explain a procedure which is effective for this purpose in the cases in which R is a sum of terms of the form e$^{ax}$φ(x), where φ(x) is an integral polynomial in x; this includes cases in which R contains terms of the form cos bx·φ(x) or sin bx·φ(x). Denote d/dx by D; it is clear that if u be any function of x, D(e$^{ax}$u) = e$^{ax}$Du + ae$^{ax}$u, or say, D(e$^{ax}$u) = e$^{ax}$(D + a)u; hence D$^2$(e$^{ax}$u), i.e. d $^2$/dx$^2$ (e$^{ax}$u), being equal to D(e$^{ax}$v), where v = (D + a)u, is equal to e$^{ax}$(D + a)v, that is to e$^{ax}$(D + a)$^2$u. In this way we find D$^n$(e$^{ax}$u) = e$^{ax}$(D + a)$^n$u, where n is any positive integer. Hence if ψ(D) be any polynomial in D with constant coefficients, ψ(D) (e$^{ax}$u) = e$^{ax}$ψ(D + a)u. Next, denoting ∫ udx by D$^{−1}$u, and any solution of the differential equation dz/dx + az = u by z = (d + a)$^{−1}$u, we have D[e$^{ax}$(D + a)$^{−1}$u] = D(e$^{ax}$z) = e$^{ax}$(D + a)z = e$^{ax}$u, so that we may write D$^{−1}$(e$^{ax}$u) = e$^{ax}$(D + a)$^{−1}$u, where the meaning is that one value of the left side is equal to one value of the right side; from this, the expression D$^{-2}$(e$^{ax}$u), which means D$^{−1}$[D$^{−1}$(e$^{ax}$u)], is equal to D$^{−1}$(e$^{ax}$z) and hence to e$^{ax}$(D + a)$^{−1}$z, which we write e$^{ax}$(D + a)$^{-2}$u; proceeding thus we obtain
D$^{-n}$(e$^{ax}$u)＝e$^{ax}$(D + a)$^{-n}$u,
where n is any positive integer, and the meaning, as before, is that one value of the first expression is equal to one value of the second. More generally, if ψ(D) be any polynomial in D with constant coefficients, and we agree to denote by [1/ψ(D)]u any solution z of the differential equation ψ(D)z = u, we have, if v = [1/ψ(D + a)]u, the identity ψ(D)(e$^{ax}$v) = e$^{ax}$ψ(D + a)v = e$^{ax}$u, which we write in the form
${\frac {1}{\psi (\mathrm {D} )}}(e^{ax}u)=e^{ax}{\frac {1}{\psi (\mathrm {D} +a)}}u.$
This gives us the first step in the method we are explaining, namely that a solution of the differential equation ψ(D)y = e$^{ax}$u + e$^{bx}$v + ... where u, v, ... are any functions of x, is any function denoted by the expression
$e^{ax}{\frac {1}{\psi (\mathrm {D} +a)}}u+e^{bx}{\frac {1}{\psi (\mathrm {D} +b)}}v+\ldots .$
It is now to be shown how to obtain one value of ${\frac {1}{\psi (\mathrm {D} +a)}}u$, when u is a polynomial in x, namely one solution of the differential equation ψ(D + a)z = u. Let the highest power of x entering in u be x$^m$; if t were a variable quantity, the rational fraction in t, ${\frac {1}{\psi (t+a)}}u$ by first writing it as a sum of partial fractions, or otherwise, could be identically written in the form
K$_r$t$^{-r}$ + K$_{r−1}$t$^{-r+1}$ + ... + K$_1$t$^{−1}$ + H + H$_1$t + ... + H$_m$t$^m$ + t$^{m+1}$φ(t)/ψ(t + a),
where φ(t) is a polynomial in t; this shows that there exists an identity of the form
1＝ψ(t + a)(K$_r$t$^{−r}$ + ... + K$_1$t$^{−1}$ + H + H$_1$t + ... + H$_m$t$^m$) + φ(t)t$^{m+1}$,
and hence an identity
u＝ψ(D + a) [K$_r$D$^{−r}$ + ... + K$_1$D$^{−1}$ + H + H$_1$D + ... + H$_m$D$^m$] u + φ(D) D$^{m+1}$u;
in this, since u contains no power of x higher than x$^m$, the second term on the right may be omitted. We thus reach the conclusion that a solution of the differential equation ψ(D + a)z = u is given by
z＝(K$_r$D$^{−r}$ + ... + K$_1$D$^{−1}$ + H + H$_1$D + ... + H$_m$D$^m$)u,
of which the operator on the right is obtained simply by expanding 1/ψ(D + a) in ascending powers of D, as if D were a numerical quantity, the expansion being carried as far as the highest power of D which, operating upon u, does not give zero. In this form every term in z is capable of immediate calculation.
Example.—For the equation
${\frac {d^{4}v}{dx^{4}}}+2{\frac {d^{2}y}{dx^{3}}}+y=x^{3}\cos x$ or $\displaystyle (\mathrm {D} ^{2}+1)^{2}=x^{3}\cos x,$
the roots of the associated algebraic equation (θ$^2$ + 1)$^2$ = 0 are θ = ±i, each repeated; the complementary function is thus
(A + Bx)e$^{ix}$ + (C + Dx)e$^{−ix}$,
where A, B, C, D are arbitrary constants; this is the same as
(H + Kx) cos x + (M + Nx) sin x,
where H, K, M, N are arbitrary constants. To obtain a particular integral we must find a value of (1 + D$^2$)−$^2$x$^3$ cos x; this is the real part of (1 + D$^2$)−$^2$ e$^{ix}$x$^3$ and hence of e$^{ix}$ [1 + (D + i)$^2$]−$^2$ x$^3$ or
e$^{ix}$ [2iD(1 + 1/2iD)]−$^2$ x$^3$,
or
−1/4e$^{ix}$ D−$^2$ (1 + iD − 3/4D$^2$ − 1/2iD$^3$ + 5/16D$^4$ + 3/16iD$^5$ ...)x$^3$,
or
−1/4e$^{ix}$ (1/20x$^5$ + 1/4ix$^4$ − 3/4x$^3$ − 3/2 ix$^2$ + 15/8 x + 9/8 i);
the real part of this is
−1/4 (1/20 x$^5$ − 3/4x$^2$ + 15/8x) cos x + 1/4 (1/4x$^4$ − 3/2x$^2$ + 9/8) sin x.
This expression added to the complementary function found above gives the complete integral; and no generality is lost by omitting from the particular integral the terms −15/32x cos x + 9/32 sin x, which are of the types of terms already occurring in the complementary function.
The symbolical method which has been explained has wider applications than that to which we have, for simplicity of explanation, restricted it. For example, if ψ(x) be any function of x, and a$_1$, a$_2$, ... a$_n$ be different constants, and [(t + a$_1$) (t + a$_2$) ... (t + a$_n$)]$^{−1}$ when expressed in partial fractions be written Σc$_m$(t + a$_m$)$^{−1}$, a particular integral of the differential equation (D + a$_1$)(D + a$_2$) ... (D + a$_n$)y = ψ(x) is given by
y＝Σc$_m$(D + a$_m$)$^{−1}$ ψ(x)＝Σc$_m$ (D + a$_m$)$^{−1}$ e$^{−a}$m$^x$e$^a$m$^x$ ψ(x)＝Σc$_m$e$^{−a}$m$^x$D$^{−1}$ (e$^a$m$^x$ψ(x) )＝Σc$_m$e$^{−a}$m$^x$ ∫ e$^a$m$^x$ψ(x)dx.
The particular integral is thus expressed as a sum of n integrals. A linear differential equation of which the left side has the form
$x^{n}{\frac {d^{n}y}{dx^{n}}}+\mathrm {P} _{1}x^{n-1}{\frac {d^{n-1}y}{dx^{n-1}}}+\ldots +\mathrm {P} _{n-1}x{\frac {dy}{dx}}+\mathrm {P} _{n}y,$
where P$_1$, ... P$_n$ are constants, can be reduced to the case considered above. Writing x = e$^t$ we have the identity
$x^{m}{\frac {d^{m}u}{dx^{m}}}=\theta (\theta -1)(\theta -2)\ldots (\theta -m+1)u,$ where θ＝d/dt
When the linear differential equation, which we take to be of the second order, has variable coefficients, though there is no general rule for obtaining a solution in finite terms, there are some results which it is of advantage to have in mind. We have seen that if one solution of the equation obtained by putting the right side zero, say y$_1$, be known, the equation can be solved. If y$_2$ be another solution of
${\frac {d^{2}y}{dx^{2}}}+\mathrm {P} {\frac {dy}{dx}}+\mathrm {Q} =0,$
there being no relation of the form my$_1$ + ny$_2$ = k, where m, n, k are constants, it is easy to see that
${\frac {d}{dx}}(y_{1}'y_{2}-y_{1}y_{2}')=\mathrm {P} (y_{1}'y_{2}-y_{1}y_{2}'),$
so that we have
y$_1$′y$_2$ − y$_1$y$_2$′＝A exp. (∫ Pdx),
where A is a suitably chosen constant, and exp. z denotes e$^z$. In terms of the two solutions y$_1$, y$_2$ of the differential equation having zero on the right side, the general solution of the equation with R = φ(x) on the right side can at once be verified to be Ay$_1$ + By$_2$ + y$_1$u − y$_2$v, where u, v respectively denote the integrals
u＝∫ y$_2$φ(x) (y$_1$′y$_2$ − y$_2$′y$_1$)$^{−1}$dx, v＝∫ y$_1$φ(x) (y$_1$′y$_2$ − y$_2$′y$_1$)$^{−1}$dx.
The equation
${\frac {d^{2}y}{dx^{2}}}+\mathrm {P} {\frac {dy}{dx}}+\mathrm {Q} =0,$
by writing y = v exp. (−1/2 ∫ Pdx), is at once seen to be reduced to d $^2$v/dx$^2$ + Iv = 0, where I = Q − 1/2dP/dx − 1/4P$^2$. If η = − 1/v dv/dx, the equation d $^2$v/dx$^2$ + Iv = 0 becomes dη/dx = I + η$^2$, a non-linear equation of the first order.
More generally the equation
${\frac {d\eta }{dx}}=\mathrm {A} +\mathrm {B} \eta +\mathrm {C} \eta ^{2}$
where A, B, C are functions of x, is, by the substitution
$\eta =-{\frac {1}{\mathrm {C} y}}{\frac {dy}{dx}}$
reduced to the linear equation
${\frac {d^{2}y}{dx^{2}}}-\left(\mathrm {B} +{\frac {1}{\mathrm {C} }}{\frac {d\mathrm {C} }{dx}}\right){\frac {dy}{dx}}+\mathrm {AC} =0.$
The equation
${\frac {d\eta }{dx}}=\mathrm {A} +\mathrm {B} \eta +\mathrm {C} \eta ^{2}$
known as Riccati’s equation, is transformed into an equation of the same form by a substitution of the form η = (aY + b)/(cY + d), where a, b, c, d are any functions of x, and this fact may be utilized to obtain a solution when A, B, C have special forms; in particular if any particular solution of the equation be known, say η$_0$, the ​substitution η = η$_0$ − 1/Y enables us at once to obtain the general solution; for instance, when
$2\mathrm {B} {=}{\frac {d}{dx}}\log \left({\frac {\mathrm {A} }{\mathrm {C} }}\right),$
a particular solution is η$_0$ = √(-A/C). This is a case of the remark, often useful in practice, that the linear equation
$\phi (x){\frac {d^{2}y}{dx^{2}}}+{\tfrac {1}{2}}{\frac {d\phi }{dx}}{\frac {dy}{dx}}+\mu y{=}0,$
where μ is a constant, is reducible to a standard form by taking a new independent variable z = ∫ dx[φ(x)]$^{-½}$.
We pass to other types of equations of which the solution can be obtained by rule. We may have cases in which there are two dependent variables, x and y, and one independent variable t, the differential coefficients dx/dt, dy/dt being given as functions of x, y and t. Of such equations a simple case is expressed by the pair
${\frac {dx}{dt}}{=}ax+by+c,{\frac {dy}{dt}}{=}a'x+b'y+c',$
wherein the coefficients a, b, c, a′, b′, c′, are constants. To integrate these, form with the constant λ the differential coefficient of z = x + λy, that is dz/dt = (a + λa′)x + (b + λb′)y + c + λc′, the quantity λ being so chosen that b + λb′ = λ(a + λa′), so that we have dz/dt = (a + λa′)z + c + λc′; this last equation is at once integrable in the form z(a + λa′) + c + λc′ = Ae$^{(a+λa′)t}$, where A is an arbitrary constant. In general, the condition b + λb′ = λ(a + λa′) is satisfied by two different values of λ, say λ$_1$, λ$_2$; the solutions corresponding to these give the values of x +λ$_1$y and x + λ$_2$y, from which x and y can be found as functions of t, involving two arbitrary constants. If, however, the two roots of the quadratic equation for λ are equal, that is, if (a − b′)² + 4a′b = 0, the method described gives only one equation, expressing x + λy in terms of t; by means of this equation y can be eliminated from dx/dt = ax + by + c, leading to an equation of the form dx/dt = Px + Q + Re$^{(a+λa′)t}$, where P, Q, R are constants. The integration of this gives x, and thence y can be found.
A similar process is applicable when we have three or more dependent variables whose differential coefficients in regard to the single independent variables are given as linear functions of the dependent variables with constant coefficients.
Another method of solution of the equations
dx/dt = ax + by + c, dy/dt = a′x + b′y + c′,
consists in differentiating the first equation, thereby obtaining
${\frac {d^{2}x}{dt^{2}}}{=}a{\frac {dx}{dt}}+b{\frac {dy}{dx}};$
from the two given equations, by elimination of y, we can express dy/dt as a linear function of x and dx/dt; we can thus form an equation of the shape d²x/dt² = P + Qx + Rdx/dt, where P, Q, R are constants; this can be integrated by methods previously explained, and the integral, involving two arbitrary constants, gives, by the equation dx/dt = ax + by + c, the corresponding value of y. Conversely it should be noticed that any single linear differential equation
${\frac {d^{2}x}{dt^{2}}}{=}u+vx+w{\frac {dx}{dt}},$
where u, v, w are functions of t, by writing y for dx/dt, is equivalent with the two equations dx/dt = y, dy/dt = u + vx + wy. In fact a similar reduction is possible for any system of differential equations with one independent variable.
Equations occur to be integrated of the form
Xdx + Ydy + Zdz = 0,
where X, Y, Z are functions of x, y, z. We consider only the case in which there exists an equation φ(x, y, z) = C whose differential
${\frac {\partial \phi }{\partial x}}dx+{\frac {\partial \phi }{\partial y}}dy+{\frac {\partial \phi }{\partial z}}dz{=}0$
is equivalent with the given differential equation; that is, μ being a proper function of x, y, z, we assume that there exist equations
${\frac {\partial \phi }{\partial x}}{=}\mu \mathrm {X} ,{\frac {\partial \phi }{\partial y}}{=}\mu \mathrm {Y} ,{\frac {\partial \phi }{\partial z}}{=}\mu \mathrm {Z}$
these equations require
${\frac {\partial }{\partial z}}(\mu \mathrm {Y} ){=}{\frac {\partial }{\partial y}}(\mu \mathrm {Z} ),$ &c.,
and hence
$\mathrm {X} \left({\frac {\partial \mathrm {Z} }{\partial y}}-{\frac {\partial \mathrm {Y} }{\partial z}}\right)+\mathrm {Y} \left({\frac {\partial \mathrm {X} }{\partial z}}-{\frac {\partial \mathrm {Z} }{\partial x}}\right)+\mathrm {Z} \left({\frac {\partial \mathrm {Y} }{\partial x}}-{\frac {\partial \mathrm {X} }{\partial y}}\right){=}0;$
conversely it can be proved that this is sufficient in order that μ may exist to render μ(Xdx + Ydy + Zdz) a perfect differential; in particular it may be satisfied in virtue of the three equations such as
${\frac {\partial \mathrm {Z} }{\partial y}}-{\frac {\partial \mathrm {Y} }{\partial z}}{=}0$
in which case we may take μ = 1. Assuming the condition in its general form, take in the given differential equation a plane section of the surface φ = C parallel to the plane z, viz. put z constant, and consider the resulting differential equation in the two variables x, y, namely Xdx + Ydy = 0; let ψ(x, y, z) = constant, be its integral, the constant z entering, as a rule, in ψ because it enters in X and Y. Now differentiate the relation ψ(x, y, z) = ƒ(z), where ƒ is a function to be determined, so obtaining
${\frac {\partial \psi }{\partial x}}dx+{\frac {\partial \psi }{\partial y}}dy+\left({\frac {\partial \psi }{\partial z}}-{\frac {df}{dz}}\right)dz{=}0;$
there exists a function σ of x, y, z such that
${\frac {\partial \psi }{\partial x}}{=}\sigma \mathrm {X} {\frac {\partial \psi }{\partial y}}{=}\sigma \mathrm {Y}$
because ψ = constant, is the integral of Xdx + Ydy = 0; we desire to prove that ƒ can be chosen so that also, in virtue of ψ(x, y, z) = ƒ(z), we have
${\frac {\partial \psi }{\partial z}}-{\frac {df}{dz}}{=}\sigma \mathrm {Z} ,$ namely ${\frac {df}{dz}}{=}{\frac {\partial \psi }{\partial z}}-\sigma \mathrm {Z} ;$
if this can be proved the relation ψ(x, y, z) − ƒ(z) = constant, will be the integral of the given differential equation. To prove this it is enough to show that, in virtue of ψ(x, y, z) = ƒ(z), the function ∂ψ/∂x − σZ can be expressed in terms of z only. Now in consequence of the originally assumed relations,
${\frac {\partial \psi }{\partial x}}{=}\mu \mathrm {X} ,{\frac {\partial \phi }{\partial y}}{=}\mu \mathrm {Y} ,{\frac {\partial \phi }{\partial z}}{=}\mu \mathrm {Z}$
we have
${\frac {\partial \psi }{\partial x}}\left/{\frac {\partial \phi }{\partial x}}\right.{=}{\frac {\sigma }{\mu }}{=}{\frac {\partial \psi }{\partial y}}\left/{\frac {\partial \phi }{\partial y}}\right.,$
and hence
${\frac {\partial \psi }{\partial x}}{\frac {\partial \phi }{\partial y}}-{\frac {\partial \psi }{\partial y}}{\frac {\partial \phi }{\partial x}}{=}0;$
this shows that, as functions of x and y, ψ is a function of φ (see the note at the end of part i. of this article, on Jacobian determinants), so that we may write ψ = F(z, φ), from which
${\frac {\sigma }{\mu }}{=}{\frac {\partial \mathrm {F} }{\partial \phi }};$ then ${\frac {\partial \psi }{\partial z}}{=}{\frac {\partial \mathrm {F} }{\partial z}}+{\frac {\partial \mathrm {F} }{\partial \phi }}{\frac {\partial \phi }{\partial z}}{=}{\frac {\partial \mathrm {F} }{\partial z}}+{\frac {\sigma }{\mu }}\cdot \mu \mathrm {Z} {=}{\frac {\partial \mathrm {F} }{\partial z}}+\sigma \mathrm {Z}$ or ${\frac {\partial \psi }{\partial z}}-\sigma \mathrm {Z} {=}{\frac {\partial \mathrm {F} }{\partial z}};$
in virtue of ψ(x, y, z) = ƒ(z), and ψ = F(z, φ), the function φ can be written in terms of z only, thus ∂F/∂z can be written in terms of z only, and what we required to prove is proved.
Consider lastly a simple type of differential equation containing two independent variables, say x and y, and one dependent variable z, namely the equation
$\mathrm {P} {\frac {\partial z}{\partial x}}+\mathrm {Q} {\frac {\partial z}{\partial y}}{=}\mathrm {R} ,$
where P, Q, R are functions of x, y, z. This is known as Lagrange’s linear partial differential equation of the first order. To integrate this, consider first the ordinary differential equations dx/dz = P/R, dy/dz = Q/R, and suppose that two functions u, v, of x, y, z can be determined, independent of one another, such that the equations u = a, v = b, where a, b are arbitrary constants, lead to these ordinary differential equations, namely such that
$\mathrm {P} {\frac {\partial u}{\partial x}}+\mathrm {Q} {\frac {\partial u}{\partial y}}+\mathrm {R} {\frac {\partial u}{\partial z}}{=}0$ and $\mathrm {P} {\frac {\partial v}{\partial x}}+\mathrm {Q} {\frac {\partial v}{\partial y}}+\mathrm {R} {\frac {\partial v}{\partial z}}{=}0.$
Then if F(x, y, z) = 0 be a relation satisfying the original differential equations, this relation giving rise to
${\frac {\partial \mathrm {F} }{\partial x}}+{\frac {\partial \mathrm {F} }{\partial z}}{\frac {\partial z}{\partial x}}{=}0$ and ${\frac {\partial \mathrm {F} }{\partial y}}+{\frac {\partial \mathrm {F} }{\partial z}}{\frac {\partial z}{\partial y}}{=}0,$ we have $\mathrm {P} {\frac {\partial \mathrm {F} }{\partial x}}+\mathrm {Q} {\frac {\partial \mathrm {F} }{\partial y}}+\mathrm {R} {\frac {\partial \mathrm {F} }{\partial z}}{=}0.$
It follows that the determinant of three rows and columns vanishes whose first row consists of the three quantities ∂F/∂x, ∂F/∂y, ∂F/∂z, whose second row consists of the three quantities ∂u/∂x, ∂u/∂y, ∂u/∂z, whose third row consists similarly of the partial derivatives of v. The vanishing of this so-called Jacobian determinant is known to imply that F is expressible as a function of u and v, unless these are themselves functionally related, which is contrary to hypothesis (see the note below on Jacobian determinants). Conversely, any relation φ(u, v) = 0 can easily be proved, in virtue of the equations satisfied by u and v, to lead to
$\mathrm {P} {\frac {dz}{dx}}+\mathrm {Q} {\frac {dz}{dy}}{=}\mathrm {R} .$
The solution of this partial equation is thus reduced to the solution of the two ordinary differential equations expressed by dx/P = dy/Q = dz/R. In regard to this problem one remark may be made which is often of use in practice: when one equation u = a has been found to satisfy the differential equations, we may utilize this to obtain the second equation v = b; for instance, we may, by means of u = a, eliminate z—when then from the resulting equations in x and y a relation v = b has been found containing x and y and a, the substitution a = u will give a relation involving x, y, z.
Note on Jacobian Determinants.—The fact assumed above that the vanishing of the Jacobian determinant whose elements are the partial derivatives of three functions F, u, v, of three variables x, y, z, ​involves that there exists a functional relation connecting the three functions F, u, v, may be proved somewhat roughly as follows:—
The corresponding theorem is true for any number of variables. Consider first the case of two functions p, q, of two variables x, y. The function p, not being constant, must contain one of the variables, say x; we can then suppose x expressed in terms of y and the function p; thus the function q can be expressed in terms of y and the function p, say q = Q(p, y). This is clear enough in the simplest cases which arise, when the functions are rational. Hence we have
${\frac {\partial q}{\partial x}}{=}{\frac {\partial \mathrm {Q} }{\partial p}}{\frac {\partial p}{\partial x}}$ and ${\frac {\partial q}{\partial y}}{=}{\frac {\partial \mathrm {Q} }{\partial p}}{\frac {\partial p}{\partial y}}+{\frac {\partial \mathrm {Q} }{\partial y}};$
these give
${\frac {\partial p}{\partial x}}{\frac {\partial q}{\partial y}}-{\frac {\partial p}{\partial y}}{\frac {\partial q}{\partial x}}{=}{\frac {\partial p}{\partial x}}{\frac {\partial \mathrm {Q} }{\partial y}};$
by hypothesis ∂p/∂x is not identically zero; therefore if the Jacobian determinant of p and q in regard to x and y is zero identically, so is ∂Q/∂y, or Q does not contain y, so that q is expressible as a function of p only. Conversely, such an expression can be seen at once to make the Jacobian of p and q vanish identically.
Passing now to the case of three variables, suppose that the Jacobian determinant of the three functions F, u, v in regard to x, y, z is identically zero. We prove that if u, v are not themselves functionally connected, F is expressible as a function of u and v. Suppose first that the minors of the elements of ∂F/∂x, ∂F/∂y, ∂F/∂z in the determinant are all identically zero, namely the three determinants such as
${\frac {\partial u}{\partial y}}{\frac {\partial v}{\partial z}}-{\frac {\partial u}{\partial z}}{\frac {\partial v}{\partial y}};$
then by the case of two variables considered above there exist three functional relations. ψ$_1$(u, v, x) = 0, ψ$_2$(u, v, y) = 0, ψ$_3$(u, v, z) = 0, of which the first, for example, follows from the vanishing of
${\frac {\partial u}{\partial y}}{\frac {\partial v}{\partial z}}-{\frac {\partial u}{\partial z}}{\frac {\partial v}{\partial y}}.$
We cannot assume that x is absent from ψ$_1$, or y from ψ$_2$, or z from ψ$_3$; but conversely we cannot simultaneously have x entering in ψ$_1$, and y in ψ$_2$, and z in ψ$_3$, or else by elimination of u and v from the three equations ψ$_1$ = 0, ψ$_2$ = 0, ψ$_3$ = 0, we should find a necessary relation connecting the three independent quantities x, y, z; which is absurd. Thus when the three minors of ∂F/∂x, ∂F/∂y, ∂F/∂z in the Jacobian determinant are all zero, there exists a functional relation connecting u and v only. Suppose no such relation to exist; we can then suppose, for example, that
${\frac {\partial u}{\partial y}}{\frac {\partial v}{\partial z}}-{\frac {\partial u}{\partial z}}{\frac {\partial v}{\partial y}}$
is not zero. Then from the equations u(x, y, z) = u, v(x, y, z) = v we can express y and z in terms of u, v, and x (the attempt to do this could only fail by leading to a relation connecting u, v and x, and the existence of such a relation would involve that the determinant
${\frac {\partial u}{\partial y}}{\frac {\partial v}{\partial z}}-{\frac {\partial u}{\partial z}}{\frac {\partial v}{\partial y}}$
was zero), and so write F in the form F(x, y, z) = Φ(u, v, x). We then have
${\frac {\partial \mathrm {F} }{\partial x}}{=}{\frac {\partial \Phi }{\partial u}}{\frac {\partial u}{\partial x}}+{\frac {\partial \Phi }{\partial v}}{\frac {\partial v}{\partial x}}+{\frac {\partial \Phi }{\partial x}},{\frac {\partial \mathrm {F} }{\partial y}}{=}{\frac {\partial \Phi }{\partial u}}{\frac {\partial u}{\partial y}}+{\frac {\partial \Phi }{\partial v}}{\frac {\partial v}{\partial y}},{\frac {\partial \mathrm {F} }{\partial z}}{=}{\frac {\partial \Phi }{\partial u}}{\frac {\partial u}{\partial z}}+{\frac {\partial \Phi }{\partial v}}{\frac {\partial v}{\partial z}};$
thereby the Jacobian determinant of F, u, v is reduced to
${\frac {\partial \Phi }{\partial x}}\left({\frac {\partial u}{\partial y}}{\frac {\partial v}{\partial z}}-{\frac {\partial u}{\partial z}}{\frac {\partial v}{\partial y}}\right);$
by hypothesis the second factor of this does not vanish identically; hence ∂Φ/∂x = 0 identically, and Φ does not contain x; so that F is expressible in terms of u, v only; as was to be proved.
Part II.—General Theory.
Differential equations arise in the expression of the relations between quantities by the elimination of details, either unknown or regarded as unessential to the formulation of the relations in question. They give rise, therefore, to the two closely connected problems of determining what arrangement of details is consistent with them, and of developing, apart from these details, the general properties expressed by them. Very roughly, two methods of study can be distinguished, with the names Transformation-theories, Function-theories; the former is concerned with the reduction of the algebraical relations to the fewest and simplest forms, eventually with the hope of obtaining explicit expressions of the dependent variables in terms of the independent variables; the latter is concerned with the determination of the general descriptive relations among the quantities which are involved by the differential equations, with as little use of algebraical calculations as may be possible. Under the former heading we may, with the assumption of a few theorems belonging to the latter, arrange the theory of partial differential equations and Pfaff’s problem, with their geometrical interpretations, as at present developed, and the applications of Lie’s theory of transformation-groups to partial and to ordinary equations; under the latter, the study of linear differential equations in the manner initiated by Riemann, the applications of discontinuous groups, the theory of the singularities of integrals, and the study of potential equations with existence-theorems arising therefrom. In order to be clear we shall enter into some detail in regard to partial differential equations of the first order, both those which are linear in any number of variables and those not linear in two independent variables, and also in regard to the function-theory of linear differential equations of the second order. Space renders impossible anything further than the briefest account of many other matters; in particular, the theories of partial equations of higher than the first order, the function-theory of the singularities of ordinary equations not linear and the applications to differential geometry, are taken account of only in the bibliography. It is believed that on the whole the article will be more useful to the reader than if explanations of method had been further curtailed to include more facts.
When we speak of a function without qualification, it is to be understood that in the immediate neighbourhood of a particular set x$_0$, y$_0$, ... of values of the independent variables x, y, ... of the function, at whatever point of the range of values for x, y, ... under consideration x$_0$, y$_0$, ... may be chosen, the function can be expressed as a series of positive integral powers of the differences x − x$_0$, y − y$_0$, ..., convergent when these are sufficiently small (see Function: Functions of Complex Variables). Without this condition, which we express by saying that the function is developable about x$_0$, y$_0$, ..., many results provisionally stated in the transformation theories would be unmeaning or incorrect. If, then, we have a set of k functions, ƒ$_1$ ... ƒ$_k$ of n independent variables x$_1$ ... x$_n$, we say that they are independent when n ≥ k and not every determinant of k rows and columns vanishes of the matrix of k rows and n columns whose r-th row has the constituents dƒ$_r$/dx$_1$, ... dƒ$_r$/dx$_n$; the justification being in the theorem, which we assume, that if the determinant involving, for instance, the first k columns be not zero for x$_1$ = xº$_1$ ... x$_n$ = xº$_n$, and the functions be developable about this point, then from the equations ƒ$_1$ = c$_1$, ... ƒ$_k$ = c$_k$ we can express x$_1$, ... x$_k$ by convergent power series in the differences x$_{k+1}$ − xº$_{k+1}$, ... x$_n$ − x$_{nº}$, and so regard x$_1$, ... x$_k$ as functions of the remaining variables. This we often express by saying that the equations ƒ$_1$ = c$_1$, ... ƒ$_k$ = c$_k$ can be solved for x$_1$, ... x$_k$. The explanation is given as a type of explanation often understood in what follows.
We may conveniently begin by stating the theorem: If each of the n functions φ$_1$, ... φ$_n$ of the (n + 1) variables x$_1$, ... x$_n$t be developable Ordinary equations of the first order. about the values xº$_1$, ... x$_n$$^0$t$^0$, the n differential equations of the form dx$_1$/dt = φ$_1$(tx$_1$, ... x$_n$) are satisfied by convergent power series
x$_r$ = xº$_r$ + (t − t$^0$) A$_{r1}$ + (t − t$_0$)² A$_{r2}$ + ...
reducing respectively to xº$_1$, ... xº$_n$ when t = t$^0$; and the only functions satisfying the equations and reducing respectively to xº$_1$, ... xº$_n$ when t = t$^0$, are those determined by continuation of these series. If the result of solving these n equations for xº$_1$, ... xº$_n$ be written in the form ω$_1$(x$_1$, ... x$_n$t) = xº$_1$, ... ω$_n$(x$_1$, ... x$_n$t) = xº$_n$, Single homogeneous partial equation of the first order. it is at once evident that the differential equation
dƒ/dt + φ$_1$dƒ/dx$_1$ + ... + φ$_n$dƒ/dx$_n$ = 0
possesses n integrals, namely, the functions ω$_1$, ... ω$_n$, which are developable about the values (xº$_1$ ... x$_n$$^0$t$^0$) and reduce respectively to x$_1$, ... x$_n$ when t = t$^0$. And in fact it has no other integrals so reducing. Thus this equation also possesses a unique integral reducing when t = t$^0$ to an arbitrary function ψ(x$_1$, ... x$_n$), this integral being. ψ(ω$_1$, ... ω$_n$). Conversely the existence of these principal integrals ω$_1$, ... ω$_n$ of the partial equation establishes the existence of the specified solutions of the ordinary equations dx$_i$/dt = φ$_i$. The following sketch of the proof of the existence of these principal integrals for the case n = 2 will show the character of more general investigations. Put x for x − x$^0$, &c., and consider the equation a(xyt) dƒ/dx + b(xyt) dƒ/dy = dƒ/dt, wherein the functions a, b are developable about x = 0, y = 0, t = 0; say
a(xyt) = a$_0$ + ta$_1$ + t²a$_2$/2! + ..., b(xyt) = b$_0$ + tb$_1$ + t²b$_2$/2! + ...,
so that
ad/dx + bd/dy = δ$_0$ + tδ$_1$ + ½t²δ$_2$ + ...,
where δ = a$_r$d/dx + b$_r$d/dy. In order that
ƒ = p$_0$ + tp$_1$ + t²p$_2$/2! + ...
p$_1$ = δ$_0$p$_0$, p$_2$ = δ$_0$p$_1$ + δ$_1$p$_0$, &c.
and in generalProof of the existence of integrals.
p$_{s+1}$ = δ$_0$p$_s$ + s$_1$δ$_1$p$_{s−1}$ + s$_2$δ$_2$p$_{s−2}$ +... + δ$_s$p$_0$,
where
s$_r$ = (s!)/(r!) (s − r)!
Now compare with the given equation another equation
A(xyt)dF/dx + B(xyt)dF/dy = dF/dt,
wherein each coefficient in the expansion of either A or B is real and positive, and not less than the absolute value of the corresponding coefficient in the expansion of a or b. In the second equation let us substitute a series
F = P$_0$ + tP$_1$ + t²P$_2$/2! + ...,
wherein the coefficients in P$_0$ are real and positive, and each not less than the absolute value of the corresponding coefficient in p$_0$; then putting Δ$_r$ = A$_r$d/dx + B$_r$d/dy we obtain necessary equations of the same form as before, namely,
P$_1$ = Δ$_0$P$_0$, P$_2$ = Δ$_0$P$_1$ + Δ$_1$P$_0$, ...
and in general P$_{s+1}$ = Δ$_0$P$_s$, + s$_1$Δ$_1$P$_{s−1}$ + ... + Δ$_s$P$_0$. These give for every coefficient in P$_{s+1}$ an integral aggregate with real positive coefficients of the coefficients in P$_s$, P$_{s−1}$, ..., P$_0$ and the coefficients in A and B; and they are the same aggregates as would be given by the previously obtained equations for the corresponding coefficients in p$_{s+1}$ in terms of the coefficients in p$_s$, p$_{s−1}$, ..., p$_0$ and the coefficients in a and b. Hence as the coefficients in P$_0$ and also in A, B are real and positive, it follows that the values obtained in succession for the coefficients in P$_1$, P$_2$, ... are real and positive; and further, taking account of the fact that the absolute value of a sum of terms is not greater than the sum of the absolute values of the terms, it follows, for each value of s, that every coefficient in p$_{s+1}$ is, in absolute value, not greater than the corresponding coefficient in P$_{s+1}$. Thus if the series for F be convergent, the series for ƒ will also be; and we are thus reduced to (1), specifying functions A, B with real positive coefficients, each in absolute value not less than the corresponding coefficient in a, b; (2) proving that the equation
AdF/dx + BdF/dy = dF/dt
possesses an integral P$_0$ + tP$_1$ + t²P$_2$/2! + ... in which the coefficients in P$_0$ are real and positive, and each not less than the absolute value of the corresponding coefficient in p$_0$. If a, b be developable for x, y both in absolute value less than r and for t less in absolute value than R, and for such values a, b be both less in absolute value than the real positive constant M, it is not difficult to verify that we may take A = B = M[1 − (x + y)/r]$^{−1}$ (1 − t/R)$^{−1}$, and obtain
$\mathrm {F} =r-(r-x-y)\left[1-{\frac {4\mathrm {MR} }{r}}\left(1-{\frac {x+y}{r}}\right)^{2}\log \left(1-{\frac {t}{\mathrm {R} }}\right)^{-1}\right]^{\frac {1}{2}},$
and that this solves the problem when x, y, t are sufficiently small for the two cases p$_0$ = x, p$_0$ = y. One obvious application of the general theorem is to the proof of the existence of an integral of an ordinary linear differential equation given by the n equations dy/dx = y$_1$, dy$_1$/dx = y$_2$, ...,
dy$_{n−1}$/dx = p − p$_1$y$_{n−1}$ − ... − p$_n$y;
but in fact any simultaneous system of ordinary equations is reducible to a system of the form
dx$_i$/dt = φ$_i$(tx$_1$, ... x$_n$).
Suppose we have k homogeneous linear partial equations of the first order in n independent variables, the general equation being a$_{σ1}$dƒ/dx$_1$ + ... + a$_{σn}$dƒ/dx$_n$ = 0, where σ = 1, ... k, and that Simultaneous linear partial equations. we desire to know whether the equations have common solutions, and if so, how many. It is to be understood that the equations are linearly independent, which implies that k ≤ n and not every determinant of k rows and columns is identically zero in the matrix in which the i-th element of the σ-th row is a$_{σi}$}(i = 1, ... n, σ = 1, ... k). Denoting the left side of the σ-th equation by Pσƒ, it is clear that every common solution of the two equations P$_σ$ƒ = 0, P$_ρ$ƒ = 0, is also a solution of the equation P$_ρ$(p$_σ$ƒ), P$_σ$(p$_ρ$ƒ), We immediately find, however, that this is also a linear equation, namely, ΣH$_i$dƒ/dx$_i$ = 0 where H$_i$ = P$_ρ$a$_σ$ − P$_σ$a$_ρ$, and if it be not already contained among the given equations, or be linearly deducible from them, it may be added to them, as not introducing any additional limitation of the possibility of their having common solutions. Proceeding thus with every pair of the original equations, and then with every pair of the possibly augmented system so obtained, and so on continually, we shall arrive at a system of equations, linearly independent of each other and therefore not more than n in number, such that the combination, in the way described, of every pair of them, leads to an equation which is linearly deducible from them. If the number of this so-called complete system is n, the equations give dƒ/dx$_1$ = 0 ... dƒ/dx$_n$ = 0, leading to the nugatory result ƒ = a constant. Suppose, then, the number of this system to be r < n; suppose, further, that from the Complete systems of linear partial equations. matrix of the coefficients a determinant of r rows and columns not vanishing identically is that formed by the coefficients of the differential coefficients of ƒ in regard to x$_1$ ... x$_r$; also that the coefficients are all developable about the values x$_1$ = xº$_1$, ... x$_n$= xº$_n$, and that for these values the determinant just spoken of is not zero. Then the main theorem is that the complete system of r equations, and therefore the originally given set of k equations, have in common n − r solutions, say ω$_{r+1}$, ... ω$_n$, which reduce respectively to x$_{r+1}$, ... x$_n$ when in them for x$_1$, ... x$_r$ are respectively put xº$_1$, ... xº$_r$; so that also the equations have in common a solution reducing when x$_1$ = xº$_1$, ... x$_r$ = xº$_r$ to an arbitrary function ψ(x$_{r+1}$, ... x$_n$) which is developable about xº$_{r+1}$, ... xº$_n$, namely, this common solution is ψ(ω$_{r+1}$, ... ω$_n$). It is seen at once that this result is a generalization of the theorem for r = 1, and its proof is conveniently given by induction from that case. It can be verified without difficulty (1) that if from the r equations of the complete system we form r independent linear aggregates, with coefficients not necessarily constants, the new system is also a complete system; (2) that if in place of the independent variables x$_1$, ... x$_n$ we introduce any other variables which are independent functions of the former, the new equations also form a complete system. It is convenient, then, from the complete system of r equations to form r new equations by solving separately for dƒ/dx$_1$, ..., dƒ/dx$_r$; suppose the general equation of the new system to be
Q$_σ$ƒ = dƒ/dx$_σ$ + c$_{σjr+1}$dƒ/dx$_{r+1}$ + ... + c$_{σn}$dƒ/dx$_n$ = 0 (σ = 1, ... r).
Then it is easily obvious that the equation Q$_ρ$Q$_σ$ƒ − Q$_σ$Q$_ρ$ƒ = 0 contains only the differential coefficients of ƒ in regard to x$_{r+1}$ ... x$_n$; as it is at most a linear function of Q$_1$ƒ, ... Q$_r$ƒ, it must be identically zero. So reduced the system is called a Jacobian system. Of this system Q$_1$ƒ=0 has n − 1 principal solutions reducing respectively Jacobian systems. to x$_2$, ... x$_n$ when
x$_1$ = xº$_1$,
and its form shows that of these the first r − 1 are exactly x$_2$ ... x$_r$. Let these n − 1 functions together with x$_1$ be introduced as n new independent variables in all the r equations. Since the first equation is satisfied by n − 1 of the new independent variables, it will contain no differential coefficients in regard to them, and will reduce therefore simply to dƒ/dx$_1$ = 0, expressing that any common solution of the r equations is a function only of the n − 1 remaining variables. Thereby the investigation of the common solutions is reduced to the same problem for r − 1 equations in n − 1 variables. Proceeding thus, we reach at length one equation in n − r + 1 variables, from which, by retracing the analysis, the proposition stated is seen to follow.
The analogy with the case of one equation is, however, still closer. With the coefficients c$_{σj}$, of the equations Q$_σ$ƒ = 0 in transposed array (σ = 1, ... r, j = r + 1, ... n) we can put down the (n − r) equations, dx$_j$ = c$_{1j}$dx$_1$ + ... + c$_{rj}$dx$_r$, equivalent to System of total differential equations. the r(n − r) equations dx$_j$/dx$_σ$ = c$_{σr}$. That consistent with them we may be able to regard x$_{r+1}$, ... x$_n$ as functions of x$_1$, ... x$_r$, these being regarded as independent variables, it is clearly necessary that when we differentiate c$_{σj}$ in regard to x$_ρ$ on this hypothesis the result should be the same as when we differentiate c$_{ρj}$, in regard to x$_σ$ on this hypothesis. The differential coefficient of a function ƒ of x$_1$, ... x$_n$ on this hypothesis, in regard to x$_{ρj}$ is, however,
dƒ/dx$_ρ$ + c$_{ρjr+1}$dƒ/dx$_{r+1}$ + ... + c$_{ρn}$dƒ/dx$_n$,
namely, is Q$_ρ$ƒ. Thus the consistence of the n − r total equations requires the conditions Q$_ρ$c$_{σj}$ − Q$_σ$c$_{ρj}$ = 0, which are, however, verified in virtue of Q$_ρ$(Q$_σ$ƒ) − Q$_σ$(Q$_ρ$ƒ) = 0. And it can in fact be easily verified that if ω$_{r+1}$, ... ω$_n$ be the principal solutions of the Jacobian system, Q$_σ$ƒ = 0, reducing respectively to x$_{r+1}$, ... x$_n$ when x$_1$ = xº$_1$, ... x$_r$ = xº$_r$, and the equations ω$_{r+1}$ = x$^0$$_{r+1}$, ... ω$_n$ = xº$_n$ be solved for x$_{r+1}$, ... x$_n$ to give x$_j$ = ψ$_j$(x$_1$, ... x$_r$, x$^0$$_{r+1}$, ... xº$_n$), these values solve the total equations and reduce respectively to x$^0$$_{r+1}$, ... xº$_n$ when x$_1$ = xº$_1$ ... x$_r$ = xº$_r$. And the total equations have no other solutions with these initial values. Conversely, the existence of these solutions of the total equations can be deduced a priori and the theory of the Jacobian system based upon them. The theory of such total equations, in general, finds its natural place under the heading Pfaffian Expressions, below.
A practical method of reducing the solution of the r equations of a Jacobian system to that of a single equation in n − r + 1 variables may be explained in connexion with a geometrical interpretation which will perhaps be clearer in a particular Geometrical interpretation and solution. case, say n = 3, r = 2. There is then only one total equation, say dz = adz + bdy; if we do not take account of the condition of integrability, which is in this case da/dy + bda/dz = db/dx + adb/dz, this equation may be regarded as defining through an arbitrary point (x$_0$, y$_0$, z$_0$) of three-dimensioned space (about which a, b are developable) a plane, namely, z − z$_0$ = a$_0$(x − x$_0$) + b$_0$(y − y$_0$), and therefore, through this arbitrary point ∞$^2$ directions, namely, all those in the plane. If now there be a surface z = ψ(x, y), satisfying dz = adz + bdy and passing through (x$_0$, y$_0$, z$_0$), this plane will touch the surface, and the operations of passing along the surface from (x$_0$, y$_0$, z$_0$) to
(x$_0$ + dx$_0$, y$_0$, z$_0$ + dz$_0$)
and then to (x$_0$ + dx$_0$, y$_0$ + dy$_0$, Z$_0$ + d$^1$z$_0$), ought to lead to the same value of d$^1$z$_0$ as do the operations of passing along the surface from (x$_0$, y$_0$, z$_0$) to (x$_0$, y$_0$ + dy$_0$, z$_0$ + δz$_0$), and then to
(x$_0$ + dx$_0$, y$_0$ + dy$_0$, z$_0$ + δ$^1$z$_0$),
namely, δ$^1$z$_0$ ought to be equal to d$^1$z$_0$. But we find
$a_{0}dx_{0}+b_{0}dy_{0}+dx_{0}dy_{0}\left({\frac {db}{dx_{0}}}+a_{0}{\frac {db}{dz_{0}}}\right),$
and so at once reach the condition of integrability. If now we put ​x = x$_o$ + t, y = y$_o$ + mt, and regard m as constant, we shall in fact be considering the section of the surface by a fixed plane y−y$_o$ = m(x−x$_o$); along this section dz = dt(a + bm); if we then integrate the equation dx/dt = a + bm, where a, b are expressed as functions of m and t, with m kept constant, finding the solution which reduces to z$_o$ for t = 0, and in the result again replace m by (y−y$_o$)/(x−x$_o$), we shall have the surface in question. In the general case the equations
dx$_j$＝c$_{1j}$ dx$_1$ + . .c$_{rj}$dx$_r$
similarly determine through an arbitrary point x$_1$ᵒ, . . . x$_n$ᵒ Mayer’s method of integration. a planar manifold of r dimensions in space of n dimensions, and when the conditions of integrability are satisfied, every direction in this manifold through this point is tangent to the manifold of r dimensions, expressed by ω$_{r+1}$ = x$^0$$_{r+1}$, . . . ω$_n$ = x$_n$ᵒ, which satisfies the equations and passes through this point. If we put x$_1$−x$_1$ᵒ = t, x$_2$−x$_2$ᵒ = m$_2$t, ... x$_r$−x$_r$ᵒ = m$_r$t, and regard m$_2$, ... m$_r$ as fixed, the (n−r) total equations take the form dx$_j$ /dt = c$_{1j}$ + m$_2$c$_{2j}$ + ... + m$_r$c$_{rj}$, and their integration is equivalent to that of the single partial equation
${df}/{dt}+\sum _{j=r+1}^{n}(c_{1j}+m_{2}c_{2j}+\ldots +m_{r}c_{rj}){df}/{dx}_{j}=0$
in the n−r + 1 variables t, x$_{r+1}$, ... x$_n$. Determining the solutions Ω$_{r+1}$, ... Ω$_n$ which reduce to respectively x$_{r+1}$, ... x$_n$ when t = 0, and substituting t = x$_1$−x$_1$ᵒ, m$_2$ = (x$_2$−x$_2$ᵒ)/(x$_1$−x$_1$ᵒ), ... m$_r$ = (x$_r$−x$_r$ᵒ)/(x$_1$−x$_1$ᵒ), we obtain the solutions of the original system of partial equations previously denoted by ω$_{r+1}$, . . . ω$_n$. It is to be remarked, however, that the presence of the fixed parameters m$_2$, ... m$_r$ in the single integration may frequently render it more difficult than if they were assigned numerical quantities.
We have above considered the integration of an equation
dz＝adz + bdy
on the hypothesis that the condition
da/dy + bda/dz＝db/dz + adb/dz.
It is natural to inquire what relations among x, y, z, if any, are implied by, or are consistent with, a differential relation adx + bdy + cdx = 0, when a, b, c are unrestricted functions of x, y, z. This problem leads to the consideration of the so-called Pfaffian Expression adx + bdy + cdz. It can be shown (1) if each of the Pfaffian Expressions.quantities db/dz−dc/dy, dc/dx−da/dz, da/dy−db/dz, which we shall denote respectively by u$_{23}$, u$_{31}$, u$_{12}$, be identically zero, the expression is the differential of a function of x, y, z, equal to dt say; (2) that if the quantity au$_{23}$ + bu$_{31}$ + cu$_{12}$ is identically zero, the expression is of the form udt, i.e. it can be made a perfect differential by multiplication by the factor 1/u; (3) that in general the expression is of the form dt + u$_1$dt$_1$. Consider the matrix of four rows and three columns, in which the elements of the first row are a, b, c, and the elements of the (r + 1)-th row, for r = 1, 2, 3, are the quantities u$_{r1}$, u$_{r2}$, u$_{r3}$, where u$_{11}$ = u$_{22}$ = u$_{33}$ = 0. Then it is easily seen that the cases (1), (2), (3) above correspond respectively to the cases when (1) every determinant of this matrix of two rows and columns is zero, (2) every determinant of three rows and columns is zero, (3) when no condition is assumed. This result can be generalized as follows: if a$_1$, ... a$_n$ be any functions of x$_1$, ... x$_n$, the so-called Pfaffian expression a$_1$dx$_1$ + ... + a$_n$dx$_n$ can be reduced to one or other of the two forms
u$_1$dt$_1$ + ... + u$_k$dt$_k$, dt + u$_1$dt$_1$ + ... + u$_{k−1}$dt$_{k−1}$,
wherein t, u$_1$ ..., t$_1$, ... are independent functions of x$_1$, ... x$_n$, and k is such that in these two cases respectively 2k or 2k−1 is the rank of a certain matrix of n + 1 rows and n columns, that is, the greatest number of rows and columns in a non-vanishing determinant of the matrix; the matrix is that whose first row is constituted by the quantities a$_1$, ... a$_n$, whose s-th element in the (r + 1)-th row is the quantity da$_r$/dx$_s$−da$_s$/dx$_r$. The proof of such a reduced form can be obtained from the two results: (1) If t be any given function of the 2m independent variables u$_1$, ... u$_m$, t$_1$, ... t$_m$, the expression dt + u$_1$dt$_1$ + ... + u$_m$dt$_m$ can be put into the form u′$_1$dt′$_1$ + ... + u′$_m$dt′$_m$. (2) If the quantities u$_1$, ..., u$_1$, t$_1$, ... t$_m$ be connected by a relation, the expression n$_1$dt$_1$ + ... + u$_m$dt$_m$ can be put into the format dt′ + u′$_1$dt′$_1$ + ... + u′$_{m−1}$dt′$_{m−1}$; and if the relation connecting u$_1$, u$_m$, t$_1$, ... t$_m$ be homogeneous in u$_1$, ... u$_m$, then t′ can be taken to be zero. These two results are deductions from the theory of contact transformations (see below), and their demonstration requires, beside elementary algebraical considerations, only the theory of complete systems of linear homogeneous partial differential equations of the first order. When the existence of the reduced form of the Pfaffian expression containing only independent quantities is thus once assured, the identification of the number k with that defined by the specified matrix may, with some difficulty, be made a posteriori.
In all cases of a single Pfaffian equation we are thus led to consider what is implied by a relation dt−u$_1$dt$_1$−. . .−u$_m$dt$_m$ = 0, in which t, u$_1$, . . . u$_m$, t$_1$ . . ., t$_m$ are, except for this equation, independent variables. This is to be satisfied in virtue of Single linear Pfaffian equation. one or several relations connecting the variables; these must involve relations connecting t, t$_1$, . . . t$_m$ only, and in one of these at least t must actually enter. We can then suppose that in one actual system of relations in virtue of which the Pfaffian equation is satisfied, all the relations connecting t, t$_1$ . . . t$_m$ only are given by
t＝ψ(t$_{s+1}$ ... t$_m$), t$_1$＝ψ$_1$(t$_{s+1}$ ... t$_m$), ... t$_s$＝ψ$_s$(t$_{s+1}$ ... t$_m$);
so that the equation
dψ−u$_1$dψ$_1$−...−u$_s$dψ$_s$−u$_{s+1}$dt$_{s+1}$−...−u$_m$dt$_m$＝0
is identically true in regard to u$_1$, ... u$_m$, t$_{s+1}$ ..., t$_m$; equating to zero the coefficients of the differentials of these variables, we thus obtain m−s relations of the form
dψ/dt$_j$−u$_1$dψ$_1$/dt$_j$−. . .−u$_s$dψ$_s$/dt$_j$−u$_j$＝0;
these m−s relations, with the previous s + 1 relations, constitute a set of m + 1 relations connecting the 2m + 1 variables in virtue of which the Pfaffian equation is satisfied independently of the form of the functions ψ,ψ$_1$, ... ψ$_s$. There is clearly such a set for each of the values s = 0, s = 1, . . ., s = m−1, s = m. And for any value of s there may exist relations additional to the specified m + 1 relations, provided they do not involve any relation connecting t, t$_1$, . . . t$_m$ only, and are consistent with the m−s relations connecting u$_1$, ... u$_m$. It is now evident that, essentially, the integration of a Pfaffian equation
a$_1$dx$_1$ + ... + a$_n$dx$_n$＝0,
wherein a$_1$, ... a$_n$ are functions of x$_1$, ... x$_n$, is effected by the processes necessary to bring it to its reduced form, involving only independent variables. And it is easy to see that if we suppose this reduction to be carried out in all possible ways, there is no need to distinguish the classes of integrals corresponding to the various values of s; for it can be verified without difficulty that by putting t′ = t−u$_1$t$_1$−...−u$_s$t$_s$, t′$_1$ = u$_1$, ... t′$_s$ = u$_s$, u′$_1$ = −t$_1$, ..., u′$_s$ = −t$_s$, t′$_{s+1}$ = t$_{s+1}$, ... t′$_m$ = t$_m$, u′$_{s+1}$ = u$_{s+1}$, ... u′$_m$ = u$_m$, the reduced equation becomes changed to dt′−u′$_1$dt′$_1$− ...−u′$_m$dt′$_m$ = 0, and the general relations changed to
t′＝ψ(t′$_{s+1}$, ... t′$_m$)−t′$_1$ψ$_1$(t′$_{s+1}$, ... t′$_m$)− ...−t′$_s$ψ$_s$(t′$_{s+1}$, ... t′$_m$),＝φ,
say, together with u′$_1$ = dφ/dt′$_1$, ..., u′$_m$ = dφ/dt′$_m$, which contain only one relation connecting the variables t′, t′$_1$, ... t′$_m$ only.
This method for a single Pfaffian equation can, strictly speaking, be generalized to a simultaneous system of (n−r) Pfaffian equations dx$_j$ = c$_{1j}$dx$_1$ + ... + c$_{rj}$dx$_r$ only in the case already treated, Simultaneous Pfaffian equations. when this system is satisfied by regarding x$_{r+1}$, ... x$_n$ as suitable functions of the independent variables x$_1$, . . . x$_r$; in that case the integral manifolds are of r dimensions. When these are non-existent, there may be integral manifolds of higher dimensions; for if
dφ＝φ$_1$dx$_1$ + ... + φ$_r$dx$_r$ + φ$_{r+1}$(c$_{1,r+1}$dx$_1$ + ... + c$_{r,r+1}$dx$_r$) + φ$_{r+2}$( ) + ...
be identically zero, then φσ + cσ,$_{r+1}$φ$_{r+1}$ + ... + cσ,$_n$φ$_n$ ≈ 0, or φ satisfies the r partial differential equations previously associated with the total equations; when these are not a complete system, but included in a complete system of r−μ equations, having therefore n−r−μ independent integrals, the total equations are satisfied over a manifold of r + μ dimensions (see E. v. Weber, Math. Annal. 1v. (1901), p. 386).
It seems desirable to add here certain results, largely of algebraic character, which naturally arise in connexion with the theory of contact transformations. For any two functions of the 2n Contact transformations. independent variables x$_1$, ... x$_n$, p$_1$, ... p$_n$ we denote by (φψ) the sum of the n terms such as dφdψ/dp$_i$dx$_i$−dψdφ/dp$_i$dx$_i$ For two functions of the (2n + 1) independent variables z, x$_1$, ... x$_n$, p$_1$, ... p$_n$ we denote by φψ the sum of the n terms such as
${\frac {d\phi }{dp_{i}}}\left({\frac {d\psi }{dx_{i}}}+p_{i}{\frac {d\psi }{dz}}\right)-{\frac {d\psi }{dp_{i}}}\left({\frac {d\phi }{dx_{i}}}+p_{i}{\frac {d\phi }{dz}}\right).$
It can at once be verified that for any three functions [ƒ[φψ]] + [φ[ψƒ]] + [ψ[ƒφ]] = dƒ/dz [φψ] + dφ/dz [ψƒ] + dψ/dz [ƒφ], which when ƒ, φ,ψ do not contain z becomes the identity (ƒ(φψ)) + (φ(ψƒ)) + (ψ(ƒφ)) = 0.Then, if X$_1$, ... X$_n$, P$_1$, ... P$_n$ be such functions Of x$_1$, ... x$_n$, p$_1$ ... p$_n$ that P$_1$dX$_1$ + ... + P$_n$dX$_n$ is identically equal to p$_1$dx$_1$ + ... + p$_n$dx$_n$, it can be shown by elementary algebra, after equating coefficients of independent differentials, (1) that the functions X$_1$, ... P$_n$ are independent functions of the 2n variables x$_1$, ... p$_n$, so that the equations x′$_i$ = X$_i$, p′$_i$ = P$_i$ can be solved for x$_1$, ... x$_n$, p$_1$, ... p$_n$, and represent therefore a transformation, which we call a homogeneous contact transformation; (2) that the X$_1$, ... X$_n$ are homogeneous functions of p$_1$, ... p$_n$ of zero dimensions, the P$_1$, ... P$_n$ are homogeneous functions of p$_1$, ... p$_n$ of dimension one, and the 1/2n(n−1) relations (X$_i$X$_j$) = 0 are verified. So also are the n² relations (P$_i$X$_i$ = 1, (P$_i$X$_j$) = 0, (P$_i$P$_j$) = 0. Conversely, if X$_1$, ... X$_n$ be independent functions, each homogeneous of zero dimension in p$_1$, ... p$_n$ satisfying the 1/2n(n−1) relations (X$_i$X$_j$) = 0, then P$_1$, ... P$_n$ can be uniquely determined, by solving linear algebraic equations, such that P$_1$dX$_1$ + ... + P$_n$dX$_n$ = p$_1$dx$_1$ + ... + p$_n$dx$_n$. If now we put n + 1 for n, put z for x$_{n+1}$, Z for X$_{n+1}$, Q$_i$ for -P$_i$/P$_{n+1}$, for i = 1, ... n, put q$_i$ for -p$_i$/p$_{n+1}$ and σ for q$_{n+1}$/Q$_{n+1}$, and then finally write P$_1$, ... P$_n$, p$_1$, ... p$_n$ for Q$_1$, ... Q$_n$, q$_1$, ... q$_n$, we obtain the following results: If ZX$_1$ ... X$_n$, P$_1$, ... P$_n$ be functions of z, x$_1$, ... x$_n$, p$_1$, ... p$_n$, such that the expression dZ−P$_1$dX$_1$−...−P$_n$dX$_n$ is identically equal to σ(dz−p$_1$dx$_1$−...−p$_n$dx$_n$), and σ not zero, then (1) the functions Z, X$_1$, ... X$_n$, P$_1$, ... P$_n$ are independent functions of z, x$_1$, ... x$_n$, p$_1$, ... p$_n$, so that the equations z′ = Z, x′$_i$ = X$_i$, p′$_i$ = P$_i$ can be solved for z, x$_1$, ... x$_n$, p$_1$, ... p$_n$ and determine a transformation which we call a (non-homogeneous) contact transformation; (2) the Z, X$_1$, ... X$_n$ verify the 1/2n(n + 1) ​identities [ZX$_i$] = 0, [X$_i$X$_j$] = 0. And the further identities
$[\mathrm {Z} \sigma ]{=}\sigma {\frac {d\mathrm {Z} }{dz}}-\sigma ^{2},[\mathrm {X} _{i}\sigma ]{=}\sigma {\frac {d\mathrm {X} _{i}}{dz}},[\mathrm {P} _{i}\sigma ]{=}{\frac {d\mathrm {P} _{i}}{dz}}$
are also verified. Conversely, if Z, x$_1$, ... X$_n$ be independent functions satisfying the identities [ZX$_i$] = 0, [X$_i$X$_j$] = 0, then σ, other than zero, and P$_1$, ... P$_n$ can be uniquely determined, by solution of algebraic equations, such that
dZ − P$_1$dX$_1$ − ... − P$_n$dX$_n$＝σ(dz − p$_1$dx$_1$ − ... − p$_n$dx$_n$).
Finally, there is a particular case of great importance arising when σ = 1, which gives the results: (1) If U, X$_1$, ... X$_n$, P$_1$, ... P$_n$ be 2n + 1 functions of the 2n independent variables x$_1$, ... x$_n$, p$_1$, ... p$_n$, satisfying the identity
dU + P$_1$dx$_1$ + ... + P$_n$dX$_n$＝p$_1$dx$_1$ + ... + p$_n$dx$_n$,
then the 2n functions P$_1$, ... P$_n$, X$_1$, ... X$_n$ are independent, and we have
(X$_i$X$_j$)＝0, (X$_i$U)＝δX$_i$, (P$_i$X$_i$)＝1, (P$_i$X$_j$)＝0, (P$_i$P$_j$)＝0, (P$_i$U) + P$_i$＝δP$_i$,
where δ denotes the operator p$_1$d/dp$_1$ + ... + p$_n$d/dp$_n$; (2) If X$_1$, ... X$_n$ be independent functions of x$_1$, ... x$_n$, p$_1$, ... p$_n$, such that (X$_i$X$_j$) = 0, then U can be found by a quadrature, such that
(X$_i$U)＝δX$_i$;
and when X$_i$, ... X$_n$, U satisfy these 1/2n(n + 1) conditions, then P$_1$, ... P$_n$ can be found, by solution of linear algebraic equations, to render true the identity dU + P$_1$dX$_1$ + ... + P$_n$dX$_n$ = p$_1$dx$_1$ + ... + p$_n$dx$_n$; (3) Functions X$_1$, ... X$_n$, P$_1$, ... P$_n$ can be found to satisfy this differential identity when U is an arbitrary given function of x$_1$, ... x$_n$, p$_1$, ... p$_n$; but this requires integrations. In order to see what integrations, it is only necessary to verify the statement that if U be an arbitrary given function of x$_1$, ... x$_n$, p$_1$, ... p$_n$, and, for r < n, X$_1$, ... X$_r$ be independent functions of these variables, such that (XσU) = δXσ, (XρXσ) = 0, for ρ, σ = 1 ... r, then the r + 1 homogeneous linear partial differential equations of the first order (Uƒ) + δƒ = 0, (Xρƒ) = 0, form a complete system. It will be seen that the assumptions above made for the reduction of Pfaffian expressions follow from the results here enunciated for contact transformations.
We pass on now to consider the solution of any partial differential equation of the first order; we attempt to explain certain ideas relatively to a single equation with any number of independent variables (in particular, an ordinary equation of the first order with one independent Partial differential equation of the first order.variable) by speaking of a single equation with two independent variables x, y, and one dependent variable z. It will be seen that we are naturally led to consider systems of such simultaneous equations, which we consider below. The central discovery of the transformation theory of the solution of an equation F(x, y, z, dz/dx, dz/dy) = 0 is that its solution can always be reduced to the solution of partial equations which are linear. For this, however, we must regard dz/dx, dz/dy, during the process of integration, not as the differential coefficients of a function z in regard to x and y, but as variables independent of x, y, z, the too great indefiniteness that might thus appear to be introduced being provided for in another way. We notice that if z = ψ(x, y) be a solution of the differential equation, then dz = dxdψ/dx + dydψ/dy; thus if we denote the equation by F(x, y, z, p, q,) = 0, and prescribe the condition dz = pdx + qdy for every solution, any solution such as z = ψ(x, y) will necessarily be associated with the equations p = dz/dx, q = dz/dy, and z will satisfy the equation in its original form. We have previously seen (under Pfaffian Expressions) that if five variables x, y, z, p, q, otherwise independent, be subject to dz − pdx − qdy = 0, they must in fact be subject to at least three mutual relations. If we associate with a point (x, y, z) the plane
Z − z＝p(X − x) + q(Y − y)
passing through it, where X, Y, Z are current co-ordinates, and call this association a surface-element; and if two consecutive elements of which the point(x + dx, y + dy, z + dz) of one lies on the plane of the other, for which, that is, the condition dz = pdx + qdy is satisfied, be said to be connected, and an infinity of connected elements following one another continuously be called a connectivity, then our statement is that a connectivity consists of not more than ∞² elements, the whole number of elements (x, y, z, p, q) that are possible being called ∞$^5$. The solution of an equation F(x, y, z, dz/dx, dz/dy) = 0 is then to be understood to mean finding in all possible ways, from the ∞$^4$ elements (x, y, z, p, q) which satisfy F(x, y, z, p, q) = 0 a set of ∞² elements forming a connectivity; or, more analytically, finding in all possible ways two relations G = 0, H = 0 connecting x, y, z, p, q and independent of F = 0, so that the three relations together may involve
dz＝pdx + qdy.
Such a set of three relations may, for example, be of the form z = ψ(x, y), p = dψ/dx, q = dψ/dy; but it may also, as another case, involve two relations z = ψ(y), x = ψ$_1$(y) connecting x, y, z, the third relation being
ψ′(y)＝pψ′$_1$(y) + q,
the connectivity consisting in that case, geometrically, of a curve in space taken with ∞¹ of its tangent planes; or, finally, a connectivity is constituted by a fixed point and all the planes passing through that point. This generalized view of the meaning of a solution of F = 0 is of advantage, moreover, in view of anomalies otherwise arising from special forms of the equation Meaning of a solution of the equation. itself. For instance, we may include the case, sometimes arising when the equation to be solved is obtained by transformation from another equation, in which F does not contain either p or q. Then the equation has ∞² solutions, each consisting of an arbitrary point of the surface F = 0 and all the ∞² planes passing through this point; it also has ∞² solutions, each consisting of a curve drawn on the surface F = 0 and all the tangent planes of this curve, the whole consisting of ∞² elements; finally, it has also an isolated (or singular) solution consisting of the points of the surface, each associated with the tangent plane of the surface thereat, also ∞² elements in all. Or again, a linear equation F = Pp + Qq − R = 0, wherein P, Q, R are functions of x, y, z only, has ∞² solutions, each consisting of one of the curves defined by
dx/P＝dy/Q＝dz/R
taken with all the tangent planes of this curve; and the same equation has ∞² solutions, each consisting of the points of a surface containing ∞¹ of these curves and the tangent planes of this surface. And for the case of n variables there is similarly the possibility of n + 1 kinds of solution of an equation F(x$_1$, ... x$_n$, z, p$_1$, ... p$_n$) = 0; these can, however, by a simple contact transformation be reduced to one kind, in which there is only one relation z′ = ψ(x′$_1$, ... x′$_n$) connecting the new variables x′$_1$, ... x′$_n$, z′ (see under Pfaffian Expressions); just as in the case of the solution
z＝ψ(y), x＝ψ$_1$(y), ψ′(y)＝pψ′$_1$(y) + q
of the equation Pp + Qq = R the transformation z′ = z − px, x′ = p, p′ = −x, y′ = y, q′ = q gives the solution
z′＝ψ(y′) + x′ψ$_1$(y′), p′＝dz′/dx′, q′＝dz′/dy′
of the transformed equation. These explanations take no account of the possibility of p and q being infinite; this can be dealt with by writing p = −u/w, q = −v/w, and considering homogeneous equations in u, v, w, with udx + vdy + wdz = 0 as the differential relation necessary for a connectivity; in practice we use the ideas associated with such a procedure more often without the appropriate notation.
In utilizing these general notions we shall first consider the theory of characteristic chains, initiated by Cauchy, which shows well the nature of the relations implied by the given differential equation; the alternative ways of carrying out the necessary integrations are suggested by considering Order of the ideas.the method of Jacobi and Mayer, while a good summary is obtained by the formulation in terms of a Pfaffian expression.
Consider a solution of F = 0 expressed by the three independent equations F = 0, G = 0, H = 0. If it be a solution in which there is more than one relation connecting x, y, z, let new variables x′, y′, z′, p′, q′ be introduced, as before explained under Pfaffian Expressions, in which z′ is of the formCharacteristic chains.
z′＝z−p$_1$x$_1$−... −p$_s$x$_s$ (s＝1 or 2),
so that the solution becomes of a form z′ = ψ(x′y′), p′ = dψ/dx′, q′ = dψ/dy′, which then will identically satisfy the transformed equations F′ = 0, G′ = 0, H′ = 0. The equation F′ = 0, if x′, y′, z′ be regarded as fixed, states that the plane Z − z′ = p′(X − x′) + q′(Y − y′) is tangent to a certain cone whose vertex is (x′, y′, z′), the consecutive point (x′ + dx′, y′ + dy′, z′ + dz′) of the generator of contact being such that
$dx'\left/{\frac {d\mathrm {F} '}{dp'}}\right.=dy'\left/{\frac {d\mathrm {F} '}{dq'}}\right.=dz'\left/\left(p'{\frac {d\mathrm {F} '}{dp'}}+q'{\frac {d\mathrm {F} '}{dq'}}\right)\right..$
Passing in this direction on the surface z′ = ψ(x′, y′) the tangent ​plane of the surface at this consecutive point is (p′ + dp′, q′ + dq′), where, since F′(x′, y′, ψ, dψ/dx′, dψ/dy′)＝0 is identical, we have dx′ (dF′/dx′ + p′dF′/dz′) + dp′dF′/dp′＝0. Thus the equations, which we shall call the characteristic equations,
$dx'\left/{\frac {d\mathrm {F} '}{dp'}}\right.=dy'\left/{\frac {d\mathrm {F} '}{dq'}}\right.=dz'\left/\left(p'{\frac {d\mathrm {F} '}{dp'}}+q'{\frac {d\mathrm {F} '}{dq'}}\right)\right.=dp'\left/\left(-{\frac {d\mathrm {F} '}{dx'}}-p'{\frac {d\mathrm {F} '}{dz'}}\right)\right.=dq'\left/\left(-{\frac {d\mathrm {F} '}{dy'}}-q'{\frac {d\mathrm {F} '}{dz'}}\right)\right.$
are satisfied along a connectivity of ∞¹ elements consisting of a curve on z′＝ψ(x′, y′) and the tangent planes of the surface along this curve. The equation F′＝0, when p′, q′ are fixed, represents a curve in the plane Z − z′＝p′(X − x′) + q′(Y − y′) passing through (x′, y′, z′); if (x′ + δx′, y′ + δy′, z′ + δz′) be a consecutive point of this curve, we find at once
$\delta x'\left({\frac {d\mathrm {F} '}{dx'}}+p'{\frac {d\mathrm {F} '}{dz'}}\right)+\delta y'\left({\frac {d\mathrm {F} '}{dy'}}q'{\frac {d\mathrm {F} '}{dz'}}\right)=0;$
thus the equations above give δx′dp′ + δy′dq′＝0, or the tangent line of the plane curve, is, on the surface z′＝ψ(x′, y′), in a direction conjugate to that of the generator of the cone. Putting each of the fractions in the characteristic equations equal to dt, the equations enable us, starting from an arbitrary element x′$_0$, y′$_0$, z′$_0$, p′$_0$, q′$_0$, about which all the quantities F′, dF′/dp′, &c., occurring in the denominators, are developable, to define, from the differential equation F′＝0 alone, a connectivity of ∞¹ elements, which we call a characteristic chain; and it is remarkable that when we transform again to the original variables (x, y, z, p, q), the form of the differential equations for the chain is unaltered, so that they can be written down at once from the equation F＝0. Thus we have proved that the characteristic chain starting from any ordinary element of any integral of this equation F＝0 consists only of elements belonging to this integral. For instance, if the equation do not contain p, q, the characteristic chain, starting from an arbitrary plane through an arbitrary point of the surface F＝0, consists of a pencil of planes whose axis is a tangent line of the surface F＝0. Or if F＝0 be of the form Pp + Qq＝R, the chain consists of a curve satisfying dx/P＝dy/Q＝dz/R and a single infinity of tangent planes of this curve, determined by the tangent plane chosen at the initial point. In all cases there are ∞³ characteristic chains, whose aggregate may therefore be expected to exhaust the ∞$^4$ elements satisfying F＝0.
Consider, in fact, a single infinity of connected elements each satisfying F＝0, say a chain connectivity T, consisting of elements specified by x$_0$, y$_0$, z$_0$, p$_0$, q$_0$, which we suppose expressed as functions of a parameter u, so that Complete integral constructed with characteristic chains.
U$_0$＝dz$_0$/du − p$_0$dx$_0$/du − q$_0$dy$_0$/du
is everywhere zero on this chain; further, suppose that each of F, dF/dp, ... , dF/dx + pdF/dz is developable about each element of this chain T, and that T is not a characteristic chain. Then consider the aggregate of the characteristic chains issuing from all the elements of T. The ∞² elements, consisting of the aggregate of these characteristic chains, satisfy F＝0, provided the chain connectivity T consists of elements satisfying F＝0; for each characteristic chain satisfies dF＝0. It can be shown that these chains are connected; in other words, that if x, y, z, p, q, be any element of one of these characteristic chains, not only is
dz/dt − pdx/dt − qdy/dt＝0,
as we know, but also U＝dz/du − pdx/du − qdy/du is also zero. For we have
${\frac {d\mathrm {U} }{dt}}={\frac {d}{dt}}\left({\frac {dz}{du}}-p{\frac {dx}{du}}-q{\frac {dy}{du}}\right)-{\frac {d}{du}}\left({\frac {dz}{dt}}-p{\frac {dx}{dt}}-q{\frac {dy}{dt}}\right)$
$={\frac {dp}{du}}{\frac {dx}{dt}}-{\frac {dp}{dt}}{\frac {dx}{du}}+{\frac {dq}{du}}{\frac {dy}{dt}}-{\frac {dq}{dt}}{\frac {dy}{du}}$
which is equal to
${\frac {dp}{du}}{\frac {d\mathrm {F} }{dp}}+{\frac {dx}{du}}\left({\frac {d\mathrm {F} }{dx}}+p{\frac {d\mathrm {F} }{dz}}\right)+{\frac {dq}{du}}{\frac {d\mathrm {F} }{dq}}+{\frac {dy}{du}}\left({\frac {d\mathrm {F} }{dy}}+q{\frac {d\mathrm {F} }{dz}}\right)=-{\frac {d\mathrm {F} }{dz}}\mathrm {U} .$
As dF/dz is a developable function of t, this, giving
$\mathrm {U} =\mathrm {U} _{0}exp\left(-\int _{t_{0}}^{t}{\frac {d\mathrm {F} }{dz}}dt\right),$
shows that U is everywhere zero. Thus integrals of F＝0 are obtainable by considering the aggregate of characteristic chains issuing from arbitrary chain connectivities T satisfying F＝0; and such connectivities T are, it is seen at once, determinable without integration. Conversely, as such a chain connectivity T can be taken out from the elements of any given integral all possible integrals are obtainable in this way. For instance, an arbitrary curve in space, given by x$_0$＝θ(u), y$_0$＝φ(u), z$_0$＝ψ(u), determines by the two equations F(x$_0$, y$_0$, z$_0$, p$_0$, q$_0$)＝0, ψ′(u)＝p$_0$θ′(u) + q$_0$φ′(u), such a chain connectivity T, through which there passes a perfectly definite integral of the equation F＝0. By taking ∞² initial chain connectivities T, as for instance by taking the curves x$_0$＝θ, y$_0$ ＝φ, z$_0$ ＝ψ to be the ∞² curves upon an arbitrary surface, we thus obtain ∞² integrals, and so ∞$^4$ elements satisfying F＝0. In general, if functions G, H, independent of F, be obtained, such that the equations F＝0, G＝b, H＝c represent an integral for all values of the constants b, c, these equations are said to constitute a complete integral. Then ∞$^4$ elements satisfying F＝0 are known, and in fact every other form of integral can be obtained without further integrations.
In the foregoing discussion of the differential equations of a characteristic chain, the denominators dF/dp, ... may be supposed to be modified in form by means of F＝0 in any way conducive to a simple integration. In the immediately following explanation of ideas, however, we consider indifferently all equations F＝constant; when a function of x, y, z, p, q is said to be zero, it is meant that this is so identically, not in virtue of F＝0; in other words, we consider the integration of F＝a, where a is an arbitrary constant. In the theory of linear partial equations we have seen that the integration Operations necessary for integration of F＝a. of the equations of the characteristic chains, from which, as has just been seen, that of the equation F＝a follows at once, would be involved in completely integrating the single linear homogeneous partial differential equation of the first order [Fƒ]＝0 where the notation is that explained above under Contact Transformations. One obvious integral is ƒ＝F. Putting F＝a, where a is arbitrary, and eliminating one of the independent variables, we can reduce this equation [Fƒ]＝0 to one in four variables; and so on. Calling, then, the determination of a single integral of a single homogeneous partial differential equation of the first order in n independent variables, an operation of order n − 1, the characteristic chains, and therefore the most general integral of F＝a, can be obtained by successive operations of orders 3, 2, 1. If, however, an integral of F＝a be represented by F＝a, G＝b, H＝c, where b and c are arbitrary constants, the expression of the fact that a characteristic chain of F＝a satisfies dG＝0, gives [FG]＝0; similarly, [FH]＝0 and [GH]＝0, these three relations being identically true. Conversely, suppose that an integral G, independent of F, has been obtained of the equation [Fƒ]＝0, which is an operation of order three. Then it follows from the identity [ƒ[φψ]] + [φ[ψƒ]] + [ψ[ƒφ]]＝dƒ/dz [ψφ] + dφ/dz [ψƒ] + dψ/dz [ƒφ] before remarked, by putting φ＝F, ψ＝G, and then [Fƒ]＝A(ƒ), [Gƒ]＝B(ƒ), that AB(ƒ) − BA(ƒ)＝dF/dz B(ƒ) − dG/dz A(ƒ), so that the two linear equations [Fƒ]＝0, [Gƒ]＝0 form a complete system; as two integrals F, G are known, they have a common integral H, independent of F, G, determinable by an operation of order one only. The three functions F, G, H thus identically satisfy the relations [FG]＝[GH]＝[FH]＝0. The ∞² elements satisfying F＝a, G＝b, H＝c, wherein a, b, c are assigned constants, can then be seen to constitute an integral of F＝a. For the conditions that a characteristic chain of G＝b issuing from an element satisfying F＝a, G＝b, H＝c should consist only of elements satisfying these three equations are simply [FG]＝0, [GH]＝0. Thus, starting from an arbitrary element of (F＝a, G＝b, H＝c), we can single out a connectivity of elements of (F＝a, G＝b, H＝c) forming a characteristic chain of G＝b; then the aggregate of the characteristic chains of F＝a issuing from the elements of this characteristic chain of G＝b will be a connectivity consisting only of elements of
(F＝a, G＝b, H＝c),
and will therefore constitute an integral of F＝a; further, it will include all elements of (F＝a, G＝b, H＝c). This result follows also from a theorem given under Contact Transformations, which shows, moreover, that though the characteristic chains of F＝a are not determined by the three equations F＝a, G＝b, H＝c, no further integration is now necessary to find them. By this theorem, since identically [FG]＝[GH]＝[FH]＝0, we can find, by the solution of linear algebraic equations only, a non-vanishing function σ and two functions A, C, such that
dG − AdF − CdH＝σ(dz − pdz − qdy);
thus all the elements satisfying F＝a, G＝b, H＝c, satisfy dz＝pdx + qdy and constitute a connectivity, which is therefore an integral of F＝a. While, further, from the associated theorems, F, G, H, A, C are independent functions and [FC]＝0. Thus C may be taken to be the remaining integral independent of G, H, of the equation [Fƒ]＝0, whereby the characteristic chains are entirely determined.
When we consider the particular equation F＝0, neglecting the case when neither p nor q enters, and supposing p to enter, we may express p from F＝0 in terms of x, y, z, q, and then eliminate it from all other equations. Then instead of the equation [Fƒ]＝0, we have, if F＝0 give p＝ψ(x, y, z, q), the equation
$\Omega f=-\left({\frac {df}{dx}}+\psi {\frac {df}{dz}}\right)+{\frac {d\psi }{dq}}\left({\frac {df}{dy}}+q{\frac {df}{dz}}\right)-\left({\frac {d\psi }{dy}}+q{\frac {d\psi }{dz}}\right){\frac {df}{dq}}=0,$
moreover obtainable by omitting the term in dƒ/dp in [p − ψ, ƒ]＝0. Let x$_0$, y$_0$, z$_0$, q$_0$, be values about which the coefficients in this equation are developable, and let ζ, η, ω be the principal solutions reducing respectively to z, y and q when x＝x$_0$. Then the equations p＝ψ, ζ＝z$_0$, η＝y$_0$, ω＝q$_0$ The single equation F＝0 and Pfaffian formulations.represent a characteristic chain issuing from the element x$_0$, y$_0$, z$_0$, ψ$_0$, q$_0$; we have seen that the aggregate of such chains issuing from the elements of an arbitrary chain satisfying
dz$_0$＝p$_0$dx$_0$ − q$_0$dy$_0$＝0
constitute an integral of the equation p＝ψ. Let this arbitrary ​chain be taken so that x$_0$ is constant; then the condition for initial values is only
dz$_0$ − q$_0$dy$_0$＝0,
and the elements of the integral constituted by the characteristic chains issuing therefrom satisfy
dζ − ωdη＝0.
Hence this equation involves dz − ψdx − qdy = 0, or we have
dz − ψdx − qdy＝σ(dζ − ωdη),
where σ is not zero. Conversely, the integration of p = ψ is, essentially, the problem of writing the expression dz − ψdx − qdy in the form σ(dζ − ωdη), as must be possible (from what was said under Pfaffian Expressions).
To integrate a system of simultaneous equations of the first order X$_1$ = a$_1$, ... X$_r$ = a$_r$ in n independent variables x$_1$, ... x$_n$ and one dependent variable z, we write p$_1$ for dz/dx$_1$, &c., and attempt to find n + 1 − r further functions Z, X$_{r+1}$ ... X$_n$, such that the equations Z = a, X$_i$ = a$_i$,(i = 1, ... n) System of equations of the first order.involve dz − p$_1$dx$_1$ − ... − p$_n$dx$_n$ = 0. By an argument already given, the common integral, if existent, must be satisfied by the equations of the characteristic chains of any one equation X$_i$ = a$_i$; thus each of the expressions [X$_i$X$_j$] must vanish in virtue of the equations expressing the integral, and we may without loss of generality assume that each of the corresponding 1/2r(r − 1) expressions formed from the r given differential equations vanishes in virtue of these equations. The determination of the remaining n + 1 − r functions may, as before, be made to depend on characteristic chains, which in this case, however, are manifolds of r dimensions obtained by integrating the equations [X$_1$ƒ] = 0, ... [X$_r$ƒ] = 0; or having obtained one integral of this system other than X$_1$, ... X$_r$, say X$_{r+1}$, we may consider the system [X$_1$ƒ] = 0, ... [X$_{r+1}$ƒ] = 0, for which, again, we have a choice; and at any stage we may use Mayer’s method and reduce the simultaneous linear equations to one equation involving parameters; while if at any stage of the process we find some but not all of the integrals of the simultaneous system, they can be used to simplify the remaining work; this can only be clearly explained in connexion with the theory of so-called function groups for which we have no space. One result arising is that the simultaneous system p$_1$ = φ$_1$, ... p$_r$ = φ$_r$, wherein p$_1$, . . . p$_r$ are not involved in φ$_1$, . . . φ$_r$, if it satisfies the 1/2r(r − 1) relations [p$_i$ − φ$_i$, p$_j$ − φ$_j$] = 0, has a solution z = ψ(x$_1$, ... x$_n$), p$_1$ = dψ/dx$_1$, ... p$_n$ = dψ/dx$_n$, reducing to an arbitrary function of x$_{r+1}$, ... x$_n$ only, when x$_1$ = xº$_1$, ... x$_r$ = xº$_r$ under certain conditions as to developability; a generalization of the theorem for linear equations. The problem of integration of this system is, as before, to put
dz − φ$_1$dx$_1$ − ... − φ$_r$dx$_r$ − p$_{r+1}$dx$_{r+1}$ − ... − p$_n$dx$_n$
into the form σ(dζ − ω$_{r+1}$ + dξ$_{r+1}$ − ... − ω$_n$dξ$_n$); and here ζ, ξ$_{r+1}$, ... ξ$_n$, ω$_{r+1}$, ... ω$_n$ may be taken, as before, to be principal integrals of a certain complete system of linear equations; those, namely, determining the characteristic chains.
If L be a function of t and of the 2n quantities x$_1$, ... x$_n$, ẋ$_1$, ... ẋ$_n$, where ẋ$_i$, denotes dx$_i$/dt, &c., and if in the n equations
${\frac {d}{dt}}\left({\frac {d\mathrm {L} }{d{\dot {x}}_{i}}}\right)={\frac {d\mathrm {L} }{dx_{i}}}$
we put p$_i$ = dL/dẋ$_i$, and so express ẋ$_i$, ... ẋ$_n$ in terms of t, x$_i$, ... x$_n$, p$_1$, ... p$_n$, assuming that the determinant of the quantities d²L/dx$_i$dẋ$_j$ is not zero; if, further, H denote the function of t, x$_1$, ... x$_n$, p$_1$, ... p$_n$, numerically equal to p$_1$ẋ$_1$ + ... + p$_n$ẋ$_n$ − L, it is easy Equations of dynamics. to prove that dp$_i$/dt = −dH/dx$_i$, dx$_i$/dt = dH/dp$_i$. These so-called canonical equations form part of those for the characteristic chains of the single partial equation dz/dt + H(t, x$_1$, ... x$_n$, dz/dx$_1$, ..., dz/dx$_n$) = 0, to which then the solution of the original equations for x$_1$ ... x$_n$ can be reduced. It may be shown (1) that if z = ψ(t, x$_1$, ... x$_n$, c$_1$, .. c$_n$) + c be a complete integral of this equation, then p$_i$ = dψ/dx$_i$, dψ/dc$_i$ = e$_i$ are 2n equations giving the solution of the canonical equations referred to, where c$_1$ ... c$_n$ and e$_1$, ... e$_n$ are arbitrary constants; (2) that if x$_i$ = X$_i$(t, x$^0$$_1$, ... pº$_n$), p$_i$ = P$_i$(t, xº$_1$, ... p$^0$$_n$) be the principal solutions of the canonical equations for t = t$^0$, and ω denote the result of substituting these values in p$_1$dH/dp$_1$ + ... + p$_n$dH/dp$_n$ − H, and Ω = ∫t$_{t0}$ ωdt, where, after integration, Ω is to be expressed as a function of t, x$_1$, ... x$_n$, xº$_1$, ... xº$_n$, then z = Ω + z$^0$ is a complete integral of the partial equation.
A system of differential equations is said to allow a certain continuous group of transformations (see Groups, Theory of) when the introduction for the variables in the differential equations of the new variables given by the equations of the group leads, for all values of the Application of theory of continuous groups to formal theories.parameters of the group, to the same differential equations in the new variables. It would be interesting to verify in examples that this is the case in at least the majority of the differential equations which are known to be integrable in finite terms. We give a theorem of very general application for the case of a simultaneous complete system of linear partial homogeneous differential equations of the first order, to the solution of which the various differential equations discussed have been reduced. It will be enough to consider whether the given differential equations allow the infinitesimal transformations of the group.
It can be shown easily that sufficient conditions in order that a complete system Π$_1$ƒ = 0 ... Π$_k$ƒ = 0, in n independent variables, should allow the infinitesimal transformation Pƒ = 0 are expressed by k equations Π$_i$Pƒ − PΠ$_i$ƒ = λ$_{i1}$Π$_1$ƒ + ... + λ$_{ik}$Π$_k$ƒ. Suppose now a complete system of n − r equations in n variables to allow a group of r infinitesimal transformations (P$_1$f, ..., P$_r$ƒ) which has an invariant subgroup of r − 1 parameters (P$_1$ƒ, ..., P$_{r−1}$ƒ), it being supposed that the n quantities Π$_1$ƒ, ..., Π$_{n-r}$ƒ, P$_1$ƒ, ..., P$_r$ƒ are not connected by an identical linear equation (with coefficients even depending on the independent variables). Then it can be shown that one solution of the complete system is determinable by a quadrature. For each of Π$_i$Pσƒ − PσΠ$_i$f is a linear function of Π$_1$ƒ, ..., Π$_{n-r}$ƒ and the simultaneous system of independent equations Π$_1$ƒ = 0, ... Π$_{n-r}$ƒ = 0, P$_1$ƒ = 0, ... P$_{r−1}$ƒ = 0 is therefore a complete system, allowing the infinitesimal transformation P$_r$ƒ. This complete system of n − 1 equations has therefore one common solution ω, and P$_r$(ω) is a function of ω. By choosing ω suitably, we can then make P$_r$(ω) = 1. From this equation and the n − 1 equations Π$_i$ω = 0, P$_{σω}$ = 0, we can determine ω by a quadrature only. Hence can be deduced a much more general result, that if the group of r parameters be integrable, the complete system can be entirety solved by quadratures; it is only necessary to introduce the solution found by the first quadrature as an independent variable, whereby we obtain a complete system of n − r equations in n − 1 variables, subject to an integrable group of r − 1 parameters, and to continue this process. We give some examples of the application of the theorem. (1) If an equation of the first order y′ = ψ(x, y) allow the infinitesimal transformation ξdƒ/dx + ηdƒ/dy, the integral curves ω(x, y) = y$^0$, wherein ω(x, y) is the solution of dƒ/dx + ψ(x, y) dƒ/dy = 0 reducing to y for x = x$^0$, are interchanged among themselves by the infinitesimal transformation, or ω(x, y) can be chosen to make ξd$_ω$/dx + ηd$_ω$/dy = 1; this, with dω/dx + ψdω/dy = 0, determines ω as the integral of the complete differential (dy − ψdx)/(η − ψξ). This result itself shows that every ordinary differential equation of the first order is subject to an infinite number of infinitesimal transformations. But every infinitesimal transformation ξdƒ/dx + ηdƒ/dy can by change of variables (after integration) be brought to the form dƒ/dy, and all differential equations of the first order allowing this group can then be reduced to the form F(x, dy/dx) = 0. (2) In an ordinary equation of the second order y ″= ψ(x, y, y′), equivalent to dy/dx = y$_1$, dy$_1$/dx = ψ(x, y, y$_1$), if H, H$_1$ be the solutions for y and y$_1$ chosen to reduce to y$^o$ and y$^o$$_1$ when x = x$^o$, and the equations H = y, H$_1$= y$_1$ be equivalent to ω = y$^o$, ω$_1$ = y$^o$$_1$, then ω, ω$_1$ are the principal solutions of Πƒ = dƒ/dx + y$_1$dƒ/dy + ψdƒ/dy$_1$ = 0. If the original equation allow an infinitesimal transformation whose first extended form (see Groups) is Pƒ = ξdƒ/dx + ηdƒ/dy + η$_1$dƒ/dy$_1$, where η$_1$δt is the increment of dy/dx when ξδt, ηδt are the increments of x, y, and is to be expressed in terms of x, y, y$_1$, then each of Pω and Pω$_1$ must be functions of ω and ω$_1$, or the partial differential equation Πƒ must allow the group Pƒ. Thus by our general theorem, if the differential equation allow a group of two parameters (and such a group is always integrable), it can be solved by quadratures, our explanation sufficing, however, only provided the form Πƒ and the two infinitesimal transformations are not linearly connected. It can be shown, from the fact that η$_1$ is a quadratic polynomial in y$_1$, that no differential equation of the second order can allow more than 8 really independent infinitesimal transformations, and that every homogeneous linear differential equation of the second order allows just 8, being in fact reducible to d²y/dx² = 0. Since every group of more than two parameters has subgroups of two parameters, a differential equation of the second order allowing a group of more than two parameters can, as a rule, be solved by quadratures. By transforming the group we see that if a differential equation of the second order allows a single infinitesimal transformation, it can be transformed to the form F(x, dγ/dx, d$^2$γ/dx$^2$); this is not the case for every differential equation of the second order. (3) For an ordinary differential equation of the third order, allowing an integrable group of three parameters whose infinitesimal transformations are not linearly connected with the partial equation to which the solution of the given ordinary equation is reducible, the similar result follows that it can be integrated by quadratures. But if the group of three parameters be simple, this result must be replaced by the statement that the integration is reducible to quadratures and that of a so-called Riccati equation of the first order, of the form dy/dx = A + By + Cy², where A, B, C are functions of x. (4) Similarly for the integration by quadratures of an ordinary equation y$_n$ = ψ(x, y, y$_1$, ... y$_{n−1}$) of any order. Moreover, the group allowed by the equation may quite well consist of extended contact transformations. An important application is to the case where the differential equation is the resolvent equation defining the group of ​transformations or rationality group of another differential equation (see below); in particular, when the rationality group of an ordinary linear differential equation is integrable, the equation can be solved by quadratures.
Following the practical and provisional division of theories of differential equations, to which we alluded at starting, into transformation theories and function theories, we pass now to give some account of the latter. These are both a necessary logical complement of the former, and the Consideration of function theories of differential equations.only remaining resource when the expedients of the former have been exhausted. While in the former investigations we have dealt only with values of the independent variables about which the functions are developable, the leading idea now becomes, as was long ago remarked by G. Green, the consideration of the neighbourhood of the values of the variables for which this developable character ceases. Beginning, as before, with existence theorems applicable for ordinary values of the variables, we are to consider the cases of failure of such theorems.
When in a given set of differential equations the number of equations is greater than the number of dependent variables, the equations cannot be expected to have common solutions unless certain conditions of compatibility, obtainable by equating different forms of the same differential coefficients deducible from the equations, are satisfied. We have had examples in systems of linear equations, and in the case of a set of equations p$_1$ = φ$_1$, . . . , p$_r$ = φ$_r$ . For the case when the number of equations is the same as that of dependent variables, the following is a general theorem which should be referred to: Let there be r equations in r dependent variables z$_1$, . . . z$_r$ and n independent A general existence theorem. variables x$_1$, . . . x$_n$; let the differential coefficient of z$_σ$ of highest order which enters be of order h$_σ$, and suppose d $^{hσ}$z$_σ$ / dx$_1$$^{hσ}$ to enter, so that the equations can be written d $^{hσ}$z$_σ$ / dx$_1$$^{hσ}$ = Φ$_σ$, where in the general differential coefficient of z$_ρ$ which enters in Φ$_σ$, say
d $^{k1+...+kn}$ z$_ρ$ / dx$_1$$^{k1}$ . . . dx$_n$$^{kn}$,
we have k$_1$ < h$_ρ$ and k$_1$ + . . . + k$_n$ ≤ h$_ρ$. Let a$_1$, . . . a$_n$, b$_1$, . . . b$_r$, and b$^ρ$$_{k1...kn}$ be a set of values of
x$_1$, . . . x$_n$, z$_1$, . . . z$_r$
and of the differential coefficients entering in Φ$_σ$ about which all the functions Φ$_1$, . . . Φ$_r$, are developable. Corresponding to each dependent variable z$_σ$, we take now a set of h$_σ$ functions of x$_2$, . . . x$_n$, say φ$_σ$, φ$_σ$;$^{(1)}$, . . . ,φ$_σ$$^{h−1}$ arbitrary save that they must be developable about a$_2$, a$_3$, . . . a$_n$, and such that for these values of x$_2$, . . . x$_n$, the function φ$_ρ$ reduces to b$_ρ$, and the differential coefficient
d $^{k2+...+kn}$ φ$_ρ$$^{k1}$ / dx$_2$$^{k2}$ . . . dx$_n$$^{kn}$
reduces to b$^ρ$$_{k1...kn}$. Then the theorem is that there exists one, and only one, set of functions z$_1$, . . . z$_r$, of x$_2$, . . . x$_n$ developable about a$_1$, . . . a$_n$ satisfying the given differential equations, and such that for x$_1$ = a$_1$ we have
z$_σ$＝φ$_σ$, dz$_σ$ / dx$_1$＝φ$_σ$$^{(1)}$, . . . d $^{hσ−1}$z$_σ$ / d $^{hσ−1}$x$_1$＝φ$_σ$$^{hσ−1}$.
And, moreover, if the arbitrary functions φ$_σ$, φ$_σ$$^{(1)}$ . . . contain a certain number of arbitrary variables t$_1$, . . . t$_m$, and be developable about the values tº$_1$, . . . tº$_m$ of these variables, the solutions z$_1$, . . . z$_r$ will contain t$_1$, . . . t$_m$, and be developable about tº$_1$, . . . tº$_m$.
The proof of this theorem may be given by showing that if ordinary power series in x$_1$ − a$_1$, . . . x$_n$ − a$_n$, t$_1$ − tº$_1$, . . . t$_m$ − tº$_m$ be substituted in the equations wherein in z$_σ$ the coefficients of (x$_1$ − a$_1$)º, x$_1$ − a$_1$, . . ., (x$_1$ − a$_1$)$^{hσ−1}$ are the arbitrary functions φ$_σ$, φ$_σ$$^{(1)}$, . . . φ$_σ$$^{h−1}$, divided respectively by 1, 1!, 2!, &c., then the differential equations determine uniquely all the other coefficients, and that the resulting series are convergent. We rely, in fact, upon the theory of monogenic analytical functions (see Function), a function being determined entirely by its development in the neighbourhood of one set of values of the independent variables, from which all its other values arise by continuation; it being of course understood that the coefficients in the differential equations are to be continued at the same time. But it is to be remarked that there is no ground for believing, if this method of continuation be utilized, that the function is single-valued; we may quite well return to the same values of the independent variables with a different Singular points of solutions. value of the function; belonging, as we say, to a different branch of the function; and there is even no reason for assuming that the number of branches is finite, or that different branches have the same singular points and regions of existence. Moreover, and this is the most difficult consideration of all, all these circumstances may be dependent upon the values supposed given to the arbitrary constants of the integral; in other words, the singular points may be either fixed, being determined by the differential equations themselves, or they may be movable with the variation of the arbitrary constants of integration. Such difficulties arise even in establishing the reversion of an elliptic integral, in solving the equation
(dx/ds)2＝(x − a$_1$)(x − a$_2$)(x − a$_3$)(x − a$_4$);
about an ordinary value the right side is developable; if we put x − a$_1$ = t$_1$$^2$, the right side becomes developable about t$_1$ = 0; if we put x = 1/t, the right side of the changed equation is developable about t = 0; it is quite easy to show that the integral reducing to a definite value x$_0$ for a value s$_0$ is obtainable by a series in integral powers; this, however, must be supplemented by showing that for no value of s does the value of x become entirely undetermined.
These remarks will show the place of the theory now to be sketched of a particular class of ordinary linear homogeneous differential equations whose importance arises from the completeness and generality with which they can be discussed. We have seen that if in the equationsLinear differential equations with rational coefficients.
dy/dx＝y$_1$, dy$_1$/dx＝y$_2$, . . ., dy$_{n−2}$/dx＝y$_{n−1}$, dy$_{n−1}$/dx＝a$_n$y + a$_{n−1}$y$_1$ + . . . + a$_1$y$_{n−1}$,
where a$_1$, a$_2$, . . ., a$_n$ are now to be taken to be rational functions of x, the value x = xº be one for which no one of these rational functions is infinite, and yº, yº$_1$, . . ., yº$_{n−1}$ be quite arbitrary finite values, then the equations are satisfied by
y＝yºu + yº$_1$u$_1$ + . . . + yº$_{n−1}$u$_{n−1}$,
where u, u$_1$, . . ., u$_{n−1}$ are functions of x, independent of yº, . . . yº$_{n−1}$, developable about x = xº; this value of y is such that for x = xº the functions y, y$_1$ . . . y$_{n−1}$ reduce respectively to yº, yº$_1$, . . . yº$_{n−1}$; it can be proved that the region of existence of these series extends within a circle centre xº and radius equal to the distance from xº of the nearest point at which one of a$_1$, . . . a$_n$ becomes infinite. Now consider a region enclosing xº and only one of the places, say Σ, at which one of a$_1$, . . . a$_n$ becomes infinite. When x is made to describe a closed curve in this region, including this point Σ in its interior, it may well happen that the continuations of the functions u, u$_1$, . . ., u$_{n−1}$ give, when we have returned to the point x, values v, v$_1$, . . ., v$_{n−1}$, so that the integral under consideration becomes changed to yº + yº$_1$v$_1$ + . . . + yº$_{n−1}$v$_{n−1}$. At xº let this branch and the corresponding values of y$_1$, . . . y$_{n−1}$ be ηº, ηº$_1$, . . . ηº$_{n−1}$; then, as there is only one series satisfying the equation and reducing to (ηº, ηº$_1$, . . . ηº$_{n−1}$) for x = xº and the coefficients in the differential equation are single-valued functions, we must have ηºu + ηº$_1$u$_1$ + . . . + ηº$_{n−1}$u$_{n−1}$ = yºv + yº$_1$v$_1$ + . . . + yº$_{n−1}$v$_{n−1}$; as this holds for arbitrary values of yº . . . yº$_{n−1}$, upon which u, . . . u$_{n−1}$ and v, . . . v$_{n−1}$ do not depend, it follows that each of v, . . . v$_{n−1}$ is a linear function of u, . . . u$_{n−1}$ with constant coefficients, say v$_i$ = A$_{i1}$u + . . . + A$_{in}$u$_{n−1}$. Then
yºv + . . . + yº$_{n−1}$v$_{n−1}$＝(Σ$_i$ A$_{i1}$ yº$_i$)u + . . . + (Σ$_i$ A$_{in}$ yº$_i$) u$_{n−1}$;
this is equal to μ(yºu + . . . + yº$_{n−1}$u$_{n−1}$) if Σ$_i$ A$_{ir}$ yº$_i$ = μyº$_{r−1}$; eliminating yº . . . yº$_{n−1}$ from these linear equations, we have a determinantal equation of order n for μ; let μ$_1$ be one of its roots; determining the ratios of yº, y$_1$º, . . . yº$_{n−1}$ to satisfy the linear equations, we have thus proved that there exists an integral, H, of the equation, which when continued round the point Σ and back to the starting-point, becomes changed to H$_1$ = μ$_1$H. Let now ξ be the value of x at Σ and r$_1$ one of the values of (1/2πi) log μ$_1$; consider the function (x − ξ)$^{−r1}$H; when x makes a circuit round x = ξ, this becomes changed to
exp (−2πir$_1$) (x − ξ)$^{−r1}$ μH,
that is, is unchanged; thus we may put H = (x − ξ)$^{r1}$φ$_1$, φ$_1$ being a function single-valued for paths in the region considered described about Σ, and therefore, by Laurent’s Theorem (see Function), capable of expression in the annular region about this point by a series of positive and negative integral powers of x − ξ, which in general may contain an infinite number of negative powers; there is, however, no reason to suppose r$_1$ to be an integer, or even real. Thus, if all the roots of the determinantal equation in μ are different, we obtain n integrals of the forms (x − ξ)$^{r1}$φ$_1$, . . ., (x − ξ)$^{rn}$φ$_n$. In general we obtain as many integrals of this form as there are really different roots; and the problem arises to discover, in case a root be k times repeated, k − 1 equations of as simple a form as possible to replace the k − 1 equations of the form yº + . . . + yº$_{n−1}$v$_{n−1}$ = μ(yº + . . . + yº$_{n−1}$u$_{n−1}$) which would have existed had the roots been different. The most natural method of obtaining a suggestion lies probably in remarking that if r$_2$ = r$_1$ + h, there is an integral [(x − ξ)$^{r1+h}$φ$_2$ − (x − ξ)$^{r1}$φ$_1$] / h, where the coefficients in φ$_2$ are ​the same functions of r$_1$ + h as are the coefficients in φ$_1$ of r$_1$; when h vanishes, this integral takes the form
(x − ξ)$^{r1}$ [dφ$_1$/dr$_1$ + φ$_1$ log (x − ξ)],
or say
(x − ξ)$^{r1}$ [φ$_1$ + ψ$_1$ log (x − ξ)];
denoting this by 2πiμ$_1$K, and (x − ξ)$^{r1}$ φ$_1$ by H, a circuit of the point ξ changes K into
$\mathrm {K} {=}{\frac {1}{2\pi i\mu _{1}}}\left[e^{2\pi ir_{1}}(x-\xi )^{r_{1}}\psi _{1}+e^{2\pi ir_{1}}(x-\xi )^{r_{1}}\phi _{1}(2\pi i+\log(x-\xi ))\right]{=}\mu _{1}\mathrm {K} +\mathrm {H} .$
A similar artifice suggests itself when three of the roots of the determinantal equation are the same, and so on. We are thus led to the result, which is justified by an examination of the algebraic conditions, that whatever may be the circumstances as to the roots of the determinantal equation, n integrals exist, breaking up into batches, the values of the constituents H$_1$, H$_2$, ... of a batch after circuit about x = ξ being H$_1$′ = μ$_1$H$_1$, H$_2$′ = μ$_1$H$_2$ + H$_1$, H$_3$′ = μ$_1$H$_3$ + H$_2$, and so on. And this is found to lead to the forms (x − ξ)$^{r1}$φ$_1$, (x − ξ)$^{r1}$ [ψ$_1$ + φ$_1$ log (x − ξ)], (x − ξ)$^{r1}$ [χ$_1$ + χ$_2$ log (x − ξ) + φ$_1$(log(x − ξ) )²], and so on. Here each of φ$_1$, ψ$_1$, χ$_1$, χ$_2$, ... is a series of positive and negative integral powers of x − ξ in which the number of negative powers may be infinite.
It appears natural enough now to inquire whether, under proper conditions for the forms of the rational functions a$_1$, ... a$_n$, it may be possible to ensure that in each of the series φ$_1$, ψ$_1$, [χ]$_1$, ... the number of negative powers shall be finite. Herein Regular equations. lies, in fact, the limitation which experience has shown to be justified by the completeness of the results obtained. Assuming n integrals in which in each of φ$_1$, ψ$_1$, χ$_1$ ... the number of negative powers is finite, there is a definite homogeneous linear differential equation having these integrals; this is found by forming it to have the form
y′ $^n$ = (x − ξ)$^{−1}$ b$_1$y′ $^{(n−1)}$ + (x − ξ)$^{−2}$ b$_2$y′ $^{(n−2)}$ + ... + (x − ξ)$^{−n}$ b$_n$y,
where b$_1$, ... b$_n$ are finite for x = ξ. Conversely, assume the equation to have this form. Then on substituting a series of the form (x − ξ)$^r$ [1 + A$_1$(x − ξ) + A$_2$(x − ξ)² + ... ] and equating the coefficients of like powers of x − ξ, it is found that r must be a root of an algebraic equation of order n; this equation, which we shall call the index equation, can be obtained at once by substituting for y only (x − ξ)$^r$ and replacing each of b$_1$, ... b$_n$ by their values at x = ξ; arrange the roots r$_1$, r$_2$, ... of this equation so that the real part of r$_i$ is equal to, or greater than, the real part of r$_{i+1}$, and take r equal to r$_1$; it is found that the coefficients A$_1$, A$_2$ ... are uniquely determinate, and that the series converges within a circle about x = ξ which includes no other of the points at which the rational functions a$_1$ ... a$_n$ become infinite. We have thus a solution H$_1$ = (x − ξ)$^{r1}$φ$_1$ of the differential equation. If we now substitute in the equation y = H$_1$∫ηdx, it is found to reduce to an equation of order n − 1 for η of the form
η′ $^{(n−1)}$ = (x − ξ)$^{−1}$ c$_1$η′ $^{(n−2)}$ + ... + (x − ξ)$^{(n−1)}$ c$_{n−1}$η,
where c$_1$, ... c$_{n−1}$ are not infinite at x = ξ. To this equation precisely similar reasoning can then be applied; its index equation has in fact the roots r$_2$ − r$_1$ − 1, ..., r$_n$ − r$_1$ − 1; if r$_2$ − r$_1$ be zero, the integral (x − ξ)$^{−1}$ψ$_1$ of the η equation will give an integral of the original equation containing log (x − ξ); if r$_2$ − r$_1$ be an integer, and therefore a negative integer, the same will be true, unless in ψ$_1$ the term in (x − ξ)$^{r1−r2}$ be absent; if neither of these arise, the original equation will have an integral (x − ξ)$^{r2}$φ$_2$. The η equation can now, by means of the one integral of it belonging to the index r$_2$ − r$_1$ − 1, be similarly reduced to one of order n − 2, and so on. The result will be that stated above. We shall say that an equation of the form in question is regular about x = ξ.
We may examine in this way the behaviour of the integrals at all the points at which any one of the rational functions a$_1$ ... a$_n$ becomes infinite; in general we must expect that beside these the value x = ∞ will be a singular point for the Fuchsian equations. solutions of the differential equation. To test this we put x = 1/t throughout, and examine as before at t = 0. For instance, the ordinary linear equation with constant coefficients has no singular point for finite values of x; at x = ∞ it has a singular point and is not regular; or again, Bessel’s equation x²y″ + xy′ + (x² − n²)y = 0 is regular about x = 0, but not about x = ∞. An equation regular at all the finite singularities and also at x = ∞ is called a Fuchsian equation. We proceed to examine particularly the case of an equation of the second order
y″ + ay′ + by = 0.
Putting x = 1/t, it becomes
d²y/dt² + (2t$^{−1}$ − at$^{−2}$) dy/dt + bt$^{−4}$ y = 0,
which is not regular about t = 0 unless 2 − at$^{−1}$ and bt$^{−2}$, that is, unless ax and bx² are finite at x = ∞; which we thus assume; putting y = t$^r$(1 + A$_1$t + ... ), we find for the index equation at x = ∞ the equation r(r − 1) + r(2 − ax)$_0$ + (bx²)$_0$ = 0. If there be Equation of the second order. finite singular points at ξ$_1$, ... ξ$_m$, where we assume m > 1, the cases m = 0, m = 1 being easily dealt with, and if φ(x) = (x − ξ$_1$) ... (x − ξ$_m$), we must have a·φ(x) and b·[φ(x)]² finite for all finite values of x, equal say to the respective polynomials ψ(x) and θ(x), of which by the conditions at x = ∞ the highest respective orders possible are m − 1 and 2(m − 1). The index equation at x = ξ$_1$ is r(r − 1) + rψ(ξ$_1$) / φ′ (ξ$_1$) + θ(ξ)$_1$ / [φ′(ξ$_1$)]² = 0, and if α$_1$, β$_1$ be its roots, we have α$_1$ + β$_1$ = 1 − ψ(ξ$_1$) / φ′ (ξ$_1$) and α$_1$β$_1$ = θ(ξ)$_1$ / [φ′(ξ$_1$)]². Thus by an elementary theorem of algebra, the sum Σ(1 − α$_i$ − β$_i$) / (x − ξ$_i$), extended to the m finite singular points, is equal to ψ(x) / φ(x), and the sum Σ(1 − α$_i$ − β$_i$) is equal to the ratio of the coefficients of the highest powers of x in ψ(x) and φ(x), and therefore equal to 1 + α + β, where α, β are the indices at x = ∞. Further, if (x, 1)$_{m−2}$ denote the integral part of the quotient θ(x) / φ(x), we have Σ α$_i$β$_i$φ′ (ξ$_i$) / (x = ξ$_i$) equal to −(x, 1)$_{m−2}$ + θ(x)/φ(x), and the coefficient of x$^{m−2}$ in (x, 1)$_{m−2}$ is αβ. Thus the differential equation has the form
y″ + y′Σ (1 − α$_i$ − β$_i$) / (x − ξ$_i$) + y[(x, 1)$_{m-2}$ + Σ α$_i$β$_i$φ′(ξ$_i$) / (x − ξ$_i$)]/φ(x) = 0.
If, however, we make a change in the dependent variable, putting y = (x − ξ$_1$)$^{α1}$ ... (x − ξ$_m$)$^{αmη}$, it is easy to see that the equation changes into one having the same singular points about each of which it is regular, and that the indices at x = ξ$_i$ become 0 and β$_i$ − α$_i$, which we shall denote by λ$_i$, for (x − ξ$_i$)$^{αj}$ can be developed in positive integral powers of x − ξ$_i$ about x = ξ$_i$; by this transformation the indices at x = ∞ are changed to
α + α$_1$ + ... + α$_m$, β + β$_1$ + ... + β$_m$
which we shall denote by λ, μ. If we suppose this change to have been introduced, and still denote the independent variable by y, the equation has the form
y″ + y′Σ (1 − λ$_i$) / (x − ξ$_i$) + y(x, 1)$_{m−2}$ / φ(x) = 0,
while λ + μ + λ$_1$ + ... + λ$_m$ = m − 1. Conversely, it is easy to verify that if λμ be the coefficient of x$^{m−2}$ in (x, 1)$_{m−2}$, this equation has the specified singular points and indices whatever be the other coefficients in (x, 1)$_{m−2}$.
Thus we see that (beside the cases m = 0, m = 1) the “Fuchsian equation” of the second order with two finite singular points is distinguished by the fact that it has a definite form when the singular points and the indices are assigned. Hypergeometric equation. In that case, putting (x − ξ$_1$) / (x − ξ$_2$) = t / (t − 1), the singular points are transformed to 0, 1, ∞, and, as is clear, without change of indices. Still denoting the independent variable by x, the equation then has the form
x(1 − x)y″ + y′[1 − λ$_1$ − x(1 + λ + μ)] − λμy = 0,
which is the ordinary hypergeometric equation. Provided none of λ$_1$, λ$_2$, λ − μ be zero or integral about x = 0, it has the solutions
F(λ, μ, 1 − λ$_1$, x), x$^{λ1}$ F(λ + λ$_1$, μ + λ$_1$, 1 + λ$_1$, x);
about x = 1 it has the solutions
F(λ, μ, 1 − λ$_2$, 1 − x), (1 − x)$^{λ2}$ F(λ + λ$_2$, μ + λ$_2$, 1 + λ$_2$, 1 − x),
where λ + μ + λ$_1$ + λ$_2$ = 1; about x = ∞ it has the solutions
x$^{−λ}$ F(λ, λ + λ$_1$, λ − μ + 1, x$^{−1}$), x$^{−μ}$ F(μ, μ + λ$_1$, μ − λ + 1, x$^{−1}$),
where F(α, β, γ, x) is the series
$1+{\frac {\alpha \beta x}{\gamma }}+{\frac {\alpha (\alpha +1)\beta (\beta +1)x^{2}}{1\cdot 2\cdot \gamma (\gamma +1)}}\ldots ,$
which converges when |x| < 1, whatever α, β, γ may be, converges for all values of x for which |x| = 1 provided the real part of γ − α − β < 0 algebraically, and converges for all these values except x = 1 provided the real part of γ − α − β > −1 algebraically.
In accordance with our general theory, logarithms are to be expected in the solution when one of λ$_1$, λ$_2$, λ − μ is zero or integral. Indeed when λ$_1$ is a negative integer, not zero, the second solution about x = 0 would contain vanishing factors in the denominators of its coefficients; in case λ or μ be one of the positive integers 1, 2, ... (−λ$_1$), vanishing factors occur also in the numerators; and then, in fact, the second solution about x = 0 becomes x$^{λ1}$ times an integral polynomial of degree (−λ$_1$) − λ or of degree (−λ$_1$) − μ. But when λ$_1$ is a negative integer including zero, and neither λ nor μ is one of the positive integers 1, 2 ... (−λ$_1$), the second solution about x = 0 involves a term having the factor log x. When λ$_1$ is a positive integer, not zero, the second solution about x = 0 persists as a solution, in accordance with the order of arrangement of the roots of the index equation in our theory; the first solution is then replaced by an integral polynomial of degree -λ or −μ$_1$, when λ or μ is one of the negative integers 0, −1, −2, ..., 1 − λ$_1$, but otherwise contains a logarithm. Similarly for the solutions about x = 1 or x = ∞; it will be seen below how the results are deducible from those for x = 0.
Denote now the solutions about x = 0 by u$_1$, u$_2$; those about x = 1 by v$_1$, v$_2$; and those about x = ∞ by w$_1$, w$_2$; in the region (S$_0$S$_1$) common to the circles S$_0$, S$_1$ of radius 1 whose centres are the points x = 0, x = 1, all the first four are valid, March of the Integral. and there exist equations u$_1$ =Av$_1$ + Bv$_2$, u$_2$ = Cv$_1$ + Dv$_2$ where A, B, C, D are constants; in the region (S$_1$S) lying inside the circle S$_1$ and outside the circle S$_0$, those that are valid are v$_1$, v$_2$, w$_1$, w$_2$, and there exist equations v$_1$ = Pw$_1$ + Qw$_2$, v$_2$ = Rw$_1$ + Tw$_2$, where P, Q, R, T are constants; thus considering any integral whose expression within the circle S$_0$ is au$_1$ + bu$_2$, where a, b are constants, the same integral will be represented within the circle S$_1$ by (aA + bC)v$_1$ + (aB + bD)v$_2$, and outside these circles will be represented by
[aA + bC)P + (aB + bD)R]w$_1$ + [(aA + bC)Q + (aB + bD)T]w$_2$.
A single-valued branch of such integral can be obtained by making a barrier in the plane joining ∞ to 0 and 1 to ∞; for instance, by excluding the consideration of real negative values of x and of real ​positive values greater than 1, and defining the phase of x and x−1 for real values between 0 and 1 as respectively 0 and π.
We can form the Fuchsian equation of the second order with three arbitrary singular points ξ$_1$, ξ$_2$, ξ$_3$, and no singular point at x＝∞, and with respective indices α$_1$, β$_1$, α$_2$, β$_2$, α$_3$, β$_3$ such that α$_1$+β$_1$+α$_2$+β$_2$+α$_3$+β$_3$＝1. This equation can then be Transformation of the equation into itself.transformed into the hypergeometric equation in 24 ways; for out of ξ$_1$, ξ$_2$, ξ$_3$ we can in six ways choose two, say ξ$_1$, ξ$_2$, which are to be transformed respectively into 0 and 1, by (x−ξ$_1$)/(x−ξ$_2$)＝t(t−1); and then there are four possible transformations of the dependent variable which will reduce one of the indices at t＝0 to zero and one of the indices at t＝1 also to zero, namely, we may reduce either α$_1$ or β$_1$ at t＝0, and simultaneously either α$_2$ or β$_2$ at t＝1. Thus the hypergeometric equation itself can be transformed into itself in 24 ways, and from the expression F(λ, μ, 1−λ$_1$, x) which satisfies it follow 23 other forms of solution; they involve four series in each of the arguments, x, x−1, 1/x, 1/(1−x), (x−1)/x, x/(x−1). Five of the 23 solutions agree with the fundamental solutions already described about x＝0, x＝1, x＝∞; and from the principles by which these were obtained it is immediately clear that the 24 forms are, in value, equal in fours.
The quarter periods K, K′ of Jacobi’s theory of elliptic functions, of which K＝∫π/20 (1−h sin $^2$θ)$^{−1/2}$dθ, and K′ is the same function of 1–h, can easily be proved to be the solutions of a hypergeometric equation of which h is the independent variable. When K, K′ are Inversion. Modular functions.regarded as defined in terms of h by the differential equation, the ratio K′/K is an infinitely many valued function of h. But it is remarkable that Jacobi’s own theory of theta functions leads to an expression for h in terms of K′/K (see Function) in terms of single-valued functions. We may then attempt to investigate, in general, in what cases the independent variable x of a hypergeometric equation is a single-valued function of the ratio ς of two independent integrals of the equation. The same inquiry is suggested by the problem of ascertaining in what cases the hypergeometric series F(α, β, γ, x) is the expansion of an algebraic (irrational) function of x. In order to explain the meaning of the question, suppose that the plane of x is divided along the real axis from −∞ to 0 and from 1 to +∞, and, supposing logarithms not to enter about x＝0, choose two quite definite integrals y$_1$, y$_2$ of the equation, say
y$_1$＝F(λ, μ, 1−λ$_1$, x), y$_2$＝x$^{λ1}$ F(λ+λ$_1$, μ+λ$_1$, 1+λ$_1$, x),
with the condition that the phase of x is zero when x is real and between 0 and 1. Then the value of ς＝y$_2$/y$_1$ is definite for all values of x in the divided plane, ς being a single-valued monogenic branch of an analytical function existing and without singularities all over this region. If, now, the values of ς that so arise be plotted on to another plane, a value p+iq of ς being represented by a point (p, q) of this ς-plane, and the value of x from which it arose being mentally associated with this point of the σ-plane, these points will fill a connected region therein, with a continuous boundary formed of four portions corresponding to the two sides of the two barriers of the x-plane. The question is then, firstly, whether the same value of ς can arise for two different values of x, that is, whether the same point (p, q) of the ς-plane can arise twice, or in other words, whether the region of the ς-plane overlaps itself or not. Supposing this is not so, a second part of the question presents itself. If in the x-plane the barrier joining −∞ to 0 be momentarily removed, and x describe a small circle with centre at x＝0 starting from a point x＝−h−ik, where h, k are small, real, and positive and coming back to this point, the original value ς at this point will be changed to a value σ, which in the original case did not arise for this value of x, and possibly not at all. If, now, after restoring the barrier the values arising by continuation from σ be similarly plotted on the ς-plane, we shall again obtain a region which, while not overlapping itself, may quite possibly overlap the former region. In that case two values of x would arise for the same value or values of the quotient y$_2$/y$_1$, arising from two different branches of this quotient. We shall understand then, by the condition that x is to be a single-valued function of x, that the region in the ς-plane corresponding to any branch is not to overlap itself, and that no two of the regions corresponding to the different branches are to overlap. Now in describing the circle about x＝0 from x＝−h−ik to −h+ik, where h is small and k evanescent,
ς＝x$^{λ1}$ F(λ＋λ$_1$, μ＋λ$_1$, 1+λ$_1$, x)/F(λ, μ, 1−λ$_1$, x)
is changed to σ＝ςe$^{2πiλ1}$. Thus the two portions of boundary of the ς-region corresponding to the two sides of the barrier (−∞, 0) meet (at ς＝0 if the real part of λ$_1$ be positive) at an angle 2πL$_1$, where L$_1$ is the absolute value of the real part of λ$_1$; the same is true for the σ-region representing the branch σ. The condition that the ς-region shall not overlap itself requires, then, L$_1$＝1. But, further, we may form an infinite number of branches σ＝ςe$^{2πiλ1}$, σ$_1$＝e$^{2πiλ1}$, . . . in the same way, and the corresponding regions in the plane upon which y$_2$/y$_1$ is represented will have a common point and each have an angle 2πL$_1$; if neither overlaps the preceding, it will happen, if L$_1$ is not zero, that at length one is reached overlapping the first, unless for some positive integer α we have 2παL$_1$＝2π, in other words L$_1$＝1/α. If this be so, the branch σ$_{α−1}$＝ςe$^{2πiαλ1}$ will be represented by a region having the angle at the common point common with the region for the branch ς; but not altogether coinciding with this last region unless λ$_1$ be real, and therefore＝±1/α; then there is only a finite number, α, of branches obtainable in this way by crossing the barrier (−∞, 0). In precisely the same way, if we had begun by taking the quotient
ς′＝(x−1)$^{λ2}$ F(λ＋λ$_2$, μ＋λ$_2$, 1+λ$_2$, 1−x) /F(λ, μ, 1−λ$_2$, 1−x)
of the two solutions about x＝1, we should have found that x is not a single-valued function of ς′ unless λ$_2$ is the inverse of an integer, or is zero; as ς′ is of the form (Aς+B)/(Cς+D), A, B, C, D constants, the same is true in our case; equally, by considering the integrals about x＝∞ we find, as a third condition necessary in order that x may be a single-valued function of ς, that λ−μ must be the inverse of an integer or be zero. These three differences of the indices, namely, λ$_1$, λ$_2$, λ−μ, are the quantities which enter in the differential equation satisfied by x as a function of ς, which is easily found to be
| − | x$_{111}$ | + | 3x$^2$$_{11}$ | ＝1/2(h−h$_1$−h$_2$)x$^{−1}$(x−1)$^{−1}$+1/2h$_1$x$^{−2}$+1/2h$_2$(x−1)$^{−2}$, |
| − | x$_1$$^3$ | + | 2x$_1$$^4$ | ＝1/2(h−h$_1$−h$_2$)x$^{−1}$(x−1)$^{−1}$+1/2h$_1$x$^{−2}$+1/2h$_2$(x−1)$^{−2}$, |
where x$_1$＝dx/dς, &c.; and h$_1$＝1−y$_1$$^2$, h$_2$＝1−λ$_2$$^2$, h$_3$＝1−(λ−μ)$^2$. Into the converse question whether the three conditions are sufficient to ensure (1) that the ς region corresponding to any branch does not overlap itself, (2) that no two such regions overlap, we have no space to enter. The second question clearly requires the inquiry whether the group (that is, the monodromy group) of the differential equation is properly discontinuous. (See Groups, Theory of.)
The foregoing account will give an idea of the nature of the function theories of differential equations; it appears essential not to exclude some explanation of a theory intimately related both to such theories and to transformation theories, which is a generalization of Galois’s theory of algebraic equations. We deal only with the application to homogeneous linear differential equations.
In general a function of variables x$_1$, x$_2$ . . . is said to be rational when it can be formed from them and the integers 1, 2, 3, . . . by a finite number of additions, subtractions, multiplications and divisions. We generalize this definition. Assume that we have assigned a fundamental series of quantities and functions Rationality group of a linear equation.of x, in which x itself is included, such that all quantities formed by a finite number of additions, subtractions, multiplications, divisions and differentiations in regard to x, of the terms of this series, are themselves members of this series. Then the quantities of this series, and only these, are called rational. By a rational function of quantities p, q, r, . . . is meant a function formed from them and any of the fundamental rational quantities by a finite number of the five fundamental operations. Thus it is a function which would be called, simply, rational if the fundamental series were widened by the addition to it of the quantities p, q, r, . . . and those derivable from them by the five fundamental operations. A rational ordinary differential equation, with x as independent and y as dependent variable, is then one which equates to zero a rational function of y, the order k of the differential equation being that of the highest differential coefficient y$^{(k)}$ which enters; only such equations are here discussed. Such an equation P＝0 is called irreducible when, firstly, being arranged as an integral polynomial in y$^{(k)}$, this polynomial Irreducibility of a rational equation.is not the product of other polynomials in y$^{(k)}$ also of rational form; and, secondly, the equation has no solution satisfying also a rational equation of lower order. From this it follows that if an irreducible equation P＝0 have one solution satisfying another rational equation Q＝0 of the same or higher order, then all the solutions of P＝0 also satisfy Q＝0. For from the equation P＝0 we can by differentiation express y$^{(k+1)}$, y$^{(k+2)}$, . . . in terms of x, y, y$^{(1)}$, . . . , y$^{(k)}$, and so put the function Q rationally in terms of these quantities only. It is sufficient, then, to prove the result when the equation Q＝0 is of the same order as P＝0. Let both the equations be arranged as integral polynomials in y$^{(k)}$; their algebraic eliminant in regard to y$^{(k)}$ must then vanish identically, for they are known to have one common solution not satisfying an equation of lower order; thus the equation P＝0 involves Q＝0 for all solutions of P＝0.
Now let y$^{(n)}$＝a$_1$y$^{(n−1)}$+ . . . +a$_n$y be a given rational homogeneous linear differential equation; let y$_1$, . . . y$_n$ be n particular functions of x, unconnected by any equation with constant coefficients of the form c$_1$y$_1$+ . . . +c$_n$y$_n$＝0, all satisfying the differential equation; let η$_1$, . . . η$_n$ be linear functions The variant function for a linear equation.of y$_1$, . . . y$_n$, say η$_i$＝A$_{i1}$y$_1$+ . . . +A$_{in}$y$_n$, where the constant coefficients A$_{ij}$ have a non-vanishing determinant; write (η)＝A(y), these being the equations of a general linear homogeneous group whose transformations may be denoted by A, B, . . . . We desire to form a rational function φ(η), or say φ(A(y)), of η$_1$, . . . η, in which the η$^2$ constants A$_{ij}$ shall all be essential, and not reduce effectively to a fewer number, as they would, for instance, if the y$_1$, . . . y$_n$ were connected by a linear equation with constant coefficients. Such a function is in fact given, if the solutions y$_1$, . . . y$_n$ be developable ​in positive integral powers about x = a, by φ(η) = η$_1$ + (x − a)$^n$ η$_2$ + . . . + (x − a)$^{(n−1)n}$ η$_n$. Such a function, V, we call a variant.
Then differentiating V in regard to x, and replacing η$_i$$^{(n)}$ by its value a$_1$η$^{(n−1)}$ + . . . + a$_n$η, we can arrange dV/dx, and similarly each of d $^2$/dx$^2$ . . . d $^N$V/dx$^N$, where N = n$^2$, as a linear function of the N quantities η$_1$, . . . η$_n$, . . . η$_1$$^{(n−1)}$, . . . η$_n$$^{(n−1)}$, and thence by elimination obtain a linear differential equation The resolvent equation.for V of order N with rational coefficients. This we denote by F = 0. Further, each of η$_1$ . . . η$_n$ is expressible as a linear function of V, dV/dx, . . . d $^{N−1}$V / dx$^{N−1}$, with rational coefficients not involving any of the n$^2$ coefficients A$_{ij}$, since otherwise V would satisfy a linear equation of order less than N, which is impossible, as it involves (linearly) the n$^2$ arbitrary coefficients A$_{ij}$, which would not enter into the coefficients of the supposed equation. In particular, y$_1$,.. y$_n$ are expressible rationally as linear functions of ω, dω/dx, . . . d $^{N−1}$ω / dx$^{N−1}$, where ω is the particular function φ(y). Any solution W of the equation F = 0 is derivable from functions ζ$_1$, . . . ζ$_n$, which are linear functions of y$_1$, . . . y$_n$, just as V was derived from η$_1$, . . . η$_n$; but it does not follow that these functions ζ$_i$, . . . ζ$_n$ are obtained from y$_1$, . . . y$_n$ by a transformation of the linear group A, B, . . . ; for it may happen that the determinant d(ζ$_1$, . . . ζ$_n$) / (dy$_1$, . . . y$_n$) is zero. In that case ζ$_1$, . . . ζ$_n$ may be called a singular set, and W a singular solution; it satisfies an equation of lower than the N-th order. But every solution V, W, ordinary or singular, of the equation F = 0, is expressible rationally in terms of ω, dω / dx, . . . d $^{N−1}$ω / dx$^{N−1}$; we shall write, simply, V = r(ω). Consider now the rational irreducible equation of lowest order, not necessarily a linear equation, which is satisfied by ω; as y$_1$, . . . y$_n$ are particular functions, it may quite well be of order less than N; we call it the resolvent equation, suppose it of order p, and denote it by γ(v). Upon it the whole theory turns. In the first place, as γ(v) = 0 is satisfied by the solution ω of F = 0, all the solutions of γ(v) are solutions F = 0, and are therefore rationally expressible by ω; any one may then be denoted by r(ω). If this solution of F = 0 be not singular, it corresponds to a transformation A of the linear group (A, B, ...), effected upon y$_1$, . . . y$_n$. The coefficients A$_{ij}$ of this transformation follow from the expressions before mentioned for η$_1$ . . . η$_n$ in terms of V, dV/dx, d $^2$V/dx$^2$, . . . by substituting V = r(ω); thus they depend on the p arbitrary parameters which enter into the general expression for the integral of the equation γ(v) = 0. Without going into further details, it is then clear enough that the resolvent equation, being irreducible and such that any solution is expressible rationally, with p parameters, in terms of the solution ω, enables us to define a linear homogeneous group of transformations of y$_1$ . . . y$_n$ depending on p parameters; and every operation of this (continuous) group corresponds to a rational transformation of the solution of the resolvent equation. This is the group called the rationality group, or the group of transformations of the original homogeneous linear differential equation.
The group must not be confounded with a subgroup of itself, the monodromy group of the equation, often called simply the group of the equation, which is a set of transformations, not depending on arbitrary variable parameters, arising for one particular fundamental set of solutions of the linear equation (see Groups, Theory of).
The importance of the rationality group consists in three propositions. (1) Any rational function of y$_1$, . . . y$_n$ which is unaltered in value by the transformations of the group can be written in rational form. (2) If any rational function be changed in form, becoming a rational function of y$_1$, . . . y$_n$, a The fundamental theorem in regard to the rationality group.transformation of the group applied to its new form will leave its value unaltered. (3) Any homogeneous linear transformation leaving unaltered the value of every rational function of y$_1$, . . . y$_n$ which has a rational value, belongs to the group. It follows from these that any group of linear homogeneous transformations having the properties (1) (2) is identical with the group in question. It is clear that with these properties the group must be of the greatest importance in attempting to discover what functions of x must be regarded as rational in order that the values of y$_1$ . . . y$_n$ may be expressed. And this is the problem of solving the equation from another point of view.
Literature.—(α) Formal or Transformation Theories for Equations of the First Order:—E. Goursat, Leçons sur l’intégration des équations aux dérivées partielles du premier ordre (Paris, 1891); E. v. Weber, Vorlesungen über das Pfaff’sche Problem und die Theorie der partiellen Differentialgleichungen erster Ordnung (Leipzig, 1900); S. Lie und G. Scheffers, Geometrie der Berührungstransformationen, Bd. i. (Leipzig, 1896); Forsyth, Theory of Differential Equations, Part i., Exact Equations and Pfaff’s Problem (Cambridge, 1890); S. Lie, “Allgemeine Untersuchungen über Differentialgleichungen, die eine continuirliche endliche Gruppe gestatten” (Memoir), Mathem. Annal.xxv. (1885), pp. 71-151; S. Lie und G. Scheffers, Vorlesungen über Differentialgleichungen mit bekannten infinitesimalen Transformationen (Leipzig, 1891). A very full bibliography is given in the book of E. v. Weber referred to; those here named are perhaps sufficiently representative of modern works. Of classical works may be named: Jacobi, Vorlesungen über Dynamik (von A. Clebsch, Berlin, 1866); Werke, Supplementband; G Monge, Application de l’analyse à la géométrie (par M. Liouville, Paris, 1850); J. L. Lagrange, Leçons sur le calcul des fonctions (Paris, 1806), and Théorie des fonctions analytiques (Paris, Prairial, an V); G. Boole, A Treatise on Differential Equations (London, 1859); and Supplementary Volume (London, 1865); Darboux, Leçons sur la théorie générale des surfaces, tt. i.-iv. (Paris, 1887–1896); S. Lie, Théorie der transformationsgruppen ii. (on Contact Transformations) (Leipzig, 1890).
(β) Quantitative or Function Theories for Linear Equations:—C. Jordan, Cours d’analyse, t. iii. (Paris, 1896); E. Picard, Traité d’analyse, tt. ii. and iii. (Paris, 1893, 1896); Fuchs, Various Memoirs, beginning with that in Crelle’s Journal, Bd. lxvi. p. 121; Riemann, Werke, 2$^r$ Aufl. (1892); Schlesinger, Handbuch der Theorie der linearen Differentialgleichungen, Bde. i.-ii. (Leipzig, 1895–1898); Heffter, Einleitung in die Theorie der linearen Differentialgleichungen mit einer unabhängigen Variablen (Leipzig, 1894); Klein, Vorlesungen über lineare Differentialgleichungen der zweiten Ordnung (Autographed, Göttingen, 1894); and Vorlesungen über die hypergeometrische Function (Autographed, Göttingen, 1894); Forsyth, Theory of Differential Equations, Linear Equations.
(γ) Rationality Group (of Linear Differential Equations):—Picard, Traité d’Analyse, as above, t. iii.; Vessiot, Annales de l’École Normale, série III. t. ix. p. 199 (Memoir); S. Lie, Transformationsgruppen, as above, iii. A connected account is given in Schlesinger, as above, Bd. ii., erstes Theil.
(δ) Function Theories of Non-Linear Ordinary Equations:—Painlevé, Leçons sur la théorie analytique des équations différentielles (Paris, 1897, Autographed); Forsyth, Theory of Differential Equations, Part ii., Ordinary Equations not Linear (two volumes, ii. and iii.) (Cambridge, 1900); Königsberger, Lehrbuch der Theorie der Differentialgleichungen (Leipzig, 1889); Painlevé, Leçons sur l’intégration des équations différentielles de la mécanique et applications (Paris, 1895).
(ε) Formal Theories of Partial Equations of the Second and Higher Orders:—E. Goursat, Leçons sur l’intégration des équations aux dérivées partielles du second ordre, tt. i. and ii. (Paris, 1896, 1898); Forsyth, Treatise on Differential Equations (London, 1889); and Phil. Trans. Roy. Soc. (A.), vol. cxci. (1898), pp. 1-86.
(ζ) See also the six extensive articles in the second volume of the German Encyclopaedia of Mathematics. (H. F. Ba.)