Bourbaki

§ 2. Linear differential equations

1. Existence of integrals of a linear differential equation

Let E be a complete normed space over the field R\mathbf{R}, and J an interval in R\mathbf{R}, not reducing to a point. One says that the differential equation

dxdt=f(t,x) \frac{d\mathbf{x}}{dt} = \mathbf{f}(t, \mathbf{x})

(1) where ff is defined on J×EJ \times E, is a linear equation if, for every t∈Jt \in J, the map x↦f(t,x)x \mapsto f(t, x) is a continuous affine linear map [^1] from EE into itself; if one puts b(t)=f(t,0)b(t) = f(t, 0) the map x↦f(t,x)−f(t,0)=f(t,x)−b(t)x \mapsto f(t, x) - f(t, 0) = f(t, x) - b(t) is then a continuous linear map from EE to itself; from now on we shall denote this map by A(t)A(t) and write A(t).xA(t).x, (or simply A(t)xA(t)x) for its value at the point x∈Ex \in E; thus the linear differential equation (1) may be written

dxdt=A(t).x+b(t) \frac{dx}{dt} = A(t).x + b(t)

(2)

where bb is a map from JJ into EE; when b=0b = 0 one says that the linear differential equation (2) is homogeneous.

Example 1

When EE is of finite dimension nn over R\mathbf{R} one can identify the endomorphism A(t)A(t) with its matrix (aij(t))\left(a_{ij}(t)\right) with respect to any basis of EE (Alg., II, p. 343); when one identifies a vector x∈Ex \in E with the column matrix (xj)(x_j) of its components with respect to the basis of EE under consideration the expression A(t).xA(t).x conforms to the general conventions of Algebra (Alg., II, p. 343, prop 2). In this case, equation (2) is equivalent to the system of scalar differential equations

dxidt=∑j=1naij(t)xj+bi(t)(1≤i≤n). \frac{dx_i}{dt} = \sum_{j=1}^n a_{ij}(t)x_j + b_i(t) \qquad (1 \leq i \leq n).

(3)

Example 2

Let GG be a complete normed algebra over R\mathbf{R}, and a(t)a(t), b(t)b(t) and c(t)c(t) three maps from JJ into GG; the equation

dxdt=a(t)x+xb(t)+c(t) \frac{dx}{dt} = a(t)x + x b(t) + c(t)

is a linear differential equation; here A(t)A(t) is the linear map x↦a(t)x+xb(t)x \mapsto a(t)x + x b(t) of GG to itself.

For every t∈Jt \in J, A(t)A(t) is an element of the set L(E)\mathscr{L}(E) of continuous linear maps from EE to itself (continuous endomorphisms of EE); one knows (Gen. Top., X, p. 298) that L(E)\mathscr{L}(E), endowed with the norm ∥U∥=sup⁡∥x∥≤1∥Ux∥\|U\| = \sup_{\|x\| \leq 1} \|Ux\| is a complete normed algebra over the field R\mathbf{R}, and that ∥UV∥≤∥U∥ ∥V∥\|UV\| \leq \|U\|\ \|V\|.

Throughout this section we shall assume that the following conditions are satisfied:

a) The map t↦A(t)t \mapsto A(t) of JJ into L(E)\mathscr{L}(E) is regulated.

b) The map t↦b(t)t \mapsto b(t) of JJ into EE is regulated.

When EE has dimension nn, L(E)\mathscr{L}(E) is isomorphic to Rn2\mathbf{R}^{n^2} (as a topological vector space) and condition a) means that each of the elements aij(t)a_{ij}(t) of the matrix A(t)A(t) is a regulated function on JJ.

Since ∥A(t′)x−A(t)x∥≤∥A(t′)−A(t)∥∥x∥\|A(t')\mathbf{x} - A(t)\mathbf{x}\| \leq \|A(t') - A(t)\|\|\mathbf{x}\|, the map

t↦A(t).x+b(t) t \mapsto A(t).\mathbf{x} + \mathbf{b}(t)
is regulated for every x∈E\mathbf{x} \in E; further,
∥A(t)x1−A(t)x2∥=∥A(t)(x1−x2)∥≤∥A(t)∥∥x1−x2∥ \|A(t)\mathbf{x}_1 - A(t)\mathbf{x}_2\| = \|A(t)(\mathbf{x}_1 - \mathbf{x}_2)\| \leq \|A(t)\|\|\mathbf{x}_1 - \mathbf{x}_2\|
for any t∈Jt \in J and x1,x2\mathbf{x}_1, \mathbf{x}_2 in EE; in other words, the right-hand side of (2) satisfies the conditions of lemma 1 of IV, p. 165 and is Lipschitz with respect to the regulated function ∥A(t)∥\|A(t)\| on J×EJ \times E. In consequence (IV, p. 173, cor. 2):

Theorem 1

Let* t↦A(t)t \mapsto A(t) be a regulated map of JJ into L(E)\mathscr{L}(E), and t↦b(t)t \mapsto \mathbf{b}(t) be a regulated map of JJ into EE. For every point (t0,x0)(t_0, \mathbf{x}_0) of J×EJ \times E the linear equation (2) admits one and only one solution defined on all of JJ and equal to x0\mathbf{x}_0 at the point t0t_0.

2. Linearity of the integrals of a linear differential equation

Solving a linear differential equation (2) is a linear problem (Alg., II, p. 240); the homogeneous linear equation

dxdt=A(t).x \frac{d\mathbf{x}}{dt} = A(t).\mathbf{x}
is said to be associated with the inhomogeneous equation (2); and one knows (Alg., II, p. 241, prop. 14) that if u1\mathbf{u}_1 is an integral of the inhomogeneous equation (2) then every integral of this equation is of the form u+u1\mathbf{u} + \mathbf{u}_1 where u1\mathbf{u}_1 is a solution of the associated homogeneous equation (4), and conversely. We shall first study in this subsection the integrals of a homogeneous equation (4).

Proposition 1

The set* I\mathscr{I} of integrals of the homogeneous linear equation (4), defined on JJ, is a vector subspace of the space C(J;E)C(J; E) of continuous maps from JJ into EE.

The proof is immediate.

Theorem 2

For every point* (t0,x0)(t_0, \mathbf{x}_0) of J×EJ \times E let u(t,t0,x0)\mathbf{u}(t, t_0, \mathbf{x}_0) be the integral of the homogeneous equation (4) defined on JJ and equal to x0\mathbf{x}_0 at t0t_0. 1. For every point t∈Jt \in J the map x0↦u(t,t0,x0)\mathbf{x}_0 \mapsto \mathbf{u}(t, t_0, \mathbf{x}_0) is a bijective bicontinuous linear map C(t,t0)C(t, t_0) of EE to itself. 2. The map t↦C(t,t0)t \mapsto C(t, t_0) of JJ into L(E)\mathscr{L}(E) is identical to the integral of the homogeneous linear differential equation

dUdt=A(t) U \frac{dU}{dt} = A(t)\ U
which takes the value II **(the identity map of* EE to itself) at the point t0t_0. 3. For any points s,t,us, t, u of JJ
C(s,u)=C(s,t)C(t,u),C(s,t)=(C(t,s))−1. C(s, u) = C(s, t)C(t, u), \qquad C(s, t) = (C(t, s))^{-1}.

By prop. 1, u(t,t0,x1)+u(t,t0,x2)\mathbf{u}(t, t_0, \mathbf{x}_1) + \mathbf{u}(t, t_0, \mathbf{x}_2) (resp. λu(t,t0,x0)\lambda \mathbf{u}(t, t_0, \mathbf{x}_0)) is an integral of (4) and takes the value x1+x2\mathbf{x}_1 + \mathbf{x}_2 (resp. λx0\lambda \mathbf{x}_0) at t0t_0, so, by th. 1 of IV, p. 179, is identical to u(t,t0,x1+x2)\mathbf{u}(t, t_0, \mathbf{x}_1 + \mathbf{x}_2) (resp. u(t,t0,λx0)\mathbf{u}(t, t_0, \lambda \mathbf{x}_0)); the map x0↦u(t,t0,x0)\mathbf{x}_0 \mapsto \mathbf{u}(t, t_0, \mathbf{x}_0) is thus a linear map C(t,t0)C(t, t_0) of E into itself, and one can write u(t,t0,x0)=C(t,t0).x0\mathbf{u}(t, t_0, \mathbf{x}_0) = C(t, t_0). \mathbf{x}_0.

Since the map (X,Y)↦XY(X, Y) \mapsto XY of L(E)×L(E)\mathscr{L}(E) \times \mathscr{L}(E) into L(E)\mathscr{L}(E) is continuous (Gen. Top., X, p. 298, prop. 8), the map t↦A(t)Ut \mapsto A(t)U of J into L(E)\mathscr{L}(E) is regulated for all U∈L(E)U \in \mathscr{L}(E); further (Gen. Top., X, p. 296)

∥A(t)X−A(t)Y∥=∥A(t)(X−Y)∥≤∥A(t)∥∥X−Y∥, \| A(t)X - A(t)Y \| = \| A(t)(X - Y) \| \leq \| A(t) \| \| X - Y \|,

so one can apply th. I of IV, p. 179 to the homogeneous linear equation (5); let V(t)V(t) be the integral of this equation defined on J and equal to II at t0t_0. One has (I, p. 6, prop. 3)

ddt(V(t)x0)=dV(t)dtx0=A(t)(V(t)x0) \frac{d}{dt} (V(t)\mathbf{x}_0) = \frac{dV(t)}{dt} \mathbf{x}_0 = A(t)(V(t)\mathbf{x}_0)

