First variation. (one-variable problem) January 14, 2013

Size: px
Start display at page:

Download "First variation. (one-variable problem) January 14, 2013"

Transcription

1 First vrition (one-vrible problem) Jnury 14, 2013 Contents 1 Sttionrity of n integrl functionl Euler eqution (Optimlity conditions) First integrls: Three specil cses Vritionl problem s limit of finite-dimensionl problem Sttionrity of boundry terms Vrition of boundry conditions Broken extreml nd the Weierstrss-Erdmn condition Functionl dependent on higher derivtives 12 4 Non-fixed intervl Trnsverslity condition Extreml broken t n unknown point Severl minimizers Euler equtions nd first integrls Vritionl boundry conditions

2 Since, however, the rules for isoperimetric curves or, in modern terms, extreml problems were not sufficiently generl, the fmous Euler undertook the tsk of reducing ll such investigtions to generl method which he gve in the work Essy on new method of determining the mxim nd minim of indefinite integrl formuls ; n originl work in which the profound science of the clculus shines through. Even so, while the method is ingenious nd rich, one must dmit tht it is not s simple s one might hope in work of pure nlysis. In Essy on new method of determining the mxim nd minim of indefinite integrl formuls, by Lgrnge, Sttionrity of n integrl functionl The technique ws developed by Euler, who lso introduced the nme Clculus of vritions in The method is bsed on n nlysis of infinitesiml vritions of minimizing curve. The min scheme of the vritionl method is s follows: Assume tht the optiml curve u(x) exists mong smooth (twice-differentible curves). Compre the optiml curve with close-by trjectories u(x) + δu(x), where δu(x) is smll in some sense. Using the smllness of δu, we simplify the comprison, deriving necessry conditions for the optiml trjectory u(x) Vritionl methods yield to only necessry conditions of optimlity becuse it is ssumed tht the compred trjectories re close to ech other; on the other hnd, they re pplicble to gret vriety of extreml problems clled vritionl problems. 1.1 Euler eqution (Optimlity conditions) Consider the problem clled the simplest problem of the clculus of vritions min u I(u), I(u) = F (x, u, u )dx, u() = u, u(b) = u b, (1) Here integrnt F clled the Lgrngin, is twice differentible function of its three rguments, I(u) is clled the cost functionl. It is not known priori whether the minimizer u 0 (x) is smooth, but let us ssume tht it is twice differentible function of x. For exmple, consider the re of the surfce of revolution. According to the clculus, the re J of the surfce is A(r) = π r(x) 1 + r (x) 2 dx, where r(x) is the vrible distnce from the xes OX of rottion. The problem of miniml re of such surfce I = min r(x) A(u), r() = R, r(b) = R b 2

