# 1911 Encyclopædia Britannica/Differences, Calculus of

See also Finite difference on Wikipedia; and our 1911 Encyclopædia Britannica disclaimer.
​DIFFERENCES, CALCULUS OF (Theory of Finite Differences), that branch of mathematics which deals with the successive differences of the terms of a series.
1. The most important of the cases to which mathematical methods can be applied are those in which the terms of the series are the values, taken at stated intervals (regular or irregular), of a continuously varying quantity. In these cases the formulae of finite differences enable certain quantities, whose exact value depends on the law of variation (i.e. the law which governs the relative magnitude of these terms) to be calculated, often with great accuracy, from the given terms of the series, without explicit reference to the law of variation itself. The methods used may be extended to cases where the series is a double series (series of double entry), i.e. where the value of each term depends on the values of a pair of other quantities.
2. The first differences of a series are obtained by subtracting from each term the term immediately preceding it. If these are treated as terms of a new series, the first differences of this series are the second differences of the original series; and so on. The successive differences are also called differences of the first, second, ... order. The differences of successive orders are most conveniently arranged in successive columns of a table thus:—
Algebra of Differences and Sums.
* Fig. 1.
3. The formal relations between the terms of the series and the differences may be seen by comparing the arrangements (A) and (B) in fig. 1. In (A) the various terms and differences are the same as in § 2, but placed differently. In (B) we take a new series of terms α, β, γ, δ, commencing with the same term α, and take the successive sums of pairs of terms, instead of the successive differences, but place them to the left instead of to the right. It will be seen, in the first place, that the successive terms in (A), reading downwards to the right, and the successive terms in (B), reading downwards to the left, consist each of a series of terms whose coefficients follow the binomial law; i.e. the coefficients in b − a, c − 2b + a, d − 3c + 3b − a, ... and in α + β, α + 2β + γ, α + 3β + 3γ + δ, ... are respectively the same as in y − x, (y − x)², (y − x)³, ... and in x + y, (x + y)², (x + y)³,.... In the second place, it will be seen that the relations between the various terms in (A) are identical with the relations between the similarly placed terms in (B); e.g. β + γ is the difference of α + 2β + γ and α + β, just as c − b is the difference of c and b: and d − c is the sum of c − b and d − 2c + b, just as β + 2γ + δ is the sum of β + γ and γ + δ. Hence if we take β, γ, δ, ... of (B) as being the same as b − a, c − 2b + a, d − 3c + 3b − a, ... of (A), all corresponding terms in the two diagrams will be the same.
Thus we obtain the two principal formulae connecting terms and differences. If we provisionally describe b − a, c − 2b + a, ... as the first, second, ... differences of the particular term a (§ 7), then (i.) the nth difference of a is
$l-nk+\ldots +(-1)^{n-2}{\frac {n\cdot n-1}{1\cdot 2}}c+(-1)^{n-1}nb+(-1)^{n}a,$
where l, k ... are the (n + 1)th, nth, ... terms of the series a, b, c, ...; the coefficients being those of the terms in the expansion of (y − x)$^n$: and (ii.) the (n + 1)th term of the series, i.e. the nth term after a, is
$a+n\beta +{\frac {n\cdot n-1}{1\cdot 2}}\gamma +\ldots$
where β, γ, ... are the first, second, ... differences of a; the coefficients being those of the terms in the expansion of (x + y)$^n$.
4. Now suppose we treat the terms a, b, c, ... as being themselves the first differences of another series. Then, if the first term of this series is N, the subsequent terms are N + a, N + a + b, N + a + b + c, ...; i.e. the difference between the (n + 1)th term and the first term is the sum of the first n terms of the original series. The term N, in the diagram (A), will come above and to the left of a; and we see, by (ii.) of § 3, that the sum of the first n terms of the original series is
$\left(N+na+{\frac {n\cdot n-1}{1\cdot 2}}\beta +\ldots \right)-N{=}na+{\frac {n\cdot n-1}{1\cdot 2}}\beta +{\frac {n\cdot n-1\cdot n-2}{1\cdot 2\cdot 3}}\gamma +\ldots$
5. As an example, take the arithmetical series
a, a + p, a + 2p, ...
The first differences are p, p, p, ... and the differences of any higher order are zero. Hence, by (ii.) of § 3, the (n + 1)th term is a + np, and, by § 4, the sum of the first n terms is na + ½n(n − 1)p = ½n{2a + (n − 1)p}.
6. As another example, take the series 1, 8, 27, ... the terms of which are the cubes of 1, 2, 3, ... The first, second and third differences of the first term are 7, 12 and 6, and it may be shown (§ 14 (i.)) that all differences of a higher order are zero. Hence the sum of the first n terms is
$n+7{\frac {n\cdot n-1}{1\cdot 2}}+12{\frac {n\cdot n-1\cdot n-2}{1\cdot 2\cdot 3}}+6{\frac {n\cdot n-1\cdot n-2\cdot n-3}{1\cdot 2\cdot 3\cdot 4}}{=}{\tfrac {1}{4}}n^{4}+{\tfrac {1}{2}}n^{3}+{\tfrac {1}{4}}n^{2}{=}{{\tfrac {1}{2}}n(n+1)}^{2}$
7. In § 3 we have described b − a, c − 2b + a, ... as the first, second, ... differences of a. This ascription of the differences to particular terms of the series is quite arbitrary. If we read the differences in the table of § 2 upwards to the right instead of downwards to the right, we might describe e − d, e − 2d + c, ... as the first, second, ... differences of e. On the other hand, the term of greatest weight in c − 2b + a, i.e. the term which has the numerically greatest coefficient, is b, and therefore c − 2b + a might properly be regarded as the second difference of b, and similarly e − 4d + 6c − 4b + a might be regarded as the fourth difference of c. These three methods of regarding the differences lead to three different systems of notation, which are described in §§ 9, 10 and 11.
Notation of Differences and Sums.
8. It is convenient to denote the terms a, b, c, ... of the series by u$_0$, u$_1$, u$_2$, u$_3$, ... If we merely have the terms of the series, u$_n$ may be regarded as meaning the (n + 1)th term. Usually, however, the terms are the values of a quantity u, which is a function of another quantity x, and the values of x, to which a, b, c, ... correspond, proceed by a constant difference h. If x$_0$ and u$_0$ are a pair of corresponding values of x and u, and if any other value x$_0$ + mh of x and the corresponding value of u are denoted by x$_m$ and u$_m$, then the terms of the series will be ... u$_{n-2}$, u$_{n−1}$, u$_n$, u$_{n+1}$, u$_{n+2}$ ..., corresponding to values of x denoted by ... x$_{n-2}$, x$_{n−1}$, x$_n$, x$_{n+1}$, x$_{n+2}$....
9. In the advancing-difference notation u$_{n+1}$ − u$_n$ is denoted by Δu$_n$. The differences Δu$_0$, Δu$_1$, Δu$_2$ ... may then be regarded as values of a function Δu corresponding to values of x proceeding by constant difference h; and therefore Δu$_{n+1}$ − Δu$_n$ denoted by ΔΔu$_n$, or, more briefly, Δ²u$_n$; and so on. Hence the table of differences in § 2, with the corresponding values of x and of u placed opposite each other in the ordinary manner of mathematical tables, becomes
| x | u | 1st Diff. | 2nd Diff. | 3rd Diff. | 4th Diff. |
| x$_{n-2}$ | u$_{n-2}$ | | Δ²u$_{n-3}$ | | Δ$^4$u$_{n-4}$ ... |
| | | Δu$_{n-2}$ | | Δ³u$_{n-3}$ | |
| x$_{n−1}$ | u$_{n−1}$ | | Δ²u$_{n-2}$ | | Δ$^4$u$_{n-3}$ ... |
| | | Δu$_{n−1}$ | | Δ³u$_{n-2}$ | |
| x$_n$ | u$_n$ | | Δ²u$_{n−1}$ | | Δ$^4$u$_{n-2}$ ... |
| | | Δu$_n$ | | Δ³u$_{n−1}$ | |
| x$_{n+1}$ | u$_{n+1}$ | | Δ²u$_n$ | | Δ$^4$u$_{n−1}$ ... |
| | | Δu$_{n+1}$ | | Δ³u$_n$ | |
| x$_{n+2}$ | u$_{n+2}$ | | Δ²u$_{n+1}$ | | Δ$^4$u$_n$  ... |
The terms of the series of which ... u$_{n−1}$, u$_n$, u$_{n+1}$, ... are the first differences are denoted by Σu, with proper suffixes, so ​that this series is ... Σu$_{n-1}$, Σu$_n$, Σu$_{n+1}$.... The suffixes are chosen so that we may have ΔΣu$_n$ = u$_n$, whatever n may be; and therefore (§ 4) Σu$_n$ may be regarded as being the sum of the terms of the series up to and including u$_{n-1}$. Thus if we write Σu$_{n-1}$ = C + u$_{n-2}$, where C is any constant, we shall have
* Σu$_n$ = Σu$_{n-1}$ + ΔΣu$_{n-1}$ = C + u$_{n-2}$ + u$_{n-1}$,
* Σu$_{n+1}$ = C + u$_{n-2}$ + u$_{n-1}$ + u$_n$,
and so on. This is true whatever C may be, so that the knowledge of ... u$_{n-1}$, u$_n$, ... gives us no knowledge of the exact value of Σu$_n$; in other words, C is an arbitrary constant, the value of which must be supposed to be the same throughout any operations in which we are concerned with values of Σu corresponding to different suffixes.
There is another symbol E, used in conjunction with u to denote the next term in the series. Thus Eu$_n$ means u$_{n+1}$, so that Eu$_n$ = u$_n$ + Δu$_n$.
10. Corresponding to the advancing-difference notation there is a receding-difference notation, in which u$_{n+1}$ − u$_n$ is regarded as a difference of u$_{n+1}$, and may be denoted by Δ′u$_{n+1}$, and similarly u$_{n+1}$ − 2u$_n$ + u$_{n-1}$ may be denoted by Δ′²u$_{n+1}$. This notation is only required for certain special purposes, and the usage is not settled (§ 19 (ii.)).
11. The central-difference notation depends on treating u$_{n+1}$ − 2u$_n$ − u$_{n-1}$ as the second difference of u$_n$, and therefore as corresponding to the value x$_n$; but there is no settled system of notation. The following seems to be the most convenient. Since u$_n$ is a function of x$_n$, and the second difference u$_{n+2}$ − 2u$_{n+1}$ + u$_n$ is a function of x$_{n+1}$, the first difference u$_{n+1}$ − u$_n$ must be regarded as a function of x$_{n+1/2}$, i.e. of ½(x$_n$ + x$_{n+1}$). We therefore write u$_{n+1}$ − u$_n$ = δu$_{n+1/2}$, and each difference in the table in § 9 will have the same suffix as the value of x in the same horizontal line; or, if the difference is of an odd order, its suffix will be the means of those of the two nearest values of x. This is shown in the table below.
In this notation, instead of using the symbol E, we use a symbol μ to denote the mean of two consecutive values of u, or of two consecutive differences of the same order, the suffixes being assigned on the same principle as in the case of the differences. Thus
μu$_{n+1/2}$ = ½(u$_n$ + u$_{n+1}$, μδu$_n$ = ½(δu$_{n-1/2}$ + δu$_{n+1/2}$, &c.
If we take the means of the differences of odd order immediately above and below the horizontal line through any value of x, these means, with the differences of even order in that line, constitute the central differences of the corresponding value of u. Thus the table of central differences is as follows, the values obtained as means being placed in brackets to distinguish them from the actual differences:—
| x | u | 1st Diff. | 2nd Diff. | 3rd Diff. | 4th Diff. |
| x$_{n-2}$ | u$_{n-2}$ | (μδu$_{n-2}$) | δ²u$_{n-2}$ | (μδ³u$_{n-2}$) | δ$^4$u$_{n-2}$ ... |
| | | δu$_{n-3/2}$ | | δ³u$_{n-3/2}$ | |
| x$_{n-1}$ | u$_{n-1}$ | (μδu$_{n-1}$) | δ²u$_{n-1}$ | (μδ³u$_{n-1}$) | δ$^4$u$_{n-1}$ ... |
| | | δu$_{n-1/2}$ | | δ³u$_{n-2}$ | |
| x$_n$ | u$_n$ | (μδu$_n$) | δ²u$_n$ | (μδ³u$_n$) | δ$^4$u$_n$  ... |
| | | δu$_{n+1/2}$ | | δ³u$_{n+1/2}$ | |
| x$_{n+1}$ | u$_{n+1}$ | (μδu$_{n+1}$) | δ²u$_{n+1}$ | (μδ³u$_{n+1}$) | δ$^4$u$_{n+1}$ ... |
| | | δu$_{n+3/2}$ | | δ³u$_{n+3/2}$ | |
| x$_{n+2}$ | u$_{n+2}$ | (μδu$_{n+2}$) | δ²u$_{n+2}$ | (μδ³u$_{n+2}$) | δ$^4$u$_{n+2}$ ... |
Similarly, by taking the means of consecutive values of u and also of consecutive differences of even order, we should get a series of terms and differences central to the intervals x$_{n-2}$ to x$_{n-1}$, x$_{n-1}$ to x$_n$, ....
The terms of the series of which the values of u are the first differences are denoted by σu, with suffixes on the same principle; the suffixes being chosen so that δσu$_n$ shall be equal to u$_n$. Thus, if
σu$_{n-3/2}$ = C + u$_{n-2}$,
then
σu$_{n-1/2}$ = C + u$_{n-2}$ + u$_{n-1}$, σ$_{n+1/2}$ = C + u$_{n-2}$ + u$_{n-1}$ + u$_n$, &c.,
and also
μσu$_{n-1}$ = C + u$_{n-2}$ + ½u$_{n-1}$, μσu$_n$ = C + u$_{n-2}$ + u$_{n-1}$ + ½u$_n$, &c.,
C being an arbitrary constant which must remain the same throughout any series of operations.
Operators and Symbolic Methods.
12. There are two further stages in the use of the symbols Δ, Σ, δ, σ, &c., which are not essential for elementary treatment but lead to powerful methods of deduction.
(i.) Instead of treating Δu as a function of x, so that Δu$_n$ means (Δu)$_n$, we may regard Δ as denoting an operation performed on u, and take Δu$_n$ as meaning Δ.u$_n$. This applies to the other symbols E, δ, &c., whether taken simply or in combination. Thus ΔEu$_n$ means that we first replace u$_n$ by u$_{n+1}$, and then replace this by u$_{n+2}$ − u$_{n+1}$.
(ii.) The operations Δ, E, δ, and μ, whether performed separately or in combination, or in combination also with numerical multipliers and with the operation of differentiation denoted by D (≡ d/dx), follow the ordinary rules of algebra: e.g. Δ(u$_n$ + v$_n$) = Δu$_n$ + Δv$_n$, ΔDu$_n$ = DΔu$_n$, &c. Hence the symbols can be separated from the functions on which the operations are performed, and treated as if they were algebraical quantities. For instance, we have
E·u$_n$ = u$_{n+1}$ = u$_n$ + Δu$_n$ = 1·u$_n$ + Δ·u$_n$,
so that we may write E = 1 + Δ, or Δ = E − 1. The first of these is nothing more than a statement, in concise form, that if we take two quantities, subtract the first from the second, and add the result to the first, we get the second. This seems almost a truism. But, if we deduce E$^n$ = (1 + Δ)$^n$, Δ$^n$ = (E-1)$^n$, and expand by the binomial theorem and then operate on u$_0$, we get the general formulae
$u_{n}{=}u_{0}+n\Delta u_{0}+{\frac {n\cdot n-1}{1\cdot 2}}\Delta ^{2}u_{0}+\ldots +\Delta ^{n}u_{0},$
$\Delta ^{n}u_{0}{=}u_{n}-nu_{n-1}+{\frac {n\cdot n-1}{1\cdot 2}}u_{n-2}+\ldots +(-1)^{n}u_{0},$
which are identical with the formulae in (ii.) and (i.) of § 3.
(iii.) What has been said under (ii.) applies, with certain reservations, to the operations Σ and σ, and to the operation which represents integration. The latter is sometimes denoted by D$^{-1}$; and, since ΔΣu$_n$ = u$_n$, and δσu$_n$ = u$_n$, we might similarly replace Σ and σ by Δ$^{-1}$ and δ$^{-1}$. These symbols can be combined with Δ, E, &c. according to the ordinary laws of algebra, provided that proper account is taken of the arbitrary constants introduced by the operations D$^{-1}$, Δ$^{-1}$, δ$^{-1}$.
Applications to Algebraical Series.
13. Summation of Series.—If u$_r$, denotes the (r + 1)th term of a series, and if v$_r$ is a function of r such that Δv$_r$ = u$_r$ for all integral values of r, then the sum of the terms u$_m$, u$_{m+1}$, ... u$_n$ is v$_{n+1}$ − v$_m$. Thus the sum of a number of terms of a series may often be found by inspection, in the same kind of way that an integral is found.
14. Rational Integral Functions.—(i.) If u$_r$ is a rational integral function of r of degree p, then Δu$_r$, is a rational integral function of r of degree p − 1.
(ii.) A particular case is that of a factorial, i.e. a product of the form (r + a + 1) (r + a + 2) ... (r + b), each factor exceeding the preceding factor by 1. We have
Δ · (r + a + 1) (r + a + 2) ... (r + b) = (b − a)·(r + a + 2) ... (r + b),
whence, changing a into a-1,
Σ(r + a + 1) (r + a + 2) ... (r + b) = const. + (r + a)(r + a + 1) ... (r + b)/(b − a + 1).
A similar method can be applied to the series whose (r + 1)th term is of the form 1/(r + a + 1) (r + a + 2) ... (r + b).
(iii.) Any rational integral function can be converted into the sum of a number of factorials; and thus the sum of a series of which such a function is the general term can be found. For example, it may be shown in this way that the sum of the pth powers of the first n natural numbers is a rational integral function of n of degree p + 1, the coefficient of n$^{p+1}$ being 1/(p + 1).
15. Difference-equations.—The summation of the series ... + u$_{n+2}$ + u$_{n-1}$ + u$_n$ is a solution of the difference-equation Δv$_n$ = u$_{n+1}$, which may also be written (E − 1)v$_n$ = u$_{n+1}$. This is a simple form of difference-equation. There are several forms which have been investigated; a simple form, more general than the above, is the linear equation with constant coefficients—
v$_{n+m}$ + a$_1$v$_{n+m-1}$ + a$_2$v$_{n+m-2}$ + ... + a$_m$v$_n$ = N,
where a$_1$, a$_2$, ... a$_m$ are constants, and N is a given function of n. This may be written
(E$^m$ + a$_1$E$^{m-1}$ + ... + a$_m$)v$_n$ = N
or
(E − p$_1$)(E − p$_2$) ... (E − p$_m$)v$_n$ = N.
The solution, if p$_1$, p$_2$, ... p$_m$ are all different, is v$_n$ = C$_1$p$_1$$^n$ + C$_2$p$_2$$^n$ + ... + C$_m$p$_m$$^n$ + V$_n$, where C$_1$, C$_2$ ... are constants, and v$_n$ = V$_n$ is any one solution of the equation. The method of finding a value for V$_n$ depends on the form of N. Certain modifications are required when two or more of the p’s are equal.
It should be observed, in all cases of this kind, that, in describing C$_1$, C$_2$ as “constants,” it is meant that the value of any one, as C$_1$, is the same for all values of n occurring in the series. A “constant” may, however, be a periodic function of n.
Applications to Continuous Functions.
16. The cases of greatest practical importance are those in which u is a continuous function of x. The terms u$_1$, u$_2$ ... of the series then represent the successive values of u corresponding to x = x$_1$, x$_2$.... The important applications of the theory in these cases are to (i.) relations between differences and differential coefficients, (ii.) ​interpolation, or the determination of intermediate values of u, and (iii.) relations between sums and integrals.
17. Starting from any pair of values x$_0$ and u$_0$, we may suppose the interval h from x$_0$ to x$_1$ to be divided into q equal portions. If we suppose the corresponding values of u to be obtained, and their differences taken, the successive advancing differences of u$_0$ being denoted by ∂u$_0$, ∂²u$_0$ ..., we have (§ 3 (ii.))
$u_{1}{=}u_{0}+q\partial u_{0}+{\frac {q\cdot q-1}{1\cdot 2}}\partial ^{2}u_{0}+\ldots .$
When q is made indefinitely great, this (writing ƒ(x) for u) becomes Taylor’s Theorem (Infinitesimal Calculus)
$f(x+h){=}f(x)+hf'(x)+{\frac {h^{2}}{1\cdot 2}}f''(x)+\ldots ,$
which, expressed in terms of operators, is
$\mathrm {E} {=}1+h\mathrm {D} +{\frac {h^{2}}{1\cdot 2}}\mathrm {D} ^{2}+{\frac {h^{3}}{1\cdot 2\cdot 3}}\mathrm {D} ^{3}+\ldots {=}e^{h\mathrm {D} }.$
This gives the relation between Δ and D. Also we have
* $u_{2}{=}u_{0}+2q\partial u_{0}+{\frac {2q\cdot 2q-1}{1\cdot 2}}\partial ^{2}u_{0}+\ldots$
* $u_{3}{=}u_{0}+3q\partial u_{0}+{\frac {3q\cdot 3q-1}{1\cdot 2}}\partial ^{2}u_{0}+\ldots$
* $\vdots \qquad \qquad \quad \vdots$
and, if p is any integer,
$u_{\frac {p}{q}}{=}u_{0}+p\partial u_{0}+{\frac {p\cdot p-1}{1\cdot 2}}\partial ^{2}u_{0}+\ldots .$
From these equations u$_{p/q}$ could be expressed in terms of u$_0$, u$_1$, u$_2$, ...; this is a particular case of interpolation (q.v.).
18. Differences and Differential Coefficients.—The various formulae are most quickly obtained by symbolical methods; i.e. by dealing with the operators Δ, E, D, ... as if they were algebraical quantities. Thus the relation E = e$^{hD}$ (§ 17) gives
hD = log$_e$ (1 + Δ) = Δ − 1/2Δ² + $^1$⁄$_3$Δ³ ...
or
h(du/dx)$_0$ = Δu$_0$ − 1/2Δ²u$_0$ + $^1$⁄$_3$Δ³u$_0$ ....
The formulae connecting central differences with differential coefficients are based on the relations μ = cosh 1/2hD = 1/2(e$^{1/2hD}$ + e$^{-1/2hD}$), δ = 2 sinh 1/2hD − e$^{1/2hD}$ − e$^{-1/2hD}$, and may be grouped as follows:—
| u$_0$ | = u$_0$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| μδu$_0$ | = (hD + $^1$⁄$_6$ h$^3$D$^3$ + $^1$⁄$_{120}$ h$^5$D$^5$ + ...)u$_0$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| δ$^2$u$_0$ | = (h$^2$D$^2$ + $^1$⁄$_{12}$ h$^4$D$^4$ + $^1$⁄$_{360}$ h$^6$D$^6$ + ...)u$_0$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| μδ$^3$u$_0$ | = (h$^3$D$^3$ + $^1$⁄$_4$ h$^5$D$^5$ + ...)u$_0$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| δ$^4$u$_0$ | = (h$^4$D$^4$ + $^1$⁄$_6$ h$^6$D$^6$ + ...)u$_0$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| μu$_{1/2}$ | = (1 + $^1$⁄$_8$ h$^2$D$^2$ + $^1$⁄$_{384}$ h$^4$D$^4$ + $^1$⁄$_{46080}$ h$^6$D$^6$ + ...)u$_{1/2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| δu$_{1/2}$ | = (hD + $^1$⁄$_{24}$ h$^3$D$^3$ + $^1$⁄$_{1920}$ h$^5$D$^5$ + ...)u$_{1/2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| μδ$^2$u$_{1/2}$ | = (h$^2$D$^2$ + $^5$⁄$_{24}$ h$^4$D$^4$ + $^{91}$⁄$_{5760}$ h$^6$D$^6$ + ...)u$_{1/2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| δ$^3$u$_{1/2}$ | = (h$^3$D$^3$ + $^1$⁄$_8$ h$^5$D$^5$ + ...)u$_{1/2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| μδ$^4$ u$_{1/2}$ | = (h$^4$D$^4$ + $^7$⁄$_{24}$ h$^6$D$^6$ + ...)u$_{1/2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| u$_0$ | = u$_0$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| hDu$_0$ | = (μδ − $^1$⁄$_6$ μδ$^3$ + $^1$⁄$_{30}$ μδ$^5$ − ...)u$_0$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| h$^2$D$^2$u$_0$ | = (δ$^2$ − $^1$⁄$_{12}$ δ$^4$ + $^1$⁄$_{90}$ δ$^6$ − ...)u$_0$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| h$^3$D$^3$u$_0$ | = (μδ$^3$ − $^1$⁄$_4$ μδ$^5$ + ...)u$_0$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| h$^4$D$^4$u$_0$ | = (δ$^4$ − $^1$⁄$_6$ δ$^6$ + ...)u$_0$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| u$_{1/2}$ | = (μ − $^1$⁄$_8$ μδ$^2$ + $^3$⁄$_{128}$ μδ$^4$ − $^5$⁄$_{1024}$ μδ$^6$ + ...)u$_{1/2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| hDu$_{1/2}$ | = (δ − $^1$⁄$_{24}$ δ$^3$ + $^3$⁄$_{640}$ δ$^5$ − ...)u$_{1/2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| h$^2$D$^2$u$_{1/2}$ | = (μδ$^2$ − $^5$⁄$_{24}$ μδ$^4$ + $^{259}$⁄$_{5760}$ μδ$^6$ − ...)u$_{1/2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| h$^3$D$^3$u$_{1/2}$ | = (δ$^3$ − $^1$⁄$_8$ δ$^5$ + ...)u$_{1/2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
| h$^4$D$^4$ u$_{1/2}$ | = (μδ$^4$ − $^7$⁄$_{24}$ μδ$^6$ + ...)u$_{1/2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ |
When u is a rational integral function of x, each of the above series is a terminating series. In other cases the series will be an infinite one, and may be divergent; but it may be used for purposes of approximation up to a certain point, and there will be a “remainder,” the limits of whose magnitude will be determinate.
19. Sums and Integrals.—The relation between a sum and an integral is usually expressed by the Euler-Maclaurin formula. The principle of this formula is that, if u$_m$ and u$_{m+1}$, are ordinates of a curve, distant h from one another, then for a first approximation to the area of the curve between u$_m$ and u$_{m+1}$ we have 1/2h(u$_m$ + u$_{m+1}$), and the difference between this and the true value of the area can be expressed as the difference of two expressions, one of which is a function of x$_m$, and the other is the same function of x$_{m+1}$. Denoting these by φ(x$_m$) and φ(x$_{m+1}$), we have
$\int _{x_{m}}^{x_{m+1}}udx{=}{\tfrac {1}{2}}h(u_{m}+u_{m+1})+\phi (x_{m+1})-\phi (x_{m}).$
Adding a series of similar expressions, we find
$\int _{f}^{x_{n}}rac{x_{m}}udx{=}h{{\tfrac {1}{2}}u_{m}+u_{m+1}+u_{m+2}+\ldots +u_{n-1}+{\tfrac {1}{2}}u_{n}}+\phi (x_{n})-\phi (x_{m}).$
The function φ(x) can be expressed in terms either of differential coefficients of u or of advancing or central differences; thus there are three formulae.
(i.) The Euler-Maclaurin formula, properly so called, (due independently to Euler and Maclaurin) is
$\int ^{x_{n}}udx{=}h\cdot \mu \sigma u_{n}-{\frac {1}{12}}h^{2}{\frac {du_{n}}{dx}}{\tfrac {1}{2}}+{\tfrac {1}{720}}h^{4}{\frac {d^{3}u_{n}}{dx^{3}}}-{\tfrac {1}{30240}}h^{6}{\frac {d^{5}u_{n}}{dx^{5}}}+\ldots {=}h\cdot \mu \sigma u_{n}-{\frac {\mathrm {B} _{1}}{2!}}h_{2}{\frac {du_{n}}{dx}}+{\frac {\mathrm {B} _{2}}{4!}}h^{4}{\frac {d^{3}u_{n}}{dx^{3}}}-{\frac {\mathrm {B} _{3}}{6!}}h^{6}{\frac {d^{5}u_{n}}{dx^{5}}}+\ldots$
where B$_1$, B$_2$, B$_3$ ... are Bernoulli’s numbers.
(ii.) If we express differential coefficients in terms of advancing differences, we get a theorem which is due to Laplace:—
${\frac {1}{h}}\int _{x_{0}}^{x_{n}}udx{=}\mu \sigma (u_{n}-u_{0})-{\tfrac {1}{12}}(\Delta u_{n}-\Delta u_{0})+{\tfrac {1}{24}}(\Delta ^{2}u_{n}-\Delta ^{2}u_{0})-{\tfrac {19}{720}}(\Delta ^{3}u_{n}-\Delta ^{3}u_{0})+{\tfrac {3}{160}}(\Delta ^{4}u_{n}-\Delta ^{4}u_{0})-\ldots$
For practical calculations this may more conveniently be written
${\frac {1}{h}}\int _{x_{0}}^{x_{n}}udx{=}\mu \sigma (u_{n}-u_{0})+{\tfrac {1}{12}}(\Delta u_{0}-{\tfrac {1}{2}}\Delta ^{2}u_{0}+{\tfrac {19}{60}}\Delta ^{3}u_{0}-\ldots )+{\tfrac {1}{12}}(\Delta 'u_{n}-{\tfrac {1}{2}}\Delta '^{2}u_{n}+{\tfrac {19}{60}}\Delta '^{3}u_{n}-\ldots ),$
where accented differences denote that the values of u are read backwards from u$_n$; i.e. Δ′u$_n$ denotes u$_{n-1}$ − u$_n$, not (as in § 10) u$_n$ − u$_{n-1}$.
(iii.) Expressed in terms of central differences this becomes
${\frac {1}{h}}\int _{x_{0}}^{x_{n}}udx{=}\mu \sigma (u_{n}-u_{0})-{\tfrac {1}{12}}\mu \delta u_{n}+{\tfrac {11}{720}}\mu \delta ^{3}u_{n}-\ldots +{\tfrac {1}{12}}\mu \delta u_{0}-{\tfrac {11}{720}}\mu \delta ^{3}u_{0}+\ldots {=}\mu (\sigma -{\tfrac {1}{12}}\delta +{\tfrac {11}{720}}\delta ^{3}-{\tfrac {191}{60480}}\delta ^{5}+{\tfrac {2497}{3628800}}\delta ^{7}-\ldots )(u_{n}-u_{0}).$
(iv.) There are variants of these formulae, due to taking hu$_{m+1/2}$ as the first approximation to the area of the curve between u$_m$ and u$_{m+1}$; the formulae involve the sum u$_{1/2}$ + u$_{3/2}$ + ... + u$_{n-1/2}$ ≡ σ(u$_n$ − u$_0$) (see Mensuration).
20. The formulae in the last section can be obtained by symbolical methods from the relation
${\frac {1}{h}}\int udx{=}{\frac {1}{h}}\mathrm {D} ^{-1}u{=}{\frac {1}{h\mathrm {D} }}\cdot u.$
Thus for central differences, if we write θ ≡ 1/2hD, we have μ = cosh θ, δ = 2 sinh θ, σ = δ$^{-1}$, and the result in (iii.) corresponds to the formula
sinh θ = θ cosh θ/(1 + 1/3 sinh² θ − 2/3·5 sinh$^4$ θ + 2·4/3·5·7 sinh$^6$ θ − . . .).
References.—There is no recent English work on the theory of finite differences as a whole. G. Boole’s Finite Differences (1st ed., 1860, 2nd ed., edited by J. F. Moulton, 1872) is a comprehensive treatise, in which symbolical methods are employed very early. A. A. Markoff’s Differenzenrechnung (German trans., 1896) contains general formulae. (Both these works ignore central differences.) Encycl. der math. Wiss. vol. i. pt. 2, pp. 919-935, may also be consulted. An elementary treatment of the subject will be found in many text-books, e.g. G. Chrystal’s Algebra (pt. 2, ch. xxxi.). A. W. Sunderland, Notes on Finite Differences (1885), is intended for actuarial students. Various central-difference formulae with references are given in Proc. Lond. Math. Soc. xxxi. pp. 449-488. For other references see Interpolation. (W. F. Sh.)