Go to National Library of New Zealand Te Puna Mātauranga o Aotearoa
Volume 35, 1902
This text is also available in PDF
(254 KB) Opens in new window
– 408 –

Art. LI.—On the Construction of a Table of Natural Sines by Means of a New Relation between the Leading Differences.

[Read before the Wellington Philosophical Society, 18th November, 1902.]

Part I.

The art of calculating tables of the numerical values of the trigonometrical ratios seems to have fallen into disuse for over a century, as the greater portion of this work was done in the seventeenth and eighteenth centuries, and modern tables are in almost every instance but reprints of the earlier ones.

It appears from the report of the British Association for Advancement of Science, 1873, on mathematical tables, that the most extensive table of natural sines is that given by François Callet in his “Tables Portatives de Logarithmes,” Paris, 1795 (Tirage, 1860). In this work the natural sines are given to fifteen places of decimals for every 0·001 of the quadrant—that is, for every 5′ 24″. In the introduction the process of calculating the table is described, and from it the following summary and extracts are made.

– 409 –

Expanding the cosine and sine by Maclaurin's theorem, we have the usual series,—

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

cos z = 1 − z2/2! + z4/4!-z6/6!+z8/8! − &c.

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

sin z = zz8/3!+z5/5!-z7/7!+z9/9! − &c.

where z is expressed in radians.

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

These series are not in convenient form for numerical calculation, so put z = m/n · π/2

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

then cos z = 1 − m2/n2 · ½!(π/2)2 + m4/n4 · ¼!(π/2)4m6/n6 · ⅙!(π/2)6+&c.

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

and sin z = m/n · π/2-m8/n8 · ⅓!(π/2)8 + m5/n5 · ⅕!(π/2)5-&c.

For convenience these may be written

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

cos z = 1 − m2/n2 · A + m4/n4 · B − m6/n6 · C + m8/n8 · D − &c. (P)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

and sin z = m/n · am3/n3 · b + m5n5 · cm7/n7 · d + &c. (Q)

where (Ćallet, pages 27-and 28)—

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

A = ½! · (π/2)2 = 1·23370, 05501, 36169, 82735, 43

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

B = ¼! · (π/2)4 = 0·25366, 95079, 01048, 01363, 66

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

C = ⅙! · (π/2)6 = 0·02086, 34807, 63352, 96087, 31

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

D = ⅛! · (π/2)8 = 0·00091, 92602, 74839, 42658, 02

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

E = 1/10! · (π/2)10 = 0·00002, 52020, 42373, 06060, 55

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

F = 1/12! · (π/2)12 = 0·00000, 04710, 87477, 88181, 72

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

G = 1/14! · (π/2)14 = 0·00000, 00063, 86603, 08379, 19

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

H = 1/16! · (π/2)16 = 0·00000, 00000, 65659, 63114, 98

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

I = 1/18! · (π/2)18 = 0·00000, 00000, 00529, 44002, 01

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

J = ½0! · (π/2)20 = 0·00000, 00000, 00003, 43773, 92

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

K = ½2! · (π/2)22 = 0·00000, 00000, 00000, 00008, 21

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

M = ½6! · (π/2)26 = 0·00000, 00000, 00000, 00000, 03 &c., and

– 410 –

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

a = π/2 = 1·57079, 63267, 94896, 61923, 13

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

b = ⅓! · (π/2)3 = 0·64596, 40975, 06246, 25365, 58

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

c = ⅕! · (π/2)5 = 0·07969, 26262, 46167, 04512, 05

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

d = 1/7! · (π/2)7 = 0·00468, 17541, 35318, 68810, 07

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

e = 1/9! · (π/2)9 = 0·00016, 04411, 84787, 35982, 19

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

f = 1/11! · (π/2)11 = 0·00000, 35988, 43235, 21208, 53

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

g = 1/13! · (π/2)13 = 0·00000, 00569, 21729, 21967, 93

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

