# Calculus of Variations/CHAPTER XIII

CHAPTER XIII: STATEMENT OF THE PROBLEM. DERIVATION OF THE NECESSARY CONDITIONS.
* 179 The general problem stated.
* 180 Existence of substitutions by which one integral remains unchanged while the other is caused to vary. An exceptional case.
* 181 Case of two variables. Convergence of the series that appear.
* 182 The nature of the substitutions that have been introduced.
* 183 Formation of certain quotients which depend only upon the nature of the curve.
* 184 Generalization, in which several integrals are to retain fixed values.
* 185 The quotient of two definite integrals being denoted by $\lambda$, it is shown that $\lambda$ has the same constant value for the whole curve.
* 186 The differential equation $G^{(0)}-\lambda G^{(1)}=0$.
* 187 Extension of the theorem of Article 97.
* 188 Discontinuities, etc.
* 189 The second variation: the three conditions formulated in Article 135 are also necessary here.
Article 179. The nature of many problems whicli arise in the Calculus of Variations presents subsidiary conditions which limit the arbitrariness that we have hitherto employed in the indefinitely small variations of the analytical structure. Such problems are the most difficult and at the same time the most interesting that occur. These last conditions which enter into the requirement for a maximum or a minimum are in general of a double nature. On the one hand, it may be proposed that among the variables there are to exist equations of condition, as indicated in Arts. 176 and 177. On the other hand, we may require that the maximum or the minimum in question satisfy a further condition, viz., it must cause another given integral to have a prescribed value. Such cases are usually called Relative Maxima and Minima.
If we limit our discussion to the region of two variables, then the problem which we have to consider may be expressed as follows (cf. Art. 17):
Let $F^{(0)}(x,y,x',y')$ and $F^{(1)}(x,y,x',y')$ be two functions of the same nature as the function $F(x,y,x',y')$ hitherto treated. The variables $x$ and $y$ are to be so determined as one-valued functions of $t$ that the curve defined through the equations $x=x(t),y=y(t)$ will cause the integral
    $1)\qquad I^{(0)}=\int _{t_{0}}^{t_{1}}F^{(0)}(x,y,x',y'){\text{d}}t$
to be a maximum or a minimum, while at the same time for the same equations the integral
    $2)\qquad I^{(1)}=\int _{t_{0}}^{t_{1}}F^{(1)}(x,y,x',y'){\text{d}}t$
will have a prescribed value; that is, for every indefinitely small variation of the curve for which the second integral retains its sign unaltered, the first integral, according as a maximum or a minimum is to enter, must be continuously smaller or continuously greater than it is for the curve $x=x(t),y=y(t)$.
Article 180. We must first show that it is possible to represent analytically the variations of a curve for which the integral $I^{(1)}$ retains a constant value.
In the place of the variables $x,y$ let us make the substitution $x+{\bar {\xi }},y+{\bar {\eta }}$. The variation of the second integral is accordingly
    $3)\qquad \Delta I^{(1)}=\int _{t_{0}}^{t_{1}}G^{(1)}{\bar {w}}{\text{d}}t+\int _{t_{0}}^{t_{1}}({\bar {\xi }},{\bar {\eta }},{\bar {\xi }}',{\bar {\eta }}')_{2}{\text{d}}t$
where $({\bar {\xi }},{\bar {\eta }},{\bar {\xi }}',{\bar {\eta }}')_{2}$ denotes that the terms within the brackets are of the second and higher dimensions in ${\bar {\xi }},{\bar {\eta }},{\bar {\xi }}',{\bar {\eta }}'$.
We have so to determine ${\bar {\xi }}$ and ${\bar {\eta }}$; that $\Delta I^{(1)}=0$. For this purpose we write
    $4)\qquad {\bar {\xi }}=\epsilon \xi +\epsilon _{1}\xi _{1}+\epsilon _{2}\xi _{2}+cdots\qquad {\bar {\eta }}=\epsilon \eta +\epsilon _{1}\eta _{1}+\epsilon _{2}\eta _{2}$
