OU blog

Personal Blogs

Valentin Fadeev

All around the Jacobi equation

Visible to anyone in the world
Edited by Valentin Fadeev, Sunday, 18 Sept 2011, 23:25

The Jacobi equation for a functional serving as a test for weak local extremum can be derived in quite different ways. The following geometrical approach was found in the book "Calculus of Variations" by L.E. Elsgoltz, 1958. It hangs upon the question of whether or not the stationary path can be included in a part of a one-parametric family of stationary paths (field of paths). There can be 2 options. Either no 2 paths of the family intersect, or all patghs of the family share 1 common point (but not more) in the given interval.

For example y equals cap c times sine of x will form a family of the first type in open delta comma a close , where delta greater than zero , a less than pi , the family of the second type in open zero comma a close , a less than pi . In the interval open zero comma a close , a greater than or equals pi no such family can be constructed.

Suppose we have a one-parametric family of stationary paths y equals y of x comma cap c . For example, we can fix one of the boundary points and use the gradient of the paths in this point as parameter C.

paths

The envelope of this family is found by eliminating C from the following system of equations:

equation sequence y equals y of x comma cap c times normal partial differential times y of x comma cap c divided by normal partial differential times cap c equals zero times open asterisk operator close

Along each path of the family normal partial differential times y of x comma cap c divided by normal partial differential times cap c is a function of only x. Denote this function as u of x for some given C. Then equation left hand side u sub x super prime equals right hand side normal partial differential squared times y of x comma cap c divided by normal partial differential times cap c times normal partial differential times x .

y equals y of x comma cap c are solutions of the E-L by assumption. Therefore:

cap f sub y of x comma y of x comma cap c comma y sub x super prime of x comma cap c minus d divided by d times x times cap f sub y super prime of x comma y of x comma cap c comma y sub x super prime of x comma cap c equals zero

Differentiating this equality by C and letting u equals normal partial differential times y of x comma cap c divided by normal partial differential times cap c we obtain:

cap f sub y times y times u plus cap f sub y times y super prime times u super prime minus d divided by d times x times open cap f sub y times y super prime times u plus cap f sub y super prime times y super prime times u super prime close equals zero

Rearranging we get:

open cap f sub y times y minus d divided by d times x times cap f sub y times y super prime close times u minus d divided by d times x times open cap f sub y super prime times y super prime times u super prime close equals zero

which is obviously a Jacobi equation.

Thus, if u has a zero somewhere in the interval, it follows from (*) above that this is a common point of the stationary path and the envelope. This is a point, conjugate to the left end of the interval.

It seemed to me at first that this proof serves only for theoretical purpose, as another way of deriving the Jacobi equation. However, the idea behind can be used to find the solution of the Jacobi equation, without actually solving the equation itself!

Consider the following example.

equation sequence cap s of y equals integral over zero under one y super prime two divided by y super four d x times y of zero equals one semicolon y of one equals one divided by two

equation left hand side y super prime times two times y super prime divided by y super four minus y super prime two divided by y super four equals right hand side cap c squared

y super prime divided by y squared equals cap c

equation left hand side negative one divided by y equals right hand side cap c times x plus cap c sub zero

(we can suppress plus minus , assuming it is absorbed by the constant)

y equals one divided by cap c sub one times x plus cap c sub two

Now apply boundary conditions:

equation sequence y of zero equals one times y of one equals one divided by two

y equals one divided by x plus one

If we only use the first condition, that is fix the left boundary, we get the one-parametric family:

y of x comma cap c equals one divided by cap c times x plus one

C being the gradient at x equals zero with a reverse sign. Now following the idea described above we can find u of x :

equation sequence u of x equals normal partial differential times y of x comma cap c divided by normal partial differential times cap c equals negative x divided by open cap c times x plus one close squared

Finally setting cap c equals one by virtue of the boundary conditions:

u of x equals negative x divided by open x plus one close squared

Now we move on to dervie the Jacobi equation through the coefficients of the second variation:

equation sequence cap p of x equals normal partial differential squared times cap f divided by normal partial differential times y super prime two equals two divided by y super four equals two left parenthesis x plus one times right parenthesis super four

equation sequence cap q of x equals normal partial differential squared times cap f divided by normal partial differential times y squared minus d divided by d times x times open normal partial differential squared times cap f divided by normal partial differential times y times normal partial differential times y super prime close equals 20 times y super prime two divided by y super six plus eight times d divided by d times x times open y super prime divided by y super five close equals 20 times y super prime two divided by y super six plus eight times y super double prime divided by y super five minus 40 times y super prime two divided by y super six equals eight times y super double prime divided by y super five minus 20 times y super prime two divided by y super six

equation sequence cap q of x equals 16 times open x plus one close super five divided by open x plus one close cubed minus 20 times open x plus one close super six divided by open x plus one close super four equals negative four times open x plus one close squared

Inserting the above results into the equation:

left parenthesis cap p times u super prime times right parenthesis super prime minus cap q times u equals zero

two times open open x plus one close super four times u super prime close super prime plus four times open x plus one close squared times u equals zero

open open x plus one close super four times u super prime close super prime plus two times open x plus one close squared times u equals zero

sum with, 3 , summands open x plus one close super four times u super double prime plus four times open x plus one close cubed times u super prime plus two left parenthesis x plus one times right parenthesis squared times u equals zero

sum with, 3 , summands u super double prime plus four divided by x plus one times u super prime plus two divided by open x plus one close squared times u equals zero

Instead of solving the equation which can be technically demanding, we shall verify if the expression for u of x found above is a solution. We do not need a general solution in this case. All non-trivial solutions of a homogeneous equation of the second order satisfying the condition u of x sub zero equals zero differ from each other only by a constant multiplier and thus have the same zeros.

u equals negative x divided by open x plus one close squared

equation sequence u super prime equals negative one divided by open x plus one close squared plus two times x divided by open x plus one close cubed equals negative x minus one plus two times x divided by open x plus one close cubed equals x minus one divided by open x plus one close cubed

equation sequence u super double prime equals one divided by open x plus one close cubed minus three times x minus one divided by open x plus one close super four equals x plus one minus three times x plus three divided by open x plus one close super four equals negative two times x minus two divided by open x plus one close super four

negative two times x plus two divided by open x plus one close super four times open x plus one close super four plus four times x minus one divided by open x plus one close cubed times open x plus one close cubed minus two times x divided by open x plus one close squared times open x plus one close squared equals zero

sum with, 3 , summands negative two times x plus four plus four times x minus four minus two times x equals zero

So we indeed have a solution.

For the integrand does not depend explicitly on x we can simplify the Jacobi equation by exchanging the roles of the variables:

equation sequence cap s of y equals integral over zero under one y super prime two divided by y super four d x equals integral over zero under one open d times y close squared divided by open d times x close squared times one divided by y super four d x equals integral over one under one divided by two one divided by x super prime times y super four d y equals negative integral over one divided by two under one one divided by x super prime times y super four d y equals cap s of x

cap g of y comma x comma x super prime equals negative one divided by x super prime times y super four

Euler-Lagrange equation has the first integral:L

cap g sub x super prime equals cap c

equation left hand side negative one divided by x super prime two times y super four equals right hand side negative c squared

x equals cap c sub one divided by y plus cap c sub two

cap c sub one plus cap c sub two equals zero

Again we leave one arbitrary constant to form a family:

x of y comma cap c equals cap c times open one divided by y minus one close

equation sequence u of y equals normal partial differential times x of y comma cap c divided by normal partial differential times cap c equals one divided by y minus one

equation sequence cap p of y equals normal partial differential squared times cap g divided by normal partial differential times x super prime two equals negative two divided by x super prime three times y super four equals negative two divided by negative one divided by y super six times y super four equals two times y squared

Since the transformed integrand does not depend explicitly on the dependent variable, Q will vanish and the Jacobi equation has the first integral:

cap p times u sub y super prime equals cap c

