VARIATIONS, Calculus of. 1. The object of this calculus is the discovery of the form of a function which shall fulfil certain conditions, not expressed in direct terms, but involved in a finite integral. The first instance of a problem of this kind appears in a scholium to Prop. 34, b. ii. of the Principia, published in 1687. (See Prob. 2. below.) Newton gave the correct result; but he supplied no demonstration, nor do we find any method of solution applicable to such cases until ten years afterwards. In the Acta Eruditorum for 1696, John Bernoulli enunciated the problem of the brachystochrone, (Prob. 1. below), and invited mathematicians to give a solution. After a considerable interval, his brother, James Bernoulli, gave in his result, which was, that the curve required is the cycloid. As the mode of demonstration which he adopted is very nearly that which was employed by all the earlier writers on the subject. It will not be uninteresting to exhibit it here.

2. The problem is the following: “To find that curve down which a body falling by the force of gravity will move from A to B in the shortest possible time.”

The principles employed in the solution are these two: 1st, that when a quantity is a maximum or minimum, a slight change in the variables will produce no variation in the value of the function; 2d, that what is true of the whole quantity is likewise true of every portion of it.

The second principle requires that the time down PQR should be less than the time down any other line, as PSR; [7:21:558]and by applying the first principle to this property, we deduce the fundamental equation, time down PQR= time down PSR ultimately.

Hence,

*****PQ , QR PS SR , . , velocity at Q v ’. at R v ’. at Q - *^ v ’. at R ' t t ' ,nate lI ’

PQ—PS SR—QR Qe Sc , .

√NS - √OR ’ ° r √NS - √OR

, . cos. ^≤rSQe cos. ^c≤QSc , . , that '* ' √NS = - √OR ≠ ,wβ "⅛"

sin. which curve makes with vertical at Q √NS ° r do. do. R - √OR ’ a well-known property of the cycloid.

3. Of the two principles employed in this solution, the first is equivalent to the theorem demonstrated in Fluxions, art. 62, and is the basis of all solutions. The second, although in the case in question it is actually true, is an assumption not warranted by the nature of the problem. It does not appear that any geometrical solution is exempt from this objection, and we cannot in consequence date the existence of the Calculus of Variations earlier than 1744, when Euler published his treatise, entitled Methodus inveniendi Lineas curvas proprietate maximi minimive gaudentes. This work was followed by a memoir from the pen of Lagrange, published in the second volume of the Miscellanea Taurinensia. By the introduction of the symbol δ to express that change which is termed a variation, Lagrange may be said to have perfected the calculus; for although certain extensions were afterwards made by himself and others, no change was afterwards introduced into the general process. We shall therefore conclude our brief sketch at this point, and refer the reader for further information to the following works: the Acta Eruditorum for 1696 and the following years; the collected works of James and John Bernoulli; Memoires de l’Acad. des Sci. 1706, 1718, &c.; Brook Taylor’s Methodus Incrementorum; the Petersburg Commentaries, vols. vi. and viii.; which, together with Euler’s tract, quoted above, and his last memoir in the New Petersburg Com. vol. x. contain Euler’s writings on the subject: the Turin Miscellanies, vols. ii. and iv. and the Theorie des Fonctions Analytiques, contain Lagrange’s perfecting of the calculus. The reader will find in Woodhouse’s Isoperimetrical Problems a complete history of the calculus, as well as an admirable digest of the different methods employed.

4. We proceed to the investigation of the theorem which is the basis of the calculus of variations.

To find the change or variation of a formula comprehended under the integral sign, when the variables x and y, on which it depends, receive the increments δ x and δ y.

Let u=∫ydx be the formula under consideration; in which y involves **** x, y, ~, &c.

**********dx dr dl∕ f Fq

Denote by p, -⅛ 3 by q, &c.; the partial differential coefficients ⅛ by M, ⅛ by N, ⅛ by P, &c.; the increment dx dy j dp of fydx by Dfydx, &c.

Then Du=fy 4.⅛γ) (dx 4- dDx)—-fydx, —fDydx +fydix, —J'Dydx -∣- yDx—fDxdy.

Now fdxDy=fdx(JΛDx 4. Nδy 4- P⅛> 4-...), and J'Dxdy-∣'Dx(lΛd.τ 4- Sdy 4. Pψ4- ·.),