where $\epsilon ,\epsilon _{1},\ldots$ are arbitrary constants and the functions $\xi ,\xi _{1},\ldots ,\eta ,\eta _{1},\ldots$ are functions similar to the quantities $$ of the preceding Chapters and vanish for $t=t_{0}$ and $t=t_{1}$. Now write
    $w_{i}=y'\xi _{i}-x'\eta _{i}$
and
    ${\bar {w}}=y'{\bar {\xi }}-x'{\bar {\eta }}=\epsilon w+\epsilon _{1}w_{1}+\epsilon _{2}w_{2}+\cdots$
Hence, from 3) we have
    $\Delta I^{(1)}=\epsilon \int _{t_{0}}^{t_{1}}G^{(1)}w{\text{d}}t+\epsilon _{1}\int _{t_{0}}^{t_{1}}G^{(1)}w_{1}{\text{d}}t+\cdots +(\epsilon ,\epsilon _{1},\cdots )_{2}$
If we write
    $5)\qquad W_{i}^{(1)}=\int _{t_{0}}^{t_{1}}G^{(1)}w_{i}{\text{d}}t$
it follows that
    $\Delta I^{(1)}=\epsilon W^{(1)}+\epsilon _{1}W_{1}^{(1)}+\epsilon _{2}W_{2}^{(1)}+\cdots +(\epsilon ,\epsilon _{1},\epsilon _{2},\ldots )_{2}$
The functions $W_{i}^{(i)}$ are completely determined as soon as definite values are given to $\xi ,\xi _{1},\ldots$; and, in order that $\Delta I^{(1)}=0$, it is necessary that
    $(A)\qquad \epsilon W^{(1)}+\epsilon _{1}W_{1}^{(1)}+\epsilon _{2}W_{2}^{(1)}+\cdots +(\epsilon ,\epsilon _{1},\epsilon _{2},\ldots )_{2}=0$
If any of the quantities $W_{i}^{(1)}$, for example $W_{\lambda }^{(1)}$, are different from zero, we are able to express $\epsilon _{\lambda }$ in a power-series of the remaining $\epsilon$'s, when these quantities have been chosen sufficiently small. The equation $\Delta I^{(1)}=0$ may consequently be satisfied for sufficiently small systems of values of the $\epsilon$'s.
Substitute one of these systems of values in 4) and it is seen that indefinitely small variations of the curve $x=x(t),y=y(t)$ exist for which the integral $I^{(1)}$ remains unaltered. These variations may be analytically represented (see the next Article).
This proof is deficient in the case where all the quantities $W^{(1)},W_{1}^{(1)},\ldots$ are zero for all values of $\xi _{i},\eta _{i}$, however $\xi _{i},\eta _{i}$ may have been chosen. When this is the case, $G^{(1)}$ must be zero along the whole curve. But this is one of the necessary conditions that the integral $I^{(1)}$have a maximum or a minimum value.
If, then, for the curve which is derived through the solution of the differential equation $G^{(0)}=0$ there also enters a maximum or a minimum value of the integral $I^{(1)}$ and consequently $G^{(1)}=0$, it is in general not possible so to vary the curve that the second integral remains unaltered.
This case is excluded from the present discussion, and is left for special investigation in each particular problem.
Article 181. Let us limit ourselves for the present to the simplest case where
    ${\bar {\xi }}=\epsilon \xi +\epsilon _{1}\xi _{1}\qquad {\bar {\eta }}=\epsilon \eta +\epsilon _{1}\eta _{1}$
and if we denote the integrals in the expansion of $\Delta I^{(1)}$ that are associated with the coefficients $e^{i}e_{1}^{j}$ by $W_{ij}^{(1)}$, the equation correresponding to (A) of the last article is
    $(A^{a})\qquad 0=W_{10}^{(1)}+\epsilon _{1}W_{01}^{(1)}+\epsilon ^{2}W_{20}^{(1)}+\epsilon \epsilon _{1}W_{11}^{(1)}+\epsilon _{1}^{2}W_{02}^{(1)}+\cdots$