equation left hand side two times y squared times u sub y super prime equals right hand side negative two times cap c sub one

u of y equals cap c sub one divided by y plus cap c sub two

cap c sub one plus cap c sub two equals zero

equation left hand side cap c sub one equals right hand side negative one

u of y equals one minus one divided by y
Thus, up to the sign we get the same expression. (I am not sure where I may be loosing the sign, but it obviously has litle effect on the argument)

Permalink Add your comment
Share post
Valentin Fadeev

Facts about integrating factors

Visible to anyone in the world
Edited by Valentin Fadeev, Sunday, 18 Sept 2011, 23:26

Even if you are faced with a plain separable ODE, the process of separation of variables itself implies multiplying both parts by some factor. Thus the integrating factor seems to be one of the most devious tricks of solving equations.

There is a general path to establish its existence. It can be found in many textbooks. I am interested in some particular cases here which give beautiful solutions.

First, for a homogeneous equation it is possible to find a closed formula for the integrating factor.

It can be shown that for equation

cap m of x comma y times d times x plus cap n of x comma y times d times y equals zero ,

where M and N are homogeneous functions of their arguments integrating factor has the form:

mu equals one divided by x times cap m of x comma y plus y times cap n of x comma y

Apply this to equation:

open a squared times x minus y close times d times x plus open x plus y close times d times y equals zero

equation sequence mu equals one divided by x times open a squared times x minus y close plus y times open x plus y close equals one divided by a squared times x squared plus y squared

Multiplying both parts by this expression we obtain:

a squared times x minus y divided by a squared times x squared plus y squared times d times x plus x plus y divided by a squared times x squared plus y squared times d times y equals zero

Rearranging:

x times d times y minus y times d times x divided by a squared times x squared plus y squared plus a squared times x times d times x plus y times d times y divided by a squared times x squared plus y squared equals zero

d of x divided by y divided by a squared postfix plus left parenthesis x divided by y times right parenthesis squared plus one divided by two times d times open a squared times x squared plus y squared close divided by a squared times x squared plus y squared equals zero

And the result becomes obvious.

For the next example it is useful to note the fact that if mu is an integrating factor for equation cap m times d times x plus cap n times d times y equals zero giving solution in the form cap u of x comma y equals cap c , then equation left hand side mu sub one equals right hand side mu times phi of cap u where phi is any differentiable function shall also be an integrating factor. Indeed

equation sequence mu sub one times open cap m times d times x plus cap n times d times y close equals phi of cap u times mu times open cap m times d times x plus cap n times d times y close equals phi of cap u times d times cap u

giving the differential for the function normal cap phi of cap u equals integral phi of cap u d cap u

This leads to the following practical trick of finding the factor. All terms of the equations are split in two groups for each of which it is easy to find the integrating factor. Then each factor is written in the most general form involving an arbitrary function as described above. Finally we try to find such functions that make both factors identical.

Consider the following equation:

open x squared times y squared minus one close times d times y plus two times x times y cubed times d times x equals zero

Rearranging the terms:

open x squared times y squared times d times y plus two times x times y cubed times d times x close minus d times y equals zero

For the second term now the integrating factor is trivial, it is 1. Hence the most general form will look like mu sub one equals phi of y .

For the first part it is easy to see that the factor should be one divided by x squared times y cubed giving solution x squared times y equals cap c , hence equation left hand side mu sub two equals right hand side one divided by x squared times y cubed times psi of x squared times y .

To make the two identical we want mu sub two to be independent of x. Setting psi of t equals t gives equation left hand side mu sub two equals right hand side one divided by y squared .

Applying this one we get:

open x squared times d times y plus two times x times y times d times x close minus d times y divided by y squared equals zero

x squared times y plus one divided by y equals cap c

Both methods were discovered in the classic book "A course on differential equations" by V.V. Stepanov

Permalink 2 comments (latest comment by Valentin Fadeev, Thursday, 15 Apr 2010, 16:02)
Share post
Valentin Fadeev

Big O-Oh