h = 1/15! · (π/2)15 = 0·00000, 00006, 68803, 51098, 11

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

i = 1/17! · (π/2)17 = 0·00000, 00000, 06066, 93573, 11

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

j = 1/19! · (π/2)19 = 0·00000, 00000, 00043, 77065, 47

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

k = 1/21! · (π/2)21 = 0·00000, 00000, 00000, 25714, 23

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

l = 1/23! · (π/2)23 = 0·00000, 00000, 00000, 00125, 39

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

m = 1/25! · (π/2)25 = 0·00000, 00000, 00000, 00000, 52, &c.

By means of these formulæ the sines and cosines of angles are readily obtained: and the calculation of the leading differences for the formation of a table of natural sines is described by Callet thus:—

“S'il est question de trouver les sinus d'une suite d'arcs en progression arithmétique; on peut à l'aide du calcul des differences finies; tirer des formules précédentes, d'autres formules qui donnent les différences premieres, secondes, troisiemes, &c., de ces quantités: pour cela, reprenons la formule Q

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

sin m/n · π/2 = m/n · am8/n8 · b + m5/n5 · cm7/n7 · d + &c.

Substituons, dans cette équation Q,m + Δm à m; il viendra une équation Q1 de laquelle ôtant l'équation Q, nous aurons

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

Δ sin m/n · π/2 = a · Δm/n

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

b · Δm/n8(3m2 + 3mΔm + δm2) + c · Δm/n5(5m4 + 10m8δm + 10m2δm2 + 5mδm8 + δm4)

– 411 –

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