which series we suppose convergent for sufficiently small values of $\epsilon$ and $\epsilon _{1}$.
Suppose next we express $\epsilon _{1}$ in terms of $\epsilon$ by the series
    $\epsilon _{1}=h_{1}\epsilon +h_{2}\epsilon ^{2}+h_{3}\epsilon ^{3}+\cdots$
Then, when this value of $\epsilon _{1}$ is substituted in $A^{a}$, by equating the coefficients of the different powers of $\epsilon$ to zero, we have
    $W_{10}^{(1)}+h_{1}W_{01}^{(1)}=0$
    $W_{20}^{(1)}+h_{1}W_{11}^{(1)}+h_{1}^{2}W_{02}^{(1)}+h_{2}W_{01}^{(1)}=0$
Hence, denoting the quotients $-{\frac {W_{ij}^{(1)}}{W_{01}^{(1)}}}$ by $V_{ij}$, where $W_{01}^{(1)}\neq 0$, we have
    $h_{1}=V_{10}$
    $h_{2}=V_{20}+h_{1}V_{11}+h_{1}^{2}V_{02}$
Further, the equation $(A^{a})$ may be written
    $(A^{b})\qquad \epsilon _{1}=\epsilon V_{10}+\epsilon ^{2}V_{20}+\epsilon \epsilon _{1}V_{11}+\epsilon _{1}^{2}V_{02}+\cdots$
Let us compare this series with the series
    $(C)\qquad \epsilon _{1}={\frac {g}{1-\left({\frac {\epsilon }{r}}+{\frac {\epsilon _{1}}{r_{1}}}\right)}}-g-g{\frac {\epsilon _{1}}{r_{1}}}$
    $\qquad =g\left[{\frac {\epsilon }{r}}+\left({\frac {\epsilon }{r}}+{\frac {\epsilon _{1}}{r_{1}}}\right)^{2}+\left({\frac {\epsilon }{r}}+{\frac {\epsilon _{1}}{r_{1}}}\right)^{3}+\cdots \right]$
Suppose from this series we have $\epsilon _{1}$ expressed in terms of $\epsilon$ in the form
    $(B^{b})\qquad \epsilon _{1}=h_{1}'\epsilon +h_{2}'\epsilon ^{2}+h_{3}'\epsilon ^{3}+\cdots$
where the $h'$s have been derived from the coefficients of powers of $\epsilon$ and $\epsilon _{1}$ as the $h$s in $(B)$ are formed from the coefficients $V$ in $(A^{b})$.
The series $(B^{b})$ is convergent for
    $\left|{\frac {\epsilon }{r}}\right|+\left|{\frac {\epsilon _{1}}{r_{1}}}\right|<1$
If, then, the coefficients $V$ of $(A^{b})$ are in absolute value less than the corresponding coefficients in $(C)$, the coefficients $h$ in $(B)$ are less in absolute value than the coefficients $h'$ in $(B^{b})$, and therefore the series $(B)$ is convergent.
Now the coefficients of $\epsilon ^{k}\epsilon _{1}^{\mu }$ in $(A^{b})$ and $(C)$ are respectively
    $V_{k,\mu }$ and ${\binom {k+\mu }{k}}{\frac {g}{r^{k}r_{1}^{\mu }}}$
where the symbol ${\binom {m}{n}}$ denotes ${\frac {m(m-1)\cdots (m-n+1)}{n!}}$. Hence for sufficiently small values of $r$ and $r_{1}$, if
    $\left|{\frac {\epsilon }{r}}\right|+\left|{\frac {\epsilon _{1}}{r_{1}}}\right|<1$
and
    $V_{k,\mu }<{\binom {k+\mu }{k}}{\frac {g}{r^{k}r_{1}^{\mu }}}$
the series $(B)$ is convergent, and when substituted in the expression for $\Delta I^{(1)}$ causes this expression to vanish.
Article 182. The expression for $\epsilon _{1}$ as a function of $\epsilon$ is had from the relation
    $\epsilon _{1}={\frac {g}{1-\left({\frac {\epsilon }{r}}+{\frac {\epsilon _{1}}{r_{1}}}\right)}}-g-g{\frac {\epsilon _{1}}{r_{1}}}$