Visible to anyone in the world
Edited by Valentin Fadeev, Sunday, 16 Jan 2011, 23:17

There are many good articles on O symbols, some more technical, some more popular. But at the end of the day you often find yourself staring at the excercise remembering all those definitions and still not knowing what to do next. This has been the case for me, until i took some freedom to play around with this device. I think i finally got it down. Here are two examples:

x1/2/(1+O(x))=x1/2(1+O(x))=x1/2+x1/2O(x)=O(x1/2)+O(x3/2)

=O(x1/2)+O(x1/2)=2O(x1/2)=O(x1/2)

Where in the first equality i used geometric expansion until the linear term and the fact that for non-zero k kO(x)=O(x), hence -O(x)=O(x). Of course "=" sign must be understood as element of through all manipulations with O.

A more demanding one:

equation sequence two times phi minus sine of two times phi divided by two times sine squared of phi equals sum with, 3 , summands two times phi minus two times phi plus eight divided by six times phi cubed plus cap o of phi super five divided by two left parenthesis phi plus cap o of phi cubed times right parenthesis squared equals four divided by three times phi cubed plus cap o of phi super five divided by two times open sum with, 3 , summands phi squared plus cap o of phi super four plus cap o of phi super six close equals

equation sequence equals four divided by three times phi cubed plus cap o of phi super five divided by two times open phi squared plus cap o of phi super four close equals two divided by three times phi times one plus cap o of phi squared divided by one plus cap o of phi squared right parenthesis equals two divided by three times phi times open one plus cap o of phi squared close times open one plus cap o of phi squared close equals

equation sequence equals two divided by three times phi times open sum with, 3 , summands one plus cap o of phi squared plus cap o of phi super four close equals two divided by three times phi times open one plus cap o of phi squared close equals two divided by three times phi plus cap o of phi cubed

Of course, it would be a bad idea in the 5th transition to cancel out the numerator and denominator simply by division. The two O-s in this case may stand for different classes of functions.

Permalink
Share post
Valentin Fadeev

Surface of revolution

Visible to anyone in the world
Edited by Valentin Fadeev, Sunday, 16 Jan 2011, 23:17

This is a simple geometric derivation of the formula for the area of the surface of revolution:

cone

If the surface area of a cone is expressed as:

equation sequence cap s sub c equals integral over zero under two times pi one divided by two times cap l times cap r d phi equals pi times cap r times cap l

Then the area of the frustum is the difference between the whole cone and the cut-away part:

equation sequence cap s sub f equals cap s minus s equals pi times cap r times cap l minus pi times r times open cap l minus l close equals pi times cap r times cap l minus pi times cap r times cap l minus l divided by cap l times open cap l minus l close equals pi times cap r times open cap l minus left parenthesis cap l minus l times right parenthesis squared divided by cap l close equals pi times cap r times l times open one plus cap l minus l divided by cap l close equals pi times open cap r plus r close times l

So, in the infinitesimal case:

equation sequence delta times cap s equals pi times open sum with, 3 , summands y plus y plus delta times y close times delta times l equals two times pi times y times delta times l plus o of delta times y equals two times pi times y times Square root of one plus y super prime two times delta times x plus o of delta times x

Permalink
Share post
Valentin Fadeev

On variable changes in definite integrals

Visible to anyone in the world
Edited by Valentin Fadeev, Sunday, 16 Jan 2011, 23:18

Some integrals yield only one type of substitution that really brings them into a convenient form. Any other method would make them more complicated. However, in some cases totally different methods can be applied with equal effect. In case of definite integrals it is of course not necessary to come back to original variable which makes things even easier. Here is one example

integral over zero under one x cubed times d times x divided by Square root of one minus x squared

The most natural way is to apply a trigonometric substitute. We will not consider this method here. Instead an algebraic trick can be employed:

equation sequence integral over zero under one x cubed times d times x divided by Square root of one minus x squared equals integral over zero under one x squared times x times d times x divided by Square root of one minus x squared equals one divided by two times integral over zero under one x squared times d of x squared divided by Square root of one minus x squared equals one divided by two times integral over zero under one t times d times t divided by Square root of one minus t equals