and for t=t0t = t_0 we have V(t)x0=Ix0=x0V(t)\mathbf{x}_0 = I\mathbf{x}_0 = \mathbf{x}_0; by th. 1 of IV, p. 179 one must have V(t).x0=C(t,t0)x0V(t).\mathbf{x}_0 = C(t, t_0)\mathbf{x}_0 for all x0∈E\mathbf{x}_0 \in E, that is, V(t)=C(t,t0)V(t) = C(t, t_0); this proves that C(t,t0)C(t, t_0) belongs to L(E)\mathscr{L}(E), in other words, that x0↦C(t,t0).x0\mathbf{x}_0 \mapsto C(t, t_0).\mathbf{x}_0 is continuous on E, and that the map t↦C(t,t0)t \mapsto C(t, t_0) is the integral of (5) which is equal to II at t0t_0.

Finally, the integral s↦C(s,u).x0s \mapsto C(s, u).\mathbf{x}_0 of (4) is equal to C(t,u).x0C(t, u).\mathbf{x}_0 at the point tt, so, by definition,

C(s,u).x0=C(s,t)(C(t,u).x0)=(C(s,t)C(t,u)).x0 C(s, u).\mathbf{x}_0 = C(s, t)(C(t, u).\mathbf{x}_0) = (C(s, t)C(t, u)).\mathbf{x}_0

for any x0∈E\mathbf{x}_0 \in E, whence the first relation (6); since C(s,s)=IC(s, s) = I one has C(s,t)C(t,s)=IC(s, t)C(t, s) = I, for any ss and tt in J; this proves (Set Theory, II, p. 86, corollary) that C(t,t0)C(t, t_0) is a bijective map of E onto itself, with inverse map C(t0,t)C(t_0, t). This completes the proof of the theorem.

One says that C(t,t0)C(t, t_0) is the resolvent of equation (2) of IV, p. 178.

Corollary 1

The map which to every point x0∈E\mathbf{x}_0 \in E associates the continuous function t↦C(t,t0).x0t \mapsto C(t, t_0).\mathbf{x}_0, defined on J, is an isomorphism of the normed space E onto the vector space I\mathscr{I} of integrals of (4), endowed with the topology of compact convergence.

It is certainly a bijective linear map of E onto I\mathscr{I}: now C(t,t0)C(t, t_0) is bounded on a compact set K⊂JK \subset J, so ∥C(t,t0).x0∥≤M∥x0∥\| C(t, t_0).\mathbf{x}_0 \| \leq M \| \mathbf{x}_0 \| for any t∈Kt \in K and x0∈E\mathbf{x}_0 \in E, which shows that this map is continuous; and since

C(t0,t0).x0=x0, C(t_0, t_0).\mathbf{x}_0 = \mathbf{x}_0,

it is clear that the inverse map is also continuous.

Corollary 2

The map (s,t)↦C(s,t)(s, t) \mapsto C(s, t) of J×JJ \times J into L(E)\mathscr{L}(E) is continuous.

By (6) we have C(s,t)=C(s,t0)(C(t,t0))−1C(s, t) = C(s, t_0) \left( C(t, t_0) \right)^{-1}; now, the map (X,Y)↦XY(X, Y) \mapsto XY of L(E)×L(E)\mathscr{L}(E) \times \mathscr{L}(E) into L(E)\mathscr{L}(E) is continuous, as is the map X↦X−1X \mapsto X^{-1} of the (open) group of invertible elements of L(E)\mathscr{L}(E) onto itself (TG, IX, p. 40, prop. 14).

One may note that the map

t↦C(t0,t)=(C(t,t0))−1 t \mapsto C(t_0, t) = (C(t, t_0))^{-1}
admits a derivative equal to −(C(t,t0))−1(dC(t,t0)/dt)(C(t,t0))−1-(C(t, t_0))^{-1} (dC(t, t_0)/dt)(C(t, t_0))^{-1} (on the complement of a countable set) (I, p. 8, prop. 4), that is to say (by IV, p. 179, formula (5)) equal to −C(t0,t)A(t)-C(t_0, t) A(t).

Corollary 3

Let K be a compact interval contained in J, and let k=sup⁡t∈K∥A(t)∥k = \sup_{t \in K} \|A(t)\|. For all t and t0t_0 in K*

∥C(t,t0)−I∥≤ek∣t−t0∣−1. \|C(t, t_0) - I\| \leq e^{k|t-t_0|} - 1 .
(7)

Indeed, ∥A(t)x0∥≤k∥x0∥\|A(t)x_0\| \leq k \|x_0\| for all t∈Kt \in K; on K the constant function equal to x0x_0 is thus an approximate integral to within k∥x0∥k \|x_0\| by equation (4) of IV, p. 18; by formula (15) of IV, p. 170, one thus has

∥C(t,t0)x0−x0∥≤∥x0∥(ek∣t−t0∣−1) \|C(t, t_0)x_0 - x_0\| \leq \|x_0\| (e^{k|t-t_0|} - 1)
for any tt and t0t_0 in K, and x0x_0 in E, which is equivalent to the inequality (7) by the definition of the norm on L(E)\mathscr{L}(E).

Proposition 2

Let B be a continuous endomorphism of E, independent of t, and commuting with A(t) for all t∈Jt \in J; then B commutes with C(t,t0)C(t, t_0) for all t and t0t_0 in J.

Indeed, by (5)

ddt(BC)=BAC=ABCandddt(CB)=ACB, \frac{d}{dt}(BC) = BAC = ABC \quad \text{and} \quad \frac{d}{dt}(CB) = ACB ,
so ddt(BC−CB)=A(BC−CB)\frac{d}{dt}(BC - CB) = A(BC - CB); but BC(t0,t0)−C(t0,t0)B=0BC(t_0, t_0) - C(t_0, t_0)B = 0, so (IV, p. 179, th. 1) BC(t,t0)−C(t,t0)B=0BC(t, t_0) - C(t, t_0)B = 0 for all t∈Jt \in J.

An important instance of prop. 2 is that where E is endowed with the structure of a normed vector space with respect to the field of complex numbers C\mathbf{C}, and where, for every t∈Jt \in J, A(t)A(t) is an endomorphism of E for this vector space structure; this means that A(t)A(t) commutes with the continuous endomorphism x↦txx \mapsto tx of E (for the vector space structure over R\mathbf{R}); then C(t,t0)C(t, t_0) commutes with this endomorphism, which means that for any tt and t0t_0 in J, the map C(t,t0)C(t, t_0) is a continuous endomorphism for the normed vector space structure of E over C\mathbf{C}.

3. Integrating the inhomogeneous linear equation

Integrating the inhomogeneous linear equation

dxdt=A(t).x+b(t) \frac{d\mathbf{x}}{dt} = A(t).\mathbf{x} + \mathbf{b}(t)

reduces to integrating the associated homogeneous equation

dxdt=A(t).x \frac{d\mathbf{x}}{dt} = A(t).\mathbf{x}

and evaluating a primitive. With the notation of th. 2 of IV, p. 179, let us put x=C(t,t0).z\mathbf{x} = C(t, t_0).\mathbf{z}, whence, from the second formula (6) of IV, p. 179, z=C(t0,t).x\mathbf{z} = C(t_0, t).\mathbf{x}; if x\mathbf{x} is an integral of (2) then z\mathbf{z} is an integral of the equation ddt(C(t,t0).z)=A(t)C(t,t0).z+b(t)\frac{d}{dt}\left(C(t, t_0).\mathbf{z}\right) = A(t)C(t, t_0).\mathbf{z} + \mathbf{b}(t); since the bilinear map

(U,y)↦U.y (U, y) \mapsto U.y

of L(E)×E\mathscr{L}(\mathbf{E}) \times \mathbf{E} into E\mathbf{E} is continuous (\emph{Gen. Top.}, X, p. 297, prop. 6), z\mathbf{z} admits a derivative (except on a countable subset of J) and one has, by the formula for differentiating a bilinear function (I, p. 6, prop. 3)

ddt(C(t,t0).z)=dC(t,t0)dt.z+C(t,t0).dzdt=A(t)C(t,t0).z+C(t,t0).dzdt \frac{d}{dt}\left(C(t, t_0).\mathbf{z}\right) = \frac{dC(t, t_0)}{dt}.\mathbf{z} + C(t, t_0).\frac{d\mathbf{z}}{dt} = A(t)C(t, t_0).\mathbf{z} + C(t, t_0).\frac{d\mathbf{z}}{dt}

(replacing dC(t,t0)/dtdC(t, t_0)/dt by A(t)C(t,t0)A(t)C(t, t_0) according to (5) (IV, p. 179)). The equation for z\mathbf{z} then reduces to C(t,t0).dz/dt=b(t)C(t, t_0).d\mathbf{z}/dt = \mathbf{b}(t), or again to

dzdt=C(t0,t).b(t) \frac{d\mathbf{z}}{dt} = C(t_0, t).\mathbf{b}(t)

by the second formula (6) of IV, p. 179. Now the right-hand side of equation (8) is a regulated function on J, having been obtained by substituting the regulated functions UU and yy in the continuous bilinear function U.yU.y (\emph{cf.} II, p. 55, cor. 2); equation (8) thus has one and only one integral taking the value x0\mathbf{x}_0 at t0t_0, given by the formula

z(t)=x0+∫t0tC(t0,s).b(s) ds. \mathbf{z}(t) = \mathbf{x}_0 + \int_{t_0}^t C(t_0, s).\mathbf{b}(s)\,ds.