Hence, it follows that
    $\left({\frac {\epsilon _{1}}{r_{1}}}\right)^{2}+\left({\frac {\epsilon }{r}}-{\frac {r_{1}}{r_{1}+g}}\right){\frac {\epsilon _{1}}{r_{1}}}+{\frac {g\epsilon }{r(r_{1}+g)}}=0$
or
    ${\frac {\epsilon _{1}}{r_{1}}}={\frac {1}{2}}\left[{\frac {r_{1}}{r_{1}+g}}-{\frac {\epsilon }{r}}\pm {\sqrt {\left({\frac {\epsilon }{r}}-{\frac {r_{1}}{r_{1}+g}}\right)^{2}-{\frac {4g\epsilon }{r(r_{1}+g)}}}}\right]$
Of the two roots we choose the one with the lower sign in order that $\epsilon _{1}$ equal zero with $\epsilon$. This root may be written
    ${\frac {\epsilon _{1}}{r_{1}}}={\frac {1}{2}}\left[{\frac {r_{1}}{r_{1}+g}}-{\frac {\epsilon }{r}}-\left({\frac {r_{1}}{r_{1}+g}}-{\frac {\epsilon }{r}}\right){\sqrt {1-\left({\frac {4g\epsilon }{r(r_{1}+g)}}\right)\left({\frac {r_{1}}{r_{1}+g}}-{\frac {\epsilon }{r}}\right)^{-2}}}\right]$
It is seen that the expression under the radical is finite, continuous and one-valued for values of $\epsilon$ such that
    ${\frac {\epsilon }{r}}<{\frac {r_{1}}{r_{1}+g}}$ and ${\frac {rg\epsilon }{r(r_{1}+g)}}<\left({\frac {r_{1}}{r_{1}+g}}-{\frac {\epsilon }{r}}\right)^{2}$
Article 183. Returning to the substitutions
    ${\bar {\xi }}=\epsilon \xi +\epsilon _{1}\xi _{1}\qquad {\bar {\eta }}=\epsilon \eta +\epsilon _{1}\eta _{1}$
we assume that the functions $\xi ,\xi _{1},\eta ,\eta _{1}$ become zero at the endpoints (or limits) of the curve and are so chosen that $W_{01}^{(1)}$ does not vanish within the limits of integration. We have then at once from $(A^{a})$ the power-series
    $\epsilon _{1}=-{\frac {W_{10}^{(1)}}{W_{01}^{(1)}}}\epsilon +\epsilon P(\epsilon )$
where the power-series $P(\epsilon )$ vanishes with $\epsilon$.
From this we have
    $6)\qquad {\bar {\xi }}=\epsilon \left(\xi -{\frac {W_{10}^{(1)}}{W_{01}^{(1)}}}\xi _{1}\right)+\epsilon \xi _{1}P(\epsilon )\qquad {\bar {\eta }}=\epsilon \left(\eta -{\frac {W_{10}^{(1)}}{W_{01}^{(1)}}}\eta _{1}\right)+\epsilon \eta _{1}P(\epsilon )$
If we subject the integral $I^{(0)}$ to the same variation, we have [cf. formula $(A^{a})$]
    $\Delta I^{(0)}=\epsilon W_{10}^{(0)}+\epsilon _{1}W_{01}^{(0)}+(\epsilon ,\epsilon _{1})_{2}$
and consequently
    $\Delta I^{(0)}=\epsilon \left(W_{10}^{(0)}-{\frac {W_{10}^{(1)}}{W_{01}^{(1)}}}W_{01}^{(0)}\right)+(\epsilon )_{2}$
If then, the integral $I^{(0)}$ is to have a maximum or a minimum value, it is necessary that
    $W_{10}^{(0)}-{\frac {W_{10}^{(1)}}{W_{01}^{(1)}}}W_{01}^{(0)}$
be equal to zero.
We have, therefore, the necessary condition
    ${\frac {W_{10}^{(0)}}{W_{10}^{(1)}}}={\frac {W_{01}^{(0)}}{W_{01}^{(1)}}}$