equation sequence negative one divided by two times integral over zero under one open one minus t minus one close times d times t divided by Square root of one minus t equals one divided by two times integral over zero under one d times t divided by Square root of one minus t minus one divided by two times integral over zero under one Square root of one minus t d t equals

negative Square root of one minus t times vertical line sub zero super one plus one divided by three left parenthesis equation sequence one minus t times right parenthesis super three solidus two times vertical line sub zero super one equals one minus one divided by three equals two divided by three

Alternatively we can use integration by parts:

equation sequence integral over zero under one x cubed times d times x divided by Square root of one minus x squared equals integral over zero under one x squared times x times d times x divided by Square root of one minus x squared equals

equation sequence equals integral over zero under one x squared times d of negative Square root of one minus x squared equals negative x squared times Square root of one minus x squared times vertical line sub zero super one plus integral over zero under one Square root of one minus x squared times two times x d x equals

equation left hand side equals right hand side negative two divided by three left parenthesis equation left hand side one minus x squared times right parenthesis super three solidus two times vertical line sub zero super one equals right hand side two divided by three

Or apply an even more exotic treatment:

let x equals one divided by t

equation left hand side integral over zero under one x cubed times d times x divided by Square root of one minus x squared equals right hand side integral over one under normal infinity d times t divided by t super four times Square root of t squared minus one

let t equals hyperbolic cosine of z

equation sequence integral over one under normal infinity d times t divided by t super four times Square root of t squared minus one equals integral over hyperbolic cosine super negative one of one under normal infinity hyperbolic sine of z times d times z divided by left parenthesis hyperbolic cosine of z times right parenthesis super four times hyperbolic sine of z equals integral over hyperbolic cosine super negative one of one under normal infinity d times z divided by left parenthesis hyperbolic cosine of z times right parenthesis super four equals

equation sequence equals integral over hyperbolic cosine super negative one of one under normal infinity left parenthesis left parenthesis hyperbolic cosine of z times right parenthesis squared minus open hyperbolic sine of z times right parenthesis squared close times d times z divided by left parenthesis hyperbolic cosine of z times right parenthesis super four equals

equation sequence equals integral over hyperbolic cosine super negative one of one under normal infinity d times z divided by left parenthesis hyperbolic cosine of z times right parenthesis squared minus integral over hyperbolic cosine super negative one of one under normal infinity left parenthesis hyperbolic tangent of z times right parenthesis squared times d times z divided by left parenthesis hyperbolic cosine of z times right parenthesis squared equals

MathJax failure: TeX parse error: Extra open brace or missing close brace

For

lim over z right arrow normal infinity of hyperbolic tangent of z equals one

and

equation sequence hyperbolic tangent of hyperbolic cosine super negative one of one equals Square root of one minus one divided by left parenthesis hyperbolic cosine of open hyperbolic cosine super negative one of one close times right parenthesis squared equals zero ,

the same result is obtained

Permalink
Share post
Valentin Fadeev

Summing multiplier

Visible to anyone in the world
Edited by Valentin Fadeev, Sunday, 16 Jan 2011, 23:19

This example illustrates the application of the method of "summing multiplier" to solving certain types in recurrencies. It can be found for instance in the book "Concrete Mathematics" by Graham, Knuth and Patashnik. In its essence it translates the idea of the integrating multiplier from the theory of ODEs.

 

Consider the recurrency:

equation left hand side cap i sub n equals right hand side negative n divided by a times cap i sub n minus one plus one divided by a times x super n times e super a times x

Rewrite it in the form

equation left hand side a times cap i sub n equals right hand side negative n times cap i sub n minus one plus x super n times e super a times x

then multiply both sides by s sub n which is to be determined:

equation left hand side s sub n times a times cap i sub n equals right hand side negative s sub n times n times cap i sub n minus one plus s sub n times x super n times e super a times x