3 is vritionl problem. To void trivil solution, vlues r() nd r(b) re fixed. To derive necessry condition of optimlity of minimizer u 0 we use the ides of clculus, computing n nlog of the derivtive of I with respect to u (clled the functionl derivtive) nd setting it to zero. We suppose tht function u 0 = u 0 (x) is minimizer nd replce u 0 with test function u 0 + δu, ssuming tht the norm δu of the vrition δu is infinitesiml. The test function u 0 +δu stisfies the sme boundry conditions s u 0. If indeed u 0 is minimizer, the increment of the cost δi(u 0 ) = I(u 0 + δu) I(u 0 ) is nonnegtive: δi(u 0 ) = 1 0 (F (x, u 0 + δu, (u 0 + δu) ) F (x, u 0, u 0))dx 0. (2) If δu is not specified, the eqution (2) is not too informtive. However, prticulr form of the vrition δu simplifies (2) nd llows for finding n eqution for the minimizer. Clculus of vritions suggests set of tests tht differ by vrious ssumed form of vritions δu. Euler Lgrnge Equtions The simplest vritionl condition (the Euler Lgrnge eqution) is derived ssuming tht the vrition δu is infinitesimlly smll nd loclized: δu = { ρ(x) if x [x0, x 0 + ε], 0 if x is outside of [x 0, x 0 + ε]. (3) Here ρ(x) is continuous function tht vnishes t points x 0 nd x 0 + ε nd is constrined s follows: ρ(x) < ε, ρ (x) < ε x. (4) The integrnd t the perturbed trjectory cn be expnded into Tylor series, F (x, u 0 + δu, (u 0 + δu) ) = F (x, u 0, u 0) + (x, u 0, u 0) δu + (x, u 0, u 0) δu + o(δu, δu ) Here, δu is derivtive of the vrition δu, δu = (δu), o(δu, δu ) denotes higher order terms which norms re smller thn δu nd δu when ε 0. Substituting this expression into (2) nd collecting liner (with respect to ε) terms, we rewrite (2) s δi(u 0 ) = ( ) (δu) + (δu) dx + o(ε) 0. (5) where The F is clculted t the exmined trjectory u 0. To simplify nottions, we omit index ( 0 ) below. 3

4 The vritions δu nd (δu) re mutully dependent nd (δu) cn be expressed in terms of δu. Integrtion by prts of the underlined term in (5) gives nd we obtin (δu) dx = 0 δi(u 0 ) = ( d dx ) δu dx + S F (x, u, u )δu dx + δu where S F denotes the functionl derivtive, δu x=b x= x=b x= + o(ε), (6) S F (x, u, u ) = d dx +. (7) The nonintegrl term in the right-hnd side of (6) is zero, becuse the boundry vlues of u re prescribed s u() = u nd u(b) = u b ; therefore their vritions δu x= nd δu x=b equl zero, δu x= = 0, δu x=b = 0 Due to the rbitrriness of δu, we rrive t the following Theorem 1.1 (Sttionrity condition) Any differentible nd bounded minimizer u 0 of the vritionl problem (1) is solution to the boundry vlue problem S F (x, u, u ) = d dx = 0 x (, b); u() = u, u(b) = u b, (8) clled the Euler Lgrnge eqution. The Euler Lgrnge eqution is lso clled the sttionry condition of optimlity becuse it expresses sttionrity of the vrition. Remrk 1.1 The sttionrity test lone does not llow to conclude whether u is true minimizer or even to conclude tht solution to (8) exists. For exmple, the function u tht mximizes I(u) stisfies the sme Euler Lgrnge eqution. The tests tht distinguish miniml trjectory from other sttionry trjectories re discussed in Chpter??. In this derivtion, it is ssumed tht the extreml u(t) is twice differentible function of x. Using the chin rule, the left-hnd side of eqution (8) cn be rewritten s S F (x, u, u ) = 2 F 2 u + 2 F u + 2 F x (9) 4

5 Exmple 1.1 Compute the Euler eqution for the problem I = min u(x) 1 We compute = u, The minimizer u 0 (x) is 0 F (x, u, u )dx u(0) = 1, u(1) =, F = 1 2 (u ) u2 = u nd the Euler eqution becomes u u = 0 in (0, 1), u(0) = 1, u(1) =. u 0 (x) = cosh(x) + cosh(1) sinh(1) sinh(x). The definition of the wek solution nturlly rises from the vritionl formultion tht does not check the behvior of the minimizer in ech point but in ech infinitesiml intervl. The minimizer cn chnge its vlues t severl points, or even t set of zero mesure without effecting the objective functionl. In mbiguous cses, one should specify in wht sense (Riemnn, Lebesgue) the integrl is defined nd chnge the definition of vrition ccordingly. 1.2 First integrls: Three specil cses In severl cses, the Euler eqution (8) cn be integrted t lest once. These re the cses when Lgrngin F (x, u, u ) does not depend on one of the rguments. Below, we investigte these cses. Assume tht F = F (x, u), nd the mini- Lgrngin is independent of u miztion problem is I(u) = min u(x) 1 0 F (x, u)dx (10) In this cse, the vrition does not involve integrtion by prts, nd the minimizer does not need to be continuous. Euler eqution (8) becomes n lgebric reltion for u = 0 (11) Curve u(x) is determined in ech point independently of neighboring points. The boundry conditions in (8) re stisfied by jumps of the extreml u(x) in the end points; these conditions do not ffect the objective functionl t ll. Exmple 1.2 Consider the problem I(u) = min u(x) 1 0 (u sin x) 2 dx, u(0) = 1; u(1) = 0. The miniml vlue J(u 0 ) = 0 corresponds to the discontinuous minimizer sin x if 0 x 1 u 0 (x) = 1 if x = 0 0 if x = 1 5

6 Formlly, the discontinuous minimizer contrdicts the ssumption posed when the Euler eqution ws derived. To be consistent, we need to repet the derivtion of the necessry condition for the problem (10) without ny ssumption on the continuity of the minimizer. This derivtion is quite obvious. Lgrngin is independent of u If Lgrngin does not depend on u, F = F (x, u ), Euler eqution (8) cn be integrted once: = constnt (12) The first order differentil eqution (12) for u is the first integrl of the problem; it defines quntity tht stys constnt everywhere long the optiml trjectory. To find the optiml trjectory, it remins to integrte the first order eqution (12) nd determine the constnts of integrtion from the boundry conditions. Exmple 1.3 Consider the problem I(u) = min u(x) 1 0 (u cos x) 2 dx, u(0) = 1; u(1) = 0. The first integrl is = u (x) cos x = C Integrting, we find the minimizer, u(x) = sin x + Cx + C 1. The constnts C nd C 1 re found from nd the boundry conditions: C 1 = 1, C = 1 sin 1, minimizer u 0 nd the cost of the problem become, respectively u 0 (x) = sin x (sin 1 + 1)x + 1 I(u 0 ) = (sin 1 + 1) 2. Notice tht the Lgrngin in the exmple (1.2) is the squre of difference between the minimizer u nd function sin x, nd the Lgrngin in the exmple (1.3) is the squre of difference of their derivtives. In the problem (1.2), the minimizer coincides with sin x, nd jumps to stisfy the prescribed boundry vlues. The minimizer u in the exmple (1.3) cnnot jump. Consider continuous pproximtion of derivtive u of discontinuous function; it increses in the proximity of the point of discontinuity nd is unbounded, such growth would increse the objective functionl, nd therefore it would be nonoptiml. We del with such problems below in Chpter??. 6

7 Lgrngin is independent of x If F = F (u, u ), eqution (8) hs the first integrl: W (u, u ) = constnt (13) where W (u, u ) = u F Indeed, compute the x-derivtive of W (u, u ) which must be equl to zero by virtue of (13): d dx W (u, u ) = [ ( 2 u + F u u + 2 F u 2 )] u u = 0 where the expression in squre brckets is the derivtive of the first term of W (u, u ). Cncelling the equl terms, we bring this eqution to the form ( u 2 F 2 u + 2 F u ) = 0 (14) The expression in prenthesis coincide with the left-hnd-side term S(x, u, u ) of the Euler eqution in the form (9), simplified for the considered cse (F is independent of x, F = F (u, u )). W is constnt t ny solution u(x) of Euler eqution. Insted of solving the Euler eqution, we my solve the first-order eqution W = 0 obtining the sme solution. Exmple 1.4 Consider the Lgrngin F = 1 2 [ (u ) 2 ω 2 u 2] The Euler eqution is The first integrl is u + ω 2 u = 0 W = 1 2 ( ω 2 u 2 + (u ) 2) = C 2 = constnt Let us check the constncy of the first integrl. The solution u of the Euler eqution is equl u = A cos(ωx) + B sin(ωx) where A nd B re constnts. Substituting the solution into the expression for the first integrl, we compute W = (u ) 2 + ω 2 u 2 = [ Aω sin(cx) + Bω cos(ωx)] 2 +ω 2 [A cos(ωx) + B sin(ωx)] 2 = ω 2 (A 2 + B 2 ) We hve shown tht W is constnt t the optiml trjectory. In mechnicl ppliction, W is the totl energy of the oscilltor. 7

8 1.3 Vritionl problem s limit of finite-dimensionl problem Here, we derive Euler eqution for finite-dimensionl problem tht pproximte the simplest vritionl problem min I(u), I(u) = u(x) F (x, u, u )dx Consider clss of piece-wise constnt discontinuous functions U N : ū(x) U N, if ū(x) = u i x [ + in ] (b ) A function ū in U N is defined by n N-dimensionl vector {u 1,... u N }. Rewriting the vritionl problem for this clss of minimizers, we replce the derivtive u (x) with finite difference Diff (u i ) Diff (u i ) = 1 (u i u i 1 ), = b N ; (15) when N, this opertor tends to the derivtive. The vritionl problem is replced by finite-dimensionl optimiztion problem: min u 1,...,u N 1 I N N I N = F i (u i, z i ), z i = Diff (z i ) = 1 (z i z i 1 ) (16) i=1 Compute the sttionry conditions for the minimum of I N (u) I N i = 0, i = 1...., N. Only two terms, F i nd F i+1, in the bove sum depend on u i : the first depends on u i directly nd lso through the opertor z i = Diff (u i ), nd the second only through z i+1 = Diff (u i+1 ): df i = i + i 1 du i i z i, df i+1 = i+1 1 du i z i+1. df k = 0, k i, k i + 1 du i Collecting the terms, we write the sttionry condition with respect to u i : I N = i + 1 ( i i i z ) i+1 = 0 z or, reclling the definition (15) of Diff -opertor, the form I N = ( ) i i+1 Diff = 0. i i z 8

9 The initil nd the finl point u 0 nd u N enter the difference scheme only once, therefore the optimlity conditions re different. They re, respectively, N+1 Diff (u N+1 ) = 0; o Diff (u 0 ) = 0. Formlly pssing to the limit N, Diff d dx, z u replcing the index ( i ) with continuous vrible x, vector of vlues {u k } of the piece-wise constnt function with the continuous function u(x), difference opertor Diff with the derivtive d dx ; then nd N F i (u i, Diff u i ) i=1 i i Diff ( i+1 z ) F (x, u, u )dx. d d x The conditions for the end points become the nturl vritionl conditions: (0) = 0, (T ) = 0, Remrk 1.2 So fr, we followed the forml scheme of necessry conditions, thereby tcitly ssuming tht ll derivtives of the Lgrngin exist, the increment of the functionl is correctly represented by the first term of its power expnsion, nd the limit of the sequence of finite-dimensionl problems exist nd does not depend on the prtition {x 1,... x N } if only x k x k 1 0 for ll k. We lso indirectly ssume tht the Euler eqution hs t lest one solution consistent with boundry conditions. If ll the mde ssumptions re correct, we obtin curve tht might be minimizer becuse it cnnot be disproved by the sttionry test. In other terms, we find tht is there is no other close-by clssicl curve tht correspond to smller vlue of the functionl. This sttement bout the optimlity seems to be rther wek but this is exctly wht the clculus of vrition cn give us. On the other hnd, the vritionl conditions re universl nd, being ppropritely used nd supplemented by other conditions, led to very detiled description of the extreml s we show lter in the course. Remrk 1.3 In the bove procedure, we ssume tht the limits of the components of the vector {u k } represent vlues of smooth function in the close-by points x 1,..., x N. At the other hnd, u k re solutions of optimiztion problems with the coefficients tht slowly vry with the number k. We need to nswer the question whether the solution of minimiztion problem tends to is differentible function of x; tht is whether the limit u k u k 1 lim k x k x k 1 exists nd this is not lwys the cse. We ddress this question lter in Chpter?? 9

10 2 Sttionrity of boundry terms 2.1 Vrition of boundry conditions Vritionl conditions nd nturl conditions The vlue of minimizer my not be specified on one or both ends of the intervl [, b]. In this cse, these vlues re clculted by the minimiztion of the gol functionl together with the minimizer. Consider vritionl problem where the boundry vlue t the right end b of the intervl is not defined nd the functionl directly depends on this vlue, min I(u), I(u) = u(x):u()=u F (x, u, u )dx + f(u(b)) (17) The Euler eqution for the problem remin the sme, S(x, u, u ) = 0, but this time it must be supplemented by vritionl boundry condition tht comes from the requirement of the sttionrity of the minimizer with respect to vrition of the boundry term. This term is ( + f ) δu(b) x=b The first term comes from the integrtion by prt in the derivtion of Euler eqution, see (6), nd the second is the vrition of the lst term in the objective functionl (17). Becuse the sign of the vrition δu(b) is rbitrry, the sttionrity condition hs the form + f x=b = 0 (18) x=b This equlity provides the missing boundry condition t the endpoint x = b for the second order Euler eqution. Similr condition cn be derived for the point x = if the vlue t this point is not prescribed. Exmple 2.1 Minimize the functionl 1 1 I(u) = min u 0 2 (u ) 2 dx + Au(1), u(0) = 0 Here, we wnt to minimize the endpoint vlue nd we do not wnt the trjectory be too steep. The Euler eqution u = 0 must be integrted with boundry conditions u(0) = 0 nd (see (18)) u (1) + A = 0 The extreml is stright line, u = Ax. The cost of the problem is I = 1 2 A2. If f = 0, the condition (18) becomes = 0 (19) x=b nd it is clled the nturl boundry condition. Exmple 2.2 Consider the problem F = (x)(u ) 2 +φ(x, u) where (x) 0 The nturl boundry condition is u x=b = 0. 10

11 2.2 Broken extreml nd the Weierstrss-Erdmn condition The clssicl derivtion of the Euler eqution requires the existence of ll second prtils of F, nd the solution u of the second-order differentil eqution is required to be twice-differentible. In some problems, F is only piece-wise twice differentible; in this cse, the extreml consists of severl curves solutions of the Euler eqution tht re computed t the intervls of smoothness. We consider the question: How to join these pieces together? The first continuity condition is continuity of the (differentible) minimizer u(x) [u] + = 0 long the optiml trjectory u(x) (20) Here [z] + = z + z denotes the jump of the vrible z. The extreml u is differentible, the first derivtive u exists t ll points of the trjectory. This derivtive does not need to be continuous. Insted, Euler eqution requests the differentibility of to ensure the existence of the term d dx in the Euler eqution. Integrting the sttionrity condition (8), we obtin sttionrity in the integrl form or If x S F (x, u, u )dx = x (x, u, u ) x = ( d (x, u, u ) dx (x, u, ) u ) dx = 0 (x, u, u ) x 0 dx + (x, u, u ) (21) x= is bounded t the optiml trjectory, the right-hnd side is continuous function of x, nd so is the left-hnd side. This requirement of continuity of n optiml trjectory is clled the Weierstrss-Erdmn condition on broken extreml. Theorem 2.1 At ny point of the optiml trjectory, the Weierstrss-Erdmn condition must be stisfied: [ ] + = 0 long the optiml trjectory u(x). (22) Exmple 2.3 (Broken extreml) Consider the Lgrngin F = 1 2 c(x)(u ) { c1 if x [, x 2 u2, c(x) = ), c 2 if x (x, b] where x is point in (, b). The Euler eqution is held everywhere in (, b) except of the point x, d dx [c 1u ] u = 0 if x [, x ) d dx [c 2u ] u = 0 if x (x, b], 11

12 At x = x, the continuity conditions hold, u(x 0) = u(x + 0), c 1 u (x 0) = c 2 u (x + 0). The derivtive u (x) itself is discontinuous; its jump is determined by the jump in the coefficients: u (x + 0) u (x 0) = c 1 c 2 These conditions together with the Euler eqution nd boundry conditions determine the optiml trjectory. 3 Functionl dependent on higher derivtives Consider more generl type vritionl problem with the Lgrngin tht depends on the minimizer nd its first nd second derivtive, J = F (x, u, u, u )dx The Euler eqution is derived similrly to the simplest cse: The vrition of the gol functionl is ( δj = δu + δu + ) δu dx Integrting by prts the second term nd twice the third term, we obtin ( δj = d ) dx + d2 dx 2 δu dx [ + δu + δu d ] x=b δu (23) dx The sttionrity condition becomes the fourth-order differentil eqution d 2 dx 2 d dx x= + = 0 (24) supplemented by two nturl boundry conditions on ech end, δu [ = 0, δu d ] dx = 0 t x = nd x = b (25) or by the correspondent min conditions posed on the minimizer u nd its derivtive u t the end points. Exmple 3.1 The equilibrium of n elstic bending bem correspond to the solution of the vritionl problem L min w(x) 0 ( 1 2 (E(x)w ) 2 q(x)w)dx (26) 12

13 where w(x) is the deflection of the point x of the bem, E(x) is the elstic stiffness of the mteril tht cn vry with x, q(x) is the lod tht bends the bem. Any of the following kinemtic boundry conditions cn be considered t ech end of the bem. (1) A clmped end: w() = 0, w () = 0 (2) simply supported end w() = 0. (3) free end (no kinemtic conditions). Let us find eqution for equilibrium nd the missing boundry conditions in the second nd third cse. The Euler eqution (24) becomes The equtions (25) become (Ew ) q = 0 (, b) δu (Eu ) = 0, δu ((Ew ) ) = 0 In the cse (2) (simply supported end), the complementry vritionl boundry condition is Eu = 0, it expresses vnishing of the bending momentum t the simply supported end. In the cse (3), the vritionl conditions re Eu = 0 nd (Ew ) = 0; the lst expresses vnishing of the bending force t the free end (the bending momentum vnishes here s well). Generliztion The Lgrngin F (x, u, u,..., u (n)) dependent on first k derivtives of udependent on higher derivtives of u is considered similrly. The sttionry condition is the 2k-order differentil eqution d dx dk ( 1)k dx k (k) = 0 supplemented t ech end x = nd x = b of the trjectory by k boundry conditions [ ] δu (k 1) (k) x=,b = 0 [ d ] δu (k 2) (k 1) dx (k) x=,b = 0... [ d dx d(k 1) ( 1)k dx (k 1) ] δu (k) x=,b = 0 If u is vector minimizer, u cn be replced by vector but the structure of the necessry conditions sty the sme. 13

14 4 Non-fixed intervl 4.1 Trnsverslity condition Free boundry Consider now the cse when the intervl [, b] is not fixed, but the end point is to be chosen so tht it minimizes the functionl. Let us compute the difference between the two functionls over two different intervls = δi = +δx (F (x, u + δu, u + δu ) F (x, u, u )) dx + The second integrl is estimted s +δx b F (x, u + δu, u + δu )dx +δx b F (x, u, u )dx F (x, u + δu, u + δu )dx F (x, u + δu, u + δu )dx = F (x, u, u ) x=b δx + o( δu, δx ) nd the first integrl is computed s before with integrtion by prts: S F (x, u, u )δu dx + δu(b) = 0 x=b 1. Suppose tht no boundry conditions re imposed t the minimizer t the point x = b. Becuse of rbitrriness of δx nd δu, we rrive t the conditions: S F (x, u, u ) = 0 x (, b), = 0, x=b nd F (x, u, u ) x=b = 0. (27) Euler eqution for the extreml stisfies n extr boundry condition (27), but hs lso n dditionl degree of freedom: unknown coordinte b. Exmple 4.1 Consider the problem s ( ) 1 min u(x),s 2 u 2 u + x dx u(0) = 0. 0 The Euler eqution u + 1 = 0 nd the condition t u(0) = 0 corresponds to the extreml u = 1 2 x2 + Ax, u = x + A where A is prmeter. The condition = u = 0 t the unknown right end x = s gives s = A. The trnsverslity condition F = 0 or ( u + x) x=a=s = 1 2 s2 s 2 + s = s (1 12 ) s = 0 We find s = 2, u = 1 2 x2 + 2x. 14

15 2. Next, consider the problem in which the boundry dt t b is prescribed, u = β, but the vlue of b is not known. In the perturbed trjectory, the boundry condition is u(b + δx) = β. The vlue of u(b + δx) is n extrpoltion of u(x) s follows u(b + δx) = u(b) + u (b)δx + o( δu, δx ) Therefore, the vlue (u + δu) x=b depends on δx, u(b) = β u (b)δx or δu(b) = u (b)δx. Combining the depending on δx terms, we obtin the condition ( F (b)) x=b u δx Becuse δx is rbitrry, the boundry conditions re: u = β nd ( F (x, u, u ) u ) x=b = 0. (28) Notice tht the condition (28) t the unknown end is identicl to the first integrl (13) of the problem in the cse when F (u, u ) is independent of x. This integrl is constnt long the trjectory. therefore condition (28) cnnot be stisfied t n isolted point. 3. Finlly, consider the problem when the rjectory ends t curve. If the boundry vlue depends on b, u(b) = φ(b), then δu = φ δx u δx. The sttionrity conditions become ( u(b) = φ(b), F (u φ ) ) x=b = 0. (29) The next exmple dels with constrint t the unknown length of the intervl nd the boundry dt. Exmple 4.2 Find the shortest pth between the origin nd curve φ(x). The pth length is given by s I = min 1 + y 2 dx, u(0) = 0 y(x),s 0 At the end point x the pth meets the curve, therefore y(s) = φ(s) or The Euler eqution δy = φ (s)δs (30) y = y = C 1 + y 2 shows tht y = constnt, therefore the pth is stright line, y = Ax s expected. At the point s, the vrition is ( u ) y F δx + y y δy = 1 δx + y δu 1 + y y 2 The sttionrity gives the reltion. δx + y δu = 0. Compring it with the constrint (30), we conclude tht y (s)φ (s) = 1, or tht the shortest pth is line orthogonl to the curve φ(x), s it is expected. 15

16 4.2 Extreml broken t n unknown point Combining the techniques, we my ddress the problem of en extreml broken in n unknown point. The position of this point is determined from the minimiztion requirement. Assume tht Lgrngin hs the form { F (x, u, u F (x, u, u ) = ) if x (, ξ) F + (x, u, u ) if x (ξ, b) where ξ is n unknown point in the intervl (, b) of the integrtion. The Euler eqution is { SF (u) if x (, ξ) S F (u) = S F+ (u) if x (ξ, b) The sttionrity conditions t the unknown point ξ consist of sttionrity of the trjectory + = 0 (31) nd sttionrity of the position of the trnsit point F + (u) u + + = F (u) u. (32) They re derived by the sme procedure s the conditions t the end point. The vrition δx of the trnsit point δx = δx + = δx increses the first prt of the trjectory nd decreses the second prt, or vise vers, which explins the structure of the sttionrity conditions. In prticulr, if the Lgrngin is independent of x, the condition (32) expresses the constncy of the first integrl (13) t the point ξ. Exmple 4.3 Consider the problem with Lgrngin { F (x, u, u + u ) = 2 + b + u 2 if x (, ξ) u 2 if x (ξ, b) nd boundry conditions The Euler eqution is u() = 0, u(b) = 1 S F (u) = { + u b + u = 0 if x (, ξ) u = 0 if x (ξ, b) The solution to this eqution tht stisfies the boundry conditions is ( ) u + (x) = C 1 sinh b+ + (x ) if x (, ξ) u (x) = C 2 (x b) + 1 if x (ξ, b) ; it depends on three constnts ξ, C 1, nd C 2 (Notice tht the coefficient does not enter the Euler equtions). These constnts re determined from three conditions 16

17 t the unknown point ξ which express (1) continuity of the extreml (2) Weierstrss-Erdmn condition (3) trnsverslity condition u + (ξ) = u (ξ), + u +(ξ) = u (ξ), + (u +(ξ)) 2 + b + u(ξ) 2 = (u (ξ)) 2. The trnsverslity condition sttes the equlity of two first integrl. It is simplified to C 2 1b + = C 2 2 From the Weierstrss-Erdmn condition, we find C 1 + b + cosh q = C 2, where q = b + + (ξ ) The first condition nd the definition of q llows for determintion of ξ: cosh q = +, ξ = + + b + cosh 1 + Finlly, we define constnts C 1 nd C 2 from the continuity C 1 sinh q = 1 + C 2 (ξ b) nd trnsverslity conditions: C 1 = sinh q b + (ξ b), C b+ 2 = sinh q b + (ξ b), 5 Severl minimizers 5.1 Euler equtions nd first integrls The Euler eqution cn be nturlly generlized to the problem with the vectorvlued minimizer I(u) = min u F (x, u, u )dx, (33) where x is point in the intervl [, b] nd u = (u 1 (x),..., u n (x)) is vector function. We suppose tht F is twice differentible function of its rguments. 17

18 Let us compute the vrition δi(u) equl to I(u + δu) I(u), ssuming tht the vrition of the extreml nd its derivtive is smll nd loclized. To compute the Lgrngin t the perturbed trjectory u+δu, we use the expnsion F (x, u + δu, u + δu ) = F (x, u, u ) + n i=1 i δu i + n δu i i=1 i We cn perform n independent vritions of ech component of vector u pplying vritions δ i u = (0,..., δu i..., 0). The increment of the objective functionl should be zero for ech of these vrition, otherwise the functionl cn be decresed by one of them. The sttionrity condition for ny of considered vritions coincides with the one-minimizer cse. δ i I(u) = ( δu i + δu i i i ) dx 0 i = 1,..., n. Proceeding s before, we obtin the system of n second-order differentil equtions, d = 0, i = 1,... n (34) dx i nd the boundry term i n x=b δu i = 0 (35) i=1 i If the vlue of u i () or u i (b) is not prescribed, the nturl boundry conditions x= or i x=b, respectively, must be stisfied. i The vector form of the system (34), S F (u) = d dx x=b = 0, δut = 0 (36) is identicl to the sclr Euler eqution. This system corresponds to n definition of differentition with respect to vector rgument u. x= Exmple 5.1 Consider the problem with the integrnd x= F = 1 2 u u 2 2 u 1 u u2 1 (37) The system of sttionrity conditions is computed to be d dx d dx = u 1 + u 2 u 1 = 0 2 = (u 2 u 1 ) = 0. If consists of two differentil equtions of second order for two unknowns u 1 (x) nd u 2 (x). 18

19 First integrls The first integrls tht re estblished for the specil cses of the sclr Euler eqution, cn lso be derived for the vector eqution. 1. If F is independent of u k, then one of the Euler equtions degenertes into lgebric reltion: = 0 k nd the one of differentil eqution in (34) becomes n lgebric one. The vrible u k (x) cn be discontinuous function of x in n optiml solution. Since the Lgrngin is independent of u k, the discontinuities of u k(x) my occur long the optiml trjectory. 2. If F is independent of u k, the first integrl exists: k = constnt For instnce, the second eqution in Exmple 5.1 cn be integrted nd replced by u 2 u 1 = constnt 3. If F is independent of x, F = F (u, u ) then first integrl exist Here W = u T T u = F = constnt (38) n i=1 u i i For the Exmple 5.1, this first integrl is computed to be ( 1 W = u u 2 (u 2 u 1) 2 u u 2 2 u 1 u ) 2 u2 1 = 1 ( u u 2 2 u 2 1) = constnt These three cses do not exhust ll possible first integrls for vector cse. For exmple, if the functionl depends only on, sy (u 1 + u 2 ), one cn hope to find new invrints by chnging the vribles. We discuss this mtter below in Sections?? nd??. Trnsverslity nd Weierstrss-Erdmn conditions These conditions re quite nlogous to the sclr cse nd their derivtion is strightforwrd. We simply list them here. The expressions, i = 1..., n remin continuous t every point of n i optiml trjectory, including the points where u i is discontinuous. If the end point of the trjectory is unknown, the condition t the end point is stisfied. u T F = 0 19

20 5.2 Vritionl boundry conditions Consider the vrition of the boundry term (35) which we rewrite here for convenience δu x=b 1 δu n = 0 (39) n x= If ll vritions δu i (), δu i (b) re free, it produces 2n boundry conditions i = 0, x = nd x = b for Euler equtions (34). In the opposite cse, when the vlues of ll minimizers re prescribed t the end points, u i () = u i, u i (b) = u b i, i = 1,..., n then the eqution (39) is stisfied, becuse ll vritions re zero. δu i () = 0, δu i (b) = 0, i = 1,..., n If the vlues of severl components of u() or u(b) re not given, the vritions of these components re free nd the corresponding nturl boundry condition supplements the boundry conditions: For ech i = 1,..., n one of the two conditions holds Either = 0 or δu i x=,b = 0. (40) x=,b i The totl number of the conditions t ech endpoint is n. The missing min boundry conditions re supplemented by the nturl conditions tht express the optimlity of the trjectory. This number grees with the number of boundry conditions needed to solve the boundry vlue problem for Euler eqution. In generl cse, p reltions (p < 2n) between boundry vlues of u re prescribed, β k (u 1 (),..., u n (), u 1 (b),..., u n (b), ) = 0 (41) t the end points x = nd x = b. In this cse, we need to find 2n p supplementry vritionl constrints t these points tht together with (41) give 2n boundry conditions for the Euler eqution (35) of the order 2n. The conditions (41) re stisfied t ll perturbed trjectories, β k (w + δw) = 0 where 2n dimensionl vector w is the direct sum of u() nd u(b) defined s: w k = u k () if k = 1,..., n w k = u k n (b) if k = n + 1,..., 2n. 20

21 The vrition δw i is constrint by liner system T β k δw = 0, w k = 1,..., p which shows tht the constrints (41) re stisfied t the vried trjectory. A mtrix form of these conditions is P δw = 0, where β 1 β n P = β p β 1... p n nd we my lso ssume tht these conditions re linerly independent. Then the solution to the system is vector δw d of the form δw d = Qv where v is n rbitrry (2n p)-dimensionl rbitrry vector nd (2n p) n mtrix Q is supplement (orthogonl mtrix) to P tht is defined s solution to the mtrix eqution P Q = 0 which leds to P δw d = P Qδv = 0 for ny v. Any dmissible vrition δw d mkes the fist vrition (36) of the objective functionl vnish; correspondingly, we hve ( ) T δw d = 0 Using the representtion of δw d nd the rbitrriness of the potentils v, we conclude tht the fist vrition vnishes is the coefficient by ech of these potentils is zero or ( ) T Q = 0 (42) This representtion provides the 2n p missing boundry conditions. Exmple 5.2 Consider gin the vritionl problem with the Lgrngin (37) ssuming tht the following boundry conditions re prescribed u 1 () = 1, β(u 1 (b), u 2 (b)) = u 2 1(b) + u 2 2(b) = 1 Find the complementry vritionl boundry conditions. At the point x =, the vrition δu 1 is zero, nd δu 2 is rbitrry. The vritionl condition is = u 2() u 1 () = 0 x= 2 Since the conditionu 1 () = 1 is prescribed, it becomes u 2() = 1 21

22 At the point x = b, the vritions δu 1 nd δu 2 re connected by the reltion β 1 δu 1 + β 2 δu 2 = 2u 1 δu 1 + 2u 2 δu 2 = 0 which implies the representtion (δu = Qδv) δu 1 = u 2 δv, δu 2 = u 1 δv where δv is n rbitrry sclr. The vritionl condition t x = b becomes ( u 2 + ) 1 u 1 δv = ( u 1u 2 + (u 2 u 1 )u 1 ) x=b δv = 0 δv 2 x=b or u 1u 2 + u 1 u 2 u 2 1 x=b = 0. We end up with four boundry conditions: u 1 () = 1, u 2() = 1, u 2 1(b) + u 2 2(b) = 1, u 1 (b)u 2(b) u 1 (b) u 2 (b) u 1 (b) 2 = 0. The second rw conditions re vritionl, they re obtined from the minimiztion requirement. Periodic boundry conditions Consider vritionl problem with periodic boundry conditions u() = u(b). The vritionl boundry conditions re = x= x=b They re obtined from the expression (39) nd equlities δu() = δu(b). 22

Chapter 4. Additional Variational Concepts

Chapter 4. Additional Variational Concepts Chpter 4 Additionl Vritionl Concepts 137 In the previous chpter we considered clculus o vrition problems which hd ixed boundry conditions. Tht is, in one dimension the end point conditions were speciied.

More information

Calculus of Variations

Calculus of Variations Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function

More information

4. Calculus of Variations

4. Calculus of Variations 4. Clculus of Vritions Introduction - Typicl Problems The clculus of vritions generlises the theory of mxim nd minim. Exmple (): Shortest distnce between two points. On given surfce (e.g. plne), nd the

More information

1.9 C 2 inner variations

1.9 C 2 inner variations 46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for

More information

Calculus of Variations: The Direct Approach

Calculus of Variations: The Direct Approach Clculus of Vritions: The Direct Approch Lecture by Andrejs Treibergs, Notes by Bryn Wilson June 7, 2010 The originl lecture slides re vilble online t: http://www.mth.uth.edu/~treiberg/directmethodslides.pdf

More information

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Partial Derivatives. Limits. For a single variable function f (x), the limit lim Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles

More information

The Basic Functional 2 1

The Basic Functional 2 1 2 The Bsic Functionl 2 1 Chpter 2: THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 2.1 Introduction..................... 2 3 2.2 The First Vrition.................. 2 3 2.3 The Euler Eqution..................

More information

Lecture 1. Functional series. Pointwise and uniform convergence.

Lecture 1. Functional series. Pointwise and uniform convergence. 1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is

More information

Problem set 1: Solutions Math 207B, Winter 2016

Problem set 1: Solutions Math 207B, Winter 2016 Problem set 1: Solutions Mth 27B, Winter 216 1. Define f : R 2 R by f(,) = nd f(x,y) = xy3 x 2 +y 6 if (x,y) (,). ()Show tht thedirectionl derivtives of f t (,)exist inevery direction. Wht is its Gâteux

More information

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007 A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus

More information

The Regulated and Riemann Integrals

The Regulated and Riemann Integrals Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004 Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when

More information

1.3 The Lemma of DuBois-Reymond

1.3 The Lemma of DuBois-Reymond 28 CHAPTER 1. INDIRECT METHODS 1.3 The Lemm of DuBois-Reymond We needed extr regulrity to integrte by prts nd obtin the Euler- Lgrnge eqution. The following result shows tht, t lest sometimes, the extr

More information

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1 Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution

More information

1 1D heat and wave equations on a finite interval

1 1D heat and wave equations on a finite interval 1 1D het nd wve equtions on finite intervl In this section we consider generl method of seprtion of vribles nd its pplictions to solving het eqution nd wve eqution on finite intervl ( 1, 2. Since by trnsltion

More information

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0) 1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this

More information

Second Variation. One-variable problem

Second Variation. One-variable problem Second Vrition. One-vrible problem Februry 11, 213 Contents 1 Locl vritions 2 1.1 Legendre Tests............................ 2 1.2 Weierstrss Tests........................... 4 1.3 Null-Lgrngins nd convexity...................

More information

MAA 4212 Improper Integrals

MAA 4212 Improper Integrals Notes by Dvid Groisser, Copyright c 1995; revised 2002, 2009, 2014 MAA 4212 Improper Integrls The Riemnn integrl, while perfectly well-defined, is too restrictive for mny purposes; there re functions which

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1 MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further

More information

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner

More information

Conservation Law. Chapter Goal. 5.2 Theory

Conservation Law. Chapter Goal. 5.2 Theory Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very

More information

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS. THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem

More information

Math& 152 Section Integration by Parts

Math& 152 Section Integration by Parts Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible

More information

Review of basic calculus

Review of basic calculus Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below

More information

Math 115 ( ) Yum-Tong Siu 1. Lagrange Multipliers and Variational Problems with Constraints. F (x,y,y )dx

Math 115 ( ) Yum-Tong Siu 1. Lagrange Multipliers and Variational Problems with Constraints. F (x,y,y )dx Mth 5 2006-2007) Yum-Tong Siu Lgrnge Multipliers nd Vritionl Problems with Constrints Integrl Constrints. Consider the vritionl problem of finding the extremls for the functionl J[y] = F x,y,y )dx with

More information

Main topics for the First Midterm

Main topics for the First Midterm Min topics for the First Midterm The Midterm will cover Section 1.8, Chpters 2-3, Sections 4.1-4.8, nd Sections 5.1-5.3 (essentilly ll of the mteril covered in clss). Be sure to know the results of the

More information

Advanced Calculus: MATH 410 Uniform Convergence of Functions Professor David Levermore 11 December 2015

Advanced Calculus: MATH 410 Uniform Convergence of Functions Professor David Levermore 11 December 2015 Advnced Clculus: MATH 410 Uniform Convergence of Functions Professor Dvid Levermore 11 December 2015 12. Sequences of Functions We now explore two notions of wht it mens for sequence of functions {f n

More information

Theoretical foundations of Gaussian quadrature

Theoretical foundations of Gaussian quadrature Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct

More information

Review of Riemann Integral

Review of Riemann Integral 1 Review of Riemnn Integrl In this chpter we review the definition of Riemnn integrl of bounded function f : [, b] R, nd point out its limittions so s to be convinced of the necessity of more generl integrl.

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties

More information

Math 8 Winter 2015 Applications of Integration

Math 8 Winter 2015 Applications of Integration Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

More information

221B Lecture Notes WKB Method

221B Lecture Notes WKB Method Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using

More information

7.2 The Definite Integral

7.2 The Definite Integral 7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where

More information

1 The Riemann Integral

1 The Riemann Integral The Riemnn Integrl. An exmple leding to the notion of integrl (res) We know how to find (i.e. define) the re of rectngle (bse height), tringle ( (sum of res of tringles). But how do we find/define n re

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

Math Advanced Calculus II

Math Advanced Calculus II Mth 452 - Advnced Clculus II Line Integrls nd Green s Theorem The min gol of this chpter is to prove Stoke s theorem, which is the multivrible version of the fundmentl theorem of clculus. We will be focused

More information

Sturm-Liouville Eigenvalue problem: Let p(x) > 0, q(x) 0, r(x) 0 in I = (a, b). Here we assume b > a. Let X C 2 1

Sturm-Liouville Eigenvalue problem: Let p(x) > 0, q(x) 0, r(x) 0 in I = (a, b). Here we assume b > a. Let X C 2 1 Ch.4. INTEGRAL EQUATIONS AND GREEN S FUNCTIONS Ronld B Guenther nd John W Lee, Prtil Differentil Equtions of Mthemticl Physics nd Integrl Equtions. Hildebrnd, Methods of Applied Mthemtics, second edition

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

Section 6.1 INTRO to LAPLACE TRANSFORMS

Section 6.1 INTRO to LAPLACE TRANSFORMS Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform

More information

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30 Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová (Mendel University) Definite integrl MENDELU / Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function

More information

Definite integral. Mathematics FRDIS MENDELU

Definite integral. Mathematics FRDIS MENDELU Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the

More information

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)

More information

Lecture 19: Continuous Least Squares Approximation

Lecture 19: Continuous Least Squares Approximation Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for

More information

LECTURE. INTEGRATION AND ANTIDERIVATIVE.

LECTURE. INTEGRATION AND ANTIDERIVATIVE. ANALYSIS FOR HIGH SCHOOL TEACHERS LECTURE. INTEGRATION AND ANTIDERIVATIVE. ROTHSCHILD CAESARIA COURSE, 2015/6 1. Integrtion Historiclly, it ws the problem of computing res nd volumes, tht triggered development

More information

PHYS 4390: GENERAL RELATIVITY LECTURE 6: TENSOR CALCULUS

PHYS 4390: GENERAL RELATIVITY LECTURE 6: TENSOR CALCULUS PHYS 4390: GENERAL RELATIVITY LECTURE 6: TENSOR CALCULUS To strt on tensor clculus, we need to define differentition on mnifold.a good question to sk is if the prtil derivtive of tensor tensor on mnifold?

More information

Generalizations of the Basic Functional

Generalizations of the Basic Functional 3 Generliztions of the Bsic Functionl 3 1 Chpter 3: GENERALIZATIONS OF THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 3.1 Functionls with Higher Order Derivtives.......... 3 3 3.2 Severl Dependent Vribles...............

More information

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS Throughout, we let [, b] be bounded intervl in R. C 2 ([, b]) denotes the spce of functions with derivtives of second order continuous up to the endpoints. Cc 2

More information

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy

More information

Introduction to Finite Element Method

Introduction to Finite Element Method Introduction to Finite Element Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pn.pl/ tzielins/ Tble of Contents 1 Introduction 1 1.1 Motivtion nd generl concepts.............

More information

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1 3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =

More information

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx Chpter 9: Green s functions for time-independent problems Introductory emples One-dimensionl het eqution Consider the one-dimensionl het eqution with boundry conditions nd initil condition We lredy know

More information

Abstract inner product spaces

Abstract inner product spaces WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the

More information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by. NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with

More information

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 UNIFORM CONVERGENCE Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 Suppose f n : Ω R or f n : Ω C is sequence of rel or complex functions, nd f n f s n in some sense. Furthermore,

More information

The Wave Equation I. MA 436 Kurt Bryan

The Wave Equation I. MA 436 Kurt Bryan 1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string

More information

MA Handout 2: Notation and Background Concepts from Analysis

MA Handout 2: Notation and Background Concepts from Analysis MA350059 Hndout 2: Nottion nd Bckground Concepts from Anlysis This hndout summrises some nottion we will use nd lso gives recp of some concepts from other units (MA20023: PDEs nd CM, MA20218: Anlysis 2A,

More information

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite

More information

Introduction to the Calculus of Variations

Introduction to the Calculus of Variations Introduction to the Clculus of Vritions Jim Fischer Mrch 20, 1999 Abstrct This is self-contined pper which introduces fundmentl problem in the clculus of vritions, the problem of finding extreme vlues

More information

Riemann is the Mann! (But Lebesgue may besgue to differ.)

Riemann is the Mann! (But Lebesgue may besgue to differ.) Riemnn is the Mnn! (But Lebesgue my besgue to differ.) Leo Livshits My 2, 2008 1 For finite intervls in R We hve seen in clss tht every continuous function f : [, b] R hs the property tht for every ɛ >

More information

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim Mth 9 Course Summry/Study Guide Fll, 2005 [1] Limits Definition of Limit: We sy tht L is the limit of f(x) s x pproches if f(x) gets closer nd closer to L s x gets closer nd closer to. We write lim f(x)

More information

10 Vector Integral Calculus

10 Vector Integral Calculus Vector Integrl lculus Vector integrl clculus extends integrls s known from clculus to integrls over curves ("line integrls"), surfces ("surfce integrls") nd solids ("volume integrls"). These integrls hve

More information

A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE. In the study of Fourier series, several questions arise naturally, such as: c n e int

A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE. In the study of Fourier series, several questions arise naturally, such as: c n e int A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE HANS RINGSTRÖM. Questions nd exmples In the study of Fourier series, severl questions rise nturlly, such s: () (2) re there conditions on c n, n Z, which ensure

More information

( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that

( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that Arc Length of Curves in Three Dimensionl Spce If the vector function r(t) f(t) i + g(t) j + h(t) k trces out the curve C s t vries, we cn mesure distnces long C using formul nerly identicl to one tht we

More information

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

1.2. Linear Variable Coefficient Equations. y + b ! = a y + b  Remark: The case b = 0 and a non-constant can be solved with the same idea as above. 1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt

More information

The Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.

The Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve. Clculus Li Vs The Fundmentl Theorem of Clculus. The Totl Chnge Theorem nd the Are Under Curve. Recll the following fct from Clculus course. If continuous function f(x) represents the rte of chnge of F

More information

5.7 Improper Integrals

5.7 Improper Integrals 458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the

More information

Lecture 4: Piecewise Cubic Interpolation

Lecture 4: Piecewise Cubic Interpolation Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 4: Piecewise Cubic Interpoltion Compiled 5 September In this lecture we consider piecewise cubic interpoltion

More information

Math 360: A primitive integral and elementary functions

Math 360: A primitive integral and elementary functions Mth 360: A primitive integrl nd elementry functions D. DeTurck University of Pennsylvni October 16, 2017 D. DeTurck Mth 360 001 2017C: Integrl/functions 1 / 32 Setup for the integrl prtitions Definition:

More information

1. On some properties of definite integrals. We prove

1. On some properties of definite integrals. We prove This short collection of notes is intended to complement the textbook Anlisi Mtemtic 2 by Crl Mdern, published by Città Studi Editore, [M]. We refer to [M] for nottion nd the logicl stremline of the rguments.

More information

Chapter 0. What is the Lebesgue integral about?

Chapter 0. What is the Lebesgue integral about? Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous

More information

Quantum Physics II (8.05) Fall 2013 Assignment 2

Quantum Physics II (8.05) Fall 2013 Assignment 2 Quntum Physics II (8.05) Fll 2013 Assignment 2 Msschusetts Institute of Technology Physics Deprtment Due Fridy September 20, 2013 September 13, 2013 3:00 pm Suggested Reding Continued from lst week: 1.

More information

MA 124 January 18, Derivatives are. Integrals are.

MA 124 January 18, Derivatives are. Integrals are. MA 124 Jnury 18, 2018 Prof PB s one-minute introduction to clculus Derivtives re. Integrls re. In Clculus 1, we lern limits, derivtives, some pplictions of derivtives, indefinite integrls, definite integrls,

More information

dt. However, we might also be curious about dy

dt. However, we might also be curious about dy Section 0. The Clculus of Prmetric Curves Even though curve defined prmetricly my not be function, we cn still consider concepts such s rtes of chnge. However, the concepts will need specil tretment. For

More information

Green s functions. f(t) =

Green s functions. f(t) = Consider the 2nd order liner inhomogeneous ODE Green s functions d 2 u 2 + k(t)du + p(t)u(t) = f(t). Of course, in prctice we ll only del with the two prticulr types of 2nd order ODEs we discussed lst

More information

Mapping the delta function and other Radon measures

Mapping the delta function and other Radon measures Mpping the delt function nd other Rdon mesures Notes for Mth583A, Fll 2008 November 25, 2008 Rdon mesures Consider continuous function f on the rel line with sclr vlues. It is sid to hve bounded support

More information

Chapter 6 Notes, Larson/Hostetler 3e

Chapter 6 Notes, Larson/Hostetler 3e Contents 6. Antiderivtives nd the Rules of Integrtion.......................... 6. Are nd the Definite Integrl.................................. 6.. Are............................................ 6. Reimnn

More information

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ), 1. Guss-Jcobi qudrture nd Legendre polynomils Simpson s rule for evluting n integrl f(t)dt gives the correct nswer with error of bout O(n 4 ) (with constnt tht depends on f, in prticulr, it depends on

More information

Best Approximation in the 2-norm

Best Approximation in the 2-norm Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion

More information

Recitation 3: More Applications of the Derivative

Recitation 3: More Applications of the Derivative Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech

More information

Line Integrals. Partitioning the Curve. Estimating the Mass

Line Integrals. Partitioning the Curve. Estimating the Mass Line Integrls Suppose we hve curve in the xy plne nd ssocite density δ(p ) = δ(x, y) t ech point on the curve. urves, of course, do not hve density or mss, but it my sometimes be convenient or useful to

More information

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), ) Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s

More information

Heat flux and total heat

Heat flux and total heat Het flux nd totl het John McCun Mrch 14, 2017 1 Introduction Yesterdy (if I remember correctly) Ms. Prsd sked me question bout the condition of insulted boundry for the 1D het eqution, nd (bsed on glnce

More information

INTRODUCTION TO INTEGRATION

INTRODUCTION TO INTEGRATION INTRODUCTION TO INTEGRATION 5.1 Ares nd Distnces Assume f(x) 0 on the intervl [, b]. Let A be the re under the grph of f(x). b We will obtin n pproximtion of A in the following three steps. STEP 1: Divide

More information

Math 100 Review Sheet

Math 100 Review Sheet Mth 100 Review Sheet Joseph H. Silvermn December 2010 This outline of Mth 100 is summry of the mteril covered in the course. It is designed to be study id, but it is only n outline nd should be used s

More information

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a). The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples

More information

Recitation 3: Applications of the Derivative. 1 Higher-Order Derivatives and their Applications

Recitation 3: Applications of the Derivative. 1 Higher-Order Derivatives and their Applications Mth 1c TA: Pdric Brtlett Recittion 3: Applictions of the Derivtive Week 3 Cltech 013 1 Higher-Order Derivtives nd their Applictions Another thing we could wnt to do with the derivtive, motivted by wht

More information

MORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)

MORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.) MORE FUNCTION GRAPHING; OPTIMIZATION FRI, OCT 25, 203 (Lst edited October 28, 203 t :09pm.) Exercise. Let n be n rbitrry positive integer. Give n exmple of function with exctly n verticl symptotes. Give

More information

MatFys. Week 2, Nov , 2005, revised Nov. 23

MatFys. Week 2, Nov , 2005, revised Nov. 23 MtFys Week 2, Nov. 21-27, 2005, revised Nov. 23 Lectures This week s lectures will be bsed on Ch.3 of the text book, VIA. Mondy Nov. 21 The fundmentls of the clculus of vritions in Eucliden spce nd its

More information

g i fφdx dx = x i i=1 is a Hilbert space. We shall, henceforth, abuse notation and write g i f(x) = f

g i fφdx dx = x i i=1 is a Hilbert space. We shall, henceforth, abuse notation and write g i f(x) = f 1. Appliction of functionl nlysis to PEs 1.1. Introduction. In this section we give little introduction to prtil differentil equtions. In prticulr we consider the problem u(x) = f(x) x, u(x) = x (1) where

More information

Integration Techniques

Integration Techniques Integrtion Techniques. Integrtion of Trigonometric Functions Exmple. Evlute cos x. Recll tht cos x = cos x. Hence, cos x Exmple. Evlute = ( + cos x) = (x + sin x) + C = x + 4 sin x + C. cos 3 x. Let u

More information

Thin-Plate Splines. Contents

Thin-Plate Splines. Contents Thin-Plte Splines Dvid Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Cretive Commons Attribution 4.0 Interntionl License. To view copy of this

More information

Plates on elastic foundation

Plates on elastic foundation Pltes on elstic foundtion Circulr elstic plte, xil-symmetric lod, Winkler soil (fter Timoshenko & Woinowsky-Krieger (1959) - Chpter 8) Prepred by Enzo Mrtinelli Drft version ( April 016) Introduction Winkler

More information

x = b a n x 2 e x dx. cdx = c(b a), where c is any constant. a b

x = b a n x 2 e x dx. cdx = c(b a), where c is any constant. a b CHAPTER 5. INTEGRALS 61 where nd x = b n x i = 1 (x i 1 + x i ) = midpoint of [x i 1, x i ]. Problem 168 (Exercise 1, pge 377). Use the Midpoint Rule with the n = 4 to pproximte 5 1 x e x dx. Some quick

More information

221A Lecture Notes WKB Method

221A Lecture Notes WKB Method A Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using ψ x, t = e

More information

Chapter 28. Fourier Series An Eigenvalue Problem.

Chapter 28. Fourier Series An Eigenvalue Problem. Chpter 28 Fourier Series Every time I close my eyes The noise inside me mplifies I cn t escpe I relive every moment of the dy Every misstep I hve mde Finds wy it cn invde My every thought And this is why

More information

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl

More information