From this it is seen that the quotient ${\frac {W_{10}^{(0)}}{W_{10}^{(1)}}}$, is independent of the arbitrary functions $\xi ,\eta$, since it does not vary if we write for $\xi ,\eta$ as functions of $t$ other functions $\xi _{1},\eta _{1}$. Consequently it follows that the value of the above quotient depends only upon the nature of the curve $x=x(t),y=y(t)$.
Article 184. We might generalize the problem treated above by requiring the curve $x=x(t),y=y(t)$ which minimizes or maximizes the integral
    $I^{(0)}=\int _{t_{0}}^{t_{1}}F^{(0)}(x,y,x',y'){\text{d}}t$
while at the same time the following integrals have a prescribed value:
    $I^{(1)}=\int _{t_{0}}^{t_{1}}F^{(1)}(x,y,x',y'){\text{d}}t$
    $I^{(2)}=\int _{t_{0}}^{t_{1}}F^{(2)}(x,y,x',y'){\text{d}}t$
    $I^{(\mu )}=\int _{t_{0}}^{t_{1}}F^{(\mu )}(x,y,x',y'){\text{d}}t$
the functions $F^{(0)},F^{(1)},\ldots ,F^{(\mu )},$ being of the same nature as the function $$ defined in Chapter I.
We must now consider the deformation of the curve caused by the variations
    ${\bar {\xi }}=\epsilon \xi +\epsilon _{1}\xi _{1}+\epsilon _{2}\xi _{2}+\cdots +\epsilon _{\mu }\xi _{\mu }\qquad {\bar {\eta }}=\epsilon \eta +\epsilon _{1}\eta _{1}+\epsilon _{2}\eta _{2}+\cdots +\epsilon _{\mu }\eta _{\mu }$
We have, then, if we write $w_{i}=y'\xi _{i}-x'\eta _{i}~~(i=1,2,\ldots ,\mu )$, and suppose that the $\xi$'s and $\eta$'s vanish for $t=t_{0}$ and $t=t_{1}$
    $\Delta I^{(0)}=\epsilon \int _{t_{0}}^{t_{1}}G^{(0)}w{\text{d}}t+\epsilon _{1}\int _{t_{0}}^{t_{1}}G^{(0)}w_{1}{\text{d}}t+\cdots +\epsilon _{\mu }\int _{t_{0}}^{t_{1}}G^{(0)}w_{\mu }{\text{d}}t+(\epsilon ,\epsilon _{1},\cdots ,\epsilon _{\mu })_{2}$
    $\Delta I^{(1)}=0=\epsilon \int _{t_{0}}^{t_{1}}G^{(1)}w{\text{d}}t+\epsilon _{1}\int _{t_{0}}^{t_{1}}G^{(1)}w_{1}{\text{d}}t+\cdots +\epsilon _{\mu }\int _{t_{0}}^{t_{1}}G^{(1)}w_{\mu }{\text{d}}t+(\epsilon ,\epsilon _{1},\cdots ,\epsilon _{\mu })_{2}$
    $\Delta I^{(2)}=0=\epsilon \int _{t_{0}}^{t_{1}}G^{(2)}w{\text{d}}t+\epsilon _{1}\int _{t_{0}}^{t_{1}}G^{(2)}w_{1}{\text{d}}t+\cdots +\epsilon _{\mu }\int _{t_{0}}^{t_{1}}G^{(2)}w_{\mu }{\text{d}}t+(\epsilon ,\epsilon _{1},\cdots ,\epsilon _{\mu })_{2}$
    $\Delta I^{(\mu )}=0=\epsilon \int _{t_{0}}^{t_{1}}G^{(\mu )}w{\text{d}}t+\epsilon _{1}\int _{t_{0}}^{t_{1}}G^{(\mu )}w_{1}{\text{d}}t+\cdots +\epsilon _{\mu }\int _{t_{0}}^{t_{1}}G^{(\mu )}w_{\mu }{\text{d}}t+(\epsilon ,\epsilon _{1},\cdots ,\epsilon _{\mu })_{2}$