The trick will be done, if we find s sub n satisfying the relation:

equation left hand side negative s sub n times n equals right hand side s sub n minus one times a

Solving for s sub n and expanding recursively we get:

equation left hand side s sub n equals right hand side left parenthesis negative one times right parenthesis super n minus one times a super n divided by n factorial

substituiting into equation:

left parenthesis negative one times right parenthesis super n minus one times a super n plus one divided by n factorial times cap i sub n postfix minus left parenthesis equation left hand side negative one times right parenthesis super n minus two times a super n divided by open n minus one close factorial times cap i sub n minus one equals right hand side left parenthesis negative one times right parenthesis super n minus one times a super n divided by n factorial times x super n times e super a times x

cap i sub zero is easily obtained from the original integral:

equation left hand side cap i sub zero equals right hand side e super a times x minus one divided by a

Summing from 1 to n we get the so called "telescopic sum" on the left side meaning that only the first and the last terms survive:

left parenthesis equation left hand side negative one times right parenthesis super n minus one times a super n plus one divided by n factorial times cap i sub n minus open negative a close times cap i sub zero equals right hand side e super a times x times n ary summation from k equals one to n over left parenthesis negative one times right parenthesis super k minus one times a super k divided by k factorial times x super k

Ultimately solve for cap i sub n :

equation left hand side cap i sub n equals right hand side minus left parenthesis negative one times right parenthesis super n times n factorial divided by a super n plus one left parenthesis equation left hand side one minus e super a times x minus e super a times x times n ary summation from k equals one to n over open negative one times right parenthesis super k times a super k divided by k factorial times x super k close equals right hand side minus left parenthesis negative one times right parenthesis super n times n factorial divided by a super n plus one left parenthesis one minus n ary summation from k equals zero to n over open negative one times right parenthesis super k times a super k divided by k factorial times x super k close equals n factorial e super a times x times n ary summation from k equals zero to n over left parenthesis negative one times right parenthesis super n minus k divided by a super n minus k plus one times k factorial minus left parenthesis negative one times right parenthesis super n times n factorial divided by a super n plus one

(Note that the second term on the right side is the solution for "homogeneous" variant of the equation, i.e. withot x super n times e super a times x term, suggesting another method borrowed from ODEs)

Permalink
Share post
Valentin Fadeev

Generalized homogeneous equations

Visible to anyone in the world
Edited by Valentin Fadeev, Sunday, 16 Jan 2011, 23:20

One of the rarely used methods of solving ODEs applies to the so-called generalized homogeneous equations. The word "generalized" means that the terms are not homogeneous in the classic sense, if all variables are assigned the same dimension. But they may be made homogeneous in a wider sense by choosing the appropriate dimension for the dependent variable. Here is one example.

equation left hand side x cubed times d times y divided by d times x equals right hand side y times open x squared plus y close

If we assign dimension 1 to x and dx and dimension m to y and dy, then the left side has dimension 3+m-1=m+2 on the right side we have m+2 and 2m. To balance things let m+2=2m, hence m=2 and we get a "generalized homogeneous equation" of the 4th order. The trick is to let:

x equals e super t comma y equals u times e super m times t

which in this case gives:

equation left hand side d times x equals right hand side e super t times d times t comma equation left hand side d times y equals right hand side e super two times t times open u super prime plus two times u close

equation left hand side d times y divided by d times x equals right hand side e super t times open u super prime plus two times u close

Hence the equation becomes:

equation left hand side e super four times t times open u super prime plus two times u close equals right hand side u times e super four times t times open one plus u close

equation left hand side u super prime plus u equals right hand side u squared

u super prime divided by u squared plus one divided by u equals one

letting z=1/y

equation left hand side d times z divided by d times x equals right hand side z minus one

z equals one plus c times x

y equals x squared divided by one plus c times x

This method can of course, be applied to higher order equations

Permalink
Share post

This blog might contain posts that are only visible to logged-in users, or where only logged-in users can comment. If you have an account on the system, please log in for full access.

Total visits to this blog: 77437