Since one has C(t,t0).∫t0tC(t0,s).b(s) ds=∫t0tC(t,t0)C(t0,s).b(s) dsC(t, t_0).\int_{t_0}^t C(t_0, s).\mathbf{b}(s)\,ds = \int_{t_0}^t C(t, t_0)C(t_0, s).\mathbf{b}(s)\,ds (II, p. 59, formula (9)), one obtains (taking account of the first formula (6) of IV, p. 179) the following result:

Proposition 3

With the notation of th. 2* (IV, p. 179), for every point (t0,x0)(t_0, \mathbf{x}_0) of J×EJ \times \mathbf{E} the integral of the linear equation (2) defined on JJ and equal to x0\mathbf{x}_0 at t0t_0 *is given by the formula

u(t)=C(t,t0).x0+∫t0tC(t,s).b(s) ds. \mathbf{u}(t) = C(t, t_0).\mathbf{x}_0 + \int_{t_0}^t C(t, s).\mathbf{b}(s)\,ds.

The method which leads to formula (10), which consists of taking the function z\mathbf{z} as a new unknown function, is often called the "method of variation of parameters".

4. Fundamental systems of integrals of a linear system of scalar differential equations

In this subsection and the next we shall consider the case where E is a vector space of finite dimension n over the field C\mathbf{C} of complex numbers (so of dimension 2n over R\mathbf{R}), and where, for every t∈Jt \in J, A(t)A(t) is an endomorphism of E with respect to the vector space structure over C\mathbf{C}. One can then identify A(t)A(t) with its matrix (aij(t))(a_{ij}(t)) with respect to a basis of E (over the field C\mathbf{C}), the aija_{ij} this time being n2n^2 complex functions defined and regulated on J; letting xjx_j (1≤j≤n1 \leq j \leq n) denote the (complex) components of a vector x∈Ex \in E with respect to the chosen basis, the linear equation

dxdt=A(t).x+b(t) \frac{d\mathbf{x}}{dt} = A(t).\mathbf{x} + \mathbf{b}(t)

is again equivalent to the system

dxidt=∑j=1naij(t)xj+bi(t)(1≤i≤n). \frac{dx_i}{dt} = \sum_{j=1}^n a_{ij}(t)x_j + b_i(t) \qquad (1 \leq i \leq n).

Theorems 1 (IV, p. 179) and 2 (IV, p. 179) and prop. 2 (IV, p. 181) then show that for every x0=(xk0)1≤k≤n\mathbf{x}_0 = (x_{k0})_{1 \leq k \leq n} in E there exists one and only one integral u=(uk)1≤k≤n\mathbf{u} = (u_k)_{1 \leq k \leq n} of the equation

dxdt=A(t).x \frac{d\mathbf{x}}{dt} = A(t).\mathbf{x}

defined on E and equal to x0\mathbf{x}_0 at the point t0t_0; this integral can be written

u(t,t0,x0)=C(t,t0).x0, \mathbf{u}(t, t_0, \mathbf{x}_0) = C(t, t_0).\mathbf{x}_0,

C(t,t0)C(t, t_0) being an invertible square matrix (cij(t,t0))(c_{ij}(t, t_0)) of order nn whose entries are continuous complex functions on J×JJ \times J and such that t↦cij(t,t0)t \mapsto c_{ij}(t, t_0) is a primitive of a regulated function on J.

In the particular case where n=1n = 1 the system (3) reduces to a single scalar equation

dxdt=a(t)x+b(t) \frac{dx}{dt} = a(t)x + b(t)

a(t)a(t) and b(t)b(t) being complex regulated functions on J); one verifies immediately that the (one element) matrix C(t,t0)C(t, t_0) is equal to exp⁡(∫t0ta(s) ds)\exp \left( \int_{t_0}^t a(s)\, ds \right); the integral of (11) equal to x0x_0 at the point t0t_0 is thus given explicitly by the formula

u(t)=x0exp⁡(∫t0ta(s) ds)+∫t0tb(s)exp⁡(∫t0τa(τ) dτ)ds. u(t) = x_0 \exp \left( \int_{t_0}^t a(s)\, ds \right) + \int_{t_0}^t b(s) \exp \left( \int_{t_0}^\tau a(\tau)\, d\tau \right) ds.

In the space C(J;E)\mathscr{C}(J; E) of continuous maps from J into E, endowed with the topology of compact convergence, the set I\mathscr{I} of integrals of equation (4) is a vector subspace (over C\mathbf{C}) isomorphic to E, therefore to Cn\mathbf{C}^n (IV, p. 180, cor. 1, and IV, p. 181, prop. 2).

A basis (uj)1≤j≤n(\mathbf{u}_j)_{1 \leq j \leq n} of this space (over the field C\mathbf{C}) is called a fundamental system of integrals of (4).

Proposition 4

For the nn integrals uj\mathbf{u}_j (1≤j≤n1 \leq j \leq n) of equation (4) to form a fundamental system it is necessary and sufficient that their values uj(t0)\mathbf{u}_j(t_0) at a point t0∈Jt_0 \in J be linearly independent vectors in EE.

Indeed, the map which to every x0∈E\mathbf{x}_0 \in E associates the integral t↦C(t,t0).x0t \mapsto C(t, t_0).\mathbf{x}_0 is an isomorphism of EE onto I\mathscr{I} (IV, p. 180, cor. 1 and IV, p. 181, prop. 2).

If (ej)1≤j≤n(\mathbf{e}_j)_{1 \leq j \leq n} is any basis of EE over C\mathbf{C}, the nn integrals

uj(t)=C(t,t0).ej(1≤j≤n) \mathbf{u}_j(t) = C(t, t_0).\mathbf{e}_j \qquad (1 \leq j \leq n)
thus form a fundamental system; if one identifies C(t,t0)C(t, t_0) with its matrix with respect to the basis (ej)(\mathbf{e}_j) the integrals uj\mathbf{u}_j are precisely the columns of the matrix C(t,t0)C(t, t_0). The integral of (4) that takes the value x0=∑j=1nλjej\mathbf{x}_0 = \sum_{j=1}^n \lambda_j \mathbf{e}_j at the point t0t_0 is then C(t,t0).x0=∑k=1nλkuk(t)C(t, t_0).\mathbf{x}_0 = \sum_{k=1}^n \lambda_k \mathbf{u}_k(t).

Given any nn integrals uj\mathbf{u}_j (1≤j≤n1 \leq j \leq n) of (4) one terms the determinant of these nn integrals at a point t∈Jt \in J with respect to a basis (ej)1≤j≤n(\mathbf{e}_j)_{1 \leq j \leq n} of EE, the determinant

Δ(t)=(u1(t), u2(t), …, un(t)) \Delta(t) = (\mathbf{u}_1(t), \ \mathbf{u}_2(t), \ \ldots, \ \mathbf{u}_n(t))
of the nn vectors uj(t)\mathbf{u}_j(t) with respect to the basis (ej)(\mathbf{e}_j) (Alg., III, p. 522). One has (Alg., III, p. 523, prop. 2)
Δ(t)=Δ(t0)det⁡(C(t,t0)). \Delta(t) = \Delta(t_0) \det \left( C(t, t_0) \right).

By prop. 4 of IV, p. 184, for (uj)1≤j≤n(\mathbf{u}_j)_{1 \leq j \leq n} to be a fundamental system of integrals of (4) it is necessary and sufficient that the determinant Δ(t)\Delta(t) of the uj\mathbf{u}_j be ≠0\neq 0 at some one point t0t_0 of JJ; the formula (14) then shows that Δ(t)≠0\Delta(t) \neq 0 at every point of JJ, in other words, that the vectors uj(t)\mathbf{u}_j(t) (1≤j≤n1 \leq j \leq n) are always linearly independent.

Proposition 5

The determinant of the matrix C(t,t0)C(t, t_0) is given by the formula*

det⁡(C(t,t0))=exp⁡(∫t0tTr⁡(A(s)) ds). \det \left( C(t, t_0) \right) = \exp \left( \int_{t_0}^t \operatorname{Tr}(A(s)) \, ds \right).
(15)

Indeed, if one puts δ(t)=det⁡(C(t,t0))\delta(t) = \det \left( C(t, t_0) \right) one has, by the formula for the derivative of a determinant (I, p. 8, formula (3))

dδdt=Tr⁡(dC(t,t0)dt(C(t,t0))−1)δ(t) \frac{d \delta}{dt} = \operatorname{Tr} \left( \frac{dC(t, t_0)}{dt} (C(t, t_0))^{-1} \right) \delta(t)

that is, by the differential equation (5) of IV, p. 179 satisfied by C(t,t0)C(t, t_0),

dδdt=Tr⁡(A(t))δ(t). \frac{d \delta}{dt} = \operatorname{Tr}(A(t)) \delta(t).

Since δ(t0)=1\delta(t_0) = 1 the formula (15) follows from the expression (12) (IV, p. 183) for the integral of a scalar linear equation.

Specifying nn linearly independent integrals of (4) determines all the integrals of this equation, as we have just seen. We shall now show that for 1≤p≤n1 \leq p \leq n the knowledge of pp linearly independent integrals uj\mathbf{u}_j (1≤j≤p1 \leq j \leq p) of equation (4) reduces the integration of this equation to that of a homogeneous linear system of n−pn - p scalar equations. Suppose that on an interval K⊂JK \subset J there are n−pn - p maps

up+k(1≤k≤n−p) \mathbf{u}_{p+k} \quad (1 \leq k \leq n - p)

of KK into EE, which are primitives of regulated functions on KK, and such that, for every t∈Kt \in K, the nn vectors uj(t)\mathbf{u}_j(t) (1≤j≤n1 \leq j \leq n) form a basis of EE.