.∙. δa=γβz4-∕[N(<⅛8y— Dxdy) 4- P(dxDp—Dxdp)+ ... j, =V 8x +f{N(Dy- p&r)4- P(⅞d— qDx) 4-. ] dx.

Also if δy— pDx be denoted by ω, we have ; 7 d(x + 6x) dx dx ’ dDu dDx dp , ^- c dχ- p -dx~di sx ' dhy d d da> = -cbt-Tx^=d- χ ^ y - p ^=dχ-- hence putting p for y and .∙. q for p, we obtain

&C. — &C.

By substitution,

βu=γftr4- f*{ N“4-P^~ ~^^^dx u ^*^'" } d jc ' I‘ d^ i ω dω dQ l’ d l Q J Q d^ ώ +∕ “dx i dx ' &C. = &C.

.∙. 8 li =yoi4-ω(P-^ 4- ^-&C.)

dω <7R

+ ώ (0_ ^+··· )

4- &c. &c.

The last line of this very elegant expression is due to Euler; the rest of the formula to Lagrange.

Let x o x 1 , y 0 y 1 , be the values of x and y at the two limits; then the whole value of δ u between the limits is, ******δu= y ,ta 1 —γ 0 8x,4-ω 1 (P 1 — &c.)—α><,(P o — & c.)

β r*' dP d 1 Q „ v ,

4- &c. + y u ,(N- — 4. — _ &c.)dx

5. To satisfy the conditions of a maximum or minimum value of u, the value of δ u, expressed only as far as to the first powers of δ x 0 , δ y 0 , &c., must be equal to zero.

But this quantity consists of two parts; the one an integrated expression, depending only on the values of δ x 0 , δ y 0 , &c.; the other an unintegrated expression, depending on the general values of δ x, δ y. Now it is evidently possible to make the latter expression assume an infinity of different values, corresponding to an assigned value of the former. It is therefore impossible to render the whole expression equal to zero for all values of ω, without making the two parts separately zero.

Our conditions are therefore,

***** κτ rfP d 2 Q o N -S + Sf- fe =»· < 1 >

γ 1 δx,—γΛr<,-∣- ω 1 (P l — &c.) — ω 0 (P<,— &C.) =0. (2.)

6. It may be remarked, that the object of our investigation is to discover a mode of satisfying certain conditions, by means of establishing a relation between y and x. To effect this, we seek to separate these quantities from the quantities δ x and δ y. This will explain why the unintegrated part of the expression has been reduced to the form in which we left it. Had ω been implicated with x and y, [7:21:559] a relation would have been assigned, not between x and y, but between x, y, and ω. We may add, that if ∫ydx be an integrable expression, the equation (1) is an identical equation.

7. In the above investigation it has been assumed that the other quantities which enter into the function γ, together with x and y, are absolute constants. If this hypothesis be not correct, our formulae require to be modified.

Suppose a, b, &c., to enter the function γ.

Let δγ=Mδ x+ Nδγ+Pδ p +&c.+Aδ a +&c..., then the addition to be made to the expression for δ u is ∫(Aδa+... ) dx. But δ a &c . depend not on the general variation, but only on the variation at the limits, (see Prob. 4); hence they are independent of δ x and δ y , and consequently form part of equation (2), which thus becomes *****r x '

γ∣8r 1 -γ o δx<,+ω 1 (P∣-&c.)-ω 0 (P.-&C.) +*α Ad⅛-)-fec.≡O.

(3.)

8. Had our object been merely to determine the form of a function which should satisfy certain conditions, we might have rested content with causing one of the quantities, as x only, to vary. By this means equation (1) would have been obtained. There are however many problems which depend on the relation between δ x and δ y, at least at the limits, (see Prob. 3.) At the commencement of the variation, for instance, we cannot remove to any point we please, from the circumstance that we are compelled to begin in a certain manner. Hence equation (2) is a most important part of the result; and by assigning a relation between the things given and sought at the extremities, shews us the position of the extremities of the curve, and thereby restricts the solution from a kind of curve (as a cycloid) to a specific curve, (a cycloid of known dimensions and position).

9. From equation (1) we might easily deduce a number of expressions, for particular cases, easy of application. Thus, if γ contains only y and p, the formula (1) gives ***** -s=°∙

Hence dy=Ndy±Pdp= dy-}-Pdp=dP .p-}-Pdp ax -d{Pp)∖ and γ=Py>-pc. (4.)

But as the reader will readily supply such formulae for himself, we proceed to exhibit the solution of a few problems.

**************Prob. 1. To find the curve of quickest descent from one given point to another.

