Main Concept - Maple Programming Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Math Apps : Calculus : Multivariate : MathApps/Brachistochrone

Brachistochrone

Main Concept

Consider the following problem, posed by Johann Bernoulli in 1696 to the most renowned mathematicians of his time:

 

Given two points A and B, of which A is the higher, what curve minimizes the time for a bead starting from rest at A to descend to B under the force of gravity alone?

 

Bernoulli had already discovered the answer (a cycloid curve), but wished to test his contemporaries. He received correct solutions from Jacob Bernoulli, Gottfried Leibniz, Guillaume de l'Hôpital, and Isaac Newton, with each arriving at that answer by vastly different methods. Today, the problem is regarded as an excellent example of the calculus of variations, and is presented as follows.

 

Variational Solution

Let A be the origin and B some lower point not directly below A. Let x,y be the horizontal and vertical distances of the bead from A, so that y>0 as the bead falls (g>0 as well). The problem stated is to minimize the time t that it takes a bead to go from A to B under gravity alone. The travel time can be expressed as the integral of the inverse of the speed v with respect to arclength s along the bead's path:

 

t=AB dsvs,

 

To express the speed in terms of arclength, notice that the potential energy of the bead is Vy = mgy where the convention V0=0 is used,  and the kinetic energy is T=12mv2 the bead starts from rest, E0 = 0 as well. Using conservation of energy E= T + V,you get:

 

 0=12mv2 - mgy.

 

Rearranging terms gives vy = 2 g y. Changing variables using:

 

ds2dx2 = 1 + dy2dx2,

 

where ds is the increment of length on the curve, the integral for minimization becomes:

t=12g0xB 1+y'2y dx,

where y' = dydx.  It remains to find the function yx which minimizes the time.

 

This problem can be approached using the calculus of variations. In particular, the function F = 1+y'2 yis varied until an extreme value for the integral is found. In most cases the Euler-Lagrange equation would be applied next:

 

FyⅆⅆxFy'=0.

 

However, F does not explicitly depend on x, and thus the simpler Beltrami identity is used instead:

 

Fy'Fy'=C,

where C is an integration constant. Applying F, you get:

 

 1+y'2yy'2 y 1+y'2 =C,

Squaring both sides yields:

 

y 1+y'2 =C2.

 

This differential equation is separable and can be integrated to give:

 

x=2ayy+a π2 arcsin1ya,

 

where a=12C2, using y0=0 and x0=0 as initial conditions. Using the following definition:

θ=π2arcsin1ya,

 

the solution can be rewritten parametrically as:

xθ = a θsin θ, yθ = a 1 cos θ.

 

This is a parameterization of a cycloid! To fully solve the problem, a and θB must be determined by satisfying the boundary conditions

xθB =xB , yθB =yB ,

 

which in general must be solved numerically.

The actual time can be then expressed in terms of θ as follows:

t

=

12g0θB 1ydx2dθ2+dy2dθ2 dθ

 

=

a2 g0θB 1 cos θ2+sin2θ1 cos θ dθ

 

=

ag θ.

Derivation of the Beltrami Identity

You can simplify the Euler-Lagrange equation for the case that F does not explicitly depend on x, that is, Fx=0. Consider the total derivative of F with respect to x:

ⅆFⅆx=Fyy'+Fy'y''+Fx.

 

The last term must be zero since F does not explicitly depend on x by assumption. Solving for Fyy' gives:

Fyy' =ⅆFⅆxFy'y''.

 

Multiplying the Euler-Lagrange equation by y' and substituting the previous equation into it gives:


ⅆFⅆxFy'y''y'ⅆⅆxFy'=0,

 

where left hand side can be written as:

ⅆⅆxFy'Fy'=0.

 

Integrating this gives the Beltrami identity:

Fy'Fy'=C,

where C is the integration constant.

 

Click or drag the endpoint and press "Play" to release the beads.

More MathApps

MathApps/Calculus