Translation:Supplementum theoriae combinationis observationum erroribus minimis obnoxiae

by Carl Friedrich Gauss

1.

In the previous Memoir, we assumed that the quantities to be determined, with the help of imperfect observations, depended on certain unknown elements, in terms of which they could be expressed: the problem then consisted in deducing from the observations, as accurately as possible, the values of these elements.

In most cases, this is indeed how the question arises. But sometimes the problem is presented differently, and at first glance it is even doubtful whether or not it can be reduced to a problem of the required form. Indeed, frequently the quantities under observation are not given in the form of a function of certain elements, and they seem to only be reducible to such a form by difficult or even ambiguous operations; whereas the nature of the problem does provide certain conditions that the observed values must rigorously satisfy.

Upon closer inspection, however, it can be seen that this case does not differ essentially from the previous one, and that it can be reduced to it. Indeed, if $\pi$ denotes the number of observed quantities, and $\sigma$ denotes the number of constraint equations, then nothing prevents us from choosing, among the former, $\pi -\sigma$ quantities, and considering them as our only unknowns, the others being regarded as functions of these, defined by the constraint equations. By this artifice, we can return to the case of the previous Memoir.

Nevertheless, although this approach often leads to the result in a quite convenient manner, one cannot deny that it is less natural, and it is therefore desirable to treat the problem in another form, which admits a very elegant solution. Moreover, since this new solution leads to faster calculations than the previous one whenever $\sigma$ is less than ${\tfrac {1}{2}}\,\pi ,$ or equivalently, whenever the multitude of elements that we have denoted by $\rho$ in the previous commentary, is greater than ${\tfrac {\pi }{2}},$ it is preferable to use the new solution, even if it would be easy, given the nature of the problem, to unambiguously eliminate the constraint equations.

2.

Let $v,$ $v',$ $v'',$ etc. denote the $\pi$ quantities whose values are provided by observation. Suppose an unknown depends on these and is expressed by a known function $u;$ let $l,$ $l',$ $l'',$ etc. be the values of the differential quotients

