Page:EB1911 - Volume 18.djvu/162

 tangential area may be expressed in terms of chordal areas. If we write Cundefined for the chordal area obtained by taking ordinates at intervals ℎ, then T1＝2Cundefined−C1. If the trapezette, as seen from above, is everywhere convex or everywhere concave, the true area lies between C1 and T1.

68. Other Rules for Trapezettes.—The extension of this method consists in dividing the trapezette into minor trapezettes, each consisting of two or more strips, and replacing each of these minor trapezettes by a new figure, whose ordinate 𝑣 is an algebraical function of 𝑥; this function being chosen so that the new figure shall coincide with the original figure so far as the given ordinates are concerned. This means that, if the minor trapezette consists of 𝑘 strips, 𝑣 will be of degree 𝑘 or 𝑘−1 in 𝑥, according as the data are the bounding ordinates or the mid-ordinates. If A denotes the true area of the original trapezette, and B the aggregate area of the substituted figures, we have A≃B, where ≃ denotes approximate equality. The value of B is found by the methods of §§ 49–55. The following are some examples.

(i) Suppose that the bounding ordinates are given, and that 𝑚 is a multiple of 2. Then we can take the strips in pairs, and treat each pair as a parabolic trapezette. Applying Simpson’s formula to each of these, we have

A≃ℎ(𝑢0 +4𝑢1 +𝑢2) + ℎ(𝑢2 + 4𝑢3 +𝑢4) + ≃ ℎ(𝑢0 +4𝑢1 + 2𝑢2 +   + 2𝑢𝑚−2 + 4𝑢𝑚−1 + 𝑢𝑚)

This is Simpson’s rule.

(ii) Similarly, if 𝑚 is a multiple of 3, the repeated application of Simpson’s second formula gives Simpson’s second rule

A≃ℎ(𝑢0 + 3𝑢1 + 3𝑢2 + 2𝑢3 + 3𝑢4 + + 3𝑢𝑚−4 + 2𝑢𝑚−3 +3𝑢𝑚−2 + 3𝑢𝑚−1 + 𝑢𝑚).

(iii) If mid-ordinates are given, and 𝑚 is a multiple of 3, the repeated application of the formula of § 55 will give

A ≃ℎ(3𝑢undefined + 2𝑢undefined; + 3𝑢undefined + 3𝑢undefined; + + 2𝑢𝑚− + 3𝑢𝑚−).

69. The formulae become complicated when the number of strips in each of the minor trapezettes is large. The method is then modified by replacing B by an expression which gives the areas of the substituted figures approximately. This introduces a further inaccuracy; but this latter may be negligible in comparison with the main inaccuracies already involved (cf. § 20 (iii)).

Suppose, for instance, that 𝑚＝6, and that we consider the trapezette as a whole; the data being the bounding ordinates. Since there are seven of these, 𝑣 will be of degree 6 in 𝑥; and we shall have (§ 54 (i))

B＝6ℎ(𝑣3 + 2𝑣3 + 4𝑣3 + 6𝑣3)＝6ℎ(𝑢3 + 2𝑢3 + 4𝑢3 + 6𝑢3).

If we replace 6𝑢3 in this expression by 6𝑢3, the method of § 68 gives

A ≃ℎ(𝑢0+ 5𝑢1 + 𝑢2 + 6𝑢3 + 𝑢4 +5𝑢5 +𝑢6);

the expression on the right-hand side being an approximate expression for B, and differing from it only by 6𝑢3. This is Weddle’s rule. If 𝑚 is a multiple of 6, we can obtain an expression for A by applying the rule to each group of six strips.

70. Some of the formulae obtained by the above methods can be expressed more simply in terms of chordal or tangential areas taken in various ways. Consider, for example, Simpson’s rule (§ 68 (i)). The expression for A can be written in the form

ℎ(𝑢0 + 𝑢1 + 𝑢2 + 𝑢3 + 𝑢𝑚−2 + 𝑢𝑚−1 +𝑢𝑚)−ℎ(𝑢0+𝑢2 + 𝑢4 +  + 𝑢𝑚−2 + 𝑢𝑚)

Now, if 𝑝 is any factor of 𝑚, there is a series of equidistant ordinates 𝑢0, 𝑢𝑝, 𝑢2𝑝, 𝑢𝑚−𝑝, 𝑢𝑚; and the chordal area as determined by these ordinates is

𝑝ℎ(𝑢0 + 𝑢𝑝 + 𝑢2𝑝 + . . . . + 𝑢𝑚−𝑝 𝑢𝑚),

which may be denoted by C𝑝. With this notation, the area as given by Simpson’s rule may be written in the form C1−C2 or C1+(C1−C2). The following are some examples of formulae of this kind, in terms of chordal areas.