For every point t1∈Jt_1 \in J there always exists an interval KK, a neighbourhood of t1t_1 in JJ, on which there are defined n−pn - p functions up+k\mathbf{u}_{p+k} (1≤k≤n−p1 \leq k \leq n - p) having the preceding properties. For, let (ei)1≤i≤n(\mathbf{e}_i)_{1 \leq i \leq n} be a basis of EE; there exist n−pn - p vectors of this basis which form with the uj(t1)\mathbf{u}_j(t_1) (1≤j≤p1 \leq j \leq p) a basis of EE (Alg., II, p. 292, th. 2); suppose for example that they are ep+1,…,en\mathbf{e}_{p+1}, \ldots, \mathbf{e}_n; since the determinant det⁡(u1(t),…,up(t),ep+1,…,en)\det(\mathbf{u}_1(t), \ldots, \mathbf{u}_p(t), \mathbf{e}_{p+1}, \ldots, \mathbf{e}_n) (with respect to the basis (ei)(\mathbf{e}_i)) is a continuous function of tt and does not vanish for t=t1t = t_1 there exists a neighbourhood KK of t1t_1 on which it does not vanish; one can then take up+k(t)=ep+k\mathbf{u}_{p+k}(t) = \mathbf{e}_{p+k} (1≤k≤n−p1 \leq k \leq n - p) for t∈Kt \in K.

There exists an invertible matrix B(t)B(t) of order nn, whose elements are primitives of regulated functions on KK, such that B(t).ej=uj(t)B(t).\mathbf{e}_j = \mathbf{u}_j(t) for 1≤j≤n1 \leq j \leq n. Let us put x=B(t).yx = B(t).y; then yy satisfies the equation dBdt.y+B(t).dydt=A(t)B(t).y\frac{dB}{dt}.y + B(t).\frac{dy}{dt} = A(t)B(t).y, which can also be written

dydt=(B(t))−1(A(t)B(t)−dBdt).y=H(t).y \frac{dy}{dt} = (B(t))^{-1} \left( A(t)B(t) - \frac{dB}{dt} \right) .y = H(t).y

where H(t)=(hjk(t))H(t) = (h_{jk}(t)) is a matrix with regulated entries on KK. By the definition of B(t)B(t) this linear equation admits the pp constant vectors ej\mathbf{e}_j (1≤j≤p1 \leq j \leq p) as integrals; one concludes immediately that necessarily hjk(t)=0h_{jk}(t) = 0 for 1≤k≤p1 \leq k \leq p; thus the components yky_k of yy (with respect to the basis (ei)(\mathbf{e}_i)) with index k≥p+1k \geq p + 1 satisfy a homogeneous linear system of n−pn - p equations; once the solutions of this system are determined, the dyj/dtdy_j/dt for indices j≤pj \leq p are linear functions of the yky_k with k≥p+1k \geq p + 1, so are known, and the primitives of these functions will give the yjy_j for indices j≤pj \leq p.

In particular, when one knows n−1n - 1 linearly independent integrals of equation (4) of IV, p. 183, integrating this equation reduces to that of a single homogeneous scalar equation, and then the evaluation of nn primitives.

Remark 1

All the above applies also to the case where E is of dimension nn over the field R\mathbf{R} and A(t)A(t) is an endomorphism of E for every t∈Jt \in J: one has only to replace C\mathbf{C} by R\mathbf{R} throughout.

Remark 2

Let A(t)=(aij(t))A(t) = (a_{ij}(t)) be a square matrix of order nn whose entries are real (resp. complex) regulated functions of tt on JJ, and let C(t,t0)=(cij(t,t0))C(t, t_0) = (c_{ij}(t, t_0)) be the resolvent matrix of the corresponding linear system (3) (IV, p. 22). Let F be an arbitrary complete normed space over R\mathbf{R} (resp. C\mathbf{C}) and let us consider the system of linear differential equations

dyidt=∑j=1naij(t)yj \frac{dy_i}{dt} = \sum_{j=1}^n a_{ij}(t) y_j

where the unknown functions yjy_j, take their values in F. It is immediate that the solution (uj)1≤j≤n(\mathbf{u}_j)_{1 \leq j \leq n} of this system such that uj(t0)=dj\mathbf{u}_j(t_0) = \mathbf{d}_j for 1≤j≤n1 \leq j \leq n (dj\mathbf{d}_j being arbitrary in F) is given by the formulae

ui(t)=∑j=1ncij(t,t0)dj(1≤i≤n). \mathbf{u}_i(t) = \sum_{j=1}^n c_{ij}(t, t_0) \mathbf{d}_j \qquad (1 \leq i \leq n).

We consider in particular the case where A(t)A(t) is an endomorphism of a vector space E of finite dimension nn over C\mathbf{C}, such that there exists a basis of E with respect to which the matrix of A(t)A(t) has real elements for every t∈Jt \in J. Then the above shows (by th. 1 of IV, p. 179) that the resolvent matrix C(t,t0)C(t, t_0) with respect to the same basis also has real elements: it suffices to consider the vector space E0E_0 over R\mathbf{R} generated by the basis of E under consideration, and to remark that the restriction of A(t)A(t) to E0E_0 is an endomorphism of this vector space.

5. Adjoint equation

Assuming always that the space E is of finite dimension nn over C\mathbf{C}, let E∗E^* be its dual (A, II, p. 40), which is a space of dimension nn over C\mathbf{C} (Alg., II, p. 299, th. 4); the canonical bilinear form ⟨x,x∗⟩\langle \mathbf{x}, \mathbf{x}^* \rangle defined on E×E∗E \times E^* (Alg., II, p. 234) is continuous on this product (being a polynomial in the components of x∈E\mathbf{x} \in E and x∗∈E∗\mathbf{x}^* \in E^*).

Given a homogeneous linear equation (4) (IV, p. 183), where t↦A(t)t \mapsto A(t) is a regulated map of J into L(E)\mathscr{L}(E), let us see if there exists a map t↦v(t)t \mapsto \mathbf{v}(t) of J into E∗E^*, a primitive of a regulated function on J, and such that the scalar function t↦⟨u(t),v(t)⟩t \mapsto \langle \mathbf{u}(t), \mathbf{v}(t) \rangle is constant on J when u\mathbf{u} is an arbitrary solution of (4); it comes to the same to write that the derivative of this function should be zero at every point where u\mathbf{u} and v\mathbf{v} are differentiable, that is, one must have

⟨dudt,v(t)⟩+⟨u(t),dvdt⟩=0 \left\langle \frac{d\mathbf{u}}{dt}, \mathbf{v}(t) \right\rangle + \left\langle \mathbf{u}(t), \frac{d\mathbf{v}}{dt} \right\rangle = 0

at such points.

Now, by (4), ⟨dudt,v(t)⟩=⟨A(t).u(t),v(t)⟩=−⟨u(t),B(t).v(t)⟩\left\langle \frac{d\mathbf{u}}{dt}, \mathbf{v}(t) \right\rangle = \langle A(t).\mathbf{u}(t), \mathbf{v}(t) \rangle = -\langle \mathbf{u}(t), B(t).\mathbf{v}(t) \rangle where −B(t)-B(t) is the transpose of A(t)A(t) (Alg., II, p. 234). The relation that v\mathbf{v} must satisfy can thus be written

⟨u(t),dvdt−B(t).v(t)⟩=0 \left\langle \mathbf{u}(t), \frac{d\mathbf{v}}{dt} - B(t).\mathbf{v}(t) \right\rangle = 0

at all points where A(t)A(t) is continuous and v(t)\mathbf{v}(t) is differentiable. Now for such a point tt and an arbitrary point x0∈E\mathbf{x}_0 \in \mathrm{E} there exists, by th. 1 of IV, p. 179, a solution u\mathbf{u} of (4) such that u(t)=x0\mathbf{u}(t) = \mathbf{x}_0; one thus must have ⟨x0,dvdt−B(t).v(t)⟩=0\left\langle \mathbf{x}_0, \frac{d\mathbf{v}}{dt} - B(t).\mathbf{v}(t) \right\rangle = 0 for all x0∈E\mathbf{x}_0 \in \mathrm{E}, which entails dvdt−B(t).v(t)=0\frac{d\mathbf{v}}{dt} - B(t).\mathbf{v}(t) = 0. Consequently:

Proposition 6

The map t↦v(t)t \mapsto \mathbf{v}(t) of J into E∗\mathrm{E}^*, a primitive of a regulated function on J, is such that ⟨u(t),v(t)⟩\langle \mathbf{u}(t), \mathbf{v}(t) \rangle is constant on J for every solution u\mathbf{u} of the equation* (4) *of IV, p. 183 if and only if v\mathbf{v} is a solution of the homogeneous linear equation

dxdt=B(t).x(16) \frac{d\mathbf{x}}{dt} = B(t).\mathbf{x} \tag{16}

where −B(t)-B(t) is the transpose of A(t)A(t).

Equation (16) is called the adjoint of (4); clearly (4) is the adjoint of (16). The elements of the matrix B(t)B(t) being regulated functions of tt on J, the results obtained above on linear equations apply to (16). In particular, the integral of (16) taking the value x0∗\mathbf{x}_0^* at the point t0t_0 can be written as H(t,t0).x0∗H(t, t_0).\mathbf{x}_0^*, where H(t,t0)H(t, t_0) is a bijective linear map of E∗\mathrm{E}^* onto itself, identical to the integral of the equation

dVdt=B(t)V(17) \frac{dV}{dt} = B(t)V \tag{17}

