tangential area may be expressed in terms of chordal areas. If we write C4 for the chordal area obtained by taking ordinates at intervals %h, then T1=2C;-C1 If the trapezette, as seen from above, is everywhere convex or everywhere concave, the truearea 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 v is an algebraical function of x; 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 k strips, v will be of degree k or k-I in x, 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 AQB, 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 rn 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&%k(u0 +4111 +110 + %h(111 +4113 -l-111) + -é.1%h(u1> + 4111 + 2111 + 4111 + 2111 + . + 21¢m 2 + 4um»-I + um)-This is Simpson's rule.

(ii) Similarly, if m is a multiple of 3, the repeated application of Simpson's second formula gives Simpsonls second rule A é %h(u11 + 3111 + 3111+ 21¢s -1- 3111 + - . . + 311.....1 + 214111-3 + Sum-2 “l” Sum-1 'l' um)(iii)

If mid-ordinates are given, and m is a multiple of 3, the repeated application of the formula of § 55 will give A -<l§ h(311§ + 211; + 311% + 311; + . . . + zum-; + Sufi#-é). 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 afurther inaccuracy; but this latter may be negligible in comparison with the main inaccuracies already involved (cf. § 20 (iii)).

Suppose, for instance, that m=6, and that we consider the trapezette as a whole; the data being the bounding ordinates. Since there are seven of these, v will be of degree 6 in x; and we shall havc (§ 54 ())

B =6h(|)3 + %5zU3 + %%§ 41J3 + ffohsilg) =6h(u3'l'f§ '52u3+'§ '§ 54u3+~§ 441°-Bsua). If we replace § *4'66°113 in this expression by § 4f56°111, the method of § 68 gives

A&f'1;h(110+ Sul 'l' u2'l'6f¢a'l' U4 +5145 +111): the expression on the right-hand side being an approximate expression for B, and differing from it only by 1;}1;H6“u3. This is'Weddle's rule. If rn 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 fhfiuo 'i' 141+ 142 'l' 143 'l* - - » Jr um-2'l'u1n-1 +%14m)-§ h(%110+1¢1 'l u4'l"- - ~ 'l'um-2'i'%um)-Now,

if p is any factor of m, there is a series of equidistant ordinates 110, 11, , uzp, um-p1 1151; and the chordal area as determined by these ordinates is

"l" up 'l" u*2p "l" ~»— 'l' urn-P 'l" Bum): which may be denoted by Cp. With this notation, the area as given by Simpson's rule may be written in the form QC1-1}C2 or C1+§ (C1-C1). The following are some examples of formulae of this kind, in terms of chordal areas.(i) m a multiple of 2 (Simpson's rule).

Aé %(4C1 " C2i'Cl + %(C1 - Cz)(ii)

m a multiple of 3 (Simpson's second rule). A& Ii(9C1 °' Cs) ¥C1 -l-%(C1 - Ca)(iii) m a multiple of 4.

A £ {5(64C1 - 2OCg-l-C4) & C1-l-§ (C1 - C1) - £5(C1 - C4). (iv) rn a multiple of 6 (Weddle's rule, or its repeated application). A &i'6(15C1-6C2+C:1) ¥ C1 +i(C1 - Cz) "'I1U(C2 “ Ca)(v) rn a multiple of 12.

A é 31F(56C1 28Cz 'l' 3Ca - C4)

&CI'i'g(C1 - Cz) ' § <Cz - Cs) 'l' nls(Cs - Ci)-There are similar formulae in terms of the tangential areas T1, Tg, T;. Thus (iii) of § 68 may be written A é § (9T1 - Ta). 7I. The general method of constructing the formulae of § 70 for chordal areas is that, if p, q, r, . . are k of the factors (including 1) of m, we take

A&PC, ,+ Qc, +Rc, + ,

where P, Q, R, . . satisfy the k equations P + Q -1- R + = 1,

PP'-l-Qqz-I-Rr2+... =0,

P114 Ll- Qq* + Rr* + .... 5 0,

ppzk-1 | Qqzk-2 ' R111-2 * =

The last k-I of these equations give

I/P11/Q21/R 3—= P2092 - q”)(P”:q2(q2-p2)(q2-r2)

r2(r2 - p2)(r2 - qz) . .:. . .

gombgiing this with the first equation, we obtain the values of Ti-ie same method applies for tangential areas, by taking A£-PTP -l-QT., --RT, -l- ..

provided that p, q, r, . are odd numbers. 72. The justification of the above methods lies in certain properties of the series of successive differences of 11. The fundamental assumption is that each group of strips of the trapezette may be replaced by a figure for which differences of 11, 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 ap roximate formula instead of by an exact formula, must be verified) in every case by actual differences.

73. Correction by means of Extreme Ordinates.-The preceding methods, though apparently simple, are open to various ob'ections 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, Ca, &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 rn; 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 11; 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 1 changes its direction. We have therefore ET' """ ' 8 in the first place to see whether the

vi

I

O.

rf). . .

reference to the

difference can be expressed in terms of the directions of the tangents.

Let KABL (fig. 10) be one of the strips, of breadth h. Draw the tangents at A and B, meeting at T; and through T draw a

line parallel to KA and LB, meeting the F 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 =§ h, and the triangles

AVD and BVE are equal.

The area of the trapezette is less (in fig. Io) than the area of the trapezium KABL by two-thirds of the area of the triangle ATB (§ 34). This latter area is ABTE - AATD = ABTG-AATD=%h“ tan GTB - Hi' tan DAT, Hence, if the angle which the tangent at the extremity of the ordinate 119 makes with the axis of x is denoted by gl/9, we have area from un to 111 =%h(110 -|- 111) -, 'gh2(tan ipl - tan (00), 111 to 112= %h(111 + 112) - {2'h2(tan 1112 - tan ¢/1), K M L

FIG. Io.

um-1 to 11.1. 2 %h(um-1 + um) - 1'§ h”(tan 1//11. - tan =P.1..1)? and thence, by summation, - M

A =C1 iI'gh2(ta.fl IP, " * tall P0)»

This, in the notation of §§ 46 and 54, may be written A=C1 + - fghzu' j§ =f'"

E 0 .

Since h = H/m, the inaccuracy in taking C1 as the area varies as I rnz. .

It might be shown in the same way that A=T1+.x1h2(1a11 -1..- 1111 1.) = T1 + [.11h=11']§§§ . 75. The above formulae apply only to a parabolic trapezette

Their generalization is given by the Euler-Maclaurin formula