In this chapter we go to work finding solutions to some important differential equations, using for this purpose the processes shown in the preceding chapters.

The beginner, who now knows how easy most of those processes are in themselves, will here begin to realize that integration is *an art*. As in all arts, so in this, facility can be acquired only by diligent and regular practice. He who would attain that facility must work out examples, and more examples, and yet more examples, such as are found abundantly in all the regular treatises on the Calculus. Our purpose here must be to afford the briefest introduction to serious work.

Transposing we have \[b \frac{dy}{dx} = -ay.\]

Now the mere inspection of this relation tells us that we have got to do with a case in which \(\dfrac{dy}{dx}\) is proportional to \(y\). If we think of the curve which will represent \(y\) as a function of \(x\), it will be such that its slope at any point will be proportional to the ordinate at that point, and will be a negative slope if \(y\) is positive. So obviously the curve will be a die-away curve (see Chapter 14), and the solution will contain \(e^{-x}\) as a factor. But, without presuming on this bit of sagacity, let us go to work.

As both \(y\) and \(dy\) occur in the equation and on opposite sides, we can do nothing until we get both \(y\) and \(dy\) to one side, and \(dx\) to the other. To do this, we must split our usually inseparable companions \(dy\) and \(dx\) from one another. \[\frac{dy}{y} = – \frac{a}{b}\, dx.\]

Having done the deed, we now can see that both sides have got into a shape that is integrable, because we recognize \(\dfrac{dy}{y}\), or \(\dfrac{1}{y}\, dy\), as a differential that we have met with (Chapter 14) when differentiating logarithms. So we may at once write down the instructions to integrate, \[\int \frac{dy}{y} = \int -\frac{a}{b}\, dx;\] and doing the two integrations, we have: \[\ln |y| = -\frac{a}{b} x + A,\] where \( A\) is the yet undetermined constant of integration. Then, delogarizing, we get

\[y = C e^{-\frac{a}{b} x},\qquad (\text{with } C=\pm e^A \text{ or }0)\]

which is *the solution* required. Now, this solution looks quite unlike the original differential equation from which it was constructed: yet to an expert mathematician they both convey the same information as to the way in which \(y\) depends on \(x\).

Now, as to the \(C\), its meaning depends on the initial value of \(y\). For if we put \(x = 0\) in order to see what value \(y\) then has, we find that this makes \(y = C e^{-0}\); and as \(e^{-0} = 1\) we see that \(C\) is nothing else than the particular value^{1} of \(y\) at starting. This we may call \(y_0\), and so write the solution as \[y = y_0 e^{-\frac{a}{b} x}.\]

Now we have done our best to get nothing but \(y\) and \(dy\) on one side, and nothing but \(dx\) on the other. But is the result on the left side integrable?

It is of the same form as the result on ; so, writing the instructions to integrate, we have: \[\int{\frac{dy}{y-\dfrac{g}{a}}} = – \int{\frac{a}{b}\, dx};\] and, doing the integration, and adding the appropriate constant,

\[\ln\left|y-\frac{g}{a}\right| = -\frac{a}{b}x + A\] \[y-\dfrac{g}{a} = Ce^{-\dfrac{a}{b}x}; \quad \text{where } C=\pm e^A \text{ or } 0\] \[y = \frac{g}{a} + Ce^{-\frac{a}{b}x},\]

which is *the solution*.

If the condition is laid down that \(y = 0\) when \(x = 0\) we can find \(C\); for then the exponential becomes \(= 1\); and we have

\[\begin{aligned} 0 &= \frac{g}{a} + C, \\C &= -\frac{g}{a}. \end{aligned}\]

Putting in this value, the solution becomes \[y = \frac{g}{a} (1-e^{-\frac{a}{b} x}).\]

But further, if \(x\) grows indefinitely, \(y\) will grow to a maximum; for when \(x=\infty\), the exponential \(= 0\), giving \(y_{\text{max.}} = \dfrac{g}{a}\). Substituting this, we get finally \[y = y_{\text{max.}}(1-e^{-\frac{a}{b} x}).\]

