dxi/dt = φi(tx1, ... xn).
Suppose we have k homogeneous linear partial equations of the first order in n independent variables, the general equation being aσ1dƒ/dx1 + ... + aσndƒ/dxn = 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, ΣHidƒ/dxi = 0 where Hi = 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ƒ/dx1 = 0 ... dƒ/dxn = 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 x1 ... xr; also that the coefficients are all developable about the values x1 = xº1, ... xn= 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 xr+1, ... xn when in them for x1, ... xr are respectively put xº1, ... xºr; so that also the equations have in common a solution reducing when x1 = xº1, ... xr = xºr to an arbitrary function ψ(xr+1, ... xn) 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 x1, ... xn 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ƒ/dx1, ..., dƒ/dxr; suppose the general equation of the new system to be
Qσƒ = dƒ/dxσ + cσjr+1dƒ/dxr+1 + ... + cσndƒ/dxn = 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 xr+1 ... xn; as it is at most a linear function of Q1ƒ, ... Qrƒ, it must be identically zero. So reduced the system is called a Jacobian system. Of this system Q1ƒ=0 has n − 1 principal solutions reducing respectively Jacobian systems. to x2, ... xn when
x1 = xº1,
and its form shows that of these the first r − 1 are exactly x2 ... xr. Let these n − 1 functions together with x1 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ƒ/dx1 = 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, dxj = c1jdx1 + ... + crjdxr, equivalent to System of total differential equations. the r(n − r) equations dxj/dxσ = cσr. That consistent with them we may be able to regard xr+1, ... xn as functions of x1, ... xr, 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 x1, ... xn on this hypothesis, in regard to xρj is, however,
dƒ/dxρ + cρjr+1dƒ/dxr+1 + ... + cρndƒ/dxn,
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 xr+1, ... xn when x1 = xº1, ... xr = xºr, and the equations ωr+1 = x0r+1, ... ωn = xºn be solved for xr+1, ... xn to give xj = ψj(x1, ... xr, x0r+1, ... xºn), these values solve the total equations and reduce respectively to x0r+1, ... xºn when x1 = xº1 ... xr = 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 (x0, y0, z0) of three-dimensioned space (about which a, b are developable) a plane, namely, z − z0 = a0(x − x0) + b0(y − y0), 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 (x0, y0, z0), this plane will touch the surface, and the operations of passing along the surface from (x0, y0, z0) to