which takes the value II at the point t0t_0. As a result one has (with the notation of IV, p. 179)

⟨C(t,t0).x0,H(t,t0).x0∗⟩=⟨x0,x0∗⟩ \left\langle C(t, t_0).\mathbf{x}_0, H(t, t_0).\mathbf{x}_0^* \right\rangle = \left\langle \mathbf{x}_0, \mathbf{x}_0^* \right\rangle

for any x0∈E\mathbf{x}_0 \in \mathrm{E} and x0∗∈E∗\mathbf{x}_0^* \in \mathrm{E}^*, which shows that

H(t,t0)=Cˇ(t,t0)(18) H(t, t_0) = \check{C}(t, t_0) \tag{18}

(the contragradient of C(t,t0)C(t, t_0)). In particular, if one knows a fundamental system of integrals of the adjoint equation (16) the matrix H(t,t0)H(t, t_0) is determined, so is C(t,t0)C(t, t_0) also, and as a result, so are all the integrals of equation (4).

Remark

Let E and F be two complete normed spaces over R\mathbf{R} (or over C\mathbf{C}), and (x,y)↦⟨x,y⟩(\mathbf{x}, \mathbf{y}) \mapsto \langle \mathbf{x}, \mathbf{y} \rangle a continuous bilinear form on E×F\mathrm{E} \times \mathrm{F}, such that the relation “⟨x,y⟩=0\langle \mathbf{x}, \mathbf{y} \rangle = 0 for all y∈F\mathbf{y} \in \mathrm{F}” (resp. ⟨x,y⟩=0\langle \mathbf{x}, \mathbf{y} \rangle = 0 for all x∈E\mathbf{x} \in \mathrm{E}”) implies x=0\mathbf{x} = 0 (resp. y=0\mathbf{y} = 0). Suppose further that for every t∈Jt \in \mathrm{J} there is a continuous linear map B(t)B(t) of F into itself, such that ⟨A(t).x,y⟩+⟨x,B(t).y⟩=0\langle A(t).\mathbf{x}, \mathbf{y} \rangle + \langle \mathbf{x}, B(t).\mathbf{y} \rangle = 0 for every (x,y)∈E×F(\mathbf{x}, \mathbf{y}) \in \mathrm{E} \times \mathrm{F}. In these circumstances one sees as before that for a map t↦v(t)t \mapsto \mathbf{v}(t) of J into F, a primitive of a regulated function, to be such that ⟨u(t),v(t)⟩\langle \mathbf{u}(t), \mathbf{v}(t) \rangle is constant for every integral u\mathbf{u} of (4), it is necessary and sufficient that v\mathbf{v} should be an integral of (16), which again we call the adjoint of (4).

6. Linear differential equations with constant coefficients

Suppose again that E is an arbitrary complete normed space over R\mathbf{R}; let A be a continuous endomorphism of A, independent of t, and consider the homogeneous linear equation

dxdt=A.x. \frac{d\mathbf{x}}{dt} = A.\mathbf{x}.

When E is of finite dimension the equation (19) is equivalent to a homogeneous system (3) (IV, p. 183) of scalar differential equations, where the coefficients aija_{ij} are constant.

By th. 1 (IV, p. 179), every integral of (19) is defined on all of R\mathbf{R}; by th. 2 (IV, p. 179) the integral of (19) taking the value x0\mathbf{x}_0 at a point t0∈Rt_0 \in \mathbf{R} can be written as C(t,t0)x0C(t, t_0)\mathbf{x}_0, where C(t,t0)C(t, t_0) is a bijective bicontinuous linear map of E onto itself that satisfies the equation

dUdt=AU \frac{dU}{dt} = AU

and is such that C(t0,t0)=IC(t_0, t_0) = I. Moreover, we have the identity

C(t+τ,t0+τ)=C(t,t0) C(t + \tau, t_0 + \tau) = C(t, t_0)

for any τ∈R\tau \in \mathbf{R}: indeed, one has dC(s,t0+τ)/ds=AC(s,t0+τ)dC(s, t_0 + \tau)/ds = AC(s, t_0 + \tau) by (20), and since A is constant one also has

dC(t+τ,t0+τ)dt=AC(t+τ,t0+τ); \frac{dC(t + \tau, t_0 + \tau)}{dt} = AC(t + \tau, t_0 + \tau);

moreover

C(t0+τ,t0+τ)=I=C(t0,t0). C(t_0 + \tau, t_0 + \tau) = I = C(t_0, t_0).

whence we have the identity (21), since the integral of (20) equal to I at the point t0t_0 is unique.

If one puts C0(t)=C(t,0)C_0(t) = C(t, 0) then C(t,t0)=C0(t−t0)C(t, t_0) = C_0(t - t_0); moreover, for every λ∈R\lambda \in \mathbf{R}, C0(λt)C_0(\lambda t) is identical to the integral of the equation

dUdt=λAU \frac{dU}{dt} = \lambda AU

which takes the value I at the point 0. We make the following definition:

Definition 1

Given a continuous endomorphism A of E we denote by eAe^A or exp⁡A\exp A the automorphism of E equal to the value at the point t=1t = 1 of the integral of the equation (20) which takes the value I at the point t=0t = 0.

With this notation the remarks preceding def. 1 show that

C(t,t0)=exp⁡(A(t−t0)). C(t, t_0) = \exp \left( A(t - t_0) \right).

The exponential notation just introduced is justified by the following properties, which are entirely analogous to those of the function exp z, for z real or complex (cf. III, p. 98 and 106):

Proposition 7

1' The map X↦eXX \mapsto e^X is a continuous map of L(E)\mathscr{L}(E) into the group of automorphisms of E (invertible elements of* L(E)\mathscr{L}(E)).

2 The map t↦eXtt \mapsto e^{Xt} of R\mathbf{R} into L(E)\mathscr{L}(E) is differentiable, and

ddt(eXt)=XeXt=eXtX. \frac{d}{dt} (e^{Xt}) = X e^{Xt} = e^{Xt} X.
(24)

3 For any X∈L(E)X \in \mathscr{L}(E) one has

eX=∑n=0∞Xnn! e^X = \sum_{n=0}^{\infty} \frac{X^n}{n!}
(25)

the right-hand side being absolutely and uniformly convergent on every bounded subset of L(E)\mathscr{L}(E); in particular, eIt=etIe^{It} = e^t I for t∈Rt \in \mathbf{R}.

4) If X and Y commute then Y and eYe^Y commute with eXe^X, and

eX+Y=eXeY. e^{X+Y} = e^X e^Y.
(26)

The relation (24) follows from the expression (23) for C(t,0)C(t, 0) and from the fact that this function is an integral of (20); by recursion on n one deduces from (24) that t↦eXtt \mapsto e^{Xt} is indefinitely differentiable on R\mathbf{R} and that

Dn(eXt)=XneXt. D^n (e^{Xt}) = X^n e^{Xt}.

By Taylor’s formula one thus can write

eX=I+X1!+X22!+⋯+Xnn!+Xn+1∫01(1−t)nn!eXt dt. e^X = I + \frac{X}{1!} + \frac{X^2}{2!} + \cdots + \frac{X^n}{n!} + X^{n+1} \int_0^1 \frac{(1-t)^n}{n!} e^{Xt} \, dt.
(27)

On the other hand, cor. 3 of IV, p. 181 shows that ∥eXt∥≤exp⁡(∥X∥∣t∣)\| e^{Xt} \| \leq \exp (\| X \| |t| ). Thus the remainder rn(X)=Xn+1∫01(1−t)nn!eXt dtr_n(X) = X^{n+1} \int_0^1 \frac{(1-t)^n}{n!} e^{Xt} \, dt in formula (27) satisfies the inequality

∥rn(X)∥≤∥X∥n+1(n+1)!e∥X∥ \| r_n(X) \| \leq \frac{\| X \|^{n+1}}{(n+1)!} e^{\| X \| }

whence one deduces the formula (25), the series on the right-hand side being absolutely and uniformly convergent on every bounded subset of L(E)\mathscr{L}(E). For every pair of elements X, TX,\ T of L(E)\mathscr{L}(E) one thus has

eX+T−eX=∑n=1∞1n!((X+T)n−Xn). e^{X+T} - e^X = \sum_{n=1}^{\infty} \frac{1}{n!} ((X+T)^n - X^n ).

Now, one can write (X+T)n−Xn=∑(Vi)V1V2…Vn(X + T)^n - X^n = \sum_{(V_i)} V_1 V_2 \ldots V_n, where the sum is taken over the 2n−12^n - 1 sequences (Vi)(V_i) of elements of L(E)\mathscr{L}(E) such that Vi=XV_i = X or Vi=TV_i = T for 1≤i≤n1 \leq i \leq n, and at least one of the ViV_i is equal to TT; the inequality

∥(X+T)n−Xn∥≤(∥X∥+∥T∥)n−∥X∥n \|(X + T)^n - X^n\| \leq (\|X\| + \|T\|)^n - \|X\|^n
follows immediately, whence
∥exp⁡(X+T)−exp⁡X∥≤exp⁡(∥X∥+∥T∥)−exp⁡∥X∥ \|\exp(X + T) - \exp X\| \leq \exp \left( \|X\| + \|T\| \right) - \exp \|X\|
which establishes the continuity of the map X↦exp⁡XX \mapsto \exp X.

Finally, if XX and YY commute, then YY commutes with eXte^{Xt} (IV, p. 181, prop. 2), so