This result is also of importance in physical science.

Let \[ ay+b\frac{dy}{dt} = g \cdot \sin 2\pi nt.\]

We shall find this much less tractable than the preceding. First divide through by \(b\). \[\frac{dy}{dt} + \frac{a}{b}y = \frac{g}{b} \sin 2\pi nt.\]

Now, as it stands, the left side is not integrable. But it can be made so by the artifice—and this is where skill and practice suggest a plan—of multiplying all the terms by \(e^{\frac{a}{b} t}\), giving us: \[\frac{dy}{dt} e^{\frac{a}{b} t} + \frac{a}{b} y e^{\frac{a}{b} t} = \frac{g}{b} e^{\frac{a}{b} t} \cdot \sin 2 \pi nt,\] which is the same as \[\frac{dy}{dt} e^{\frac{a}{b} t} + y \frac{d(e^{\frac{a}{b} t})}{dt} = \frac{g}{b} e^{\frac{a}{b} t} \cdot \sin 2 \pi nt;\] and this being a perfect differential may be integrated thus:—since, if \(u = ye^{\frac{a}{b} t}\), \(\dfrac{du}{dt} = \dfrac{dy}{dt} e^{\frac{a}{b} t} + y \dfrac{d(e^{\frac{a}{b} t})}{dt}\),

\[\begin{align} y e^{\frac{a}{b} t} &= \frac{g}{b} \int e^{\frac{a}{b} t} \cdot \sin 2 \pi nt \cdot dt + C, \\ y &= \frac{g}{b} e^{-\frac{a}{b} t} \int e^{ \frac{a}{b} t} \cdot \sin 2\pi nt \cdot dt + Ce^{-\frac{a}{b} t}. \tag{A}\end{align}\]