(i) 𝑚 a multiple of 2 (Simpson’s rule).


 * A ≃ (4C1 − C2≃C1 + (C1 − C2)

(ii) 𝑚 a multiple of 3 (Simpson’s second rule).


 * A ≃ (9C1 − C3) ≃ C1 +(C1 − C3)

(iii) 𝑚 a multiple of 4.

A ≃ (64C1 − 2OC2+C4) ≃ C1 + (C1 − C2) − (C1 − C4).

(iv) 𝑚 a multiple of 6 (Weddle’s rule, or its repeated application).

A ≃(15C1−6C2+C3) ≃ C1 +(C1 − C2) −(C2 − C3).

(v) 𝑚 a multiple of 12.

A ≃ (56C1 − 28C2 + 8C3 − C4) ≃C1+(C1 − C2) − (C2 − C3) + (C3 − C4).

There are similar formulae in terms of the tangential areas T1, T2, T3. Thus (iii) of § 68 may be written A ≃ (9T1 − T3).

71. The general method of constructing the formulae of § 70 for chordal areas is that, if 𝑝, 𝑞, 𝑟, are 𝑘 of the factors (including 1) of 𝑚, we take

A≃PC𝑝 + QC𝑞, +RC𝑟, +. . . ,

where P, Q, R,. . . satisfy the 𝑘 equations The last 𝑘−1 of these equations give

1/P : 1/Q : 1/R :. . ＝ 𝑝2(𝑝2 − 𝑞2)(𝑝2−𝑞2)(𝑞2−𝑟2). . . : 𝑞2(𝑞2−𝑝2)(𝑞2−𝑟2) :. . . 𝑟2(𝑟2 − 𝑝2)(𝑟2 − 𝑞2)

Combining this with the first equation, we obtain the values of P, Q, R,

The same method applies for tangential areas, by taking

A ≃PT𝑝 +QT𝑞, +RT𝑟, +. ..

provided that 𝑝, 𝑞, 𝑟,   are odd numbers.

72. The justification of the above methods lies in certain properties of the series of successive differences of 𝑢. The fundamental assumption is that each group of strips of the trapezette may be replaced by a figure for which differences of 𝑢, above those of a certain order, vanish (§ 54). The legitimacy of this assumption, and of the further assumption which enables the area of the new figure to be expressed by an approximate formula instead of by an exact formula, must be verified) in every case by reference to the actual differences.

73. Correction by means of Extreme Ordinates.—The preceding methods, though apparently simple, are open to various objections in practice, such as the following: (i) The assignment of different coefficients of different ordinates, and even the selection of ordinates for the purpose of finding C2, C3, &c. (§ 70), is troublesome. (ii) This assignment of different coefficients means that different weights are given to different ordinates; and the relative weights ma not agree with the relative accuracies of measurement. (iii) Different formulae have to be adopted for different values of 𝑚; the method is therefore unsuitable for the construction of a table giving successive values of the area up to successive ordinates. (iv) In order to find what formula may be applied, it is necessary to take the successive differences of 𝑢; and it is then just as easy, in most cases, to use a formula which directly involves these differences and therefore shows the degree of accuracy of the approximation.

The alternative method, therefore, consists in taking a simple formula, such as the trapezoidal rule, and correcting it to suit the mutual relations of the differences.

74. To illustrate the method, suppose that we use the chordal area C1, and that the trapezette is in fact parabolic. The difference between C1 and the true area is made up of a series of areas bounded by chords and arcs; this difference becoming less as we subdivide the figure into a greater number of strips.

The fact that C1 does not give the true area is due to the fact that in passing from one extremity of the top of any strip to the other extremity the tangent to the trapezette changes its direction. We have therefore in the first place to see whether the difference can be expressed in terms of the directions of the tangents.

Let KABL (fig. 10) be one of the strips, of breadth ℎ. Draw the tangents at A and B, meeting at T; and through T draw a line parallel to KA and LB, meeting the arc AB in C and the chord AB in V. Draw AD and BE perpendicular to this line, and DF and TG perpendicular to LB. Then AD＝EB＝ℎ, and the triangles AVD and BVE are equal.

The area of the trapezette is less (in fig. 10) than the area of the trapezium KABL by two-thirds of the area of the triangle ATB (§ 34). This latter area is

∆BTE − ∆ATD ＝ ∆BTG − ∆ATD＝ℎ2 tan GTB − ℎ2 tan DAT.

Hence, if the angle which the tangent at the extremity of the ordinate 𝑢θ makes with the axis of 𝑥 is denoted by, we have and thence, by summation,

A＝C1 − ℎ2(tan 𝑚 − 0).

This, in the notation of §§ 46 and 54, may be written

A＝C1 + − ℎ2𝑢′.

Since ℎ＝H/𝑚, the inaccuracy in taking C1 as the area varies as 1/𝑚2.

It might be shown in the same way that

A＝T1+ℎ2(tan 𝑚−0)＝T1 + ℎ2𝑢′.

75. The above formulae apply only to a parabolic trapezette Their generalization is given by the Euler-Maclaurin formula