By means of the last $\mu$ equations, if the determinant
    $\left|\int _{t_{0}}^{t_{1}}G^{(i)}w_{j}{\text{d}}t\right|_{i,j=1,2,\cdots ,\mu }$
is different from zero, we may, for sufficiently small values of $\epsilon _{1},\epsilon _{2},\ldots ,\epsilon _{\mu }$, express these quantities as convergent power-series in $\epsilon$
These power-series when substituted in $\Delta I^{(0)}$ cause it to have the form
    $\Delta I^{(0)}=\epsilon D+(\epsilon )_{2}$
where
    $D=\left|\int _{t_{0}}^{t_{1}}G^{(i)}w_{j}{\text{d}}t\right|_{i,j=0,1,\cdots ,\mu ~{\text{and}}~w_{0}=w}$
In order that the integral $I^{(1)}$ have a maximum or a minimum value, it is therefore necessary that
    $D=$
This determinant, when expanded, may be written in the form
    $\int _{t_{0}}^{t_{1}}[\lambda _{0}G^{(0)}+\lambda _{1}G^{(1)}+\cdots +\lambda _{\mu }G^{(\mu )}]w{\text{d}}t=0$
where $\lambda _{i}$ is the first minor of $\int _{t_{0}}^{t_{1}}G^{(i)}w{\text{d}}t$ in the determinant $D$.
Hence, as before (cf. Art. 79, where we had $G=0$), we have here
    $\lambda _{0}G^{(0)}+\lambda _{1}G^{(1)}+\cdots +\lambda _{\mu }G^{(\mu )}=0$
Article 185. Similarly, if in Art. 183 we denote the quotient ${\frac {W_{10}^{(0)}}{W_{10}^{(1)}}}$ by $\lambda$ and then give to $W_{01}^{(0)}$ and $W_{01}^{(1)}$ their values, we have
    $\int _{t_{0}}^{t_{1}}(G^{(0)}-\lambda G^{(1)})w{\text{d}}t=0$
From this it follows that
    $G^{(0)}-\lambda G^{(1)}$
We may prove a very important theorem regarding the constant $\lambda$, viz: -it has one and the same value for the whole curve; i. e., we always have the same value of $\lambda$, whatever part of the curve $x=x(t),y=y(t)$ we may vary. Consider the values of $t$ laid off on a straight line, and suppose that the constant $\lambda$ has a definite value for, say, the interval $t_{2}\ldots t_{3}$ which also corresponds to a certain portion of curve. This value (see Art. 183) is independent of the manner in which the portion of curve $t_{2}\ldots t_{3}$ has been varied. Next consider an interval $t'\ldots t''$ which includes the interval $t_{2}\ldots t_{3}$; then, there belongs to all the possible variations of the interval $t'\ldots t''$, also that variation by which $t'\ldots t_{2}$ and $t_{3}\ldots t''$ remain unchanged and only $t_{2}\ldots t_{3}$, varies. As $\lambda$ has a definite value for this interval and is independent of the manner in which the curve has been varied, it must have the same value for $t'\ldots t''$.
Article 186. The differential equation $G^{(0)}-\lambda G^{(1)}=0$ is the same as the one we would have if we require that the integral
    $\int _{t_{0}}^{t_{1}}f(x,y,x',y'){\text{d}}t$
have a maximum or a minimum value, where $F$ is written for the function
    $F^{(0)}-\lambda F^{(1)}$
Through this differential equation (See Art. 90) ;$x$ and $y$ are expressible in terms of $t$ and $\lambda$ and two constants of integration $\alpha$ and $\beta$ in the form
    $x=\phi (t,\alpha ,\beta ,\lambda )\qquad y=\psi (t,\alpha ,\beta ,\lambda )$