Let x and y be the co-ordinates of a point in the curve, y being measured vertically downwards from the upper point.

Then fc= -? — ∕*A +p dx . ; n w ] 1 ; c ) 1 ⅛ e quantity √2√ √; corresponding to y contains only y and p.

Ifγ= —, P≡; ÷ andformula(4)gives √,y √ 1 +p 2 ∙√.J'

¼∕ 1 + p t p" l ds _ 1

√y ~ Ji^+p 2 -Jy +c ' ° r dx ~ c Jy which is the differential equation to a cycloid. We may solve the problem by taking x as the vertical co-ordinate, in which case y will contain only p, and equation (1) will dV give - r _0, or P=c.

α⅛

But P, in this case, is ”, "' r"- . a — _ Aj∖+p τ ., y ∕x dg C/ J x the same result as before, with the exception that x and y are interchanged. The solution of this equation shews that the cusp of the cycloid is at the highest point, and its axis vertical.

Prob. 2. Required the curve which by revolution round its axis generates the solid on which the resistance produced by motion in a fluid shall be less than on any other curve having the same extreme co-ordinates.

The resistance varies ********** dx∖ ∙'. by formula (4)

W ’ _(W W ∖ ,

1+P* - V+P ’ (>+∕∕∕ p orθ= 2 ⅛ + c.

(1+P 4 ) ,ψ

This equation expresses the same property as that which Newton, without demonstration, assigned to the curve.

Prob. 3. To find the shortest distance between two given curves.

The expression for y is here √1+ p 2 where, as the line on which the distance is measured is a right line, p=c. But by formula (2)

γ 1 δx 1 —γ 0 δ x 0 +P 1 ω 1 —P 0 ω 0 =0, when δ x 1 , &c . are indefinitely small. Also, if the curves have no connexion, we may make δ x 1 δ y 1 assume what values soever we please, without affecting δ x 0 , δ y 0 ; we must therefore have separately

********γ 1 ⅛r l 4-P 1 ω∣=0, and γ, ∂r 0 +P o ω o =0.

The latter gives √l +∕√ δx 0 + - ∙ p ° (δy, t -p,,δx ο )=0, *l 1 +p: or δ∙r i 4- p ο δy c ∙=0, or ~ = ι = —, when δx a and

° ai u Pa c δy 0 are indefinitely diminished.

But, since the measurement must begin in the first curve, the limit of expresses the tangent of the angle in

***O j '.

which the tangent to the first curve cuts the axis. Our equation, then, shews that the line of shortest distance cuts the first curve at right angles. For the same reason it cuts the second curve also at right angles. Its position is consequently determined.

Prob. 4. To find the curve of quickest descent from one given curve to another.

Let h be the vertical ordinate of the point at which motion begins:

****1 ∕ ’ √ 1 -PP 2 then tχx- i =z! — ' - dx, from which (by Prob. 1.) V 2gJ t ∕y-A it is evident that the curve is a cycloid. To determine the angles at which the cycloid cuts the given curves, we have the following equation, ***(3): γ∣ &r,—y β δx 0 4-P 1 ω 1 —P,<0,4- ∙√⅛ dx=0. The problem admits of several cases.

1. If the body starts from a given point in the first curve, that curve performs no part in the problem. Here *** δx,, ω a , and hh, are all zero.

∙ ’ ∙ γ 1 Aτι + Pι(ty,-p,δar,)=O.

But γ 1 —P 1 p 1 =c, (see Prob. 1.)

1 , r , Pi

.∙. , -=c, and P 1 = —- -5=π- —; — ., = r>∣c,

1 +Pι ’ —b ∙J 1 +∕> 1 2 V y 1 —h .∙. cix l -∣- p l coy l = 0; or ⅛ for the second curve, at the [7:21:560]point in which the cycloid cuts it, is equal to **— which shews that the curve and cycloid cut each other at right angles.

2. If the motion is supposed to commence at the same horizontal line, whatever he the point at which the body reaches the first curve; we have h=d constant, and the equation gives

***γ 1 Ax 1 —γ o Ar β -∣-P∣ω∣—P o ω o =0, from which it is evident, as in Problem 3, that the cycloid cuts both curves at right angles.

<Iy

3. If motion begins at the first curve, A=y rt and = -⅞=---i(--i=÷

hence our equation becomes

7∣δ* 1 —γ 0 Ja∙,+ P 1 ω l —P.ω 0 — ~ dx=A, *. or γ 1 ⅛r 1 + P 1 ω,=0, y,ta,+ P^> 0 4-iy 0 (P 1 -P<,)=0.

From the first of these it follows that the lower curve is cut at right angles. The second gives

(γ.-P.P.)*e.+P 1 8y,=0, or c + p ι^°= 0 ∙ Also the first equation is c4-P 1 ^ l =O.

.∙.; from which it appears that the tangents to

0Λ' 0 o*c∣

the two curves at the two points of section are parallel to one another.

10. We have hitherto solved only problems of absolute maxima and minima. But a very simple consideration will enable us to apply the same formulae to the investigation of relative maxima and minima. The problem of isoperimetricals (prob. 5 below) will illustrate this class of questions. Here the integrated function is not required to be a maximum or a minimum absolutely, but only one consistent with the further condition that another integrated function shall not change its value. Although, therefore, δ ∫ydx is actually zero, yet δ x and δ y are not, as in our previous investigations, any quantities whatever, but only such as will consist with the required condition that δ ∫ √1 +p 2 dx shall also be zero. To Euler we owe the idea of substituting for the quantity u , which is made to vary, not fydx, as in other cases, but ∫ydx+∫vdx∙, ∫νdx being the quantity which is to remain unchanged. By this substitution we determine the relation between x and y , which renders γ+av a maximum or a minimum. But as that relation involves the arbitrary constant a, we restrict it by assigning to that quantity such a value as shall render ∫vdx of the required magnitude. We have therefore found a relation, which not only makes the sum of two quantities a maximum or a minimum, but which likewise reduces one of those quantities to a specified value. This relation then evidently makes the other quantity a maximum or a minimum. Thus the problem is solved in all its generality.

Prob. 5. AB is a given line, PMQ a line perpendicular to AB at the point M. The points P and Q are supposed to generate two curves by the motion of the line perpendicular to AB, having a relation such, that although both are unknown, PM is a known function of QM. It is required to determine both curves, when the length of AQB is given, and the area APBMA is a maximum.

Let QM=y, then PM= f (y ), where f is a known fonction.

Also length ABQ= ∫ √1+ p 2 dx , and area APB= ∫ f(y ) dx , the limits of both integrals being the same.

*****∙∙∙ y=∕(y)+≡^Γ+p i >

, on*

and by formula (4)∕(y) -∣- a -J 1 +p ’ = ~f==i+ c < .J 1 -pp θ r i c -∕ Ο)} √ 1 +p'=o, the differential equation to the curve AQB.

Ify=φ(ar) be the solution of this equation, j=z∕'[ψ('r); is the equation to the curve APB.