ddt(eXteYt)=XeXteYt+eXt(YeYt)=(X+Y)eXteYt. \frac{d}{dt} (e^{Xt} e^{Yt}) = X e^{Xt} e^{Yt} + e^{Xt} (Y e^{Yt}) = (X + Y) e^{Xt} e^{Yt}.
Since, on the other hand, eXteYte^{Xt} e^{Yt} is equal to II for t=0t = 0, one has eXteYt=e(X+Y)te^{Xt} e^{Yt} = e^{(X+Y)t}, whence formula (26). From this latter, one deduces in particular that for arbitrary real ss and tt one has
eX(s+t)=eXseXt e^{X(s+t)} = e^{Xs} e^{Xt}
and also that
e−X=(eX)−1. e^{-X} = (e^X)^{-1}.

On the contrary one notes that (26) need not remain valid if XX and YY are no longer assumed to commute: indeed it would imply that exp⁡X\exp X and exp⁡Y\exp Y always commute, which is not the case, as is shown by simple examples (IV, p. 204, exerc. 3).

Now let us assume that EE is a vector space of finite dimension over the field C\mathbf{C}, and AA an endomorphism of EE (for its vector space structure over C\mathbf{C}) which one can identify with its matrix with respect to a basis of EE; then, for all t∈Rt \in \mathbf{R}, eAte^{At} is an automorphism of EE for the same structure (IV, p. 181, prop. 2). Let rkr_k (1≤k≤q1 \leq k \leq q) be the distinct roots (in C\mathbf{C}) of the characteristic polynomial φ(r)=det⁡(A−rI)\varphi(r) = \det(A - rI) of the endomorphism AA (the "characteristic roots" of AA); if nkn_k is the order of multiplicity of rkr_k then ∑k=1qnk=n\sum_{k=1}^q n_k = n. One knows that (Alg., VII, 31, n° 3) to each root rkr_k there corresponds a subspace EkE_k of EE, of dimension nkn_k, such that EkE_k is stable under AA, and that EE is the direct sum of the EkE_k: EkE_k can be defined as the subspace of vectors x\mathbf{x} such that

(A−rkI)nk.x=0. (A - r_k I)^{n_k} . \mathbf{x} = 0.
Let a\mathbf{a} be any vector in EE; one can write a=∑k=1qak\mathbf{a} = \sum_{k=1}^q \mathbf{a}_k, where ak∈Ek\mathbf{a}_k \in E_k; the integral of the equation (19) of IV, p. 188, taking the value a\mathbf{a} at the point t=0t = 0 is thus given by
u(t)=eAt.a=∑k=1qeAt.ak=∑k=1qerkte(A−rkI)t.ak. \mathbf{u}(t) = e^{At} . \mathbf{a} = \sum_{k=1}^q e^{At} . \mathbf{a}_k = \sum_{k=1}^q e^{r_k t} e^{(A - r_k I)t} . \mathbf{a}_k.

But since ak∈Ek\mathbf{a}_k \in E_k one has

e(A−rkI)t⋅ak=ak+t1!(A−rkI)⋅ak+t22!(A−rkI)2⋅ak+⋯ e^{(A - r_k I)t} \cdot \mathbf{a}_k = \mathbf{a}_k + \frac{t}{1!}(A - r_k I) \cdot \mathbf{a}_k + \frac{t^2}{2!}(A - r_k I)^2 \cdot \mathbf{a}_k + \cdots
+tnk−1(nk−1)!(A−rkI)nk−1⋅ak. + \frac{t^{n_k - 1}}{(n_k - 1)!}(A - r_k I)^{n_k - 1} \cdot \mathbf{a}_k .

Thus every integral of the equation (19) of IV, p. 188, can be written as

u(t)=∑k=1qerktpk(t) \mathbf{u}(t) = \sum_{k=1}^q e^{r_k t} \mathbf{p}_k(t)

where pk(t)\mathbf{p}_k(t) is a polynomial in tt, with coefficients in the vector space EkE_k, and of degree ≤nk−1\leq n_k - 1. In particular, if all the roots of the characteristic equation of AA are simple, the spaces EkE_k (1≤k≤n1 \leq k \leq n) are all of dimension 1 over the field C\mathbf{C}, so there exist nn vectors ckc_k such that the nn functions erktcke^{r_k t} c_k (1≤k≤n1 \leq k \leq n) form a fundamental system of integrals of the equation (19) of IV, p. 188.

The characteristic roots of the endomorphism AA are also called characteristic roots of the linear equation (19) of IV, p. 188. One may observe that one obtains the characteristic equation of AA by writing that the function certc e^{r t} is an integral of (19) for a vector c≠0c \neq 0.

When one has determined the roots rkr_k (1≤k≤q1 \leq k \leq q) explicitly, and thus the order of multiplicity nkn_k of rkr_k, in practice one obtains the integrals of (19) by writing that this equation is satisfied by the expression (32) of IV, p. 191, where pk\mathbf{p}_k is an arbitrary polynomial of degree ≤nk−1\leq n_k - 1, with coefficients in EE; on identifying the coefficients of erkte^{r_k t} (for 1≤k≤q1 \leq k \leq q) in the two sides of the equation so obtained one obtains linear equations for the coefficients of the polynomials pk\mathbf{p}_k: one establishes easily that these equations determine the terms of degree >0> 0 of pk\mathbf{p}_k as functions of the constant term, and that the latter is a solution of the equation (A−rkI)nk⋅x=0(A - r_k I)^{n_k} \cdot \mathbf{x} = 0, which defines the subspace EkE_k (method of "undetermined coefficients").

Remark

When there exists a basis of EE such that the matrix of AA with respect to this basis has real elements (cf. IV, p. 186, Remark 2), the characteristic equation of AA has real coefficients. For every x=(ξk)1≤k≤n\mathbf{x} = (\xi_k)_{1 \leq k \leq n} of EE, expressed in the basis under consideration, let x‾=(ξ‾k)1≤k≤n\overline{\mathbf{x}} = (\overline{\xi}_k)_{1 \leq k \leq n}; the map x↦x‾\mathbf{x} \mapsto \overline{\mathbf{x}} is an antilinear involution of EE. One knows (Alg., VII) that, if rkr_k is a non-real root of the characteristic equation, and EkE_k is the corresponding stable subspace, then r‾k\overline{r}_k is a characteristic root having the same multiplicity nkn_k as rkr_k, and the image Ek′E'_k of EkE_k under the map x↦x‾\mathbf{x} \mapsto \overline{\mathbf{x}} is the stable subspace corresponding to r‾k\overline{r}_k. One deduces from this that if uj\mathbf{u}_j (1≤j≤nk1 \leq j \leq n_k) are nkn_k linearly independent integrals with values in EkE_k, then the 2nk2n_k integrals uj+u‾j, i(uj−u‾j)\mathbf{u}_j + \overline{\mathbf{u}}_j,\ i(\mathbf{u}_j - \overline{\mathbf{u}}_j) are linearly independent, and have, with respect to the chosen basis of EE, components which are real functions of EE. If rkr_k is a real characteristic root Remark 2 of IV, p. 186, shows that (with the same notation) there are nkn_k linearly independent integrals vj\mathbf{v}_j (1≤j≤nk1 \leq j \leq n_k) with values in EkE_k whose components are real. That way one obtains a fundamental system of integrals of (19) whose components are all real.

7. Linear equations of order nn

One calls a linear differential equation of order nn an equation of the form

Dnx−a1(t)Dn−1x−…−an−1(t)Dx−an(t)x=b(t) D^n x - a_1(t) D^{n-1} x - \ldots - a_{n-1}(t) D x - a_n(t) x = b(t)

where the aka_k (1≤k≤n1 \leq k \leq n) and bb are real (complex) functions of the real variable tt, defined on an interval JJ of R\mathbf{R}. The general procedure of IV, p. 164 shows that this equation is equivalent to the linear system of nn first order equations