d · Δm/n7 (7m6+21m5Δm + 35m4Δm2 + 35m8Δm8 + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ e · Δm/n9 (9m8 + 36m7Δm + 84m6Δm2 + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

f · Δm/n11 (11m10 + 55m9Δm + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ g · Δm/n13 (13m12 + &c.) (ΔQ)

“Nous trouverons de měme en faisant Δm constant

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

Δ2 sin m/n · π/2 = − b · Δm2/n3 (6m + 6δm)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ c · Δm2/n5 (20m3 + 60m2Δm + 70mΔm2 + 30Δm3)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

d · Δm2/n7 (42m5 + 210m4Δm + 490m3Δm2 + 630Δm2Δm3 + 434m · Δm4 + 126Δm5)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ e · Δm2/n7 (72m7 + 504m6Δm + 1764m5 Δm2 + 3780m4Δm3 + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

f · Δm2/n11 (110m9 + 990m8Δm + 4620m7 Δm2 + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

g · Δm2/n18 (156m11 + 171m10 Δm + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

h · Δm2/n15 (210m13 + &c.) + &c. (Δ2Q

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

Δ3 sin m/n · π/2 = − 6b · Δm3/n3 + c · Δm3/n5 (60m2 + 180mΔm + 150Δm2)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

d · Δm3/n7 (210m4 + 1260m8Δm + 3150m2Δm2 + 3780ΔmΔm3 + 1806 Δm4)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ e · Δm3/n9 (504m6 + 4536m5Δm + 18900m4Δm2 + 45360m3Δm3 + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

f · Δm8/n11 (990m8 + 11880m7Δm + 69300m6Δm2 + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ g · Δm8/n13 (1716m10 + 25740m9Δm + 193050m8Δm2 + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

h · Δm3/n15 (2730m12 + 49140m11Δm + &c.) + &c. (Δ8Q)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

Δ4 sin m/n π/2 = c · Δm4/n5 (120m + 240Δm)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

d · Δm4/n7 (840m3 + 5040m2Δm + 10920m2 + 8400Δm8)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ e · Δm4/n9 (3024m5 + 30240m4Δm + 131040m3Δm2 + &c.)

– 412 –

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

f · Δm4/n11 (7920m7 + 110880m6Δm + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ g · Δm4/n13 (17160m9 + 308880m8Δm + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

h · Δm4/n15 (32760m11 + &c.) + &c. (Δ4Q)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

Δ5 m/n π/2 = 120c · Δm5/n5d · Δm5/n7 (2520m2 + 12600mΔm + 16800Δm2)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ e · Δm5/n9 (15120m4 + 151200m3Δm + 604800m2Δm2 + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

f · Δm5/n11 (55440m6 + 831600m5Δm + 5544000m4Δm2 + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ g · Δm5/n13 (154440m8 + 3088800m7Δm + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

h · Δm5/n15 (360360m10 + &c.) + &c. (Δ5Q)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

Δ6 sin m/n π/2 = − d · Δm6/n7 (5040m + 15120Δm)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ e · Δm6/n9 (60480m3 + 544320m2Δm + 1723680mΔm2 + 1905120Δm3)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

f · Δm6/n11 (332640m5 + 4989600m4Δm + 31600800m3Δm2 + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

+ g · Δm6/n13 (1235520m7 + 25945920m6Δm + &c.)

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

h · Δm6/n15 (3603600m9 + &c.) + &c. (Δ6Q)

Ainsi des autres. Nota. — Ces expressions Δm2, Δm3, &c., tiennent lieu de celles-ci (Δm)2, (Δm)3, &c.” (Callet, pages 59 and 60.)

Other expressions for the differences of sin x are, if Δx = 2θ, as shown below,—

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

Δ4n sin x = 24n sin (x + 4nθ) sin4nθ

Δ4n+1 sin x = 24n+1 cos (x + 4n + 1 θ) sin4n+1θ

Δ4n+2 sin x = − 24n+2 sin (x + 4n + 2 θ) sin4n+2θ

Δ4n+3 sin x = − 24n+3 cos (x + 4n + 3 θ) sin4n+3θ

(De Morgan, “Differential and Integral Calculus.”)

These four expressions may be combined into one expression as follows (Boole, “Finite Differences”):—

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

Δm sin x = 2m sin[x + m (π/2 + θ)] sin mθ

Now, Colenso (“Plane Trigonometry”) shows that when the tabular interval is small it is possible to use the expression

– 413 –

Δ2 sin x = − 22 sin (x + 2θ) sin2θ with advantage in the calculation of a table of natural sines.

It appeared that this relation between the sines and their second differences might be a general one, and upon investigation this was found to be the case, the general relation

Δm sin x = − 22 Δm−2 sin (x + 2θ) · sin2θ (A)

being obtained.

This will now be proved; thus we have

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

Δm sin x = 2m sin [x + m (π/2 + θ)] sinmθ

and similarly

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

Δm-2 sin (x + 2θ) = 2m-2 sin [x + 2θ + (m − 2)(π/2 + θ)] sinm-2θ = 2m-2 sin [x + m (π/2 + θ) -π] sinm-2θ
= − 2m-2 sin [x + (m (π/2 + θ)] sinm-2θ

[The section below cannot be correctly rendered as it contains complex formatting. See the image of the page for a more accurate rendering.]

hence Δm sin x = − 22Δm-2 sin [x + 2θ). sin2θ

It will now be shown how it is possible by means of this relation to calculate readily the leading differences, and thus dispense with the cumbersome series for these differences given above by Callet.

For convenience (A) is preferably written

Δm sin x = − 2(1 − cos 2θ) Δm-2 sin (x + 2θ)

or Δm sin x = − 2(1 − cos Δx) Δm-2 sin (x + Δx)

or Δm sin x = − k · Δm-2 sin (x + Δx) (A1

where k = 2 (1 − cos Δx) and is a constant depending on the tabular interval only.

Now, Δm sin x = − k ·[Δm-2 sin x + Δm-1 sin x] (A2)

hence any difference is expressed in terms of the two preceding differences.

The formation of the leading differences then reduces to the very simple operation shown in (A2) above. It will only be necessary to compare this with the series given above (ΔQ, Δ2Q, &c.) to see how much simpler the method here described is.

In Part II. the application of this method to the calculation of a table of natural sines will be given.