Cor. lff(y)=y, we get α*=(c-y)*4-(c ʹ —a:)*, the equation to a circle.

Prob. 6. To find the curve which, with a given length, contains between its chord and its arc the greatest possible area.

By cor. to last problem the equation is α ’ =(c— y) ’

Let the origin be at A and AB=2r; then since y=0, both when x=S) and or=2r, we get

ÿ*— 2cy~2rx—x'.

And from the equation to the limits (γ 1 —P l p,)Ax 1 =.0, or γ 1 —P 1 p∣=0 .∙. c—0 ; and the curve required is the semicircle.

Prob. 7. Given the length of the curve contained between two points in a horizontal line; required its nature, that the centre of gravity of the arc may be the lowest possible. Here we have to make fx>J 1 +p i dx a maximum, whi1st ∕√l +p ’ dx is constant.

.∙. γ=ar√ 1 -∣-p ’-∣-αV 1 -∣∙p 1 , and^*=0, P=c, QX or x*+⅛ r .

√ 1 +p ’ the integral of which is y~c log. (ar-}-α-∣-v∖.r-∣-αy—c*) +c, the equation to the catenary.

The solution of more complicated problems is not consistent with our limits. For the investigation of formulae in cases where γ involves an integral, or where it is given by a differential equation, the reader is referred to Woodhouse’s Treatise. Considerations on the mode of distinguishing a maximum from a minimum will be found in Lagrange’s Theorie des Fonctions Analytiques. Variation of the compass. See Magnetism. Variation , in Music, means either the change induced upon a melody by embellishments given by the vocal or instrumental performer, or the new forms and more studied ornaments given by a composer to a theme of his own or to some popular air. Among instrumental variations, some of the most beautiful are to be found in Hadyn’s andante movements, and in Mozart and Beethoven’s variations to their own andantes, or to popular themes on which they have displayed their skill and genius.