{dxkdt=xk+1(1≤k≤n−1)dxndt=a1(t)xn+a2(t)xn−1+⋯+an(t)x1+b(t) \left\{ \begin{aligned} \frac{dx_k}{dt} &= x_{k+1} & (1 \leq k \leq n-1) \\ \frac{dx_n}{dt} &= a_1(t)x_n + a_2(t)x_{n-1} + \cdots + a_n(t)x_1 + b(t) \end{aligned} \right.

that is to say, to the linear equation

dxdt=A(t).x+b(t) \frac{d\mathbf{x}}{dt} = A(t).\mathbf{x} + \mathbf{b}(t)

where we have put x=(x1,x2,…,xn)∈Cn\mathbf{x} = (x_1, x_2, \ldots, x_n) \in \mathbf{C}^n, b(t)=(0,0,…,0,b(t))\mathbf{b}(t) = (0, 0, \ldots, 0, b(t)), and where the matrix A(t)A(t) is defined by

A(t)=(010…0001…0⋮⋮⋮⋱⋮000…1an(t)an−1(t)an−2(t)…a1(t)). A(t) = \begin{pmatrix} 0 & 1 & 0 & \ldots & 0 \\ 0 & 0 & 1 & \ldots & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & 0 & \ldots & 1 \\ a_n(t) & a_{n-1}(t) & a_{n-2}(t) & \ldots & a_1(t) \end{pmatrix}.

The study of the linear equation of order nn thus consists of applying the general results above to the particular linear equation (35). For every interval JJ where the functions aja_j (1≤j≤n1 \leq j \leq n) and bb are regulated there exists one and only one function uu, defined on JJ, having a continuous derivative of order n−1n-1, and a regulated derivative of order nn on this interval (except at the points of a countable set), satisfying (33) on the complement of a countable subset of JJ, and such that

u(t0)=x0,Du(t0)=x0′,…,Dn−1u(t0)=x0(n−1) u(t_0) = x_0, \qquad Du(t_0) = x_0', \ldots, \quad D^{n-1}u(t_0) = x_0^{(n-1)}

where t0t_0 is an arbitrary point of JJ, and x0,x0′,…,x0(n−1)x_0, x_0', \ldots, x_0^{(n-1)} are nn arbitrary complex numbers.

For the pp integrals uju_j (1≤j≤p1 \leq j \leq p) of the homogeneous equation

Dnx−a1(t)Dn−1x−…−an−1(t)Dx−an(t)x=0 D^n x - a_1(t) D^{n-1} x - \ldots - a_{n-1}(t) D x - a_n(t) x = 0

associated with (33) to be linearly independent (in the space C(J;C)C(J; \mathbf{C}) of continuous maps from JJ into C\mathbf{C}, considered as a vector space over C\mathbf{C}), it is necessary and sufficient that the corresponding pp integrals uj=(uj,Duj,…,Dn−1uj)\mathbf{u}_j = (u_j, Du_j, \ldots, D^{n-1}u_j) of the homogeneous equation dx/dt=A(t).xd\mathbf{x}/dt = A(t).\mathbf{x} should be linearly independent (in the space C(J;Cn)C(J; \mathbf{C}^n) of continuous maps from JJ into Cn\mathbf{C}^n). It is clear that this condition is necessary.

Conversely, if there are nn complex constants λj\lambda_j, not all zero, such that ∑j=1nλjuj(t)=0\sum_{j=1}^n \lambda_j u_j(t) = 0 identically on JJ, one deduces that ∑j=1nλjDkuj(t)=0\sum_{j=1}^n \lambda_j D^k u_j(t) = 0 on JJ for every integer kk such that 1≤k≤n−11 \leq k \leq n - 1, which means that ∑j=1nλjuj(t)=0\sum_{j=1}^n \lambda_j \mathbf{u}_j(t) = 0 on JJ.

Consequently (IV, p. 180, cor. 1)

Proposition 8

The set of integrals of the homogeneous linear equation (37), defined on JJ, is a vector space of dimension nn over the field C\mathbf{C}.

Given any nn integrals uju_j (1≤j≤n1 \leq j \leq n) of equation (37) one calls the Wronskian of this system of integrals the determinant (with respect to the canonical basis of Cn\mathbf{C}^n) of the corresponding nn integrals uj\mathbf{u}_j of the equation dx/dt=A(t).xd\mathbf{x}/dt = A(t).\mathbf{x}, that is to say, the function

W(t)=∣u1(t)u2(t)…un(t)Du1(t)Du2(t)…Dun(t)⋮⋮⋱⋮Dn−1u1(t)Dn−1u2(t)…Dn−1un(t)∣. W(t) = \begin{vmatrix} u_1(t) & u_2(t) & \ldots & u_n(t) \\ Du_1(t) & Du_2(t) & \ldots & Du_n(t) \\ \vdots & \vdots & \ddots & \vdots \\ D^{n-1}u_1(t) & D^{n-1}u_2(t) & \ldots & D^{n-1}u_n(t) \end{vmatrix}.

For the nn integrals uju_j to be linearly independent it is necessary and sufficient that W(t)≠0W(t) \neq 0 on JJ; moreover, it suffices for this that W(t0)≠0W(t_0) \neq 0 at only one t0t_0 of JJ (IV, p. 184, prop. 4); further, one has (IV, p. 184, prop. 5)

W(t)=W(t0)exp⁡(∫t0ta1(s) ds). W(t) = W(t_0) \exp \left( \int_{t_0}^t a_1(s)\, ds \right).
(38)

We identify the resolvent C(t,t0)C(t, t_0) of equation (35) with its matrix with respect to the canonical basis of Cn\mathbf{C}^n; the columns vj(t,t0)\mathbf{v}_j(t, t_0) (1≤j≤n1 \leq j \leq n) of this matrix are then nn linearly independent integrals

vj(t,t0)=(vj(t,t0),Dvj(t,t0),…,Dn−1vj(t,t0)) \mathbf{v}_j(t, t_0) = (v_j(t, t_0), Dv_j(t, t_0), \ldots, D^{n-1}v_j(t, t_0))

of the homogeneous equation dx/dt=A(t).xd\mathbf{x}/dt = A(t).\mathbf{x} which correspond to nn linearly independent integrals vj(t,t0)v_j(t, t_0) of equation (37) such that

Dk−1vj(t0,t0)=δjk D^{k-1}v_j(t_0, t_0) = \delta_{jk}

(Kronecker delta) for 1≤j≤n,1≤k≤n1 \leq j \leq n, 1 \leq k \leq n (on agreeing to put D0vj=vjD^0 v_j = v_j). It results in particular that the method of variation of parameters (IV, p. 182) applied to equation (35) here gives as a particular integral of (33), equal to 0 together with its first n−1n-1 derivatives at the point t0t_0, the function

w(t)=∫t0tvn(t,s)b(s) ds. w(t) = \int_{t_0}^t v_n(t, s)b(s)\, ds.
(39)

In the particular case of the equation Dnx=b(t)D^n x = b(t) the formula (39) again gives the formula expressing the nthn^{th} primitive of the regulated function b(t)b(t) which vanishes with its first n−1n - 1 derivatives at the point t0t_0, namely

w(t)=∫t0tb(s)(t−s)n−1(n−1)! ds w(t) = \int_{t_0}^t b(s) \frac{(t-s)^{n-1}}{(n-1)!} \, ds
(II, p. 62, formula (19)): the integral of Dnx=0D^n x = 0 which vanishes together with its first n−2n - 2 derivatives at the point t0t_0, and whose n−1n - 1th derivative is equal to 1 there, is actually the polynomial (t−t0)n−1/(n−1)!(t-t_0)^{n-1}/(n-1)!.

8. Linear equations of order nn with constant coefficients

If the coefficients aja_j in equation (33) are constant, the corresponding matrix AA is constant; the characteristic equation is obtained by writing that erte^{rt} is a solution, which gives

rn−a1rn−1−…−an−1r−an=0. r^n - a_1 r^{n-1} - \ldots - a_{n-1} r - a_n = 0.
Let rjr_j (1≤j≤q1 \leq j \leq q) be the distinct roots of this equation, and njn_j (1≤j≤q1 \leq j \leq q) the multiplicity of the root rj(∑j=1qnj=n)r_j \left( \sum_{j=1}^q n_j = n \right). By the results of IV, p. 188 to 194, to each root rjr_j there corresponds, for the homogeneous equation
Dnx−a1Dn−1x−⋯−an−1Dx−anx=0 D^n x - a_1 D^{n-1} x - \cdots - a_{n-1} D x - a_n x = 0
a system of njn_j linearly independent integrals
ujk(t)=erjtpjk(t), u_{jk}(t) = e^{r_j t} p_{jk}(t),
where pjkp_{jk} is a polynomial (with complex coefficients) of degree ≤nj−1\leq n_j - 1 (1≤k≤nj1 \leq k \leq n_j); further, the nn integrals ujku_{jk} (1≤j≤q, 1≤k≤nj1 \leq j \leq q,\ 1 \leq k \leq n_j) so obtained are linearly independent. It follows that the njn_j polynomials pjkp_{jk} (1≤k≤nj1 \leq k \leq n_j) are linearly independent in the space of polynomials in tt of degree ≤nj−1\leq n_j - 1, so form a basis (over C\mathbf{C}) of this space, since the latter is of dimension njn_j. In other words:

Proposition 9

Let rjr_j (1≤j≤q1 \leq j \leq q) be the distinct roots of the characteristic equation (40), and let njn_j be the multiplicity of the root rjr_j (1≤j≤q1 \leq j \leq q). Then the nn functions tkerjtt^k e^{r_j t} (1≤k≤nj, 1≤j≤q1 \leq k \leq n_j,\ 1 \leq j \leq q) are linearly independent integrals of the homogeneous equation (41).

One can prove this result directly in the following way. It follows from equation (41) that the nthn^{th} derivative of every integral of this equation is differentiable on R\mathbf{R}, from which one deduces immediately, by induction on the integer m>nm > n, that every integral of (41) admits a derivative of order mm, that is, is indefinitely differentiable on R\mathbf{R}. Let D\mathscr{D} be the (non-topological) vector space over C\mathbf{C} of indefinitely differentiable complex-valued functions on R\mathbf{R}; the map x↦Dxx \mapsto Dx is an endomorphism of this space, and equation (41) can be written

f(D)x=0(42) f(D)x = 0 \tag{42}

where f(D)=Dn−a1Dn−1−…−an−1D−anf(D) = D^n - a_1 D^{n-1} - \ldots - a_{n-1} D - a_n (Alg., IV, p. 4).

Proposition 10

Let g and h be two relatively prime polynomials such that f=ghf = gh. The subspace of solutions of (42) is the direct sum of the subspaces of solutions of the two equations

g(D)x=0,h(D)x=0. g(D)x = 0, \qquad h(D)x = 0.

Now, by the Bezout identity (Alg., VII. 2, th. 1) there exist two polynomials p(D)p(D) and q(D)q(D) such that p(D)g(D)+q(D)h(D)=1p(D)g(D) + q(D)h(D) = 1. For every solution xx of (42) one can write x=y+zx = y + z, where y=p(D)g(D)xy = p(D)g(D)x and z=q(D)h(D)xz = q(D)h(D)x; and then h(D)y=p(D)(f(D)x)=0h(D)y = p(D)(f(D)x) = 0 and g(D)z=q(D)(f(D)x)=0g(D)z = q(D)(f(D)x) = 0. On the other hand, if both g(D)x=0g(D)x = 0 and h(D)x=0h(D)x = 0, one deduces that

x=p(D)(g(D)x)+q(D)(h(D)x)=0, x = p(D)(g(D)x) + q(D)(h(D)x) = 0,

which completes the proof.

With the preceding notation one can write

f(D)=∏j=1q(D−rj)nj f(D) = \prod_{j=1}^q (D - r_j)^{n_j}

and prop. 10, by induction on qq, shows that the subspace of solutions of (42) is the direct sum of the subspaces of solutions of the qq equations

(D−rj)njx=0(1≤j≤q).(43) (D - r_j)^{n_j} x = 0 \qquad (1 \leq j \leq q). \tag{43}

Now, for every complex number rr one has

D(ertx)=ert(D+r)x(44) D(e^{rt} x) = e^{rt}(D + r)x \tag{44}

so (43) is equivalent to

Dnj(e−rjtx)=0 D^{n_j}(e^{-r_j t} x) = 0

and thus has as solutions the functions erjtpj(t)e^{r_j t} p_j(t), where pjp_j runs through the set of polynomials of degree ≤nj−1\leq n_j - 1; thus one recovers prop. 9 of IV, p. 194.

Assuming that the homogeneous equation (41) has been solved (that is, assuming that the characteristic roots have been found), one knows that the method of variation of parameters allows one to find the solutions of the inhomogeneous equation

Dnx−a1Dn−1x−…−an−1Dx−anx=b(t)(45) D^n x - a_1 D^{n-1} x - \ldots - a_{n-1} Dx - a_n x = b(t) \tag{45}

where b(t)b(t) is an arbitrary regulated function (IV, p. 182); we note that if b(t)b(t) is indefinitely differentiable on an interval J then all the integrals of (45) are indefinitely differentiable on J. In the particular case b(t)=eαtp(t)b(t) = e^{\alpha t} p(t), where pp is a polynomial (with complex coefficients) and α\alpha is an arbitrary complex number, one obtains an integral of (45) more simply by the following method. Put x=eαtyx = e^{\alpha t} y; the equation

f(D)x=eαtp(t) f(D)x = e^{\alpha t} p(t)

can, by (44), be written

f(α+D)y=p(t) f(\alpha + D)y = p(t)

or, by Taylor’s formula applied to the polynomial f(D)f(D),

f(n)(α)n!Dny+f(n−1)(α)(n−1)!Dn−1y+⋯+f′(α)1!Dy+f(α)y=p(t).(46) \frac{f^{(n)}(\alpha)}{n!} D^n y + \frac{f^{(n-1)}(\alpha)}{(n-1)!} D^{n-1} y + \cdots + \frac{f'(\alpha)}{1!} Dy + f(\alpha) y = p(t). \tag{46}

Let mm be the degree of the polynomial p(t)=∑k=0mλktm−kp(t) = \sum_{k=0}^m \lambda_k t^{m-k}; if f(α)≠0f(\alpha) \neq 0 (that is, if α\alpha is not a characteristic root), there exists one and only one polynomial u(t)=∑k=0mcktm−ku(t) = \sum_{k=0}^m c_k t^{m-k} of degree mm which is a solution of (46), for the coefficients ckc_k are determined by the system of linear equations

f(α)ck+(m−k+11)f′(α)ck−1+(m−k+22)f′′(α)ck−2+⋯+(mk)f(k)(α)c0=λk(0≤k≤m) f(\alpha)c_k + \binom{m-k+1}{1} f'(\alpha)c_{k-1} + \binom{m-k+2}{2} f''(\alpha)c_{k-2} + \cdots + \binom{m}{k} f^{(k)}(\alpha)c_0 = \lambda_k \qquad (0 \leq k \leq m)

which clearly admits one and only one solution. If, on the other hand, α\alpha is a characteristic root, and hh is its multiplicity, the preceding calculation shows that there is one and only one polynomial of degree mm such that every solution of Dhy=v(t)D^h y = v(t) is an integral; in other words, every polynomial solution of (46) is then of degree m+hm + h ("resonance").

9. Systems of linear equations with constant coefficients

With the notation of n° 8, let us consider more generally a system of mm differential equations of the form

∑k=1npjk(D)xk=bj(t)(1≤j≤m)(47) \sum_{k=1}^n p_{jk}(D) x_k = b_j(t) \qquad (1 \leq j \leq m) \tag{47}

where the unknowns xkx_k (1≤k≤n1 \leq k \leq n) and the right-hand sides bjb_j (1≤j≤m1 \leq j \leq m) are complex functions of the real variable tt, and where the pjk(D)p_{jk}(D) are polynomials (of any degree) with constant (complex) coefficients in the differentiation operator DD (1≤j≤m, 1≤k≤n1 \leq j \leq m,\ 1 \leq k \leq n).

Such systems are not of the same type as those considered in IV, p. 1164 (formula (5)), as the following example shows:

{Dx1=a(t)D2x1+Dx2+x2=b(t). \begin{cases} D x_1 = a(t) \\ D^2 x_1 + D x_2 + x_2 = b(t). \end{cases}

We shall confine ourselves to the case where the functions bj(t)b_j(t) are indefinitely differentiable on the interval J, and we shall look only for solutions (xk)1≤k≤n(x_k)_{1 \leq k \leq n} which are indefinitely differentiable on J. On putting b(t)=(b1(t),…,bm(t))\mathbf{b}(t) = (b_1(t), \ldots, b_m(t)), (a map from J into Cm\mathbf{C}^m), and x=(x1,x2,…,xn)\mathbf{x} = (x_1, x_2, \ldots, x_n), the system (47) can be written as

P(D)x=b(t) P(D) \mathbf{x} = \mathbf{b}(t)

where P(D)P(D) is the matrix (pjk(D))(p_{jk}(D)) with m rows and n columns, whose coefficients belong to the ring C[D]\mathbf{C}[D] of polynomials in D, with coefficients in C\mathbf{C}. Let fj(D)f_j(D) (1≤j≤r≤Min⁡(m,n)1 \leq j \leq r \leq \operatorname{Min}(m, n)) be the nonzero similarity invariants of the matrix P(D)P(D); one knows (Alg., VII, p. 32) that these are well-determined monic polynomials, such that fjf_j divides fj+1f_{j+1} for 1≤j≤r−11 \leq j \leq r - 1 (rr being the rank of P(D)P(D)); further, there are two square matrices U(D)U(D) and V(D)V(D) of order m and n respectively, invertible (in the rings of square matrices of order m and n respectively, with coefficients in the ring C[D]\mathbf{C}[D] of polynomials in D with complex coefficients), and such that all the entries of the matrix Q(D)=(qjk(D))=U(D)P(D)V(D)Q(D) = (q_{jk}(D)) = U(D) P(D) V(D) are zero, apart from the diagonal terms qjj(D)=fj(D)q_{jj}(D) = f_j(D) for 1≤j≤r1 \leq j \leq r. Now we put y=V−1(D)x\mathbf{y} = V^{-1}(D) \mathbf{x}; the equation (49) is equivalent to the equation

U(D)(P(D)(V(D)y))=U(D)b U(D) (P(D)(V(D) \mathbf{y})) = U(D) \mathbf{b}

that is, to

Q(D)y=U(D)b Q(D) \mathbf{y} = U(D) \mathbf{b}

since U(D)U(D) is invertible. Now, if y=(y1,y2,…,yn)\mathbf{y} = (y_1, y_2, \ldots, y_n), and if

U(D)b(t)=(c1(t),…,cm(t)), U(D) \mathbf{b}(t) = (c_1(t), \ldots, c_m(t)),

then equation (50) can be written

fj(D)yj=cj(t)for 1≤j≤r(51) f_j(D) y_j = c_j(t) \quad \text{for } 1 \leq j \leq r \tag{51}
0=cj(t)for r+1≤j≤m.(52) 0 = c_j(t) \quad \text{for } r + 1 \leq j \leq m. \tag{52}

The system thus does not admit any indefinitely differentiable solutions unless the conditions (52) are satisfied; the determination of the yjy_j for indices j≤rj \leq r then reduces to integrating the r linear differential equations with constant coefficients (51); the yjy_j for indices >r> r are arbitrary indefinitely differentiable functions. Once the solutions y\mathbf{y} of (50) have thus been determined, one deduces the solutions of (47) from the formula x=V(D)y\mathbf{x} = V(D) \mathbf{y}.

Remark 1

Some of the polynomials fj(D)f_j(D) may reduce to nonzero constants; the corresponding yjy_j are then completely determined.

Remark 2

When the bjb_j are all zero, that is, if the system (47) is homogeneous, the conditions (52) are always satisfied; if, further, r=nr = n, one sees that the set of solutions of (47) is a vector space over C\mathbf{C}, of dimension equal to the sum of the degrees of the fj(D)f_j(D), that is, to the degree of det⁡(P(D))\det(P(D)).

Remark 3

Given the polynomials pjk(D)p_{jk}(D), a system (47) which admits solutions when the right-hand sides are indefinitely differentiable (or differentiable of a certain order) may not admit them when the right hand sides are arbitrary regulated functions: this is shown by the example (48), which has no solution when a(t)a(t) is not a primitive. We shall not undertake here to seek for supplementary possibility conditions which enter when the right-hand sides are arbitrary regulated functions.

Exercises

See the exercises for § 2.

[^1]: Recall that if EE is of finite dimension then every affine linear map from EE into itself is continuous (Gen. Top., VI, p. 33 and 37).

Transcribed from the printed volume. Not translated.