The last term is obviously a term which will die out as \(t\) increases, and may be omitted. The trouble now comes in to find the integral that appears as a factor. To tackle this we resort to the device (see Chapter 20) of integration by parts, the general formula for which is \(\int u dv = uv – \int v du\). For this purpose write \[\begin{aligned} &\left\{ \begin{aligned} u &= e^{\frac{a}{b} t}; \\ dv &= \sin 2\pi nt \cdot dt. \end{aligned} \right.\end{aligned}\] We shall then have \[\begin{aligned} &\left\{ \begin{aligned} du &= e^{\frac{a}{b} t} \times \frac{a}{b}\, dt; \\ v &= – \frac{1}{2\pi n} \cos 2\pi nt. \end{aligned} \right.\end{aligned}\]

Inserting these, the integral in question becomes: \[\begin{align} \int e^{\frac{a}{b} t} &{} \cdot \sin 2 \pi n t \cdot dt \\ &= -\frac{1}{2 \pi n} \cdot e^{\frac{a}{b} t} \cdot \cos 2 \pi nt -\int -\frac{1}{2\pi n} \cos 2 \pi nt \cdot e^{\frac{a}{b} t} \cdot \frac{a}{b}\, dt \\ &= -\frac{1}{2 \pi n} e^{\frac{a}{b} t} \cos 2 \pi nt +\frac{a}{2 \pi nb} \int e^{\frac{a}{b} t} \cdot \cos 2 \pi nt \cdot dt. \tag {B}\end{align}\]

The last integral is still irreducible. To evade the difficulty, repeat the integration by parts of the left side, but treating it in the reverse way by writing:

\[\begin{aligned} &\left\{ \begin{aligned} u &= \sin 2 \pi n t ; \\ dv &= e^{\frac{a}{b} t} \cdot dt; \end{aligned} \right. \\ &\left\{ \begin{aligned} du &= 2 \pi n \cdot \cos 2 \pi n t \cdot dt; \\ v &= \frac{b}{a} e ^{\frac{a}{b} t} \end{aligned} \right. \end{aligned}\]

Inserting these, we get \[\begin{align} \int e^{\frac{a}{b} t} &{} \cdot \sin 2 \pi n t \cdot dt\\ &= \frac{b}{a}\cdot e^{\frac{a}{b} t} \cdot \sin 2 \pi n t – \frac{2 \pi n b}{a} \int e^{\frac{a}{b} t} \cdot \cos 2 \pi n t \cdot dt. \tag{C} \end{align}\]

Noting that the final intractable integral in (C) is the same as that in (B), we may eliminate it, by multiplying (B) by \(\dfrac{2 \pi nb}{a}\), and multiplying (C) by \(\dfrac{a}{2 \pi nb}\), and adding them.

The result, when cleared down, is: \[\begin{align} \int e^{\frac{a}{b} t} \cdot \sin 2 \pi n t \cdot dt &= e^{\frac{a}{b} t} \left\{\frac{ ab \cdot \sin 2 \pi nt – 2 \pi n b^2 \cdot \cos 2 \pi n t}{ a^2 + 4 \pi^2 n^2 b^2 } \right\} \tag {D} &\end{align}\] Inserting this value in (A) , we get \[\begin{aligned} y &= g \left\{\frac{ a \cdot \sin 2 \pi n t – 2 \pi n b \cdot \cos 2 \pi nt}{ a^2 + 4 \pi^2 n^2 b^2}\right\}. &\end{aligned}\]

To simplify still further, let us imagine an angle \(\phi\) such that \(\tan \phi = \dfrac{2 \pi n b}{ a}\).

\[\begin{aligned} \sin \phi &= \frac{2 \pi nb}{\sqrt{a^2 + 4 \pi^2 n^2 b^2}}, \\ \cos \phi &= \frac{a}{\sqrt{a^2 + 4 \pi^2 n^2 b^2}}. \end{aligned}\]

Substituting these, we get:\[ y = g \frac{\cos \phi \cdot \sin 2 \pi nt – \sin \phi \cdot \cos 2 \pi nt}{\sqrt{a^2 + 4 \pi^2 n^2 b^2}},\] which may be written \[ y = g \frac{\sin(2 \pi nt – \phi)}{\sqrt{a^2 + 4 \pi^2 n^2 b^2}},\] which is *the solution* desired.

This is indeed none other than the equation of an alternating electric current, where \(g\) represents the amplitude of the electromotive force, \(n\) the frequency, \(a\) the resistance, \(b\) the coefficient of self-induction of the circuit, and \(\phi\) is an angle of lag.

We could integrate this expression directly, if \(M\) were a function of \(x\) only, and \(N\) a function of \(y\) only; but, if both \(M\) and \(N\) are functions that depend on both \(x\) and \(y\), how are we to integrate it? Is it itself an exact differential? That is: have \(M\) and \(N\) each been formed by partial differentiation from some common function \(U\), or not? If they have, then \[\left\{ \begin{aligned} \frac{\partial U}{\partial x} = M, \\ \frac{\partial U}{\partial y} = N. \end{aligned} \right.\] And if such a common function exists, then \[\frac{\partial U}{\partial x}\, dx + \frac{\partial U}{\partial y}\, dy\] is an exact differential (compare Chapter 16).

Now the test of the matter is this. If the expression is an exact differential, it must be true that

\[\frac{dM}{dy} = \frac{dN}{dx}; \] for then \[\frac{d(dU)}{dx\, dy} = \frac{d(dU)}{dy\, dx},\]

which is necessarily true.

Take as an illustration the equation \[(1 + 3 xy)\, dx + x^2\, dy = 0.\]

Is this an exact differential or not? Apply the test. \[\left\{ \begin{aligned} \frac{d(1 + 3xy)}{dy}=3x, \\ {\dfrac{d(x^2)}{dx}} = 2x, \end{aligned} \right.\] which do not agree. Therefore, it is not an exact differential, and the two functions \(1+3xy\) and \(x^2\) have not come from a common original function.

It is possible in such cases to discover, however, *an integrating factor*, that is to say, a factor such that if both are multiplied by this factor, the expression will become an exact differential. There is no one rule for discovering such an integrating factor; but experience will usually suggest one. In the present instance \(2x\) will act as such. Multiplying by \(2x\), we get \[(2x + 6x^2y)\, dx + 2x^3\, dy = 0.\]

Now apply the test to this. \[\left\{ \begin{aligned} \frac{d(2x + 6x^2y)}{dy}=6x^2, \\ {\dfrac{d(2x^3)}{dx}} = 6x^2, \end{aligned} \right.\] which agrees. Hence this is an exact differential, and may be integrated. Now, if \(w = 2x^3y\),

\[\begin{aligned} dw=6x^2y\, dx + 2x^3\, dy. \\ \int 6x^2y\, dx + \int 2x^3\, dy=w=2x^3y; \\ U = x^2 + 2x^3y + C. \end{aligned}\]

In this case we have a differential equation of the second degree, in which \(y\) appears in the form of a second differential coefficient, as well as in person.

Transposing, we have \(\dfrac{d^2 y}{dt^2} = – n^2 y\).

It appears from this that we have to do with a function such that its second differential coefficient is proportional to itself, but with reversed sign. In Chapter XV. we found that there was such a function—namely, the *sine* (or the *cosine* also) which possessed this property. So, without further ado, we may infer that the solution will be of the form \(y = A \sin (pt + q)\). However, let us go to work.

Multiply both sides of the original equation by \(2\dfrac{dy}{dt}\) and integrate, giving us \[2\dfrac{d^2 y}{dt^2}\, \dfrac{dy}{dt} + 2x^2 y \dfrac{dy}{dt} = 0,\] and, as \[2 \frac{d^2y}{dt^2}\, \frac{dy}{dt} = \frac{d \left(\dfrac{dy}{dt}\right)^2}{dt},\quad \left(\frac{dy}{dt}\right)^2 + n^2 (y^2-C^2) = 0,\] \(C\) being a constant. Then, taking the square roots, \[\frac{dy}{dt} = -n \sqrt{ y^2 – C^2}\quad \text{and}\quad \frac{dy}{\sqrt{C^2 – y^2}} = n \cdot dt.\]

But it can be shown that (see Chapter 15) \[\frac{1}{\sqrt{C^2 – y^2}} = \frac{d (\arcsin \dfrac{y}{C})}{dy};\] whence, passing from angles to sines, \[\arcsin \frac{y}{C} = nt + C_1\quad \text{and}\quad y = C \sin (nt + C_1),\] where \(C_1\) is a constant angle that comes in by integration.

Or, preferably, this may be written \[y = A \sin nt + B \cos nt, \]which is the solution.

Here we have obviously to deal with a function \(y\) which is such that its second differential coefficient is proportional to itself. The only function we know that has this property is the exponential function (see Chapter 14), and we may be certain therefore that the solution of the equation will be of that form.

Proceeding as before, by multiplying through by \(2 \dfrac{dy}{dx}\), and integrating, we get \(2\dfrac{d^2 y}{dx^2}\, \dfrac{dy}{dx} – 2n^2 y \dfrac{dy}{dx}=0\), and as

\[\begin{aligned} 2\frac{d^2 y}{dx^2}\, \frac{dy}{dx} = \frac{d \left(\dfrac{dy}{dx}\right)^2}{dx},\quad \left(\frac{dy}{dx}\right)^2 – n^2 (y^2 + c^2) = 0, \\ \frac{dy}{dx} – n \sqrt{y^2 + c^2} = 0, \end{aligned}\]

where \(c\) is a constant, and \(\dfrac{dy}{\sqrt{y^2 + c^2}} = n\, dx\).

Now, if \(w = \ln ( y+ \sqrt{y^2+ c^2}) = \ln u\),

\[\begin{aligned} \frac{dw}{du} = \frac{1}{u},\quad \frac{du}{dy} = 1 + \frac{y}{\sqrt{y^2 + c^2}} = \frac{y + \sqrt{ y^2 + c^2}}{\sqrt{y^2 + c^2}} \\ \frac{dw}{dy} = \frac{1}{\sqrt{ y^2 + c^2}}.\end{aligned}\]

Hence, integrating, this gives us

\[\ln (y + \sqrt{y^2 + c^2} ) = nx + \ln C, \]\[ y + \sqrt{y^2 + c^2} = C e^{nx}. \tag{1}\] Now \[( y + \sqrt{y^2 + c^2} ) \times ( -y + \sqrt{y^2 + c^2} ) = c^2 ; \] \[ -y + \sqrt{y^2 + c^2} = \dfrac{c^2}{C} e^{-nx}. \tag{2}\]

Subtracting (2) from (1) and dividing by \(2\), we then have \[y = \frac{1}{2} C e^{nx} – \frac{1}{2}\, \frac{c^2}{C} e^{-nx},\] which is more conveniently written \[y = A e^{nx} + B e^{-nx}.\] Or, the solution, which at first sight does not look as if it had anything to do with the original equation, shows that \(y\) consists of two terms, one of which grows logarithmically as \(x\) increases, and of a second term which dies away as \(x\) increases.

Let \[ b \frac{d^2y}{dt^2} + a \frac{dy}{dt} + gy = 0.\]

Examination of this expression will show that, if \(b = 0\), it has the form of Example 1, the solution of which was a negative exponential. On the other hand, if \(a = 0\), its form becomes the same as that of Example 6, the solution of which is the sum of a positive and a negative exponential. It is therefore not very surprising to find that the solution of the present example is

\[ \begin{aligned} y &= (e^{-mt})(A e^{nt} + B e^{-nt}), \\ m &= \frac{a}{2b}\quad \text{and}\quad n = \sqrt{\frac{a^2}{4b^2}} – \frac{g}{b}. \end{aligned}\]

The steps by which this solution is reached are not given here; they may be found in advanced treatises.

It was seen (Chapter 15) that this equation was derived from the original \[y = F(x+at) + f(x-at),\] where \(F\) and \(f\) were any arbitrary functions of \(t\).

Another way of dealing with it is to transform it by a change of variables into \[\frac{d^2y}{du \cdot dv} = 0,\] where \(u = x + at\), and \(v = x – at\), leading to the same general solution. If we consider a case in which \(F\) vanishes, then we have simply \[y = f(x-at);\] and this merely states that, at the time \(t = 0\), \(y\) is a particular function of \(x\), and may be looked upon as denoting that the curve of the relation of \(y\) to \(x\) has a particular shape. Then any change in the value of \(t\) is equivalent simply to an alteration in the origin from which \(x\) is reckoned. That is to say, it indicates that, the form of the function being conserved, it is propagated along the \(x\) direction with a uniform velocity \(a\); so that whatever the value of the ordinate \(y\) at any particular time \(t_0\) at any particular point \(x_0\), the same value of \(y\) will appear at the subsequent time \(t_1\) at a point further along, the abscissa of which is \(x_0 + a(t_1 – t_0)\). In this case the simplified equation represents the propagation of a wave (of any form) at a uniform speed along the \(x\) direction.

If the differential equation had been written \[m \frac{d^2y}{dt^2} = k\, \frac{d^2y}{dx^2},\] the solution would have been the same, but the velocity of propagation would have had the value \[a = \sqrt{\frac{k}{m}}.\]

You have now been personally conducted over the frontiers into the enchanted land. And in order that you may have a handy reference to the principal results, the author, in bidding you farewell, begs to present you with a passport in the shape of a convenient collection of standard forms (see TABLE OF STANDARD FORMS.). In the middle column are set down a number of the functions which most commonly occur. The results of differentiating them are set down on the left; the results of integrating them are set down on the right. May you find them useful!