The curve represented by these equations is a solution of the problem, when indeed a solution is possible.
Article 187. We prove next a very important theorem which often gives a criterion whether a sudden change in direction can take place or not within a stretch where the variation is unrestricted (cf. Art. 97). Suppose that on a position $t=t'$, where the variation is unrestricted, a sudden change in direction is experienced. On either side of $t'$ take two points $t_{1}$ and $t_{2}$ so near to $t'$ that within the intervals $t_{1}\ldots t'$ and $t'\ldots t_{2}$ a similar discontinuity in change of direction is not had. Among the possible variations there is one such that the whole curve remains unchanged except the interval $t_{1}\ldots t_{2}$, which is, of course, varied in such a way that the integral $I^{(1)}$ retains its value. The variation of the integral $I^{(0)}$ depends then only upon the variation of the sum of integrals
    $\int _{t_{1}}^{t'}F^{(0)}(x,y,x',y'){\text{d}}t+\int _{t_{'}}^{t_{2}}F^{(0)}(x,y,x',y'){\text{d}}t$
We cause a variation in the stretch $t_{1}\ldots t_{2}$ by writing
    ${\bar {\xi }}=\epsilon \xi +\epsilon _{1}\xi _{1}\qquad {\bar {\eta }}=\epsilon \eta +\epsilon _{1}\eta _{1}$
where we assume that
(A)
    $\xi ,\xi _{1},\eta ,\eta _{1}$ are all zero for $t=t_{1}$ and $t=t_{2}$
    $\xi ,\xi _{1},\eta _{1}$ are zero for $t=t_{1}$
    $\eta \neq 0$ for $t=t_{1}$
