Appendix. M. Germain, A. Magnus, and V. van Steenberghe. December 2004.
|
|
- Stewart Goodman
- 5 years ago
- Views:
Transcription
1 Appendix to Should developing countries participate in the Clean Development Mechanism under the Kyoto Protocol? The low-hanging fruits and baseline issues by M. Germain, A. Magnus, and V. van Steenberghe. December We give here the mathematical derivation of some propositions and statements of the paper. Contents. Preliminaries: the second period imization problem The absolute baseline case a. Proof of proposition 2.a (existence, unicity and characterization of λ a ) a.. The equation for λ a : the imum condition a.2. Existence of λ a : the first order condition a.3. Unicity of the solution of the first order condition b. Proof of Proposition 2.b (behaviour of e,, and y with respect to τ ) The relative baseline case a. Proof of proposition 3.a (existence, unicity and characterization of λ r ) b. Discussion of Result 3.b (behaviour of e,, and y with respect to τ ) c. Proof of Proposition 3.c (comparison of y and under both baselines) d. Proof of Proposition 3.d (comparison of Π under both baselines)... 8
2 Appendix to low hanging fruits etc. 2. Preliminaries: the second period imization problem. We see that the full problem (7) contains several times the same subproblem: we first have to imize with respect to a variable e 0 an expression of the form Ae γ Be+C, with A,B > 0, and γ <. The derivative γae γ B has exactly one positive zero: [ γa /[ γ e =, (A.) B which leads to Be = γae γ, so that the imum is [ γae γ + C = γ A /[ γ B γ/[ γ + C, where γ = [ γγ γ/[ γ. (A.2) We apply now this to Π 2. fig A. In Π 2 (λ,τ 2 ), for a given τ 2, A and B depend on λ 2 : A = λ β 2, B = p e+τ 2 +p λ 2, and so β[p e+τ 2 fig. A.2 β[p e+τ 2 F λ 2 Π 2 λ Π 2 (λ,τ 2 ) = λ 2 λ F(λ 2 ;τ 2 ) + C, where F(λ 2 ;τ 2 ) := e2>0 eγ 2 λ 2 [p e + τ 2 + p λ 2 e 2 = γ λ β/[ γ 2 [p e + τ 2 + p λ 2 γ/[ γ, (A.3) F(λ 2 ;τ 2 ) increases between λ 2 = 0 and λ 2 = β[p e + τ 2 /[, and decreases afterwards (see fig. A.). Therefore, when λ β[p e + τ 2 /[, the imum over λ 2 λ of F(λ 2 ;τ 2 ) is the imal value of F, i.e., F ( ) β[pe + τ 2 ;τ 2 = γ [p e + τ 2 α/[ γ p β/[ γ αp, (A.4) where γ = γ β β/[ γ α α/[ γ /γ γ/[ γ = [ γβ β/[ γ α α/[ γ ; whereas, when λ β[p e + τ 2 /[, the imum over λ 2 λ of F(λ 2 ;τ 2 ) is merely F(λ ;τ 2 ). Thus: Π 2 (λ,τ 2 ) = γ λ β/[ γ [p e + τ 2 + p λ γ/[ γ + C, if λ β[p e + τ 2 /[, (A.5) = γ [p e + τ 2 α/[ γ p β/[ γ + C if λ β[p e + τ 2 /[, Remar that Π 2 (λ,τ 2 ) is a non increasing function of λ (see figure A.2). This also holds for the integral term of (7), as this term is a positive linear combination of functions Π 2 (λ,τ 2 ) for an interval of values of τ 2.
3 Appendix to low hanging fruits etc The absolute baseline case. 2.a. Proof of proposition 2.a (existence, unicity and characterization of λ a ). 2.a.. The equation for λ a : the imum condition. We now come to the full problem (4), with an absolute baseline: e 0,λ 0 eγ λβ [p e + p λ e + τ [e e τ Emphasizing the dependence in e, using (2): λ 0 e 0 f(τ 2 ) e γ 2 e 2 0 λβ 2 [p e + p λ 2 e 2 + τ 2 [ẽ 2 δ[e e e 2 λ 2 λ e γ λβ [p e + p λ + τ ρδ τ 2 e τ f(τ 2 ) λ 2 λ e 2 0 which is first solved with respect to e : Π a (λ,τ ) := F(λ ;τ ρδ τ 2 ) + ρ λ 0 dτ 2. e γ 2 λβ 2 [p e + p λ 2 + τ 2 e 2 dτ 2 + τ e + τ 2 ρ[ẽ 2 δe τ f(τ 2 )[ F(λ 2 ;τ 2 )dτ 2 + τ e + τ 2 ρ[ẽ 2 δe, λ 2 λ using the definition of F in (A.3). Recall that the first term is an increasing function of λ as long as λ λ ref, where (A.6) λ ref := β[p e + τ ρδ τ 2, (A.7) and that the second term contains only non-increasing functions of λ. Moreover, if λ is small enough, so that the condition seen above λ β[p e + τ 2 /[ holds for all the τ 2 s in the integral, i.e., if λ β[p e + /[, all the Π 2 s are constant functions of λ, and the net result is a function Π a behaving lie the first term F plus a constant. case (2) fig. A.3 case () F+const. β[p e+ Π a λ λ ref λ a λ ref β[p e+ fig. A.4 F+const. Π a λ When λ > β[p e + /[, the sum is F plus an actually decreasing function of λ, and is smaller than the continuation of the rule F plus a constant (shown by a dashed line in the figures A.3 and A.4 nearby, whereas the actual Π a [solid line is even smaller). Therefore, () If τ τ 2min + ρδ τ 2, Π a reaches its imum value at λ = λ ref = β[p e + τ ρδ τ 2 /[, as this latter value is smaller than β[p e + /[, which is the place where a more complicated behaviour occurs, but we don t have to care, as the imum has already been encountered. (2) If τ > + ρδ τ 2, the actual imum of Π a occurs between β[p e + /[ and λ ref = β[p e +τ ρδ τ 2 /[, as the λ derivative is still positive at the first bound (F is still increasing, and the integral term still contains constant values with respect to λ ), and is negative at the second bound (the integral now contains actually decreasing functions of λ ).
4 Appendix to low hanging fruits etc. 4 We have therefore established parts (i) and (iii) of Proposition 2.a, i.e., that the imum is reached at λ a = λ ref in case (); that β[p e + τ 2min < λ a < λ ref in case (2). 2.a.2. Existence of λ a : the first order condition. In order to establish more properties of λ a when τ > τ 2min + ρδ τ 2 (case (2) above) we get a closer loo at the equation for λ a from the λ derivative Π a (λ,τ ) = τ2 F(λ ;τ ρδ τ 2 ) + ρ f(τ 2 )[ F(λ 2 ;τ 2 )dτ 2 = 0, λ λ λ 2 λ or F(λ ;τ ρδ τ 2 ) αp λ /β p e + ρ f(τ 2 ) F(λ ;τ 2 ) dτ 2 = 0, λ τ 2min λ as we saw in (A.5) that Π 2 (λ,τ 2 ) is a constant function with respect to λ while λ β[p e + τ 2 /[, and leaves therefore no trace in the integral of the λ derivative, so that we only have to consider the derivative when λ > β[p e + τ 2 /[, or τ 2 < λ /β p e in the integral. Remar that the upper bound in the integral is actually larger than the lower bound, i.e., λ /β p e > τ 2min, as λ > β[p e + /[ (case (2) above). We now perform the λ derivative of F(λ;τ) = const. λ β/[ γ 2 [p e + τ + p λ 2 γ/[ γ : βλ β/[ γ αp λ /β p e [p e + τ + p λ ρδ τ 2 γ/[ γ γp λ β/[ γ [p e + τ + p λ ρδ τ 2 γ/[ γ f(τ 2 )βλ β/[ γ [p e +τ 2 +p λ γ/[ γ γp λ β/[ γ [p e +τ 2 +p λ γ/[ γ dτ 2 = 0, The condition on λ amounts to [β[p e + τ ρδ τ 2 λ [p e + τ + p λ ρδ τ 2 /[ γ fig. A.5 G β[p e+ λ + ρ αp λ /β p e f(τ 2 )[β[p e + τ 2 λ [p e + τ 2 + p λ /[ γ dτ 2 = 0. Remar that the first term is positive, while the second term is negative when β[p e + τ 2min /[ < λ < λ ref = β[p e + τ ρδ τ 2 /[. Managing to eep only positive terms: G(λ,τ ) =, where λ ref αp λ /β p e λ β[p e + τ 2 [p e + τ + p λ ρδ τ 2 /[ γ G(λ,τ ) = ρ f(τ 2 ) β[p e + τ ρδ τ 2 λ [p e + τ 2 + p λ /[ γ dτ 2. (A.8) G(λ,τ ) = 0 when λ = β[p e + /[, and G(λ,τ ) + when λ β[p e + τ ρδ τ 2 /[, so there must be at least one intermediate value of λ where G = (see fig. A.5). This is another way to ensure parts (i) and (iii) of Proposition 2.a, i.e., the existence of a solution β[p e + /[ < λ < λ ref = β[p e + τ ρδ τ 2 /[.
5 Appendix to low hanging fruits etc. 5 2.a.3. Unicity of the solution of the first order condition. Π a? fig. A.6 λ When λ imizes (7), the equation G = is certainly solved. The converse is not clear, as solutions of G =, i.e., of Π a / λ = 0, may be any ind of stationary point for (7) (see fig. A.6). The situation is non ambiguous when G is a monotonous function of λ, as G = can then have at most one root. We will loo for a condition ensuring that G(λ,τ ) is an increasing function of λ. Thus, we have to consider [ G(λ,τ ) λ αp λ /β p e = ρ f(τ 2 )[τ ρδ τ 2 τ 2 p γ αβ[ γ [ λ β[p e + τ 2 [β[p e + τ ρδ τ 2 λ [p e + τ + p λ ρδ τ 2 [p e + τ 2 + p λ λ β[p e + τ 2 [p e + τ + p λ ρδ τ 2 /[ γ β[p e + τ ρδ τ 2 λ [p e + τ 2 + p λ /[ γ dτ 2 τ ρδ τ 2 τ 2 is positive in the whole integration interval, as τ 2 < λ /β p e, so λ > β[p e + τ 2 /[, and we saw that λ < β[p e + τ ρδ τ 2 /[. αp We also see that the ratio λ β[p e + τ 2 is positive for τ β[p e + τ ρδ τ 2 λ 2 in the (open) integration interval. To be sure that G is increasing, it remains to show that the big intermediate factor is positive too, or that αβ[ γ > λ β[p e + τ 2 β[p e + τ ρδ τ 2 λ p e + τ 2 + p λ p e + τ + p λ ρδ τ 2 on the whole τ 2 interval, i.e., at τ 2 = τ 2min, as the right-hand side is a decreasing function of τ 2. And as we do not now more about λ that β[p e + τ 2min /[ < λ < β[p e + τ ρδ τ 2 /[, we tae the first fraction of the right hand side at its largest possible value, which is reached at λ = β[p e +τ ρδ τ 2 /[, and so for the second fraction, but this one at λ = β[p e + /[. The sufficient condition of unicity is then αβ[ γ > α 2 β 2 [τ ρδ τ 2 2 [α[p e + τ 2min + β[p e + τ ρδ τ 2 [β[p e + τ 2min + α[p e + τ ρδ τ 2 which depends only on nown values. The condition can be rewored as where x = τ ρδ τ 2 p e +, or γ < stronger sufficient condition γ > αβx 2 [γ + βx[γ + αx, + αβ x 2, or, as αβ is always smaller than γ 2 /4, a γ 2 + x γ < + x + x + x 2 /4. (A.9) For a given x (or, equivalently, τ ), returns to scale must be sufficiently decreasing (γ sufficiently low), and conversely. So that part (ii) of Proposition 2.a is established.
6 Appendix to low hanging fruits etc. 6 2.b. Proof of Proposition 2.b (behaviour of e,, and y with respect to τ ). We now discuss e,, and y, which are positive powers of λ divided by powers of p λ + p e + τ ρδ τ 2 : constant λ a /[p λ + p e + τ ρδ τ 2 b, where the coefficients a and b are a b e γ β γ α y γ β γ Remar that a < b. Things are simple when τ + ρδ τ 2, as exact formulas are available: λ = β[p e + τ ρδ τ 2 /[, so that λ a [p λ + p e + τ ρδ τ 2 b = constant [p e + τ ρδ τ 2 a b is obviously a decreasing function of τ. When τ > τ 2min + ρδ τ 2, things are less explicit. We discuss first variations of λ with respect to τ. As G must eep a fixed value G =, dg = [ G/ λ dλ + [ G/ τ dτ = 0, so, dλ /dτ = [ G/ τ /[ G/ λ. As seen above, G being an increasing function of λ, G/ λ > 0 (see fig. A.5). For a fixed λ, G(λ,τ ) will increase or decrease according to the sign of G/ τ which, from (A.8), is the sign of [p e + τ + p λ ρδ τ 2 /[ γ = γ[β[p e +τ ρδ τ 2 [ βp λ [p e + τ + p λ ρδ τ 2 γ/[ γ τ β[p e + τ ρδ τ 2 λ [β[p e + τ ρδ τ 2 λ 2. So, we see that G decreases when τ increases for a fixed λ, as long as λ β[p e +τ ρδ τ 2 /[p [ β. Therefore, λ will increase with τ as long as λ β[p e + τ ρδ τ 2 /[p [ β, and decreases if λ is lower than this bound. Let us derivate λ a /[p λ + p e + τ ρδ τ 2 b with respect to τ. The sign of the derivative is the sign of [[a bp λ + a[p e + τ ρδ τ 2 dλ dτ bλ. For e, a = β and b =, the sign of dλ /dτ is exactly the sign of [ βp λ β[p e +τ ρδ τ 2, as seen above. So, e is a decreasing function of τ, and this settles the first part of Proposition 2.b. For and y, it will be shown here that these functions of τ are smaller than their values at λ ref = β[p e + τ ρδ τ 2 /[, i.e., that [λ a a [p λ a + p e + τ ρδ τ 2 b < λ a ref [p λ ref + p e + τ ρδ τ 2 b, nowing that λ a < λ ref. Indeed, λ a /[p λ + p e + τ ρδ τ 2 b is an increasing function of λ up to λ = β[p e + τ ρδ τ 2 /[ if a[p λ + p e + τ ρδ τ 2 bp λ = [[a bβ/α + a[p e + τ ρδ τ 2 0, or aγ bβ, which holds for and y. Now, the upper bound λ a ref /[p λ ref +p e +τ ρδ τ 2 b = const. [p e + τ ρδ τ 2 a b is a decreasing function of τ, so that the remaining part of Proposition 2.b is established.
7 [ + τ e y /[ γ Appendix to low hanging fruits etc. 7 We rewrite (4) with a relative baseline: e 0,λ 0 eγ λβ [p e + p λ e + τ τ 3. The relative baseline case. [ e e y y Emphasizing the dependence in e, using () and (2): [ e + τ λ 0 e 0 y τ f(τ 2 ) e γ 2 e 2 0 λβ 2 [p e + p λ 2 e 2 + τ 2 [ẽ 2 δ[e e e 2 λ 2 λ e γ λβ [p e + p λ + τ ρδ τ 2 e f(τ 2 ) λ 2 λ which is first solved with respect to e : λ 0 Π r (λ,τ ) := e 2 0 dτ 2. e γ 2 λβ 2 [p e + p λ 2 + τ 2 e 2 dτ 2 + τ 2 ρ[ẽ 2 δe [ e /[ γ + τ F(λ,τ ρδ τ 2 ) y τ f(τ 2 ) F(λ 2,τ 2 ) dτ 2 + τ 2 [ẽ 2 δe (A.0) λ 2 λ where, as before, F is given by (A.3). Optimal emissions of period with respect to λ and τ are given by [ e r = λ β [ + τ e /y γ /[ γ p e + τ + p λ ρδ τ 2 Note that e and y are the solutions of the elementary problem when τ = 0: λ = β[p e ρδ τ 2 /[, e γ = γλ β /[p e + p λ ρδ τ 2 = βλ β /p, or e/y = e γ λ β = β/[p λ, or e = α [ β/[ γ β β/[ γ [p e ρδ τ 2 [β /[ γ p β/[ γ, (A.) e y = α p e ρδ τ 2. 3.a. Proof of proposition 3.a (existence, unicity and characterization of λ r ). (A.2) As seen before, the integral does not depend on λ if τ τ 2min + ρδ τ 2, and the imum is reached at λ = λ ref = β[p e + τ ρδ τ 2 /[, as in the absolute case. Fig. A.7 G λ r λ a λ ref When τ > + ρδ τ 2, we consider the λ derivative which is [ e /[ γ + τ [β[p e + τ ρδ τ 2 λ [p e + τ + p λ ρδ τ 2 /[ γ λ y αp λ /β p e + ρ f(τ 2 )[β[p e + τ 2 λ [p e + τ 2 + p λ /[ γ dτ 2 = 0,
8 Appendix to low hanging fruits etc. 8 or G(λ,τ ) = [ e /[ γ + τ, with the same G as above, in (A.8). As before, G(λ,τ ) = 0 when λ = β[p e + /[, and G(λ,τ ) + when λ β[p e + τ ρδ τ 2 /[, so there must be at least one intermediate value of λ where the equation for λ is satisfied: β[p e + /[ < λ < β[p e +τ ρδ τ 2 /[ (see fig. A.7). As the right-hand side is larger than, and if G is an increasing function of λ (remember the sufficient condition (A.9) γ [ + x/[/x + x 2 /4, which is still valid here), we have β[p e + τ 2min /[ < λ a < λ r < β[p e + τ ρδ τ 2 /[, and the whole Proposition 3.a is established. 3.b. Discussion of Result 3.b (behaviour of e,, and y with respect to τ ). e, = λ e, and y = e γ λβ now contain a power of + τ e y too, so e,, and y write lie constant λ a [p λ + p e + τ ρδ τ 2 [ b + τ e c, y where the coeficients a, b, c are a b c e γ β γ α y γ β γ γ We loo at the variation with respect to τ assuming λ = β[p e + τ ρδ τ 2 /[, which is exact when τ τ 2min + ρδ τ 2, and found numerically to be very close for larger τ s. Then, the above expresion for e,, y becomes: constant [p e + τ ρδ τ 2 a b [ + τ e /y c, = constant [p e + τ ρδ τ 2 a b [p e + ατ ρδ τ 2 c. The τ derivative at τ = 0 has the sign of a b+αc = γ,0, and α[γ for e,, and y. Whereas for large τ, the behaviour is τ a b+c always increasing in the long run (β, α, and β), which are indications of a U shape behaviour of e,, and y with respect to τ. This U shape behaviour is indeed confirmed numerically. Finally, e and y, computed assuming λ r = λ ref, are minimal when [a b + αc[p e ρδ τ 2 + α[a b + cτ = 0, so argmin y = α argmin e. This closes the discussion of result 3.b. 3.c. Proof of Proposition 3.c (comparison of y and under both baselines). We already saw that and y are increasing functions of λ, so r a and y r y a hold, as λ a λ r λ ref = β[p e + τ ρδ τ 2 /[, and as there is a further factor + τ e /y > in the relative case. Note that for e, the expected result is not established, but the factor + τ e /y present in the formula for e r maes the inequality e r e a valid in most cases. Indeed, the inequality is verified numerically. 3.d. Proof of Proposition 3.d (comparison of Π under both baselines). Π r versus Π a : We show that Π r Π a while τ is smaller than the smallest positive root of p e ρδ τ 2 X γ [ [ + αx /[ γ [ + X α/[ γ = 0. (A.3) α From (A.6), (A.0), and (A.3), we may write for a given τ : Π a (λ ) = H a (λ ) + K(λ ), where H a (λ ) = F(λ ;τ ρδ τ 2 ) + τ e ; y
9 Appendix to low hanging fruits etc. 9 [ and Π r (λ ) = H r (λ ) + K(λ ), where H r (λ ) = + τ e /[ γ y F(λ ;τ ρδ τ 2 ) ; and where K(λ ) is the integral term K(λ ) = ρ τ f(τ 2 ) λ2 λ F(λ 2,τ 2 )dτ 2 + τ 2 ρ[ẽ 2 δe. Then, Π a (λ a ) Π r (λ r ) [Π a (λ ref ) Π r (λ ref ) = = 0, λ a Π a(λ )dλ + λr λ a Π a(λ )dλ λ r Π r(λ )dλ λ r [H a(λ ) + K (λ )dλ + λ r Π a(λ )dλ + λ r Π r(λ )dλ λ r [H r(λ ) + K (λ )dλ = λ r [H r(λ ) H a(λ )dλ as Π a is negative between λ a and λ r, Π a = H a + K being a decreasing function of λ when λr λ > λ a (see fig. A.4), so that Π a(λ )dλ 0. λ a Finally, H r and H a are positive as H r and H a, both of the form const. F+ const. (see section ), reach their imum at λ = λ ref, and H r = [ + τ e /y /[ γ H a > H a. So, Π a (λ a ) Π r (λ r ) Π a (λ ref ) Π r (λ ref ). We even have an equality while τ τ 2min + ρδ τ 2, as λ a = λ r = λ ref in this case. We compute now the right-hand side in order to estimate the τ interval where this lower bound is positive. It is [ [ e /[ γ Π a (λ ref ) Π r (λ ref ) = τ e + τ γ λ β/[ γ [p e + τ + p λ ρδ τ 2 γ/[ γ y = τ α [ β/[ γ β β/[ γ [p e ρδ τ 2 [β /[ γ p β/[ γ [ [ ατ /[ γ + γ λ β/[ γ p e ρδ τ ref [p e + τ + p λ ref ρδ τ 2 γ/[ γ 2 [ [ τ = const. γ ατ /[ γ [ τ α/[ γ + + p e ρδ τ 2 α p e ρδ τ 2 p e ρδ τ 2 = const. X γ [ [ + αx /[ γ [ + X α/[ γ, (A.4) α from (A.), (A.2), and λ = β[p e ρδ τ 2 τ, and where X =. p e ρδ τ 2 The lower bound Π a (λ ref ) Π r (λ ref ) is positive for small positive τ, as the Taylor expansion of α[2 γ (A.4) with respect to X about X = 0 starts with const. 2[ γ X2 > 0, but the bound behaves lie X [ α/[ γ < 0 for large X, so there is a finite lowest positive root, however comfortably large if γ is not too close to. An empirical formula for a valid interval is 0 < X < 6 γ.
3.2 Configuration model
3.2 Configuration model 3.2.1 Definition. Basic properties Assume that the vector d = (d 1,..., d n ) is graphical, i.e., there exits a graph on n vertices such that vertex 1 has degree d 1, vertex 2 has
More informationWAITING FOR A BAT TO FLY BY (IN POLYNOMIAL TIME)
WAITING FOR A BAT TO FLY BY (IN POLYNOMIAL TIME ITAI BENJAMINI, GADY KOZMA, LÁSZLÓ LOVÁSZ, DAN ROMIK, AND GÁBOR TARDOS Abstract. We observe returns of a simple random wal on a finite graph to a fixed node,
More informationNotes on the Z-transform, part 4 1
Notes on the Z-transform, part 4 5 December 2002 Solving difference equations At the end of note 3 we saw how to solve a difference equation using Z-transforms. Here is a similar example, solve x k+2 6x
More informationDependence and independence
Roberto s Notes on Linear Algebra Chapter 7: Subspaces Section 1 Dependence and independence What you need to now already: Basic facts and operations involving Euclidean vectors. Matrices determinants
More informationUniversity of Utrecht
University of Utrecht Bachelor thesis The bifurcation diagram of the second nontrivial normal form of an axially symmetric perturbation of the isotropic harmonic oscillator Author: T. Welker Supervisor:
More informationThe fine-grained Gibbs entropy
Chapter 12 The fine-grained Gibbs entropy 12.1 Introduction and definition The standard counterpart of thermodynamic entropy within Gibbsian SM is the socalled fine-grained entropy, or Gibbs entropy. It
More informationGraphing Linear Inequalities
Graphing Linear Inequalities Linear Inequalities in Two Variables: A linear inequality in two variables is an inequality that can be written in the general form Ax + By < C, where A, B, and C are real
More informationSection 2.4: Add and Subtract Rational Expressions
CHAPTER Section.: Add and Subtract Rational Expressions Section.: Add and Subtract Rational Expressions Objective: Add and subtract rational expressions with like and different denominators. You will recall
More informationNOTES ON EXISTENCE AND UNIQUENESS THEOREMS FOR ODES
NOTES ON EXISTENCE AND UNIQUENESS THEOREMS FOR ODES JONATHAN LUK These notes discuss theorems on the existence, uniqueness and extension of solutions for ODEs. None of these results are original. The proofs
More informationMechanics Physics 151
Mechanics Physics 151 Lecture 15 Special Relativity (Chapter 7) What We Did Last Time Defined Lorentz transformation Linear transformation of 4-vectors that conserve the length in Minkowski space Derived
More informationFor those who want to skip this chapter and carry on, that s fine, all you really need to know is that for the scalar expression: 2 H
1 Matrices are rectangular arrays of numbers. hey are usually written in terms of a capital bold letter, for example A. In previous chapters we ve looed at matrix algebra and matrix arithmetic. Where things
More informationOpen quantum systems
Chapter 4 Open quantum systems 4.1 Quantum operations Let s go bac for a second to the basic postulates of quantum mechanics. Recall that when we first establish the theory, we begin by postulating that
More informationSeveral Generating Functions for Second-Order Recurrence Sequences
47 6 Journal of Integer Sequences, Vol. 009), Article 09..7 Several Generating Functions for Second-Order Recurrence Sequences István Mező Institute of Mathematics University of Debrecen Hungary imezo@math.lte.hu
More informationMathematical modelling and numerical analysis of thin elastic shells
Mathematical modelling and numerical analysis of thin elastic shells Zhang Sheng Wayne State University, Michigan Zhang Sheng (szhang@wayne.edu) Thin elastic shell May, 202 / 9 Models for mechanics of
More informationCournot and Bertrand Competition in a Differentiated Duopoly with Endogenous Technology Adoption *
ANNALS OF ECONOMICS AND FINANCE 16-1, 231 253 (2015) Cournot and Bertrand Competition in a Differentiated Duopoly with Endogenous Technology Adoption * Hongkun Ma School of Economics, Shandong University,
More informationCollected trivialities on algebra derivations
Collected trivialities on algebra derivations Darij Grinberg December 4, 2017 Contents 1. Derivations in general 1 1.1. Definitions and conventions....................... 1 1.2. Basic properties..............................
More informationand the compositional inverse when it exists is A.
Lecture B jacques@ucsd.edu Notation: R denotes a ring, N denotes the set of sequences of natural numbers with finite support, is a generic element of N, is the infinite zero sequence, n 0 R[[ X]] denotes
More informationAP Calculus Chapter 9: Infinite Series
AP Calculus Chapter 9: Infinite Series 9. Sequences a, a 2, a 3, a 4, a 5,... Sequence: A function whose domain is the set of positive integers n = 2 3 4 a n = a a 2 a 3 a 4 terms of the sequence Begin
More informationTheorem (Special Case of Ramsey s Theorem) R(k, l) is finite. Furthermore, it satisfies,
Math 16A Notes, Wee 6 Scribe: Jesse Benavides Disclaimer: These notes are not nearly as polished (and quite possibly not nearly as correct) as a published paper. Please use them at your own ris. 1. Ramsey
More informationLECTURE 21: THE HESSIAN, LAPLACE AND TOPOROGOV COMPARISON THEOREMS. 1. The Hessian Comparison Theorem. We recall from last lecture that
LECTURE 21: THE HESSIAN, LAPLACE AND TOPOROGOV COMPARISON THEOREMS We recall from last lecture that 1. The Hessian Comparison Theorem K t) = min{kπ γt) ) γt) Π γt) }, K + t) = max{k Π γt) ) γt) Π γt) }.
More informationAn Introduction to Ergodic Theory
An Introduction to Ergodic Theory Normal Numbers: We Can t See Them, But They re Everywhere! Joseph Horan Department of Mathematics and Statistics University of Victoria Victoria, BC December 5, 2013 An
More informationREAL ANALYSIS LECTURE NOTES: 1.4 OUTER MEASURE
REAL ANALYSIS LECTURE NOTES: 1.4 OUTER MEASURE CHRISTOPHER HEIL 1.4.1 Introduction We will expand on Section 1.4 of Folland s text, which covers abstract outer measures also called exterior measures).
More informationSection 12.6: Non-homogeneous Problems
Section 12.6: Non-homogeneous Problems 1 Introduction Up to this point all the problems we have considered are we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationA Propagating Wave Packet The Group Velocity
Lecture 7 A Propagating Wave Pacet The Group Velocity Phys 375 Overview and Motivation: Last time we looed at a solution to the Schrödinger equation (SE) with an initial condition (,) that corresponds
More informationFACTORING IN QUADRATIC FIELDS. 1. Introduction
FACTORING IN QUADRATIC FIELDS KEITH CONRAD For a squarefree integer d other than 1, let 1. Introduction K = Q[ d] = {x + y d : x, y Q}. This is closed under addition, subtraction, multiplication, and also
More informationLecture 6 Random walks - advanced methods
Lecture 6: Random wals - advanced methods 1 of 11 Course: M362K Intro to Stochastic Processes Term: Fall 2014 Instructor: Gordan Zitovic Lecture 6 Random wals - advanced methods STOPPING TIMES Our last
More informationSection 3.1: Definition and Examples (Vector Spaces), Completed
Section 3.1: Definition and Examples (Vector Spaces), Completed 1. Examples Euclidean Vector Spaces: The set of n-length vectors that we denoted by R n is a vector space. For simplicity, let s consider
More informationNow, the above series has all non-negative terms, and hence is an upper bound for any fixed term in the series. That is to say, for fixed n 0 N,
l p IS COMPLETE Let 1 p, and recall the definition of the metric space l p : { } For 1 p
More informationTime series models in the Frequency domain. The power spectrum, Spectral analysis
ime series models in the Frequency domain he power spectrum, Spectral analysis Relationship between the periodogram and the autocorrelations = + = ( ) ( ˆ α ˆ ) β I Yt cos t + Yt sin t t= t= ( ( ) ) cosλ
More informationNUCLEAR NORM PENALIZED ESTIMATION OF INTERACTIVE FIXED EFFECT MODELS. Incomplete and Work in Progress. 1. Introduction
NUCLEAR NORM PENALIZED ESTIMATION OF IERACTIVE FIXED EFFECT MODELS HYUNGSIK ROGER MOON AND MARTIN WEIDNER Incomplete and Work in Progress. Introduction Interactive fixed effects panel regression models
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationKiddie Talk - The Diamond Lemma and its applications
Kiddie Tal - The Diamond Lemma and its applications April 22, 2013 1 Intro Start with an example: consider the following game, a solitaire. Tae a finite graph G with n vertices and a function e : V (G)
More informationContest Number Theory
Contest Number Theory Andre Kessler December 7, 2008 Introduction Number theory is one of the core subject areas of mathematics. It can be somewhat loosely defined as the study of the integers. Unfortunately,
More informationMT101-03, Area Notes. We briefly mentioned the following triumvirate of integration that we ll be studying this semester:
We briefly mentioned the following triumvirate of integration that we ll be studying this semester: AD AUG MT101-03, Area Notes NC where AUG stands for Area Under the Graph of a positive elementary function,
More informationAppendix A: Matrices
Appendix A: Matrices A matrix is a rectangular array of numbers Such arrays have rows and columns The numbers of rows and columns are referred to as the dimensions of a matrix A matrix with, say, 5 rows
More information(V.C) Complex eigenvalues and other topics
V.C Complex eigenvalues and other topics matrix Let s take a closer look at the characteristic polynomial for a 2 2 namely A = a c f A λ = detλi A = λ aλ d bc b d = λ 2 a + dλ + ad bc = λ 2 traλ + det
More informationare the q-versions of n, n! and . The falling factorial is (x) k = x(x 1)(x 2)... (x k + 1).
Lecture A jacques@ucsd.edu Notation: N, R, Z, F, C naturals, reals, integers, a field, complex numbers. p(n), S n,, b(n), s n, partition numbers, Stirling of the second ind, Bell numbers, Stirling of the
More informationPhysics Oct A Quantum Harmonic Oscillator
Physics 301 5-Oct-2005 9-1 A Quantum Harmonic Oscillator The quantum harmonic oscillator (the only kind there is, really) has energy levels given by E n = (n + 1/2) hω, where n 0 is an integer and the
More informationRandomization in the First Hitting Time Problem
Randomization in the First Hitting Time Problem Ken Jacson Alexander Kreinin Wanhe Zhang February 1, 29 Abstract In this paper we consider the following inverse problem for the first hitting time distribution:
More information0.1. Linear transformations
Suggestions for midterm review #3 The repetitoria are usually not complete; I am merely bringing up the points that many people didn t now on the recitations Linear transformations The following mostly
More informationProperties of a Taylor Polynomial
3.4.4: Still Better Approximations: Taylor Polynomials Properties of a Taylor Polynomial Constant: f (x) f (a) Linear: f (x) f (a) + f (a)(x a) Quadratic: f (x) f (a) + f (a)(x a) + 1 2 f (a)(x a) 2 3.4.4:
More informationThe Probabilistic Method
Lecture 3: Tail bounds, Probabilistic Method Today we will see what is nown as the probabilistic method for showing the existence of combinatorial objects. We also review basic concentration inequalities.
More information6.1 Moment Generating and Characteristic Functions
Chapter 6 Limit Theorems The power statistics can mostly be seen when there is a large collection of data points and we are interested in understanding the macro state of the system, e.g., the average,
More information9 Instantons in Quantum Mechanics
9 Instantons in Quantum Mechanics 9.1 General discussion It has already been briefly mentioned that in quantum mechanics certain aspects of a problem can be overlooed in a perturbative treatment. One example
More informationWKB Approximation in 3D
1 WKB Approximation in 3D We see solutions ψr of the stationary Schrodinger equations for a spinless particle of energy E: 2 2m 2 ψ + V rψ = Eψ At rst, we just rewrite the Schrodinger equation in the following
More informationRETRACTED On construction of a complex finite Jacobi matrix from two spectra
Electronic Journal of Linear Algebra Volume 26 Volume 26 (203) Article 8 203 On construction of a complex finite Jacobi matrix from two spectra Gusein Sh. Guseinov guseinov@ati.edu.tr Follow this and additional
More informationBohr & Wheeler Fission Theory Calculation 4 March 2009
Bohr & Wheeler Fission Theory Calculation 4 March 9 () Introduction The goal here is to reproduce the calculation of the limiting Z /A against spontaneous fission Z A lim a S. (.) a C as first done by
More information2 Generating Functions
2 Generating Functions In this part of the course, we re going to introduce algebraic methods for counting and proving combinatorial identities. This is often greatly advantageous over the method of finding
More informationMath 200 University of Connecticut
RELATIVISTIC ADDITION AND REAL ADDITION KEITH CONRAD Math 200 University of Connecticut Date: Aug. 31, 2005. RELATIVISTIC ADDITION AND REAL ADDITION 1 1. Introduction For three particles P, Q, R travelling
More informationA second look at waves
A second loo at waves ravelling waves A first loo at Amplitude Modulation (AM) Stationary and reflected waves Lossy waves: dispersion & evanescence I thin this is the MOS IMPORAN of my eight lectures,
More information1. Kinematics, cross-sections etc
1. Kinematics, cross-sections etc A study of kinematics is of great importance to any experiment on particle scattering. It is necessary to interpret your measurements, but at an earlier stage to determine
More informationCLOSE-TO-CLEAN REGULARIZATION RELATES
Worshop trac - ICLR 016 CLOSE-TO-CLEAN REGULARIZATION RELATES VIRTUAL ADVERSARIAL TRAINING, LADDER NETWORKS AND OTHERS Mudassar Abbas, Jyri Kivinen, Tapani Raio Department of Computer Science, School of
More information33. SOLVING LINEAR INEQUALITIES IN ONE VARIABLE
get the complete book: http://wwwonemathematicalcatorg/getfulltextfullbookhtm 33 SOLVING LINEAR INEQUALITIES IN ONE VARIABLE linear inequalities in one variable DEFINITION linear inequality in one variable
More informationFactor Substitution, Income Distribution, and Growth in a Generalized Neoclassical Model
University of Heidelberg Department of Economics Discussion Paper Series No. 453 Factor Substitution, Income Distribution, and Growth in a Generalized Neoclassical Model Andreas Irmen and Rainer Klump
More informationPositive and null recurrent-branching Process
December 15, 2011 In last discussion we studied the transience and recurrence of Markov chains There are 2 other closely related issues about Markov chains that we address Is there an invariant distribution?
More information10-725/36-725: Convex Optimization Spring Lecture 21: April 6
10-725/36-725: Conve Optimization Spring 2015 Lecturer: Ryan Tibshirani Lecture 21: April 6 Scribes: Chiqun Zhang, Hanqi Cheng, Waleed Ammar Note: LaTeX template courtesy of UC Berkeley EECS dept. Disclaimer:
More informationConstrained optimization.
ams/econ 11b supplementary notes ucsc Constrained optimization. c 2016, Yonatan Katznelson 1. Constraints In many of the optimization problems that arise in economics, there are restrictions on the values
More informationSign of derivative test: when does a function decreases or increases:
Sign of derivative test: when does a function decreases or increases: If for all then is increasing on. If for all then is decreasing on. If then the function is not increasing or decreasing at. We defined
More informationInequalities for the perimeter of an ellipse. G.J.O. Jameson, Math. Gazette 98 (2014)
Inequalities for the perimeter of an ellipse G.J.O. Jameson, Math. Gazette 98 (14) integral The perimeter of the ellipse x /a + y /b = 1 is 4J(a, b), where J(a, b) is the elliptic J(a, b) = (a cos θ +
More informationLecture Notes on Compositional Reasoning
15-414: Bug Catching: Automated Program Verification Lecture Notes on Compositional Reasoning Matt Fredrikson Ruben Martins Carnegie Mellon University Lecture 4 1 Introduction This lecture will focus on
More informationPh.D. Katarína Bellová Page 1 Mathematics 1 (10-PHY-BIPMA1) EXAM SOLUTIONS, 20 February 2018
Ph.D. Katarína Bellová Page Mathematics 0-PHY-BIPMA) EXAM SOLUTIONS, 0 February 08 Problem [4 points]: For which positive integers n does the following inequality hold? n! 3 n Solution: Trying first few
More informationGenerating Function Notes , Fall 2005, Prof. Peter Shor
Counting Change Generating Function Notes 80, Fall 00, Prof Peter Shor In this lecture, I m going to talk about generating functions We ve already seen an example of generating functions Recall when we
More informationLine adsorption in a mean-field density functional model
Chemistry Physical & Theoretical Chemistry fields Okayama University Year 2006 Line adsorption in a mean-field density functional model Kenichiro Koga Benjamin Widom Okayama University, koga@cc.okayama-u.ac.jp
More information1 Matrices and Systems of Linear Equations. a 1n a 2n
March 31, 2013 16-1 16. Systems of Linear Equations 1 Matrices and Systems of Linear Equations An m n matrix is an array A = (a ij ) of the form a 11 a 21 a m1 a 1n a 2n... a mn where each a ij is a real
More informationPhysics 411 Lecture 8. Parametrized Motion. Lecture 8. Physics 411 Classical Mechanics II
Physics 411 Lecture 8 Parametrized Motion Lecture 8 Physics 411 Classical Mechanics II September 14th 2007 We have our fancy new derivative, but what to do with it? In particular, how can we interpret
More informationThe line element for the hyperbolic plane was given in problem 12, of chapter 8 in Hartle
Physics 4445 Solution for homework Fall 20 Cornell University (46 points) I. HARTLE CHAPTER 2, PROBLEM (8 POINTS) The line element for the hyperbolic plane was given in problem 2, of chapter 8 in Hartle
More informationThe Inverse Function Theorem via Newton s Method. Michael Taylor
The Inverse Function Theorem via Newton s Method Michael Taylor We aim to prove the following result, known as the inverse function theorem. Theorem 1. Let F be a C k map (k 1) from a neighborhood of p
More informationSlope Fields: Graphing Solutions Without the Solutions
8 Slope Fields: Graphing Solutions Without the Solutions Up to now, our efforts have been directed mainly towards finding formulas or equations describing solutions to given differential equations. Then,
More informationNonlinear Regression
Nonlinear Regression 28.11.2012 Goals of Today s Lecture Understand the difference between linear and nonlinear regression models. See that not all functions are linearizable. Get an understanding of the
More informationNOTES ON DIOPHANTINE APPROXIMATION
NOTES ON DIOPHANTINE APPROXIMATION Jan-Hendrik Evertse January 29, 200 9 p-adic Numbers Literature: N. Koblitz, p-adic Numbers, p-adic Analysis, and Zeta-Functions, 2nd edition, Graduate Texts in Mathematics
More informationMath Precalculus I University of Hawai i at Mānoa Spring
Math 135 - Precalculus I University of Hawai i at Mānoa Spring - 2014 Created for Math 135, Spring 2008 by Lukasz Grabarek and Michael Joyce Send comments and corrections to lukasz@math.hawaii.edu Contents
More information1 Matrices and Systems of Linear Equations
March 3, 203 6-6. Systems of Linear Equations Matrices and Systems of Linear Equations An m n matrix is an array A = a ij of the form a a n a 2 a 2n... a m a mn where each a ij is a real or complex number.
More information8.324 Quantum Field Theory II. Problem Set 7 Solutions. d λ dλ
8.4 Quantum Field Theory II Problem Set 7 Solutions. a Changing the variable λ = λλ we can find β using the chain rule of differentiation. We find, β λ = µ d λ dµ = µ dλ d λ dµ dλ = βλ d λ dλ. Thus β transforms
More informationChapter 1 Review of Equations and Inequalities
Chapter 1 Review of Equations and Inequalities Part I Review of Basic Equations Recall that an equation is an expression with an equal sign in the middle. Also recall that, if a question asks you to solve
More informationStochastic Shortest Path Problems
Chapter 8 Stochastic Shortest Path Problems 1 In this chapter, we study a stochastic version of the shortest path problem of chapter 2, where only probabilities of transitions along different arcs can
More informationRelevant sections from AMATH 351 Course Notes (Wainwright): 1.3 Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls): 1.1.
Lecture 8 Qualitative Behaviour of Solutions to ODEs Relevant sections from AMATH 351 Course Notes (Wainwright): 1.3 Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls): 1.1.1 The last few
More informationKUMMER S LEMMA KEITH CONRAD
KUMMER S LEMMA KEITH CONRAD Let p be an odd prime and ζ ζ p be a primitive pth root of unity In the ring Z[ζ], the pth power of every element is congruent to a rational integer mod p, since (c 0 + c 1
More information29. GREATEST COMMON FACTOR
29. GREATEST COMMON FACTOR Don t ever forget what factoring is all about! greatest common factor a motivating example: cutting three boards of different lengths into same-length pieces solving the problem:
More informationElementary Linear Transform Theory
6 Whether they are spaces of arrows in space, functions or even matrices, vector spaces quicly become boring if we don t do things with their elements move them around, differentiate or integrate them,
More information6.5 Systems of Inequalities
6.5 Systems of Inequalities Linear Inequalities in Two Variables: A linear inequality in two variables is an inequality that can be written in the general form Ax + By < C, where A, B, and C are real numbers
More informationMini-Course on Limits and Sequences. Peter Kwadwo Asante. B.S., Kwame Nkrumah University of Science and Technology, Ghana, 2014 A REPORT
Mini-Course on Limits and Sequences by Peter Kwadwo Asante B.S., Kwame Nkrumah University of Science and Technology, Ghana, 204 A REPORT submitted in partial fulfillment of the requirements for the degree
More informationn if n is even. f (n)=
6 2. PROBABILITY 4. Countable and uncountable Definition 32. An set Ω is said to be finite if there is an n N and a bijection from Ω onto [n]. An infinite set Ω is said to be countable if there is a bijection
More informationSOME OBSERVATIONS ON THE CRYPTOGRAPHIC HASH FUNCTIONS
SOME OBSERVATIONS ON THE CRYPTOGRAPHIC HASH FUNCTIONS by Lavinia Ciungu Abstract In this paper we will mae a discussion on the conditions when a strongly collisionfree hash function is also one-way hash
More informationLesson 3: Solving Equations A Balancing Act
Opening Exercise Let s look back at the puzzle in Lesson 1 with the t-shape and the 100-chart. Jennie came up with a sum of 380 and through the lesson we found that the expression to represent the sum
More informationExercises for Chap. 2
1 Exercises for Chap..1 For the Riemann-Liouville fractional integral of the first kind of order α, 1,(a,x) f, evaluate the fractional integral if (i): f(t)= tν, (ii): f(t)= (t c) ν for some constant c,
More informationQuasi-Stationary Distributions in Linear Birth and Death Processes
Quasi-Stationary Distributions in Linear Birth and Death Processes Wenbo Liu & Hanjun Zhang School of Mathematics and Computational Science, Xiangtan University Hunan 411105, China E-mail: liuwenbo10111011@yahoo.com.cn
More informationLimits and Continuity
Chapter Limits and Continuity. Limits of Sequences.. The Concept of Limit and Its Properties A sequence { } is an ordered infinite list x,x,...,,... The n-th term of the sequence is, and n is the index
More informationMeasures. 1 Introduction. These preliminary lecture notes are partly based on textbooks by Athreya and Lahiri, Capinski and Kopp, and Folland.
Measures These preliminary lecture notes are partly based on textbooks by Athreya and Lahiri, Capinski and Kopp, and Folland. 1 Introduction Our motivation for studying measure theory is to lay a foundation
More informationAPPROXIMATING CONTINUOUS FUNCTIONS: WEIERSTRASS, BERNSTEIN, AND RUNGE
APPROXIMATING CONTINUOUS FUNCTIONS: WEIERSTRASS, BERNSTEIN, AND RUNGE WILLIE WAI-YEUNG WONG. Introduction This set of notes is meant to describe some aspects of polynomial approximations to continuous
More information1 I (x)) 1/2 I. A fairly immediate corollary of the techniques discussed in the last lecture is Theorem 1.1. For all 1 < p <
1. Lecture 4 1.1. Square functions, paraproducts, Khintchin s inequality. The dyadic Littlewood Paley square function of a function f is defined as Sf(x) := ( f, h 2 1 (x)) 1/2 where the summation goes
More informationExact formulae for the prime counting function
Notes on Number Theory and Discrete Mathematics Vol. 19, 013, No. 4, 77 85 Exact formulae for the prime counting function Mladen Vassilev Missana 5 V. Hugo Str, 114 Sofia, Bulgaria e-mail: missana@abv.bg
More informationTHE P-ADIC NUMBERS AND FINITE FIELD EXTENSIONS OF Q p
THE P-ADIC NUMBERS AND FINITE FIELD EXTENSIONS OF Q p EVAN TURNER Abstract. This paper will focus on the p-adic numbers and their properties. First, we will examine the p-adic norm and look at some of
More informationComputational Techniques Prof. Dr. Niket Kaisare Department of Chemical Engineering Indian Institute of Technology, Madras
Computational Techniques Prof. Dr. Niket Kaisare Department of Chemical Engineering Indian Institute of Technology, Madras Module No. # 07 Lecture No. # 04 Ordinary Differential Equations (Initial Value
More informationAbel Summation MOP 2007, Black Group
Abel Summation MOP 007, Blac Group Zachary Abel June 5, 007 This lecture focuses on the Abel Summation formula, which is most often useful as a way to tae advantage of unusual given conditions such as
More informationLecture 23: November 19
10-725/36-725: Conve Optimization Fall 2018 Lecturer: Ryan Tibshirani Lecture 23: November 19 Scribes: Charvi Rastogi, George Stoica, Shuo Li Charvi Rastogi: 23.1-23.4.2, George Stoica: 23.4.3-23.8, Shuo
More informationRiemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E,
Tel Aviv University, 26 Analysis-III 9 9 Improper integral 9a Introduction....................... 9 9b Positive integrands................... 9c Special functions gamma and beta......... 4 9d Change of
More informationResonance and response
Chapter 2 Resonance and response Last updated September 20, 2008 In this section of the course we begin with a very simple system a mass hanging from a spring and see how some remarkable ideas emerge.
More information6. Age structure. for a, t IR +, subject to the boundary condition. (6.3) p(0; t) = and to the initial condition
6. Age structure In this section we introduce a dependence of the force of infection upon the chronological age of individuals participating in the epidemic. Age has been recognized as an important factor
More informationGlobal Well-posedness of the Incompressible Magnetohydrodynamics
Global Well-posedness of the Incompressible Magnetohydrodynamics arxiv:165.439v1 [math.ap] May 16 Yuan Cai Zhen Lei Abstract This paper studies the Cauchy problem of the incompressible magnetohydrodynamic
More informationMultiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions
Multiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions March 6, 2013 Contents 1 Wea second variation 2 1.1 Formulas for variation........................
More information