${\begin{aligned}{\frac {{d}u}{{d}v}}&,&{\frac {{d}u}{{d}v'}}&,&{\frac {{d}u}{{d}v''}},\ldots \end{aligned}}$

corresponding to the true values of $v,$ $v',$ $v'',$ etc.. If one substituted the true values for $v,$ $v',$ $v'',$ etc. in the function $u,$ one would also obtain the true value of $u;$ but if the observations are affected by errors $e,$ $e',$ $e'',$ etc., there resulting total error for $u$ will be represented by

$le+l'e'+l''e''+\ldots ,$

provided (when $u$ is not linear) that one can neglect the squares and products of errors $e,$ $e',$ $e'',$ etc., as we will always assume.

Although the magnitude of errors $e,$ $e',$ $e'',$ etc., is uncertain, the uncertainty attached to the value found for $u$ can generally be measured by the mean error to be feared in the adopted determination. According to the principles developed in the previous commentary, this average error is

${\sqrt {{l}^{2}{m}^{2}+{l'}^{2}{m'}^{2}+{l''}^{2}{m''}^{2}+\ldots }},$

where $m,$ $m',$ $m'',$ etc. are the average errors of the various observations. If all observations are affected by the same degree of uncertainty, this expression becomes

$m\,{\sqrt {{l}^{2}+{l'}^{2}+{l''}^{2}+\ldots }}.$

and it is clear that in this calculation, it is permissible to replace $l,$ $l',$ $l'',$ etc., with the values taken by the differential coefficients

${\begin{aligned}{\frac {{d}u}{{d}v}}&,&{\frac {{d}u}{{d}v'}}&,&{\frac {{d}u}{{d}v''}},\ldots ,\end{aligned}}$

when one replaces $v,$ $v',$ $v'',$ etc. with their observed values.

3.

When the quantities $v,$ $v',$ $v'',$ etc. are independent, the unknown can only be determined in one way, and the uncertainty attached to the result cannot be avoided or reduced. In this case, the observations provide a value of the unknown which is not arbitrary.

It is quite different when the quantities $v,$ $v',$ $v'',$ etc. are subject to certain relations, which we will suppose are expressed by $\sigma$ conditional equations

${\begin{aligned}{X}&=0,&{Y}&=0,&{Z}&=0,\ldots ,\end{aligned}}$

where ${X},$ ${Y},$ ${Z},$ etc. are given functions of the unknowns $v,$ $v',$ $v'',$ etc. In this case, there are infinitely many ways to represent our unknown as a combination of the quantities $v,$ $v',$ $v'',$ etc., since it is clear that the function $u$ can be replaced by any other function ${U}$ such that the difference ${U}-u$ vanishes identically upon setting

${\begin{aligned}{X}&=0,&{Y}&=0,&{Z}&=0,\ldots .\end{aligned}}$

If the observations were rigorously exact, this substitution would not change the result at all; but, due to inevitable errors, each form adopted for $u$ will correspond to a different result, and the error committed, instead of being

$le+l'e'+l''e''+\ldots ,$

will be

${L}e+{L}'e'+{L}''e''+\ldots ,$

where ${L},$ ${L}',$ ${L}'',$ etc. denote the differential quotients

${\begin{aligned}{\frac {{d}{U}}{{d}v}}&,&{\frac {{d}{U}}{{d}v'}}&,&{\frac {{d}{U}}{{d}v''}},\ldots .\end{aligned}}$

Although it is impossible to assign the value of the various errors, we can, however, compare the mean errors to be feared in the various combinations. The most advantageous combination will be the one that gives the minimum value to the mean error. This error being

${\sqrt {{L}^{2}{m}^{2}+{L''}^{2}{m'}^{2}+{L''}^{2}{m''}^{2}+\,\ldots }},$

we must seek to make the sum

${L}^{2}{m}^{2}+{L''}^{2}{m'}^{2}+{L''}^{2}{m''}^{2}+\ldots ,$

as small as possible.

4.

The infinite number of functions ${U},$ by which $u$ can be replaced, will differ from each other in our investigations only by the values they provide for ${L},$ ${L}',$ ${L}'',$ etc.: therefore, we must first seek the relationships that exist between the systems of values that these coefficients can take. Let

${\begin{array}{c}{\begin{array}{lll}a,&a',&a'',\ldots ,\\b,&b',&b'',\ldots ,\\c,&c',&c'',\ldots ,\end{array}}\\\cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

denote the values taken by the coefficients

${\begin{array}{c}{\begin{array}{lll}{\dfrac {dX}{{d}v}},&{\dfrac {dX}{{d}v'}},&{\dfrac {dX}{{d}v''}},\ldots ,\\[0.75ex]{\dfrac {dY}{{d}v}},&{\dfrac {dY}{{d}v'}},&{\dfrac {dY}{{d}v''}},\ldots ,\\[0.75ex]{\dfrac {dZ}{{d}v}},&{\dfrac {dZ}{{d}v'}},&{\dfrac {dZ}{{d}v''}},\ldots ,\end{array}}\\\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

when we substitute the true values of $v,$ $v',$ $v'',$ etc. It is clear that if increments ${d}v,$ ${d}v',$ ${d}v'',$ etc. are given to $v,$ $v',$ $v'',$ etc. which do not change ${X},$ ${Y},$ ${Z},$ etc. and thus leave each of them with a value of zero, these increments, which satisfy the equations

${\begin{array}{l}{\begin{alignedat}{4}0{}={}&a\,&{d}v&{}+{}a'&{d}v'&{}+{}a''&{d}v''&{}+{}\ldots ,\\0{}={}&b\,&{d}v&{}+{}b'&{d}v'&{}+{}b''&{d}v''&{}+{}\ldots ,\\0{}={}&c\,&{d}v&{}+{}c'&{d}v'&{}+{}c''&{d}v''&{}+{}\ldots ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

will not change the value of ${U}-u,$ and consequently, we will have

$0=(l-{L})\,{d}v+(l'-{L}')\,{d}v'+(l''-{L}'')\,{d}v''+\ldots .$

We easily conclude that ${L},$ ${L}',$ ${L}'',$ etc., must have the form

${\begin{array}{l}{\begin{alignedat}{6}&{L}&{}={}&l&{}+{}&a&x{}+{}&b&y{}+{}&c&z+\ldots ,\\&{L}'&{}={}&l'&{}+{}&a'&x{}+{}&b'&y{}+{}&c'&z+\ldots ,\\&{L}''&{}={}&l''&{}+{}&a''&x{}+{}&b''&y{}+{}&c''&z+\ldots ,\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}$

where $x,$ $y,$ $z,$ etc. are determined multipliers. Conversely, it is clear that for all values of $x,$ $y,$ $z,$ etc., we can form a function ${U},$ for which the values ${L},$ ${L'},$ ${L''},$ etc., will precisely be those provided by these equations, and this function can, according to the above, be substituted for $u.$ The simplest form we can give it is

${U}=u+x\,{X}+y\,{Y}+z\,{Z}+\ldots ,$

but the most general is

${U}=u+x\,{X}+y\,{Y}+z\,{Z}+\ldots +u',$

where $u'$ denotes a function of $v,$ $v',$ $v'',$ etc. which vanishes identically when ${X},$ ${Y},$ ${Z},$ etc., are zero, and whose value, in the current case, will be maximum or minimum. However, this makes no difference for our purposes.

5.

It is now easy to assign values to $x,$ $y,$ $z,$ etc. such that the sum

${L}^{2}m^{2}+{L'}^{2}{m'}^{2}+{L''}^{2}{m''}^{2}+\ldots$

achieves its minimum value. It is clear that to achieve this goal, the knowledge of the absolute mean errors $m,$ $m',$ $m'',$ etc. is not necessary. Rather, it suffices to know only their ratios. Indeed, instead of these quantities, let us introduce the weights of the observations, $p,$ $p',$ $p'',$ etc., i.e. numbers inversely proportional to the squares $m^{2},$ ${m'}^{2},$ ${m''}^{2},$ etc. The quantities $x,$ $y,$ $z,$ etc. must then be determined in such a way that the polynomial

${\frac {\left(ax+by+cz+\ldots +l\right)^{2}}{p}}+{\frac {\left(a'x+b'y+c'z+\ldots +l'\right)^{2}}{p'}}+\ldots$

achieves its minimum value. Suppose that $x^{0},$ $y^{0},$ $z^{0},$ etc., are the determined values of $x,$ $y,$ $z,$ etc. that correspond to this minimum.

We now introduce the following notation:

${\begin{array}{l}{\begin{alignedat}{4}&{\frac {a^{2}}{p}}&{}+{}&{\frac {{a'}^{2}}{p'}}&{}+{}&{\frac {{a''}^{2}}{p''}}&{}+{}&\ldots =(aa),\\[0.5ex]&{\frac {ab}{p}}&{}+{}&{\frac {a'b'}{p'}}&{}+{}&{\frac {a''b''}{p''}}&{}+{}&\ldots =(ab),\\[0.5ex]&{\frac {ac}{p}}&{}+{}&{\frac {a'c'}{p'}}&{}+{}&{\frac {a''c''}{p''}}&{}+{}&\ldots =(ac),\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \\{\begin{alignedat}{4}&{\frac {b^{2}}{p}}\,&{}+{}&\;{\frac {{b'}^{2}}{p'}}\,&{}+{}&\;{\frac {{b''}^{2}}{p''}}\,&{}+{}&\ldots =(bb),\\[0.5ex]&{\frac {c^{2}}{p}}&{}+{}&\;{\frac {{c'}^{2}}{p'}}\,&{}+{}&\;{\frac {{c''}^{2}}{p''}}\,&{}+{}&\ldots =(cc),\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \\{\begin{alignedat}{4}&{\frac {a\,l}{p}}\,&{}+{}&{\frac {a'l'}{p'}}\,&{}+{}&{\frac {a''l''}{p''}}&{}+{}&\ldots =(al),\\[0.5ex]&{\frac {b\,l}{p}}&{}+{}&{\frac {b'l'}{p'}}&{}+{}&{\frac {b''l''}{p''}}&{}+{}&\ldots =(bl),\\[0.5ex]&{\frac {c\,l}{p}}&{}+{}&{\frac {c'l'}{p'}}&{}+{}&{\frac {c''l''}{p''}}&{}+{}&\ldots =(cl),\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

The minimality condition clearly requires that we have

$\left.{\begin{array}{l}{\begin{alignedat}{4}0{}={}&(aa)\,&x^{0}{}+{}&(ab)\,&y^{0}{}+{}&(ac)\,&z^{0}+\ldots {}+{}&(al),\\0{}={}&(ba)\,&x^{0}{}+{}&(bb)\,&y^{0}{}+{}&(bc)\,&z^{0}+\ldots {}+{}&(bl),\\0{}={}&(ca)\,&x^{0}{}+{}&(cb)\,&y^{0}{}+{}&(cc)\,&z^{0}+\ldots {}+{}&(cl),\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right\}\quad$

(1)

After these equations have provided the values of $x^{0},$ $y^{0},$ $z^{0},$ etc., we set

$\left.{\begin{array}{l}{\begin{alignedat}{5}&a&x^{0}{}+{}&b&y^{0}{}+{}&c&z^{0}+\ldots {}+{}&l&{}={}&{L},\\&a'&x^{0}{}+{}&b'&y^{0}{}+{}&c'&z^{0}+\ldots {}+{}&l'&{}={}&{L}',\\&a''&x^{0}{}+{}&b''&y^{0}{}+{}&c''&z^{0}+\ldots {}+{}&l''&{}={}&{L}'',\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}\right\}\quad$

(2)

and the most suitable function to determine our unknown, corresponding to the minimum mean error, will be one whose differential coefficients, for the considered values of the variables, are equal to ${L},$ ${L'},$ ${L''},$ etc. The weight of the determination thus obtained will be

${P}={\frac {1}{{\dfrac {{L}^{2}}{p}}+{\dfrac {{L'}^{2}}{p'}}+{\dfrac {{L''}^{2}}{p''}}+\ldots }},\quad$

(3)

meaning that ${\frac {1}{P}}$ will precisely be the value taken by the polynomial considered above, for the values of the variables $x,$ $y,$ $z,$ etc. which satisfy the equations (1).

6.

In the previous article, we showed how to determine the function ${U}$ that provides the most suitable determination of the unknown $u.$ Let us now examine what value results from this. Let's denote it by ${K}$: it will be obtained by substituting, in ${U},$ the observed values of the quantities $v,$ $v',$ $v'',$ etc. Let $k$ be the value that $u$ takes when the same substitutions are made, and let $\varkappa$ be the true value of the unknown, obtained by substituting the true values of $v,$ $v',$ $v'',$ etc., either in $u,$ or in ${U}.$ We have

${\begin{alignedat}{4}k&{}={}\varkappa {}+{}&l\,&e{}+{}&l'&e'{}+{}&l''&e''\ldots ,\\{K}&{}={}\varkappa {}+{}&{L}\,&e{}+{}&{L}'&e'{}+{}&{L}''&e''\ldots ,\end{alignedat}}$

and, consequently,

${K}=k+({L}-l)\,e+({L}'-l')\,e'+({L}''-l'')\,e''+\ldots .$

Replacing ${L}-l,$ ${L'}-l',$ ${L''}-l'',$ etc. with their values provided by (2), and setting

$\left\{{\begin{array}{l}{\begin{alignedat}{4}&a\,&e{}+{}&a'&e'{}+{}&a''&e''+\ldots &{}={}{\mathfrak {A}},\\&b\,&e{}+{}&b'&e'{}+{}&b''&e''+\ldots &{}={}{\mathfrak {B}},\\&c\,&e{}+{}&c'&e'{}+{}&c''&e''+\ldots &{}={}{\mathfrak {C}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}\right.$

(4)

we will have

${K}=k+{\mathfrak {A}}\,x^{0}+{\mathfrak {B}}\,y^{0}+{\mathfrak {C}}\,z^{0}+\ldots .$

(5)

It is not possible to calculate ${\mathfrak {A}},$ ${\mathfrak {B}},$ ${\mathfrak {C}},$ etc. by means of the formulas (4), as the errors $e,$ $e',$ $e'',$ etc. appearing in them have unknown values, but it is easy to see that the quantities ${\mathfrak {A}},$ ${\mathfrak {B}},$ ${\mathfrak {C}},$ etc., are none other than the values of ${X},$ ${Y},$ ${Z},$ etc., corresponding to the observed values of $v,$ $v',$ $v'',$ etc. In this way, the systems of equations (1), (3), (5) exhibit the complete solution to our problem. In fact, it is clear that the remark made at the end of art. 2, regarding the quantities $l,$ $l',$ $l'',$ etc., can be applied to the calculation of $a,$ $a',$ $a'',$ etc. $b,$ $b',$ $b'',$ etc. etc., or in other words, the true values of $v,$ $v',$ $v'',$ etc. can be replaced with the observed values.

7.

In place of formula (3), which represents the weight of the most probable determination, several other expressions can be given, which are worth developing.

First, observe that if the equations (2) are multiplied by ${\frac {a}{p}},$ ${\frac {a'}{p'}},$ ${\frac {a''}{p''}},$ etc. resp. and then added, the result is

$(aa)\,x^{0}+(ab)\,y^{0}+(ac)\,z^{0}+\ldots ={\frac {a\,{L}}{p}}+{\frac {a'{L}'}{p'}}+{\frac {a''{L}''}{p''}}+\ldots .$

The left hand side is $=0,$ and therefore, denoting the right hand side by $(a\,{L}),$ we have

$(a\,{L})=0,$

and similarly

${\begin{aligned}(b\,{L})&=0,&(c\,{L})&=0,\ldots .\end{aligned}}$

Next, multiplying the equations (2) by ${\frac {L}{p}},$ ${\frac {{L}'}{p'}},$ ${\frac {{L}''}{p''}},$ etc. and adding, we find that

${\frac {l\,{L}}{p}}+{\frac {l'{L}'}{p'}}+{\frac {l''{L}''}{p''}}+\ldots ={\frac {{L}^{2}}{p}}+{\frac {{L'}^{2}}{p'}}+{\frac {{L''}^{2}}{p''}}+\ldots ,$

and thus we obtain a second expression for the weight,

${P}={\frac {1}{{\dfrac {l\,{L}}{p}}+{\dfrac {l'{L}'}{p'}}+{\dfrac {l''{L}''}{p''}}+\ldots }}\cdot$

If, finally, if we multiply the same equations (2) by ${\frac {l}{p}},$ ${\frac {l'}{p'}},$ ${\frac {l''}{p''}},$ etc. and add, we obtain a third expression for the weight,

${P}={\frac {1}{(al)\,x^{0}+(bl)\,y^{0}+(cl)\,z^{0}+\ldots +(ll)}},$

where, with the same notation as before,

${\frac {{l}^{2}}{p}}+{\frac {{l'}^{2}}{p'}}+{\frac {{l''}^{2}}{p''}}+\ldots =(ll).$

One can easily derive from this a fourth expression for the weight,

${\begin{aligned}{\frac {1}{P}}&=(ll)-(aa)\,{x^{0}}^{2}-(bb)\,{y^{0}}^{2}-(cc)\,{z^{0}}^{2}-\ldots \\&{}-{}2\,(ab)\,x^{0}y^{0}-2\,(ac)\,x^{0}z^{0}-2\,(bc)\,y^{0}z^{0}-\ldots .\end{aligned}}$

8.

The general solution we have just outlined is particularly adapted to the case where there is only one unknown to determine. However, when seeking the most plausible values for several unknowns, all of which depend on the same observations, or when it is unknown which unknowns should be derived from the observations, it is necessary to proceed differently, as we will now explain.

Let us consider the quantities $x,$ $y,$ $z,$ etc. as indeterminate, and set

$\left.{\begin{array}{l}{\begin{alignedat}{4}&(aa)\,&x{}+{}&(ab)\,&y{}+{}&(ac)\,&z+\ldots &{}={}\xi ,\\&(ab)\,&x{}+{}&(bb)\,&y{}+{}&(bc)\,&z+\ldots &{}={}\eta ,\\&(ac)\,&x{}+{}&(bc)\,&y{}+{}&(cc)\,&z+\ldots &{}={}\zeta ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}\right\}$

(6)

Suppose that we deduce, by elimination, that

$\left.{\begin{array}{l}{\begin{alignedat}{4}&(\alpha \alpha )\,&\xi {}+{}&(\alpha \beta )\,&\eta {}+{}&(\alpha \gamma )\,&\zeta +\ldots &{}={}x,\\&(\beta \alpha )\,&\xi {}+{}&(\beta \beta )\,&\eta {}+{}&(\beta \gamma )\,&\zeta +\ldots &{}={}y,\\&(\gamma \alpha )\,&\xi {}+{}&(\gamma \beta )\,&\eta {}+{}&(\gamma \gamma )\,&\zeta +\ldots &{}={}z,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right\}$

(7)

Let us first of all observe that the symmetrically placed coefficients are necessarily equal, i.e.

${\begin{array}{l}{\begin{aligned}(\beta \alpha )&=(\alpha \beta ),\\(\gamma \alpha )&=(\alpha \gamma ),\\(\gamma \beta )&=(\beta \gamma ),\end{aligned}}\\\;\cdots \cdots \cdots \cdots \end{array}}$

This follows from the theory of elimination, but we will also give a direct demonstration below.

Thus we have

$\left.{\begin{array}{l}{\begin{alignedat}{4}x^{0}{}={}&{}-{}(\alpha \alpha )\,&(al){}-{}&(\alpha \beta )\,&(bl){}-{}&(\alpha \gamma )\,&(cl)-\ldots ,\\y^{0}{}={}&{}-{}(\beta \alpha )\,&(al){}-{}&(\beta \beta )\,&(bl){}-{}&(\beta \gamma )\,&(cl)-\ldots ,\\z^{0}{}={}&{}-{}(\gamma \alpha )\,&(al){}-{}&(\gamma \beta )\,&(bl){}-{}&(\gamma \gamma )\,&(cl)-\ldots ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right\}$

(8)

and so, by setting

$\left.{\begin{array}{l}{\begin{alignedat}{4}&(\alpha \alpha )\,&{\mathfrak {A}}{}+{}&(\alpha \beta )\,&{\mathfrak {B}}{}+{}&(\alpha \gamma )\,&{\mathfrak {C}}+\ldots &{}={}{A},\\&(\beta \alpha )\,&{\mathfrak {A}}{}+{}&(\beta \beta )\,&{\mathfrak {B}}{}+{}&(\beta \gamma )\,&{\mathfrak {C}}+\ldots &{}={}{B},\\&(\gamma \alpha )\,&{\mathfrak {A}}{}+{}&(\gamma \beta )\,&{\mathfrak {B}}{}+{}&(\gamma \gamma )\,&{\mathfrak {C}}+\ldots &{}={}{C},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}\right\}$

(9)

we obtain

${K}=k-{A}\,(al)-{B}\,(bl)-{C}\,(cl)-\ldots ;$

and if we further set

$\left.{\begin{array}{l}{\begin{alignedat}{5}&a&{A}{}+{}&b&{B}{}+{}&c&{C}+\ldots &{}={}p&&\varepsilon ,\\&a'&{A}{}+{}&b'&{B}{}+{}&c'&{C}+\ldots &{}={}p'&&\varepsilon ',\\&a''&{A}{}+{}&b''&{B}{}+{}&c''&{C}+\ldots &{}={}p''&&\varepsilon '',\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right\}$

(10)

it follows that

${K}=k-l\,\varepsilon -l'\varepsilon '-l''\varepsilon ''-\ldots .$

(11)

9.

Comparing equations (7) and (9) shows that the auxiliary quantities ${A},$ ${B},$ ${C},$ etc., are the values taken by the indeterminates $x,$ $y,$ $z,$ etc., when we set $\xi ={\mathfrak {A}},$ $\eta ={\mathfrak {B}},$ $\zeta ={\mathfrak {C}}$ etc., and thus we have

$\left.{\begin{array}{l}{\begin{alignedat}{4}&(aa)\,&{A}{}+{}&(ab)\,&{B}{}+{}&(ac)\,&{C}+\ldots &{}={}{\mathfrak {A}},\\&(ba)\,&{A}{}+{}&(bb)\,&{B}{}+{}&(bc)\,&{C}+\ldots &{}={}{\mathfrak {B}},\\&(ca)\,&{A}{}+{}&(cb)\,&{B}{}+{}&(cc)\,&{C}+\ldots &{}={}{\mathfrak {C}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}\right\}$

(12)

Multiplying equations (10) by ${\frac {a}{p}},$ ${\frac {a'}{p'}},$ ${\frac {a''}{p''}},$ etc. respectively, and adding, we obtain

and similarly,

$\left.{\begin{array}{l}{\begin{alignedat}{4}{\mathfrak {A}}{}&={}a\,\varepsilon {}+\;a'\varepsilon '{}&+{}&a''&\varepsilon ''+\ldots ,\\&&&\\{\mathfrak {B}}{}&={}b\,\varepsilon {}+\;b'\varepsilon '{}&+{}&b''&\varepsilon ''+\ldots ,\\{\mathfrak {C}}{}&={}c\,\varepsilon {}+\;c'\varepsilon '{}&+{}&c''&\varepsilon ''+\ldots ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right\}$

(13)

Now, since ${\mathfrak {A}}$ is the value taken by ${X}$ after the observed values are substituted for $v,$ $v',$ $v'',$ etc., it is easy to see that if we apply the corrections $-\varepsilon ,$ $-\varepsilon ',$ $-\varepsilon '',$ etc., to each of these quantities, the value of ${X}$ becomes $0$, and similarly ${Y},$ ${Z},$ etc., vanish under this assumption. In a similar way, equation (11) also shows that ${K}$ is the value taken by $u$ after the same substitutions have been made.

We will refer to the application of the corrections $-\varepsilon ,$ $-\varepsilon ',$ $-\varepsilon '',$ etc. as compensation of observations. This evidently leads us to the following very important consequence. Namely, the observations, compensated in the way we have describde, exactly satisfy all of the conditional equations, and cause any function of the observed quantities to take the value resulting from the most suitable combination of the unaltered observations. Since the conditional equations are too few to deduce the exact values of the errors, we will have found, by the above, the most plausible errors. We will henceforth refer to the quantities $\varepsilon ,$ $\varepsilon ',$ $\varepsilon '',$ etc. by this name.

10.

Since the number of observations is greater than the number of conditional equations, an infinite number of systems can be found, other than the system of the most plausible corrections, which exactly satisfy the condition equations.

It is important to examine the relationships that link these various systems. Let $-{E},$ $-{E}',$ $-{E}'',$ etc., be such a system of corrections, different from the most plausible system. Then we have

${\begin{array}{l}{\begin{alignedat}{4}&a\,&{E}&{}+{}a'&{E}'&{}+{}a''&{E}''+\ldots &{}={}{\mathfrak {A}},\\&b\,&{E}&{}+{}b'&{E}'&{}+{}b''&{E}''+\ldots &{}={}{\mathfrak {B}},\\&c\,&{E}&{}+{}c'&{E}'&{}+{}c''&{E}''+\ldots &{}={}{\mathfrak {C}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

Multiplying these equations by ${A},$ ${B},$ ${C},$ etc., and adding, we obtain, with the help of the equations (10),

$p\,\varepsilon \,{E}+p'\varepsilon '{E}'+p''\varepsilon ''{E}''+\ldots ={A}\,{\mathfrak {A}}+{B}\,{\mathfrak {B}}+{C}\,{\mathfrak {C}}+\ldots .$

But combining equations (13) in the same manner gives

$p\,{\varepsilon }^{2}+p'{\varepsilon '}^{2}+p''{\varepsilon ''}^{2}+\ldots ={A}\,{\mathfrak {A}}+{B}\,{\mathfrak {B}}+{C}\,{\mathfrak {C}}+\ldots ;$

(14)

From the combination of these results, one easily deduces

$p\,{E}^{2}+p'{E'}^{2}+p''{E''}^{2}+\ldots$
$=p\,{\varepsilon }^{2}+p'{\varepsilon '}^{2}+p''{\varepsilon ''}^{2}+\ldots +p\,\left({E}-\varepsilon \right)^{2}+p'\left({E}'-\varepsilon '\right)^{2}+p''\left({E}''-\varepsilon ''\right)^{2}+\ldots ,$

Therefore, the sum

$p\,{E}^{2}+p'{E'}^{2}+p''{E''}^{2}+\ldots$

is necessarily greater than the sum

$p\,{\varepsilon }^{2}+p'{\varepsilon '}^{2}+p''{\varepsilon ''}^{2}+\ldots ,$

which can be stated as follows:

Theorem. The sum of squares of the corrections by which the observations can be reconciled with the condition equations, being multiplied by the weights of the corresponding observations respectively, gives a minimum, if the most plausible corrections are adopted.

This is the very principle of least squares, from which equations (12) and (10) can be easily and immediately derived. The minimum value of the sum, which we will now denote by ${S},$ is equal, by equation (14), to

${A}\,{\mathfrak {A}}+{B}\,{\mathfrak {B}}+{C}\,{\mathfrak {C}}+\ldots .$

11.

The determination of the most plausible errors, being independent of $l,$ $l',$ $l'',$ etc., clearly provides the most convenient preparation, regardless of how the observations are to be subsequently used. Furthermore, it is easy to see that, to achieve this goal, it is not necessary to perform elimination indefinitely, i.e. to calculate $(\alpha \alpha ),$ $(\alpha \beta ),$ etc., it suffices to deduce from the equations (12), by a definite elimination, the auxiliary quantities ${A},$ ${B},$ ${C},$ etc., which we will call, in what follows, the correlates of the condition equations

${\begin{aligned}{X}&=0,&{Y}&=0,&{Z}&=0,\ldots ,\end{aligned}}$

and to substitute these quantities into equation (10).

This method leaves nothing to be desired when only the most plausible values of the quantities provided by the observation are required, but it is different when one also seeks the weights of each of the found values. Whichever of the four previous formulas one wants to use, it is essential to know ${L},$ ${L'},$ ${L''},$ etc., or equivalently, $x^{0},$ $y^{0},$ $z^{0},$ etc.; for this reason, it will be useful to more closely study the elimination that provides these quantities, and to obtain a more convenient method for determining the weights.

12.

The relations between the quantities we are concerned with are notably simplified by considering the indefinite second-degree function

${\begin{aligned}(aa)\,x^{2}&+2\,(ab)\,xy+2\,(ac)\,xz+\ldots \\&+(bb)\,y^{2}+2\,(bc)\,yz+\ldots +(cc)\,z^{2}+\ldots ,\end{aligned}}$

which we will denote by ${T}.$ This function is obviously equal to

${\frac {\left(ax+by+cz+\ldots \right)^{2}}{p}}+{\frac {\left(a'x+b'y+c'z+\ldots \right)^{2}}{p'}}+\ldots$

(15)

Moreover, we clearly have

${T}=\xi \,x+\eta \,y+\zeta \,z+\ldots ;$

(16)

and if, with the help of the equations (7), one expresses $x,$ $y,$ $z,$ etc. in terms of $\xi ,$ $\eta ,$ $\zeta ,$ etc.,

${\begin{aligned}{T}=(\alpha \alpha )\,\xi ^{2}&+2\,(\alpha \beta )\,\xi \eta +2\,(\alpha \gamma )\,\xi \zeta +\ldots +(\beta \beta )\,\eta ^{2}\\&+2\,(\beta \gamma )\,\eta \zeta +\ldots +(\gamma \gamma )\,\zeta ^{2}+\ldots .\end{aligned}}$

The theory developed above provides two sets of determined values for the quantities $x,$ $y,$ $z,$ etc., $\xi ,$ $\eta ,$ $\zeta ,$ etc.

The first one is

${\begin{aligned}x&=x^{0},&y&=y^{0},&z&=z^{0},\ldots ,\\\xi &=-(al),&\eta &=-(bl),&\zeta &=-(cl),\ldots ;\end{aligned}}$

whose corresponding value of $T$ is

${T}=(ll)-{\frac {1}{P}},$

as can be seen by comparing the third form of the weight ${P}$ with equation (16), or by directly considering the fourth form.

The second set of values is

${\begin{aligned}x&={A},&y&={B},&z&={C},\ldots ,\\\xi &={\mathfrak {A}},&\eta &={\mathfrak {B}},&\zeta &={\mathfrak {C}},\,\ldots ;\end{aligned}}$

and the corresponding value of ${T}$ is

${T}={S},$

as is evident from formulas (10) and (15), and also from formulas (14) and (16).

13.

We must now subject the function ${T}$ to a transformation similar to that which has been indicated in Theoria Motus Corporum Coelestium art. 182, and was further developed in Disquisitione de elementis ellipticis Palladis. To this end, let us set

${\begin{array}{l}{\begin{alignedat}{2}(bb,1)&{}={}(bb)&{}-{}&{\frac {(ab)^{2}}{(aa)}},\\[0.75ex](bc,1)&{}={}(bc)&{}-{}&{\frac {(ab)\,(ac)}{(aa)}},\\[0.75ex](bd,1)&{}={}(bd)&{}-{}&{\frac {(ab)\,(ad)}{(aa)}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \\{\begin{alignedat}{2}(cc,2)&{}={}(cc)&{}-{}&{\frac {(ac)^{2}}{(aa)}}-{\frac {(bc,1)^{2}}{(bb,1)}},\\[0.75ex](cd,2)&{}={}(cd)&{}-{}&{\frac {(ac)\,(ad)}{(aa)}}-{\frac {(bc,1)\,(bd,1)}{(bb,1)}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \\{\begin{alignedat}{2}(dd,3)&{}={}(dd)&{}-{}&{\frac {(ad)^{2}}{(aa)}}-{\frac {(bd,1)^{2}}{(bb,1)}}-{\frac {(cd,2)^{2}}{(cc,2)}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

(17)

After also setting^[1]

${\begin{array}{l}{\begin{alignedat}{4}&(bb,1)\,&y{}+{}&(bc,1)\,&z{}+{}&(bd,1)\,&w+\ldots &{}={}\eta ',\\&&&(cc,2)\,&z{}+{}&(cd,2)\,&w+\ldots &{}={}\zeta '',\\&&&&&(dd,3)\,&w+\ldots &{}={}\varphi ''',\end{alignedat}}\\\qquad \qquad \qquad \qquad \quad \;\cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}$

we will have

${T}={\frac {\xi ^{2}}{(aa)}}+{\frac {{\eta '}^{2}}{(bb,1)}}+{\frac {{\zeta ''}^{2}}{(cc,2)}}+{\frac {{\varphi '''}^{2}}{(dd,3)}}+\ldots ,$

where $\eta ',$ $\zeta '',$ $\varphi ''',$ etc. are derived from $\xi ,$ $\eta ,$ $\zeta ,$ $\varphi ,$ etc. using the following equations:

${\begin{array}{l}{\begin{alignedat}{6}&\eta '&{}={}&\eta &{}-{}&{\frac {(ab)}{(aa)}}&\,\xi &,\\[0.75ex]&\zeta ''&{}={}&\zeta &{}-{}&{\frac {(ac)}{(aa)}}&\,\xi &{}-{}{\frac {(bc,1)}{(bb,1)}}&\,\eta '&,\\[0.75ex]&\varphi '''&{}={}&\varphi &{}-{}&{\frac {(ad)}{(aa)}}&\,\xi &{}-{}{\frac {(bd,1)}{(bb,1)}}&\,\eta '&{}-{}{\frac {(cd,2)}{(cc,2)}}&\,\zeta '',\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}$

We can then easily derive all the formulas required for our purpose. Namely, to determine the correlates ${A},$ ${B},$ ${C},$ etc., we set

$\left.{\begin{array}{l}{\begin{alignedat}{6}&{\mathfrak {B}}'&{}={}&{\mathfrak {B}}&{}-{}&{\frac {(ab)}{(aa)}}&\,{\mathfrak {A}}&,\\[0.75ex]&{\mathfrak {C}}''&{}={}&{\mathfrak {C}}&{}-{}&{\frac {(ac)}{(aa)}}&\,{\mathfrak {A}}&{}-{}{\frac {(bc,1)}{(bb,1)}}&\,{\mathfrak {B}}'&,\\[0.75ex]&{\mathfrak {D}}'''&{}={}&{\mathfrak {D}}&{}-{}&{\frac {(ad)}{(aa)}}&\,{\mathfrak {A}}&{}-{}{\frac {(bd,1)}{(bb,1)}}&\,{\mathfrak {B}}'&{}-{}{\frac {(cd,2)}{(cc,2)}}&\,{\mathfrak {C}}'',\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}\right\}$

(18)

and then ${A},$ ${B},$ ${C},$ ${D},$ etc. are obtained from the following formulas, starting with the last one:

$\left.{\begin{array}{l}{\begin{aligned}(aa)\,{A}+(ab)\,{B}+(ac)\,{C}+(ad)\,{D}+\ldots &={\mathfrak {A}},\\(bb,1)\,{B}+(bc,1)\,{C}+(bd,1)\,{D}+\ldots &={\mathfrak {B}}',\\(cc,2)\,{C}+(cd,2)\,{D}+\ldots &={C}'',\\(dd,3)\,{D}+\ldots &={\mathfrak {D}}''',\end{aligned}}\\\quad \qquad \qquad \qquad \qquad \qquad \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}\right\}$

(19)

For the sum ${S},$ we have the new formula

${S}={\frac {{\mathfrak {A}}^{2}}{(aa)}}+{\frac {{{\mathfrak {B}}'}^{2}}{(bb,1)}}+{\frac {{{\mathfrak {C}}''}^{2}}{(cc,2)}}+{\frac {{{\mathfrak {D}}'''}^{2}}{(dd,3)}}+\ldots ,$

(20)

and finally, the weight ${P}$ which is to be assigned to the most plausible determination of the quantity $u,$ is given by the formula

${\frac {1}{P}}=(ll)-{\frac {{(al)}^{2}}{(aa)}}-{\frac {{(bl,1)}^{2}}{(bb,1)}}-{\frac {{(cl,2)}^{2}}{(cc,2)}}-{\frac {{(dl,3)}^{2}}{(dd,3)}}-\ldots ;$

(21)

In this formula, we have

$\left.{\begin{alignedat}{5}&(bl,1)&{}={}&(bl)&{}-{}&{\frac {(ab)\,(al)}{(aa)}},\\&(cl,2)&{}={}&(cl)&{}-{}&{\frac {(ac)\,(al)}{(aa)}}&{}-{}&{\frac {(bc,1)\,(bl,1)}{(bb,1)}},\\&(dl,3)&{}={}&(dl)&{}-{}&{\frac {(ad)\,(al)}{(aa)}}&{}-{}&{\frac {(bd,1)\,(bl,1)}{(bb,1)}}&{}-{}&{\frac {(cd,2)\,(cl,2)}{(cc,2)}}\cdot \end{alignedat}}\right\}$

(22)

Formulas (17)....(21), whose simplicity leaves nothing to be desired, provide the complete solution to our problem.

14.

Having solved our primary problem, we will now address some secondary questions, which will shed more light on the subject.

We will first investigate whether it is possible for the elimination that provides $x,$ $y,$ $z,$ etc., as functions of $\xi ,$ $\eta ,$ $\zeta ,$ etc., to become impossible in certain cases. This would obviously happen if the functions $\xi ,$ $\eta ,$ $\zeta ,$ etc., were not independent of each other. Suppose, for a moment, that this is the case, and that one of them can be expressed as a function of the others, such that we have an identity

$\alpha \xi +\beta \eta +\gamma \zeta +\ldots =0,$

where $\alpha ,$ $\beta ,$ $\gamma ,$ etc. are determined numbers. Then we will have

${\begin{array}{l}{\begin{alignedat}{5}\alpha &\,(aa)&{}+{}\beta &\,(ab)&{}+{}\gamma &\,(ac)+\ldots &{}={}&0,\\\alpha &\,(ab)&{}+{}\beta &\,(bb)&{}+{}\gamma &\,(bc)+\ldots &{}={}&0,\\\alpha &\,(ac)&{}+{}\beta &\,(bc)&{}+{}\gamma &\,(cc)+\ldots &{}={}&0,\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots ,\end{array}}$

and if we set

${\begin{array}{l}{\begin{alignedat}{5}\alpha &\,a&{}+{}\beta &\,b&{}+{}\gamma &\,c&{}+{}\ldots {}={}&p\,\theta ,\\\alpha &\,a'&{}+{}\beta &\,b'&{}+{}\gamma &\,c'&{}+{}\ldots {}={}&p'\theta ',\\\alpha &\,a''&{}+{}\beta &\,b''&{}+{}\gamma &\,c''&{}+{}\ldots {}={}&p''\theta '',\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}$

it follows automatically that

${\begin{array}{l}{\begin{alignedat}{5}a&\,\theta {}+{}&a'&\theta '{}+{}&a''&\theta ''+\ldots &{}={}&0,\\b&\,\theta {}+{}&b'&\theta '{}+{}&b''&\theta ''+\ldots &{}={}&0,\\c&\,\theta {}+{}&c'&\theta '{}+{}&c''&\theta ''+\ldots &{}={}&0,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

and consequently,

$p\,{\theta }^{2}+p'{\theta '}^{2}+p''{\theta ''}^{2}+\ldots =0\;;$

where $p,$ $p',$ $p'',$ etc. are, by their nature, all positive, this equation requires

${\begin{aligned}\theta &=0,&\theta '&=0,&\theta ''&=0,\ldots .\end{aligned}}$

Let us now consider the complete differentials ${dX},$ ${dY},$ ${dZ},$ etc. corresponding to the values of $v,$ $v',$ $v'',$ etc., immediately provided by the observations. By our previous results, these differentials

${\begin{array}{l}{\begin{alignedat}{3}a&\,{d}v{}+{}&a'&{d}v'{}+{}&a''&{d}v''+\ldots ,\\b&\,{d}v{}+{}&b'&{d}v'{}+{}&b''&{d}v''+\ldots ,\\c&\,{d}v{}+{}&c'&{d}v'{}+{}&c''&{d}v''+\ldots ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}$

will be related to each other in such a way that when they are multiplied by $\alpha ,$ $\beta ,$ $\gamma ,$ etc. respectively and added, the sum is identically zero, so that among the equations

${\begin{aligned}{X}&=0,&{Y}&=0,&{Z}&=0,\ldots ,\end{aligned}}$

at least one can be regarded as superfluous, because it will be satisfied as soon as the others are.

Upon closer examination of the question, it is seen that this conclusion is only applicable to values of the variables that differ infinitesimally from those provided by observation. In particular, there are two cases to distinguish: first, the case where one of the equations

${\begin{aligned}{X}&=0,&{Y}&=0,&{Z}&=0,\ldots ,\end{aligned}}$

is generally and absolutely contained within the others, and can therefore be eliminated; and second, the case where, for the particular values of $v,$ $v',$ $v'',$ etc. provided by the observations, one of the functions ${X},$ ${Y},$ ${Z},$ etc., ${X}$ for example, achieves its maximum or minimum value, or, more generally, a value where its differential vanishes, while the other equations remain satisfied. However, since we only consider variations of our variables whose squares are negligible, this second case (which in practice will rarely occur) can be assimilated into the first, and one of the constraint equations can be removed as redundant. If the remaining equations are independent in the sense we have indicated, it is certain, according to the above, that the elimination is possible. However, we reserve the right to return to this matter, which deserves to be examined as a theoretical subtlety rather than as a question of practical utility.

15.

In the previous commentary, arts. 37 sqq., we have shown how to approximate the weight of a determination a posteriori. If approximate values of $\pi$ quantities are provided by equally precise observations, and they are compared with the values resulting from the most plausible assumptions one can make about the $\rho$ elements on which they depend, it has been seen that one must add the squares of the obtained differences, divide the sum by $\pi -\rho ,$ and then the resulting quotient can be regarded as an approximate value of the square of the average error inherent in this kind of observation. If the observations are of unequal precision, the only modification that must be made is that one must multiply the squares of the differences by the respective weights of the corresponding observations, and the mean error obtained in this way relates to observations whose weight is taken to be unity.

In the current case, the sum of the squares we are talking about obviously coincides with the sum ${S}$, and the difference $\pi -\rho$ with the number $\sigma$ of the condition equations. Consequently, for the average error of observations whose weight is $1$, we will have the expression ${\sqrt {\tfrac {S}{\sigma }}},$ so that the determination will be increasingly reliable for larger values of $\sigma .$

However, it is worthwhile to establish this result independently of the reasoning in the first investigation. To this end, it will be helpful to introduce some new notation. Suppose that, corresponding to the values

${\begin{aligned}\xi &=a,&\eta &=b,&\zeta &=c,\ldots ,\end{aligned}}$

we have

${\begin{aligned}x&=\alpha ,&y&=\beta ,&z&=\gamma ,\ldots ,\end{aligned}}$

so that

${\begin{alignedat}{5}&\alpha &{}={}a&\,(\alpha \alpha )&{}+{}b&\,(\alpha \beta )&{}+{}c&\,(\alpha \gamma )&{}+{}\ldots ,\\&\beta &{}={}a&\,(\alpha \beta )&{}+{}b&\,(\beta \beta )&{}+{}c&\,(\beta \gamma )&{}+{}\ldots ,\\&\gamma &{}={}a&\,(\alpha \gamma )&{}+{}b&\,(\beta \gamma )&{}+{}c&\,(\gamma \gamma )&{}+{}\ldots ;\end{alignedat}}$

and furthermore, corresponding to the values

${\begin{aligned}\xi &=a',&\eta &=b',&\zeta &=c',\ldots ,\end{aligned}}$

we have

${\begin{aligned}x&=\alpha ',&y&=\beta ',&z&=\gamma ',\ldots ;\end{aligned}}$

and finally, corresponding to the values

${\begin{aligned}\xi &=a'',&\eta &=b'',&\zeta &=c'',\ldots ,\end{aligned}}$

we have

${\begin{aligned}x&=\alpha '',&y&=\beta '',&z&=\gamma '',\ldots ,\end{aligned}}$

and so on.

Combining equations (4) and (9) then yields

${\begin{alignedat}{4}{A}&{}={}\alpha &\,e{}+{}&\alpha '&e'{}+{}&\alpha ''&e''+\ldots ,\\{B}&{}={}\beta &\,e{}+{}&\beta '&e'{}+{}&\beta ''&e''+\ldots ,\\{C}&{}={}\gamma &\,e{}+{}&\gamma '&e'{}+{}&\gamma ''&e''+\ldots ,\end{alignedat}}$

and, since

${\mathfrak {S}}={\mathfrak {A}}\,{A}+{\mathfrak {B}}\,{B}+{\mathfrak {C}}\,{C}+\ldots ,$

we will clearly have

${\begin{alignedat}{7}{S}&{}={}(a&\,e{}+{}&a'&e'{}+{}&a''&e''+\ldots )\;&(\alpha &\,e{}+{}&\alpha '&e'{}+{}&\alpha ''&e''+\ldots )\\&{}+{}(b&\,e{}+{}&b'&e'{}+{}&b''&e''+\ldots )\;&(\beta &\,e{}+{}&\beta '&e'{}+{}&\beta ''&e''+\ldots )\\&{}+{}(c&\,e{}+{}&c'&e'{}+{}&c''&e''+\ldots )\;&(\gamma &\,e{}+{}&\gamma '&e'{}+{}&\gamma ''&e''+\ldots )&{}+{}\ldots .\end{alignedat}}$

16.

The series of observations that provide the quantities $v,$ $v',$ $v'',$ etc., affected by random errors $e,$ $e',$ $e'',$ etc., can be considered as a test that does not reveal the magnitude of each error, but instead, by means of the rules explained above, allows us to determine the quantity ${S}$ as a function of the errors. In such a test, some errors may be larger and others smaller; but the greater the number of errors used, the greater the probability that ${S}$ differs little from its average value. The difficulty thus boils down to finding the average value of ${S}.$

By the principles outlined in the first Memoir, which it is unnecessary to reproduce here, we find for this average value

$(a\,\alpha +b\,\beta +c\,\gamma +\ldots )\,{m}^{2}+(a'\alpha '+b'\beta '+c'\gamma '+\ldots )\,{m'}^{2}+\ldots .$

Let $\mu$ denote the average error corresponding to observations with weight 1, so that we have

$\mu ^{2}=pm^{2}=p'{m'}^{2}=p''{m''}^{2}=\ldots$

Then the previous expression can be written as follows:

${\begin{aligned}&\left({\frac {a\,\alpha }{p}}+{\frac {a'\alpha '}{p'}}+{\frac {a''\alpha ''}{p''}}+\ldots \right)\mu ^{2}+\left({\frac {b\,\beta }{p}}+{\frac {b'\beta '}{p'}}+{\frac {b''\beta ''}{p''}}+\ldots \right)\mu ^{2}\\+{}&\left({\frac {c\,\gamma }{p}}+{\frac {c'\gamma '}{p'}}+{\frac {c''\gamma ''}{p''}}+\ldots \right)\mu ^{2},\end{aligned}}$

but we have found

${\frac {a\,\alpha }{p}}+{\frac {a'\alpha '}{p'}}+{\frac {a''\alpha ''}{p''}}+\ldots =(aa)(\alpha \alpha )+(ab)(\alpha \beta )+(ac)(\alpha \gamma )+\ldots ,$

and thus the right-hand side is unity, as can be easily recognized by comparing equations (6) and (7). Similarly, we find that

${\begin{alignedat}{4}&{\frac {b\,\beta }{p}}&&+{\frac {b'\beta '}{p'}}&&+{\frac {b''\beta ''}{p''}}&&+\ldots =1,\\&{\frac {c\,\gamma }{p}}&&+{\frac {c'\gamma '}{p'}}&&+{\frac {c''\gamma ''}{p''}}&&+\ldots =1,\end{alignedat}}$

and so on.

Hence the average value of ${S}$ is $=\sigma \mu ^{2},$ and insofar as it is permissible to regard the random value of ${S}$ as equal to the average value, we conclude that $\mu ={\sqrt {\tfrac {S}{\sigma }}}.$

17.

One can assess the confidence warranted by this determination by calculating either the mean error to be feared, or its square. The latter will be the square root of the average value of the expression $\left({\frac {S}{\sigma }}-\mu ^{2}\right)^{2},$ which can be developed by reasoning similar to that which was explained in the first commentary (arts. 39 sqq.). We will omit these for the sake of brevity, and simply give the result.

The mean error to be feared in determining the square $\mu ^{2}$ is expressed by

${\sqrt {{\frac {2\,\mu ^{4}}{\sigma }}+{\frac {\nu ^{4}-3\,\mu ^{4}}{\sigma ^{2}}}\,{N}}},$

where $\nu ^{4}$ is the average value of the fourth powers of the errors with weight 1, and ${N}$ is the sum

${\begin{aligned}&(a\,\alpha +b\,\beta +c\,\gamma +\ldots )^{2}\\+&(a'\alpha '+b'\beta '+c'\gamma '+\ldots )^{2}\\+&(a''\alpha ''+b''\beta ''+c''\gamma ''+\ldots )^{2}\\+&\ldots .\end{aligned}}$

This sum cannot generally be simplified; however it can be shown, by a method similar to that used in art. 40 of the previous commentary, that its value is between $\pi$ and ${\frac {\sigma ^{2}}{\pi }}.$ With the hypothesis on which we originally established the method of least squares, the term containing this sum disappears, because $\nu ^{4}=3\,\mu ^{4},$ and the precision that must be attributed to the determination $\mu ={\sqrt {\frac {S}{\sigma }}}$ is therefore the same as if one had operated on $\sigma$ observations affected by exactly known errors, following the precepts of arts. 15, 16 of the previous commentary.

18.

For the compensation of observations, we have seen above that there are two operations to perform: first, it is necessary to determine the correlates of the conditional equations, that is, ${A},$ ${B},$ ${C},$ etc., which satisfy equations (12); secondly, we must substitute these quantities into equation (10). The compensation thus obtained can be called perfect or complete, as opposed to an imperfect or incomplete compensation. We will use the latter term to describe those resulting from the same equations (10), in which values of ${A},$ ${B},$ ${C}$ are substituted which do not satisfy the equations (12), i.e. which satisfy only some or none. We will not consider such a system of corrections here, and we will not even describe them as a compensation. When equations (10) are satisfied, systems (12) and (13) become equivalent, and the difference of which we speak can then be stated as follows: Fully compensated observations satisfy the conditional equations

${\begin{aligned}{X}&=0,&{Y}&=0,&{Z}&=0,\ldots ;\end{aligned}}$

Partially compensated observations satisfy only some of these equations, and perhaps none; a compensation for which all equations are satisfied is necessarily complete.

19.

It follows from the very definition of a compensation that a combination of two compensations is again a compensation, and it does not matter whether the rules given to obtain a perfect compensation are applied to the raw observations or to imperfectly compensated observations.

Let $-\theta ,$ $-\theta ',$ $-\theta '',$ etc., be a system of incomplete compensations resulting from the formulas

$\left\{{\begin{array}{l}{\begin{alignedat}{5}&\theta \,p&{}={}{A}^{0}&\,a&{}+{}{B}^{0}&\,b&{}+{}{C}^{0}&\,c&{}+{}\ldots ,\\&\theta 'p'&{}={}{A}^{0}&\,a'&{}+{}{B}^{0}&\,b'&{}+{}{C}^{0}&\,c'&{}+{}\ldots ,\\&\theta ''p''&{}={}{A}^{0}&\,a''&{}+{}{B}^{0}&\,b''&{}+{}{C}^{0}&\,c''&{}+{}\ldots ,\end{alignedat}}\\\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \end{array}}\right.$

(I)

Since the observations thus changed do not satisfy all of the conditional equations, let ${\mathfrak {A}}^{\ast },$ ${\mathfrak {B}}^{\ast },$ ${\mathfrak {C}}^{\ast },$ etc., be the values taken by ${X},$ ${Y},$ ${Z},$ etc., when the values obtained for $v,$ $v',$ $v'',$ etc., are substituted into them. We must then find the values ${A}^{\ast },$ ${B}^{\ast },$ ${C}^{\ast },$ etc., which satisfy the equations

$\left\{{\begin{array}{l}{\begin{alignedat}{5}{\mathfrak {A}}^{\ast }&{}={}&{A}^{\ast }&\,(aa)&{}+{}{B}^{\ast }&\,(ab)&{}+{}{C}^{\ast }&\,(ac)&{}+{}\ldots ,\\{\mathfrak {B}}^{\ast }&{}={}&{A}^{\ast }&\,(ab)&{}+{}{B}^{\ast }&\,(bb)&{}+{}{C}^{\ast }&\,(bc)&{}+{}\ldots ,\\{\mathfrak {C}}^{\ast }&{}={}&{A}^{\ast }&\,(ac)&{}+{}{B}^{\ast }&\,(bc)&{}+{}{C}^{\ast }&\,(cc)&{}+{}\ldots ,\end{alignedat}}\\\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \end{array}}\right.$

(II)

Once this is done, the complete compensation of the observations thus modified will be accomplished by the new changes $-\varkappa ,$ $-\varkappa ',$ $-\varkappa '',$ etc.; $\varkappa ,$ $\varkappa ',$ $\varkappa '',$ etc., being deduced from the formulas

$\left\{{\begin{array}{l}{\begin{alignedat}{5}&\varkappa \,p&{}={}{A}^{\ast }&\,a&{}+{}{B}^{\ast }&\,b&{}+{}{C}^{\ast }&\,c&{}+{}\ldots ,\\&\varkappa 'p'&{}={}{A}^{\ast }&\,a'&{}+{}{B}^{\ast }&\,b'&{}+{}{C}^{\ast }&\,c'&{}+{}\ldots ,\\&\varkappa ''p''&{}={}{A}^{\ast }&\,a''&{}+{}{B}^{\ast }&\,b''&{}+{}{C}^{\ast }&\,c''&{}+{}\ldots ,\end{alignedat}}\\\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \end{array}}\right.$

(III)

Let us now examine how these corrections agree with the complete compensation of the raw observations. First of all, it is clear that we have

${\begin{array}{l}{\begin{alignedat}{5}{\mathcal {A}}^{\ast }&{}={}{\mathcal {A}}&{}-{}&a&\,\theta &{}-{}a'&\theta '&{}-{}a''&\theta ''{}-{}\ldots ,\\{\mathcal {B}}^{\ast }&{}={}{\mathcal {B}}&{}-{}&b&\,\theta &{}-{}b'&\theta '&{}-{}b''&\theta ''{}-{}\ldots ,\\{\mathcal {C}}^{\ast }&{}={}{\mathcal {C}}&{}-{}&c&\,\theta &{}-{}c'&\theta '&{}-{}c''&\theta ''{}-{}\ldots ,\end{alignedat}}\\\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \end{array}}$

In these equations, if we replace $\theta ,$ $\theta ',$ $\theta '',$ etc. with their values from (I), and if we replace ${\mathcal {A}}^{\ast },$ ${\mathcal {B}}^{\ast },$ ${\mathcal {C}}^{\ast },$ etc., with their values from (II), then we obtain

${\begin{array}{l}{\begin{alignedat}{3}{\mathcal {A}}{}={}&({A}^{0}+{A}^{\ast })\,(aa)&{}+{}&({B}^{0}+{B}^{\ast })\,(ab)&{}+{}\ldots ,\\{\mathcal {B}}{}={}&({A}^{0}+{A}^{\ast })\,(ab)&{}+{}&({B}^{0}+{B}^{\ast })\,(bb)&{}+{}\ldots ,\\{\mathcal {C}}{}={}&({A}^{0}+{A}^{\ast })\,(ac)&{}+{}&({B}^{0}+{B}^{\ast })\,(bc)&{}+{}\ldots ,\end{alignedat}}\\\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \cdots \cdots \cdots \end{array}}$

It follows that the correlates of the conditional equations (12) are

${\begin{aligned}{A}&={A}^{0}+{A}^{\ast },&{B}&={B}^{0}+{B}^{\ast },&{C}&={C}^{0}+{C}^{\ast },\ldots ,\end{aligned}}$

and then equations (10), (I), and (III) show that

${\begin{aligned}\varepsilon &=\theta +\varkappa ,&\varepsilon '&=\theta '+\varkappa ',&\varepsilon ''&=\theta ''+\varkappa '',\ldots ;\end{aligned}}$

Hence the perfect compensation produces the same value for each unknown, whether calculated directly or obtained through an incomplete compensation.

20.

When there are a large number of conditional equations, determining the correlates ${A},$ ${B},$ ${C},$ etc. may require calculations so lengthy that the calculator may be deterred: in such cases, it may be advantageous to obtain a complete compensation using a series of approximations based on the theorem of the preceding article. For this purpose, the conditional equations will be divided into two or more groups, and first, a compensation will be sought that satisfies the equations of the first group. Then, the values modified by this initial calculation will be treated, and they will be corrected again, considering only the equations of the second group. This second calculation will generally yield results that no longer satisfy the equations of the first group, and if only two groups were formed, one must then return to the first group and satisfy its equations with new corrections. The observations will then undergo a fourth compensation, in which only the conditions of the second group are considered, and by alternating in this way between the two groups of equations, corrections will be formed that will necessarily become smaller and smaller. If the choice of groups has been made skillfully, one will arrive at stable values after a few iterations.

When forming more than two groups, one must proceed in the same manner, with the various groups being used successively until the last, after which one returns to the first group to repeat the process in the same order. We have only indicated this method, the success of which will depend greatly on the skill of the calculator.

21.

We have yet to provide the proof of the lemma which was assumed in art. 8. For the sake of perspicuity, let us use notation which is better adapted to this matter.

Let $x^{0},$ $x',$ $x'',$ etc., be indeterminates, and suppose that the equations

${\begin{array}{l}{\begin{alignedat}{4}&n^{00}\,&x^{0}&{}+{}n^{01}\,&x'&{}+{}n^{02}\,&x''&{}+{}\ldots ={X}^{0},\\&n^{10}\,&x^{0}&{}+{}n^{11}\,&x'&{}+{}n^{12}\,&x''&{}+{}\ldots ={X}',\\&n^{20}\,&x^{0}&{}+{}n^{21}\,&x'&{}+{}n^{22}\,&x''&{}+{}\ldots ={X}'',\end{alignedat}}\\\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \end{array}}$

have, through elimination, led to the following:

${\begin{array}{l}{\begin{alignedat}{4}&{N}^{00}\,&{X}^{0}&{}+{}{N}^{01}\,&{X}'&{}+{}{N}^{02}\,&{X}''&{}+{}\ldots =x^{0},\\&{N}^{10}\,&{X}^{0}&{}+{}{N}^{11}\,&{X}'&{}+{}{N}^{12}\,&{X}''&{}+{}\ldots =x',\\&{N}^{20}\,&{X}^{0}&{}+{}{N}^{21}\,&{X}'&{}+{}{N}^{22}\,&{X}''&{}+{}\ldots =x'',\end{alignedat}}\\\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \end{array}}$

By substituting, in the first two equations of the second system, the values of ${X}^{0},$ ${X}',$ ${X}'',$ etc., provided by the first, we obtain two identical equations:

${\begin{array}{r}{\begin{alignedat}{4}x^{0}{}={}&{N}^{00}&\,(&n^{00}\,&x^{0}&{}+{}n^{01}\,&x'&{}+{}n^{02}\,&x''&{}+{}\ldots )\\{}+{}&{N}^{01}&\,(&n^{10}\,&x^{0}&{}+{}n^{11}\,&x'&{}+{}n^{12}\,&x''&{}+{}\ldots )\\{}+{}&{N}^{02}&\,(&n^{20}\,&x^{0}&{}+{}n^{21}\,&x'&{}+{}n^{22}\,&x''&{}+{}\ldots )\\{}+{}&{N}^{03}&\,(&n^{30}\,&x^{0}&{}+{}n^{31}\,&x'&{}+{}n^{32}\,&x''&{}+{}\ldots )\end{alignedat}}\\\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\[1ex]{\begin{alignedat}{4}x'{}={}&{N}^{10}&\,(&n^{00}\,&x^{0}&{}+{}n^{01}\,&x'&{}+{}n^{02}\,&x''&{}+{}\ldots )\\{}+{}&{N}^{11}&\,(&n^{10}\,&x^{0}&{}+{}n^{11}\,&x'&{}+{}n^{12}\,&x''&{}+{}\ldots )\\{}+{}&{N}^{12}&\,(&n^{20}\,&x^{0}&{}+{}n^{21}\,&x'&{}+{}n^{22}\,&x''&{}+{}\ldots )\end{alignedat}}\\\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \end{array}}$

These equations being identical, we can substitute any quantities we wish for $x^{0},$ $x',$ $x'',$ etc. Let's choose in the first

${\begin{aligned}x^{0}&={N}^{10},&x'&={N}^{11},&x''&={N}^{12},\ldots ,\end{aligned}}$

and in the second

${\begin{aligned}x^{0}&={N}^{00},&x'&={N}^{01},&x''&={N}^{02},\ldots ,\end{aligned}}$

Subtracting the two identities element-wise yields ${\begin{aligned}{N}^{10}-{N}^{01}{}={}&({N}^{00}\,{N}^{11}-{N}^{10}\,{N}^{01})\,(n^{01}-n^{10})\\{}+{}&({N}^{00}\,{N}^{12}-{N}^{10}\,{N}^{02})\,(n^{02}-n^{20})\\{}+{}&({N}^{00}\,{N}^{13}-{N}^{10}\,{N}^{03})\,(n^{03}-n^{30})\\{}+{}&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\[0.75ex]{}+{}&({N}^{01}\,{N}^{12}-{N}^{11}\,{N}^{02})\,(n^{12}-n^{21})\\{}+{}&({N}^{01}\,{N}^{13}-{N}^{11}\,{N}^{03})\,(n^{13}-n^{31})\\{}+{}&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\[0.75ex]{}+{}&({N}^{02}\,{N}^{13}-{N}^{12}\,{N}^{03})\,(n^{23}-n^{32})\\{}+{}&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \end{aligned}}$ This can be written more succinctly as: ${N}^{10}-{N}^{01}=\sum \left({N}^{0\alpha }\,{N}^{1\beta }-{N}^{1\alpha }\,{N}^{0\beta }\right)\left(n^{\alpha \beta }-n^{\beta \alpha }\right),$ where $\alpha ,$ $\beta$ denote two arbitrarily chosen indices; from this, we conclude that the equalities ${\begin{aligned}n^{01}&=n^{10},&n^{02}&=n^{20},\end{aligned}}$ and more generally, $n^{\alpha \beta }=n^{\beta \alpha },\ldots ,$ and so ${N}^{01}={N}^{10}\,;$ Moreover, since the order of the indeterminates is arbitrary, it is evident that, under the assumed hypothesis, we have ${N}^{\alpha \beta }={N}^{\beta \alpha }$ in general.

22.

Since the method outlined in this commentary is intended primarily for useful application in high geodesy, the reader may appreciate it if we include some examples drawn from this branch of science.

The conditional equations that exist among the angles of a system of triangles can generally be divided into three categories.

I. The sum of the horizontal angles formed around the same vertex, encompassing the entire horizon, must be equal to four right angles.

II. The sum of the angles of each triangle can always be considered as known; for even when the triangle is situated on a curved surface, the excess of the sum of its angles over two right angles can be calculated with such approximation that it is permissible to consider the result as absolutely exact.

III. Finally, a third kind of relation is obtained by examining the ratios of the sides in triangles that form a closed network. If, indeed, the triangles are placed in such a way that the second triangle has a side $a$ in common with the first, and a side $b$ in common with the third; if the fourth triangle has two sides $c$ and $d$ respectively in common with the third and the fifth, and so on, until the last triangle, which has a common side $k$ with the previous one, and with the first of all a common side $l,$ then the quotients

${\begin{aligned}{\frac {a}{l}}&,&{\frac {b}{a}}&,&{\frac {c}{b}}&,&{\frac {d}{c}},\ldots &,&{\frac {l}{k}}&,\end{aligned}}$

can be calculated by means of the angles opposite to them in the triangle whose two compared sides are part of it, and since the product of these fractions is obviously unity, there will be a relation between the sines of the various measured angles (reduced by one-third of the spherical or spheroidal excess when operating on a curved surface).

In somewhat complicated networks, it often happens that there are too many equations of the second and third category, and some of them are contained in the rest. On the other hand, it may happen in rare cases that some new equations must be added to those of the second category. For example, this will occur when the network contains polygons not divided into triangles by measurements; one must then introduce equations related to figures with more than three sides. On another occasion, we will return to this and provide more details on these various circumstances, but examining them now would divert us from our purpose. However, we cannot refrain from making an essential remark here for those who would like to make the rigorous application of our theory: we always assume that the quantities denoted by $v,$ $v',$ $v'',$ etc., have been observed directly, or deduced from observations, in such a way that their determinations are independent of each other, or at least can be regarded as such. In common practice, one observes the angles that can be considered as the elements $v,$ $v',$ $v'',$ etc. themselves. But one must not forget that if the system contains, in addition, triangles whose angles have not been directly observed and have been deduced from those that are known, by additions or subtractions, then these angles should not be counted among the quantities determined by the observation, and they must be included in the calculation as functions of the elements used to form them. It will be different if one adopts the method of observations of Struve (Astronomische Nachrichten II, p. 431), which consists of determining all directions around the same vertex, by relating them all to a single arbitrary direction. The angles measured in this way will then be taken for $v,$ $v',$ $v'',$ etc., and the angles of the triangles will all appear as differences. The equations of the first category should, in this case, be suppressed as superfluous, for they will be identically satisfied. The method that I myself followed in the triangulations carried out during the last few years differs from the two previous methods; however, it can be reconciled, as to the result, with Struve's method, in the sense that, at each station, one must regard $v,$ $v',$ $v'',$ etc., as the angles formed by the directions emanating from it, with a single arbitrarily chosen line.

We will give two examples: the first relates to the first mode of operation, and the second is related to observations made according to the second method.

23.

The first example comes from the work of de Krayenhof, Précis historique des opérations trigonométriques faites en Hollande. We will seek to compensate the part of the observations related to the system of triangles that are contained between Harlingen, Sneek, Oldeholtpade, Ballum, Leeuwarden, Dockum, Drachten, Oosterwolde, and Gröningen. Between these points, nine triangles numbered 121, 122, 123, 124, 125, 127, 128, 131, 132 were formed in the cited work. The observed angles are as follows:

Triangle 121.
0. Harlingen		50° 58′ 15″ 238
1. Leeuwarden		82° 47′ 15″ 351
2. Ballum		46° 14′ 27″ 202
Triangle 122.
3. Harlingen		51° 05′ 39″ 717
4. Sneek		70° 48′ 33″ 445
5. Leeuwarden		58° 05′ 48″ 707
Triangle 123.
6. Sneek		49° 30′ 40″ 051
7. Drachten		42° 52′ 59″ 382
8. Leeuwarden		87° 36′ 21″ 057
Triangle 124.
9. Sneek		45° 36′ 07″ 492
10. Oldeholtpade		67° 52′ 00″ 048
11. Drachten		66° 31′ 56″ 513
Triangle 125.
12. Drachten		53° 55′ 24″ 745
13. Oldeholtpade		47° 48′  52″ 580
14. Oosterwolde		78° 15′  42″ 347
Triangle 127.
15. Leeuwarden		59° 24′  00″ 645
16. Dockum		76° 34′  09″ 021
17. Ballum		44° 01′  51″ 040
Triangle 128.
18. Leeuwarden		72° 06′  32″ 043
19. Drachten		46° 53′  27″ 163
20. Dockum		61° 00′  04″ 494
Triangle 131.
21. Dockum		57° 01′  55″ 292
22. Drachten		83° 33′  14″ 515
23. Gröningen		39° 24′  52″ 397
Triangle 132.
24. Oosterwolde		81° 54′  17″ 447
25. Gröningen		31° 52′  49″ 094
26. Drachten		66° 12′  57″ 246

The consideration of these triangles shows that the twenty-seven angles directly provided by observation have thirteen necessary relationships among them, namely: two of the first kind, nine of the second kind, and two of the third kind. However, it is not necessary to write down all these equations here in their final form, because for the calculation we only need the quantities which were denoted in the general theory by ${\mathcal {A}}$, $a$, $a'$, etc., ${\mathcal {B}}$, $b$, $b'$, etc. etc. Therefore, we can immediately write down the equations (13), which bring out these quantities. Instead of $\varepsilon$, $\varepsilon '$, etc., we will simply write here $(0),$ $(1),$ $(2),$ etc.

In this way, the two equations of the first kind correspond to the following:

${\begin{alignedat}{4}(1)&{}+{}&(5)&{}+{}&(8)&{}+{}&(15)+(18)=-2''{,}197,\\(7)+(11)&{}+{}&(12)&{}+{}&(19)&{}+{}&(22)+(26)=-0''{,}436.\end{alignedat}}$

Next, we find, for the spheroidal excesses of the nine triangles: $1''.749;$ $1''.147;$ $1''.243;$ $1''.698;$ $0''.873;$ $1''.167;$ $1''.104;$ $2''.161;$ $1''.403.$ Then, we have the conditional equation of the second kind:

$v^{(0)}+v^{(1)}+v^{(2)}-180^{\circ }0'1''{,}749=0,$

and so on for the others, and we have the following nine equations:

${\begin{alignedat}{3}(0)&{}+{}&(1)&{}+{}&(2)&{}={}-3''\!{,}958,\\(3)&{}+{}&(4)&{}+{}&(5)&{}={}+0\,{,}722,\\(6)&{}+{}&(7)&{}+{}&(8)&{}={}-0\,{,}753,\\(9)&{}+{}&(10)&{}+{}&(11)&{}={}+2\,{,}355,\\(12)&{}+{}&(13)&{}+{}&(14)&{}={}-1\,{,}201,\\(15)&{}+{}&(16)&{}+{}&(17)&{}={}-0\,{,}461,\\(18)&{}+{}&(19)&{}+{}&(20)&{}={}+2\,{,}596,\\(21)&{}+{}&(22)&{}+{}&(23)&{}={}+0\,{,}043,\\(24)&{}+{}&(25)&{}+{}&(26)&{}={}-0\,{,}616.\end{alignedat}}$

The conditional equations of the third kind are more easily expressed using logarithms. The first one is:

${\begin{alignedat}{3}\log \sin \,(v^{0}-0''{,}583)&{}-{}\log \sin \,(v^{2}&{}-{}0''{,}583)&{}-{}\log \sin \,(v^{3}&{}-{}0''{,}382)\\&{}+{}\log \sin \,(v^{4}&{}-{}0''{,}382)&{}-{}\log \sin \,(v^{6}&{}-{}0''{,}414)\\&{}+{}\log \sin \,(v^{7}&{}-{}0''{,}414)&{}-{}\log \sin \,(v^{16}&{}-{}0''{,}389)\\&{}+{}\log \sin \,(v^{17}&{}-{}0''{,}389)&{}-{}\log \sin \,(v^{19}&{}-{}0''{,}368)\\&{}+{}\log \sin \,(v^{20}&{}-{}0''{,}368)&={}0.\end{alignedat}}$

It seems unnecessary to develop the other one in finite form. To these two equations correspond the following, in which the coefficients refer to the seventh decimal place of Brigg's logarithms:

${\begin{aligned}&17{,}068\,(0)-20{,}174\,(2)-16{,}993\,(3)+7{,}328\,(4)-17{,}976\,(6)\\&\qquad {}+{}22{,}672\,(7)-5{,}028\,(16)+21{,}780\,(17)-19{,}710\,(19)\\&\qquad {}+{}11{,}671\,(20)=-371,\\&17{,}976\,(6)-0{,}880\,(8)-20{,}617\,(9)+8{,}564\,(10)-19{,}082\,(13)\\&\qquad {}+{}4{,}375\,(14)+6{,}798\,(18)-11{,}671\,(20)+13{,}657\,(21)\\&\qquad {}-{}25{,}620\,(23)-2{,}995\,(24)+33{,}854\,(25)=+370,\end{aligned}}$

No reason leads us to attribute unequal weights to the various observations, so we will assume

$p^{0}=p^{1}=p^{2}=\ldots =1.$

Denoting the correlates of the condition equations, in the same order as these equations were written, by

${A},$ ${B},$ ${C},$ ${D},$ ${E},$ ${F},$ ${G},$ ${H},$ ${I},$ ${K},$ ${L},$ ${M},$ ${N},$

we determine them by the following equations:

${\begin{alignedat}{2}&{}-{}&2''{,}197&{}={}5\,{A}+{C}+{D}+{E}+{H}+{I}+5{,}917\,{N},\\&{}-{}&0{,}436&{}={}6\,{B}+{E}+{F}+{G}+{I}+{K}+{L}+2{,}962\,{M},\\&{}-{}&3{,}958&{}={}{A}+3\,{C}-3{,}106\,{M},\\&{}+{}&0{,}722&{}={}{A}+3\,{D}-9{,}665\,{M},\\&{}-{}&0{,}753&{}={}{A}+{B}+3\,{E}+4{,}696\,{M}+17{,}096\,{N},\\&{}+{}&2{,}355&{}={}{B}+3\,{F}-12{,}053\,{N},\\&{}-{}&1{,}201&{}={}{B}+3\,{G}-14{,}707\,{N},\\&{}-{}&0{,}461&{}={}{A}+3\,{H}+16{,}752\,{M},\\&{}+{}&2{,}596&{}={}{A}+{B}+3\,{I}-8{,}039\,{M}-4{,}874\,{N},\\&{}+{}&0{,}043&{}={}{B}+3\,{K}-11{,}963\,{N},\\&{}-{}&0{,}616&{}={}{B}+3\,{L}+30{,}859\,{N},\\&{}-{}&371&{}={}2{,}962\,{B}-3{,}106\,{C}-9{,}665\,{D}+4{,}696\,{E}\\&&&\qquad \qquad \;\,+16{,}752\,{H}-8{,}039\,{I}+2902{,}27\,{M}\\&&&\qquad \qquad \;\,-459{,}33\,{N},\\&{}+{}&370&{}={}5{,}917\,{A}+17{,}096\,{E}-12{,}053\,{F}-14{,}707\,{G}\\&&&\qquad \qquad \;\,-4{,}874\,{I}-11{,}963\,{K}+30{,}859\,{L}\\&&&\qquad \qquad \;\,-459{,}33\,{M}+3385{,}96\,{N}\end{alignedat}}$

We then deduce by elimination:

${\begin{alignedat}{2}&{A}&{}={}&-0{,}598\\&{B}&{}={}&-0{,}255\\&{C}&{}={}&-1{,}234\\&{D}&{}={}&+0{,}086\\&{E}&{}={}&-0{,}447\\&{F}&{}={}&+1{,}351\\&{G}&{}={}&+0{,}271\end{alignedat}}$

${\begin{alignedat}{2}&{H}&{}={}&+0{,}659\\&{I}&{}={}&+1{,}050\\&{K}&{}={}&+0{,}577\\&{L}&{}={}&-1{,}351\\&{M}&{}={}&-0{,}109792\\&{N}&{}={}&+0{,}119681.\\\quad \end{alignedat}}$

Finally, the most plausible errors are given by the formulas

${\begin{array}{l}{\begin{aligned}(0)&={C}+17{,}068\,{M}\\(1)&={A}+{C}\\(2)&={C}-20{,}174\,{M}\\(3)&={D}-16{,}993\,{M},\end{aligned}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

and we obtain the following numerical values, to which we add, for comparison, the corrections adopted by de Krayenhof:

	de Kr.		de Kr.
${\begin{alignedat}{2}&(0)&{}={}-{}&3''\!{,}108\\&(1)&{}={}-{}&1\,{,}832\\&(2)&{}={}+{}&0\,{,}981\\&(3)&{}={}+{}&1\,{,}952\\&(4)&{}={}-{}&0\,{,}719\\&(5)&{}={}-{}&0\,{,}512\\&(6)&{}={}+{}&3\,{,}648\\&(7)&{}={}-{}&3\,{,}221\\&(8)&{}={}-{}&1\,{,}180\\&(9)&{}={}-{}&1\,{,}116\\&(10)&{}={}+{}&2\,{,}376\\&(11)&{}={}+{}&1\,{,}096\\&(12)&{}={}+{}&0\,{,}016\\&(13)&{}={}-{}&2\,{,}013\\\end{alignedat}}$	${\begin{aligned}{}-{}&2''\!{,}090\\[0.125ex]{}+{}&0\,{,}116\\[0.125ex]{}-{}&1\,{,}982\\[0.125ex]{}+{}&1\,{,}722\\[0.125ex]{}+{}&2\,{,}848\\[0.125ex]{}-{}&3\,{,}848\\[0.125ex]{}-{}&0\,{,}137\\[0.125ex]{}+{}&1\,{,}000\\[0.125ex]{}-{}&1\,{,}614\\[0.125ex]&\,0\\[0.125ex]{}+{}&5\,{,}928\\[0.125ex]{}-{}&3\,{,}570\\[0.125ex]{}+{}&2\,{,}414\\[0.125ex]{}-{}&6\,{,}014\\[0.125ex]\end{aligned}}$	${\begin{alignedat}{2}&(14)&{}={}+{}&0''\!{,}795\\&(15)&{}={}+{}&0\,{,}061\\&(16)&{}={}+{}&1\,{,}211\\&(17)&{}={}-{}&1\,{,}732\\&(18)&{}={}+{}&1\,{,}265\\&(19)&{}={}+{}&2\,{,}959\\&(20)&{}={}-{}&1\,{,}628\\&(21)&{}={}+{}&2\,{,}211\\&(22)&{}={}+{}&0\,{,}322\\&(23)&{}={}-{}&2\,{,}489\\&(24)&{}={}-{}&1\,{,}709\\&(25)&{}={}+{}&2\,{,}701\\&(26)&{}={}-{}&1\,{,}606\\&&&\quad \end{alignedat}}$	${\begin{aligned}{}+{}&2''\!{,}400\\[0.125ex]{}+{}&1\,{,}273\\[0.125ex]{}+{}&5\,{,}945\\[0.125ex]{}-{}&7\,{,}674\\[0.125ex]{}+{}&1\,{,}876\\[0.125ex]{}+{}&6\,{,}251\\[0.125ex]{}-{}&5\,{,}530\\[0.125ex]{}+{}&3\,{,}486\\[0.125ex]{}-{}&3\,{,}454\\[0.125ex]&\,0\\[0.125ex]{}+{}&0\,{,}400\\[0.125ex]{}+{}&2\,{,}054\\[0.125ex]{}-{}&3\,{,}077\\[0.125ex]&\quad \end{aligned}}$

|}

The sum of the squares of our corrections is 97.8845; the average error, as indicated by the 27 observed angles, is therefore

${\sqrt {\frac {97{,}8845}{13}}}=2''{,}7440.$

The sum of the squares of the corrections by de Krayenhof is 341.4201.

24.

The triangles whose vertices, in the triangulation of Hanover, are placed at Falkenberg, Breithorn, Hauselberg, Wulfsode, and Wilsede, will provide us with a second example.

The following directions have been observed:^[2]

At the station in Falkenberg
0. Wilsede		187° 47′ 30″ 311
1. Wulfsode		225° 09′ 39″ 676
2. Hauselberg		266° 13′ 56″ 239
3. Breithorn		274° 14′ 43″ 634
À la station de Breithorn.
4. Falkenberg		094° 33′ 40″ 755
5. Hauselberg		122° 51′ 23″ 054
6. Wilsede		150° 18′ 35″ 100
À la station de Hauselberg.
7. Falkenberg		086° 29′ 06″ 872
8. Wilsede		154° 37′ 09″ 624
9. Wulfsode		189° 02′ 56″ 376
10. Breithorn		302° 47′ 37″ 732
At the station in Wulfsode.
11. Hauselberg		009° 05′ 36″ 593
12. Falkenberg		045° 27′ 33″ 556
13. Wilsede		118° 44′ 13″ 159
At the station in Wilsede.
14. Falkenberg		007° 51′ 01″ 027
15. Wulfsode		298° 29′ 49″ 519
16. Breithorn		330° 03′ 07″ 392
17. Hauselberg		334° 25′ 26″ 746

These observations allow us to form seven triangles.

Triangle I.
Falkenberg		008° 00′ 47″ 395
Breithorn		028° 17′ 42″ 299
Hauselberg		143° 41′ 29″ 140
Triangle II.
Falkenberg		086° 27′ 13″ 323
Breithorn		055° 44′ 54″ 345
Wilsede		037° 47′ 53″ 635
Triangle III.
Falkenberg		041° 04′ 16″ 563
Hauselberg		102° 33′ 49″ 504
Wulfsode		036° 21′ 56″ 963
Triangle IV.
Falkenberg		078° 26′ 25″ 928
Hauselberg		068° 08′ 02″ 752
Wilsede		035° 25′ 34″ 281
Triangle V.
Falkenberg		037° 22′ 09″ 365
Wulfsode		073° 16′ 39″ 603
Wilsede		069° 21′ 11″ 508
Triangle VI.
Breithorn		027° 27′ 12″ 046
Hauselberg		148° 10′ 28″ 108
Wilsede		004° 22′ 19″ 354
Triangle VII.
Hauselberg		034° 25′ 46″ 752
Wulfsode		109° 38′ 36″ 566
Wilsede		035° 55′ 37″ 227

We have here seven conditional equations of the second kind (there is obviously no need to form any of the first); to form them, we must first find the spheroidal excesses of the seven triangles, and for this it is essential to know the length of one side. The one that connects Wilsede to Wulfsode is 22877.94 meters. We conclude that the spheroidal excesses of the various triangles are:

I...$0''202;$ II...$2''442;$ III...$1'257;$ IV...$1''919;$ V...$1''957;$ VI...$0''321;$ VII...$1''295.$

If we denote the angles that determine the directions indicated above by $v^{(0)},$ $v^{(1)},$ $v^{(2)},$ $v^{(3)},$ etc., then the angles of the first triangle, marked with the same indices, will be

${\begin{aligned}v^{(3)}-v^{(2)}&,&v^{(5)}-v^{(4)}&,&360+v^{(7)}-v^{(10)}&,\end{aligned}}$

and the first condition equation is therefore

${}-{}v^{(2)}+v^{(3)}-v^{(4)}+v^{(5)}+v^{(7)}-v^{(10)}+179^{\circ }59'59''{,}798=0.$

The remaining six triangles will yield six similar equations, but a little attention shows that these equations are not independent. Indeed, the second is identical to the sum of the first, fourth, and sixth; the sum of the third and fifth is identical to that of the fourth and seventh: for this reason we will ignore the second and fifth. Instead of writing the remaining equations in finite form, we will write here the corresponding equations (13), with symbols $(0),$ $(1),$ $(2),$ etc. taking the place of $\varepsilon ,$ $\varepsilon ',$ $\varepsilon '',$ etc.:

${\begin{alignedat}{8}&-&1''\!{,}368&{}={}&{}-{}&(2)&{}+{}&(3)&{}-{}&(4)&{}+{}&(5)&{}+{}&(7)&{}-{}&(10),\\&+&1\,{,}773&{}={}&{}-{}&(1)&{}+{}&(2)&{}-{}&(7)&{}+{}&(9)&{}-{}&(11)&{}+{}&(12),\\&+&1\,{,}042&{}={}&{}-{}&(0)&{}+{}&(2)&{}-{}&(7)&{}+{}&(8)&{}+{}&(14)&{}-{}&(17),\\&-&0\,{,}813&{}={}&{}-{}&(5)&{}+{}&(6)&{}-{}&(8)&{}+{}&(10)&{}-{}&(16)&{}+{}&(17),\\&-&0\,{,}750&{}={}&{}-{}&(8)&{}+{}&(9)&{}-{}&(11)&{}+{}&(13)&{}-{}&(15)&{}+{}&(17).\end{alignedat}}$

We can obtain eight conditional equations of the third kind from the triangles of the system, to which end one may combine three of the four triangles I, II, IV, VI, or three of the triangles III, IV, V, VII; however, a little attention shows that it suffices to consider two belonging to each of these two systems of triangles, and those will encompass all the others. Thus, for the sixth conditional equation, we have

${\begin{alignedat}{6}&\log \sin(v^{(3)}&{}-{}&v^{(2)}&{}-{}&0''{,}067)&{}-{}&\log \sin(v^{(5)}&{}-{}&v^{(4)}&{}-{}&0''{,}067)\\{}+{}&\log \sin(v^{(14)}&{}-{}&v^{(17)}&{}-{}&0''{,}640)&{}-{}&\log \sin(v^{(2)}&{}-{}&v^{(0)}&{}-{}&0''{,}640)\\{}+{}&\log \sin(v^{(6)}&{}-{}&v^{(5)}&{}-{}&0''{,}107)&{}-{}&\log \sin(v^{(17)}&{}-{}&v^{(16)}&{}-{}&0''{,}107)=0,\end{alignedat}}$

and for the seventh,

${\begin{alignedat}{6}&\log \sin(v^{(2)}&{}-{}&v^{(1)}&{}-{}&0''{,}419)&{}-{}&\log \sin(v^{(12)}&{}-{}&v^{(11)}&{}-{}&0''{,}419)\\{}+{}&\log \sin(v^{(14)}&{}-{}&v^{(17)}&{}-{}&0''{,}640)&{}-{}&\log \sin(v^{(2)}&{}-{}&v^{(0)}&{}-{}&0''{,}640)\\{}+{}&\log \sin(v^{(13)}&{}-{}&v^{(11)}&{}-{}&0''{,}432)&{}-{}&\log \sin(v^{(17)}&{}-{}&v^{(15)}&{}-{}&0''{,}432)=0,\end{alignedat}}$

to which correspond the equations (13),

${\begin{alignedat}{1}{}+{}25{}={}&{}+{}4{,}31\,(0)-153{,}88\,(2)+149{,}57\,(3)\\&{}+{}39{,}11\,(4)-79{,}64\,(5)+40{,}53\,(6)\\&{}+{}31{,}90\,(14)+275{,}39\,(16)-307{,}29\,(17),\\[0.75ex]{}-{}\;3{}={}&{}+{}4{,}31\,(0)-24{,}16\,(1)+19{,}85\,(2)\\&{}+{}36{,}11\,(11)-28{,}59\,(12)-7{,}52\,(13)\\&{}+{}31{,}90\,(14)+29{,}06\,(15)-60{,}96\,(17).\end{alignedat}}$

If we attribute the same certainty to the various directions, setting

$p^{(0)}=p^{(1)}=p^{(2)}=\ldots =1$,

and if we denote the correlates of the seven conditional equations by ${A},$ ${B},$ ${C},$ ${D},$ ${E},$ ${F},$ ${G},$ then their determination will depend on the following equations:

${\begin{aligned}-1.368=&+6\,{A}-2\,{B}-2\,{C}-2\,{D}+184.72\,{F}-19.85\,{G},\\+1.773=&-2\,{A}+6\,{B}+2\,{C}+2\,{E}-153.88\,{F}-20.69\,{G},\\+1.042=&-2\,{A}+2\,{B}+6\,{C}-2\,{D}-2\,{E}+181.00\,{F}\\&+108.40\,{G},\\-0.813=&-2\,{A}-2\,{C}+6\,{D}+2\,{E}-462.51\,{F}-60.96\,{G},\\-0.750=&+2\,{B}-2\,{C}+2\,{D}+6\,{E}-307.29\,{F}-133.65\,{G},\\+25=&+184.72\,{A}-153.88\,{B}+181.00\,{C}-462.51\,{D}\\&-307.29\,{E}+224868\,{F}+16694.1\,{G},\\-\;3=&-19.85\,{A}-20.69\,{B}+108.40\,{C}-60.96\,{D}\\&-133.65\,{E}+16694.1\,{F}+8752.39\,{G}.\end{aligned}}$

From this, by elimination, we deduce

${\begin{alignedat}{2}&{A}&{}={}&-0.225,\\&{B}&{}={}&+0.344,\\&{C}&{}={}&-0.088,\\&{D}&{}={}&-0.171,\\&{E}&{}={}&-0.323,\\&{F}&{}={}&+0.000215915,\\&{G}&{}={}&-0.005474620,\end{alignedat}}$

and the most probable errors are given by the formulas

${\begin{array}{l}{\begin{aligned}(0)&={}-{}{C}+4.31\,{F}+4.31\,{G},\\(1)&={}-{}{B}-24.16\,{G},\\(2)&={}-{}{A}+{B}+{C}-153.88\,{F}+19.85\,{G},\end{aligned}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

from which we derive the following numerical values:

${\begin{alignedat}{2}&(0)&{}={}&+0''\!{,}065\\&(1)&{}={}&-0{,}212\\&(2)&{}={}&+0{,}339\\&(3)&{}={}&-0{,}193\\&(4)&{}={}&+0{,}233\\&(5)&{}={}&-0{,}071\\&(6)&{}={}&-0{,}162\\&(7)&{}={}&-0{,}481\\&(8)&{}={}&+0{,}406\end{alignedat}}$

${\begin{alignedat}{2}&(9)&{}={}&+0''\!{,}021\\&(10)&{}={}&+0{,}054\\&(11)&{}={}&-0{,}219\\&(12)&{}={}&+0{,}501\\&(13)&{}={}&-0{,}282\\&(14)&{}={}&-0{,}256\\&(15)&{}={}&+0{,}164\\&(16)&{}={}&+0{,}230\\&(17)&{}={}&-0{,}139.\end{alignedat}}$

The sum of the squares of these errors is equal to $1{,}2288;$ hence, the average error resulting from the 18 observed directions is

$={\sqrt {\frac {1{,}2288}{7}}}=0''{,}4190.$

25.

To give an example of the last part of our theory, let us determine the precision with which the compensated observations determine the Falkenberg-Breithorn side, using the Wilsede-Wulfsode side. The function $u,$ by which it is expressed, in this case, is:

$u=22877^{\text{m}}94\times {\frac {\sin(v^{(13)}-v^{(12)}-0''\!,652)\cdot \sin(v^{(14)}-v^{(16)}-0''\!,814)}{\sin(v^{(1)}-v^{(0)}-0''\!,652)\cdot \sin(v^{(6)}-v^{(4)}-0''\!,814)}};$

Its value, deduced from the corrected observations, is

$26766^{\text{m}}68.$

Differentiation of this equation provides the following, with ${d}v^{(0)},$ ${d}v^{(1)},$ etc. expressed in seconds,

${\begin{aligned}{d}u&=\left[0^{\text{m}}{,}16991\,({d}v^{(0)})-{d}v^{(1)}\right]+0^{\text{m}}{,}08836\,({d}v^{(4)}-{d}v^{(6)})\\&{}-{}\left[0^{\text{m}}{,}03899\,({d}v^{(12)})-{d}v^{(13)}\right]+0^{\text{m}}{,}16731\,({d}v^{(14)}-{d}v^{(16)})\,;\end{aligned}}$

From this we find that:

${\begin{alignedat}{3}(al){}={}&{}-{}&0&{,}08836,\\(bl){}={}&{}+{}&0&{,}13092,\\(cl){}={}&{}-{}&0&{,}00260,\\(dl){}={}&{}+{}&0&{,}07895,\\(el){}={}&{}+{}&0&{,}03899,\\(fl){}={}&{}-{}&40&{,}1315,\\(gl){}={}&{}+{}&10&{,}9957,\\(ll){}={}&{}+{}&0&{,}13238,\end{alignedat}}$

Taking the meter as the unit of length, the methods indicated above yield,

${\frac {1}{P}}=0{,}08329\quad$ or $\quad {P}=12{,}006.$

We conclude that the mean error to be feared in the value of the Falkenberg-Breithorn side is $0{,}2886m$ (where $m$ denotes the average error to be feared in the observed directions, this error being expressed in seconds), and consequently, if we adopt the value of $m$ announced above, this mean error to be feared is $0^{\text{m}}1209.$

Moreover, inspecting the system of triangles immediately shows that one cans completely disregard the Hauselberg station without breaking the network that connects the other four. However, it would not be good practice to eliminate the operations related to this point,^[3] as they certainly contribute to increasing the precision of the whole. To show more clearly what increase in precision results from this, we will finish by again performing the calculation, after excluding all results related to the Hauselberg point. Of the eighteen directions mentioned above, eight then cease to be used, and the most plausible errors on those that remain are:

${\begin{alignedat}{2}&(0)&{}={}&+0''\!{,}327\\&(1)&{}={}&-0\,{,}206\\&(3)&{}={}&-0\,{,}121\\&(4)&{}={}&+0\,{,}121\\&(6)&{}={}&-0\,{,}121\end{alignedat}}$

${\begin{alignedat}{2}&(12)&{}={}&+0''\!{,}206\\&(13)&{}={}&-0\,{,}206\\&(14)&{}={}&+0\,{,}327\\&(15)&{}={}&+0\,{,}206\\&(16)&{}={}&+0\,{,}121\end{alignedat}}$

The value of the Falkenberg-Breithorn side then becomes $26766^{\text{m}}{,}63,$ a result little different from that obtained earlier. But calculating the weight yields

${\frac {1}{P}}=0{,}13082\quad$ or $\quad {P}=7{,}644,$

and the mean error to be feared, in meters, is

$=0{,}36169m=0^{\text{m}}1515.$

We see that by adding the operations related to Hauselberg, the weight of the determination of the Falkenberg-Breithorn side is increased in the ratio of $7{,}644$ to $12{,}006,$ that is, in the ratio of unity to $1{,}571.$

1. In the previous calculations, it was sufficient to use three letters from each series to clarify the formulas' law; it seemed necessary here to include a fourth one to make the algorithm more evident. We have followed the series of letters $a,$ $b,$ $c;$ $A,$ $B,$ $C;$ ${\mathfrak {A}},$ ${\mathfrak {B}},$ ${\mathfrak {C}};$ in the natural way, with $d,$ $D,$ ${\mathfrak {D}}.$ However, lacking the required letters, we have followed the series $x,$ $y,$ $z$ with $w,$ and the series $\xi ,$ $\eta ,$ $\zeta$ with $\varphi$.
2. The starting points to which the individual directions are referred are here regarded as arbitrary, although in reality they coincide with the meridian lines of the stations. The observations will be published in full at a later date; meanwhile, the figure can be found in "Astronomische Nachrichten" Vol. I. p. 441
3. The greater part of these observations had already been made before Breithorn's point was discovered and incorporated into the system.