We may then always determine $\epsilon _{1}$ as a power-series in $\epsilon$ so that $\Delta I^{(1)}=0$.
If by $\phi$ we denote an expression of the form $\phi ^{(0)}-\lambda \phi ^{(1)}$, we have (Art. 79)
    $\Delta I^{(0)}=\epsilon \int _{t_{1}}^{t'}Gw{\text{d}}t+epsilon\int _{t'}^{t_{2}}Gw{\text{d}}t+\epsilon \left[(\xi -\lambda \xi _{1}){\frac {\partial F}{\partial x'}}+(\eta -\lambda \eta _{1}){\frac {\partial F}{\partial y'}}\right]_{t_{1}}^{t'}+\epsilon \left[(\xi -\lambda \xi _{1}){\frac {\partial F}{\partial x'}}+(\eta -\lambda \eta _{1}){\frac {\partial F}{\partial y'}}\right]_{t'}^{t_{2}}+\epsilon (\epsilon )$
If the curve $x=x(t),y=y(t)$ minimizes or maximizes the integral $I^{(0)}$, it is necessary that the coefficient of $\epsilon$ on the right-hand side of the above expression be zero. Since $G=0$ for unrestricted variation, it follows from the assumption (A) that
    $\eta _{t'}\left[\left({\frac {\partial F}{\partial y'}}\right)_{t'}^{-}-\left({\frac {\partial F}{\partial y'}}\right)_{t'}^{+}\right]=0$
If in the assumptions (A) we assume for $t=t_{1}$ that $\eta =0$ and $\xi \neq 0$, we have an analogous equation for $x'$.
It therefore follows (cf. Art. 97) that
    $\left[{\frac {\partial (F^{(0)}-\lambda F^{(1)})}{\partial x'}}\right]_{t'}^{-}=\left[{\frac {\partial (F^{(0)}-\lambda F^{(1)})}{\partial x'}}\right]_{t'}^{+}$
    $\left[{\frac {\partial (F^{(0)}-\lambda F^{(1)})}{\partial y'}}\right]_{t'}^{-}=\left[{\frac {\partial (F^{(0)}-\lambda F^{(1)})}{\partial y'}}\right]_{t'}^{+}$
We have then the theorem : Along those positions which are free to vary of the curve which satisfies the differential equation $G=0$, the quantities ${\frac {\partial F}{\partial x'}}$ and ${\frac {\partial F}{\partial y'}}$ vary everywhere in a continuous manner, even on such positions of the curve where a sudden change in its direction takes place.
Article 188. It is obvious that these discontinuities may all be avoided, if we assume that $\xi ,\eta ,\xi _{1},\eta _{1}$ vanish at such points. This we may suppose has been done. We may also impose many other restrictions upon the curve ; for example, that it is to go through certain fixed points, or that it is to contain certain given portions of curve, or that it is to pass through a certain limited region. In all these cases there are points on the curve which cannot vary in a free manner. But whatever condition may be imposed upon the curve, the following theorem is true.
All points which are free to vary and there always exist such points must satisfy the differential equation $G^{(0)}-\lambda G^{(1)}=0$, and for all such points the constant $\lambda$ has the same value.
Article 189. The second variation. We assume that the variations at the limits and at all points of the curve where there is a discontinuity in the direction, vanish. We also suppose that the variations ${\bar {\xi }},{\bar {\eta }}$ have been so chosen that $\Delta I^{(1)}=0$.
We then have (cf. Art. 115):
    $\Delta I^{(0)}=\epsilon \delta I^{(0)}+{\frac {\epsilon ^{2}}{2}}\int _{t_{0}}^{t_{1}}\left[F_{1}^{(0)}\left({\frac {{\text{d}}w}{{\text{d}}t}}\right)^{2}+F_{2}^{(0)}w^{2}\right]{\text{d}}t+(\epsilon )_{3}$
    $0=\epsilon \delta I^{(1)}+{\frac {\epsilon ^{2}}{2}}\int _{t_{0}}^{t_{1}}\left[F_{1}^{(1)}\left({\frac {{\text{d}}w}{{\text{d}}t}}\right)^{2}+F_{2}^{(1)}w^{2}\right]{\text{d}}t+(\epsilon )_{3}$
and consequently
    $\Delta I^{(0)}=\epsilon [\delta I^{(0)}-\lambda \delta I^{(1)}]+{\frac {\epsilon ^{2}}{2}}\int _{t_{0}}^{t_{1}}\left[F_{1}\left({\frac {{\text{d}}w}{{\text{d}}t}}\right)^{2}+F_{2}w^{2}\right]{\text{d}}t+(\epsilon )_{3}$
Since
    $\delta I^{(0)}-\lambda \delta I^{(1)}=\int _{t_{0}}^{t_{1}}(G^{(0)}-\lambda G^{(1)})w{\text{d}}t=0$
it follows that
    $\Delta I^{(0)}={\frac {\epsilon ^{2}}{2}}\int _{t_{0}}^{t_{1}}\left[F_{1}\left({\frac {{\text{d}}w}{{\text{d}}t}}\right)^{2}+F_{2}w^{2}\right]{\text{d}}t+(\epsilon )_{3}$
This last integral may be written at once (Art. 119) in the form
    $\Delta I^{(0)}={\frac {\epsilon ^{2}}{2}}\int _{t_{0}}^{t_{1}}F_{1}\left({\frac {{\text{d}}w}{{\text{d}}t}}-{\frac {w}{u}}{\frac {{\text{d}}u}{{\text{d}}t}}\right)^{2}{\text{d}}t$
where $u$ is determined from the differential equation (Art. 118)
    $J=F_{1}{\frac {{\text{d}}^{2}u}{{\text{d}}t^{2}}}+{\frac {{\text{d}}F_{1}}{{\text{d}}t}}{\frac {{\text{d}}u}{{\text{d}}t}}-F_{2}u=0$
It follows here as a necessary condition for the existence of a maximum or a minimum that $F_{1}$ for all portions of the curve at which there is free variation, must in the first case be everywhere negative' and in the second case everywhere positive' and must also be different from 0 and $\infty$. In order that this transformation of the integral be possible the equation $J=0$ must admit of being integrated in such a way that $u$ is different from zero on all portions of curve, which vary freely (Art. 128).
We shall determine in Chapter XVII whether the three necessary conditions thus formulated are also sufficient for a maximum or a minimum value of the integral $I^{(1)}$. By means of the example in the next Chapter, we shall also show that if there exists a curve, for which the first integral has a maximum or a minimum value while the second integral retains a given value, then the curve is determined through the three conditions, which are the same here as those formulated in Art. 135. The behavior of the ${\mathcal {E}}$-function is then decisive regarding whether there in reality exists a maximum or a minimum.
* ↑ Cf. Lectures on the Theory of Maxima and Minima, etc., p. 20.
* ↑ Cf. Lectures on the Theory of Maxima and Minima of Functions of Several Variables, p. 21.