Symmetry Methods for Differential and Difference Equations. Peter Hydon

Size: px
Start display at page:

Download "Symmetry Methods for Differential and Difference Equations. Peter Hydon"

Transcription

1 Lecture 2: How to find Lie symmetries Symmetry Methods for Differential and Difference Equations Peter Hydon University of Kent

2 Outline 1 Reduction of order for ODEs and O Es 2 The infinitesimal generator 3 Prolongation formulae for scalar ODEs 4 How to solve the LSC 5 Beyond point symmetries 6 The LSC for a given O E 7 Inherited symmetries 8 Summary: the main results in Lecture 2

3 Reduction of order for ODEs and O Es Reduction of order for ODEs A given ODE of order p may be written in canonical coordinates as s (p) = Ω ( r, ṡ, s,..., s (p 1)), s (k) = dk s dr k. (1) The variable s is missing, because the ODE has the Lie symmetries (ˆr, ŝ) = (r, s + ε), which give ŝ (k) = s (k) for each k 1. In other words, the ODE can be written entirely in terms of differential functions that are invariant under the Lie symmetries. Consequently, (1) reduces to v (p 1) = Ω ( r, v, v,..., v (p 2)), where v = ṡ. (2) If (2) can be solved, the solution of (1) is s = v(r; c 1,... c p 1 ) dr + c p.

4 Reduction of order for ODEs and O Es Note Sometimes, it is more convenient to choose v to be a different function of r and ṡ. Any choice will do, provided that v depends locally diffeomorphically on ṡ. This is to ensure that (r, ṡ) can be obtained from (r, v) and vice versa, and that the local smooth structure is preserved. Whichever choice of v is made, the reduced ODE will always be of the form v (p 1) = Ω ( r, v, v,..., v (p 2)).

5 Reduction of order for ODEs and O Es Example The ODE y = 1 xy 2 (3) has Lie symmetries with ξ(x, y) = x 2, η(x, y) = xy. The canonical coordinates ( y (r, s) = x, 1 ) ( 1 give (ṡ, s) = x xy y, x 3 y ) (xy y) 3. With v = ṡ, the ODE reduces to v = v 3 r 2, which happens to be separable. Therefore, the general solution of (3) can be obtained by integration.

6 Reduction of order for ODEs and O Es Reduction of order for O Es A given O E of order p may be written in canonical coordinates as s p s p 1 = Ω ( n, s 1 s, s 2 s 1,..., s p 1 s p 2 ). (4) The variable s is missing, because the ODE has the Lie symmetries (ˆn, ŝ) = (n, s + ε), so ŝ k+1 ŝ k = s k+1 s k for each k 0. In other words, the O E can be written entirely in terms of functions that are invariant under the Lie symmetries. Consequently, (1) reduces to v p 1 = Ω ( n, v, v 1,..., v p 2 ), where v = s1 s. (5) If (5) can be solved, the solution of (4) is s = v(n; c 1,... c p 1 ) + c p.

7 Reduction of order for ODEs and O Es Example The O E u 2 = u 1(uu 1 + u 1 u) (u 1 ) 2 + u 1 u has Lie symmetries whose characteristic is Q(n, u) = u 2. Thus s 2 = (s 1) 2 ss 1 + s s 1 s + 1, where s = 1 u. With the substitution v = s 1 s, this reduces to v 1 = v v + 1, whose general solution is v = 2 { 1 + c 1 ( 1) n} / { (c 1 ) 2 1 }.Then s = v(n) + c 2 = 2n + c { 1 1 ( 1) n } (c 1 ) 2 + c 2, 1 which is easily inverted to obtain u.

8 The infinitesimal generator The infinitesimal generator Question How does an arbitrary locally-smooth function F (x, y) vary along the orbit of a 1-parameter local Lie group? Answer Replace (x, y) by (ˆx, ŷ) and expand in powers of ε: F (ˆx, ŷ) = F (x, y) + ε { ξ(x, y)f,x + η(x, y)f,y } + O(ε 2 ), where the subscripts indicate the partial derivative with respect to each variable that follows the comma. Equivalently, F (ˆx, ŷ) = F (x, y) + εx {F (x, y)} + O(ε 2 ), where X = ξ(x, y) x + η(x, y) y is the infinitesimal generator of the local Lie group. (We use the shorthand a for a here and from now on.)

9 The infinitesimal generator The infinitesimal generator turns up in the formal solution of the equations that define the tangent vector to the orbit at (ˆx, ŷ): dˆx dε = ξ(ˆx, ŷ), dŷ The formal solution is dε = η(ˆx, ŷ), (ˆx, ŷ) ε=0 = (x, y). ˆx = exp{εx }x, ŷ = exp{εx }y, where exp{εx } = k=0 ε k k! X k. The infinitesimal generator determines (ˆx, ŷ) locally, to all orders in ε. Indeed, wherever the above series converges, F (ˆx, ŷ) = exp{εx }F (x, y).

10 The infinitesimal generator Theorem F (x, y) is invariant if and only if X {F (x, y)} = 0. More generally, if a 1-parameter local Lie group acts on any space with coordinates x α, such that the infinitesimal generator is ˆx α = x α + ε ζ α (x) + O(ε 2 ), X = ζ α x α (summed over α). A locally-smooth function F (x) is invariant if and only if X {F (x)} = 0.

11 The infinitesimal generator The change-of-variables formula is straightforward for infinitesimal generators. Given local coordinates (a, b) on the (x, y)-plane, X = (Xa) a + (Xb) b. In particular, canonical coordinates (r, s) satisfy Xr = 0, Xs = 1, so X = s. [By the above theorem, r is invariant under the action of the 1-parameter Lie group.] More generally, in a space with coordinates x α, X = (Xx α ) x α.

12 Prolongation formulae for scalar ODEs The symmetry condition for a p th -order ODE, is y (p) = ω(x, y, y,..., y (p 1) ), (6) ŷ (p) = ω (ˆx, ŷ, ŷ (1),..., ŷ (p 1)) when (6) holds. (Remember: solutions must be mapped to solutions.) This condition involves not only ˆx and ŷ, but also the derivatives ŷ (k) for k from 1 to p. Consequently, we must prolong the infinitesimal generator to act on the p th jet space, with local coordinates (x, y, y,..., y (p) ). To deal with all p at once, write X = ξ x + η y + η (1) y + η (2) y +. We will use (6) to evaluate the output on solutions of the ODE.

13 Prolongation formulae for scalar ODEs The symmetry condition gives, to first-order in ε, ( ) X y (p) ω(x, y,..., y (p 1) ) = 0 when y (p) = ω(x, y,..., y (p 1) ). Explicitly, this linearized symmetry condition (LSC) amounts to η (p) ξ ω,x η ω,y η (1) ω,y η (p 1) ω,y (p 1) = 0, (7) when (6) holds. [Here ω is shorthand for ω(x, y,..., y (p 1) ).] The LSC holds as an identity in x, y,..., y (p 2) and y (p 1). Example For a given first-order ODE, y = ω(x, y), the LSC is η (1) ξ ω,x η ω,y = 0 when y = ω(x, y); this holds as an identity in x and y (because the ODE eliminates derivatives of y).

14 Prolongation formulae for scalar ODEs The coefficients η (k) are defined by the expansion ŷ (k) = y (k) + εη (k) + O(ε 2 ). To calculate η (1), use the total derivative operator, to expand ŷ (1) as follows: D x = x + y y + y y + y y, ŷ (1) = D xŷ D x ˆx = y + εd x η + O(ε 2 ) 1 + εd x ξ + O(ε 2 ) = y +ε ( D x η y D x ξ ) +O(ε 2 ). Similarly, each η (k) is calculated recursively from ŷ (k) = y (k) + εd x η (k 1) + O(ε 2 ) 1 + εd x ξ + O(ε 2 ) = y (k) +ε ( D x η (k 1) y (k) D x ξ ) +O(ε 2 ).

15 Prolongation formulae for scalar ODEs Examples Assuming that ξ and η depend only on x and y, the prolongation formulae yield: η (1) = η,x + ( ) η,y ξ,x y ξ,y y 2, η (2) = η,xx + (2η,xy ξ,xx )y + (η,yy 2ξ,xy )y 2 ξ,yy y 3 + {η,y 2ξ,x 3ξ,y y }y, η (3) = η,xxx + (3η,xxy ξ,xxx )y + 3(η,xyy ξ,xxy )y 2 + (η,yyy 3ξ,xyy )y 3 ξ,yyy y 4 + 3{η,xy ξ,xx + (η,yy 3ξ,xy )y 2ξ,yy y 2 }y 3ξ,y y 2 + {η,y 3ξ,x 4ξ,y y }y.

16 How to solve the LSC How to solve the LSC We are now able to write out the LSC as a set of partial differential equations (PDEs) for ξ(x, y) and η(x, y). Usually, for ODEs of order p 2, one can solve these equations. This is because the LSC can be split according to the dependence of each term on y,..., y (p 1). The easiest way to understand how this works in practice is to study some examples.

17 How to solve the LSC Example The LSC for the simplest second-order ODE, y =0, is η (2) = 0 when y = 0, which amounts to η,xx + (2η,xy ξ,xx )y + (η,yy 2ξ,xy )y 2 ξ,yy y 3 = 0. As ξ and η are independent of y, the LSC splits into η xx = 0, 2η xy ξ xx = 0, η yy 2ξ xy = 0, ξ yy = 0. (8) The general solution of the last two of (8) is ξ(x, y) = A(x)y + B(x), η(x, y) = A (x)y 2 + C(x)y + D(x), where A, B, C and D are arbitrary locally-smooth functions.

18 How to solve the LSC The remaining equations in (8) amount to A (x)y 2 +C (x)y +D (x) = 0, 3A (x)y +2C (x) B (x) = 0, which are split (by comparing powers of y) into A (x) = 0, C (x) = 0, D (x) = 0, B (x) = 2C (x). All that remains is to solve this simple system of ODEs (in terms of some arbitrary constants, c i ) and reconstruct ξ and η: ξ(x, y) = c 1 + c 3 x + c 5 y + c 7 x 2 + c 8 xy, η(x, y) = c 2 + c 4 y + c 6 x + c 7 xy + c 8 y 2. Therefore, every symmetry generator X is a linear combination of X 1 = x, X 2 = y, X 3 = x x, X 4 = y y, X 5 = y x, X 6 = x y, X 7 = x 2 x + xy y, X 8 = xy x + y 2 y.

19 How to solve the LSC Example The ODE y = y 3 arises in models of thin-film flow. The LSC is which amounts to η (3) + 3y 4 η = 0 when y = y 3, 0 = η,xxx + (3η,xxy ξ,xxx )y + 3(η,xyy ξ,xxy )y 2 + (η,yyy 3ξ,xyy )y 3 ξ,yyy y 4 + 3{η,xy ξ,xx + (η,yy 3ξ,xy )y 2ξ,yy y 2 }y 3ξ,y y 2 + {η,y 3ξ,x 4ξ,y y }y 3 + 3y 4 η. Although this looks horrible, only the red term is quadratic in y. Consequently, ξ = A(x).

20 How to solve the LSC Taking this into account, the LSC reduces to 0 = η,xxx + ( 3η,xxy A (x) ) y + 3η,xyy y 2 + η,yyy y 3 +3{η,xy A (x) + η,yy y }y + {η,y 3a (x)}y 3 + 3y 4 η. Only the blue terms are linear in y, so (equating powers of y in such terms) η,xy A (x) = 0, η,yy = 0 η = ( A (x) + c 1 ) y + B(x). Substitute this into the remaining terms of the LSC and compare powers of y, then y, to obtain the general solution of the LSC: ξ(x, y) = 4c 1 x + c 2, η(x, y) = 3c 1 y. Thus every symmetry generator X is a linear combination of X 1 = 4x x 3y y, X 2 = x.

21 How to solve the LSC For a given first-order ODE, the LSC has infinitely many linearly independent solutions. However, the LSC cannot be split with respect to derivatives of y, so one cannot find its general solution. The best that can be done is to restrict ξ and η and hope that this yields some symmetries. Example The LSC for the ODE y = (1 y 2 )/(xy) + 1 is ( 1 y 2 0 =η,x + xy ( 1 y 2 + x 2 y ( (η,y ) 1 y 2 + 1) ξ,x xy ( 1 + y 2 ) ξ + xy 2 ) η. + 1) 2 ξ,y Let us see whether there are any symmetry generators for which ξ = α(x), η = β(x)y + γ(x). (This is a common restriction.)

22 How to solve the LSC The LSC amounts to ( 1 y 0 = β y+γ +(β α 2 ) xy ) + 1 +α ( 1 y 2 x 2 y ) +(βy +γ) ( 1 + y 2 xy 2 ). Now one can split the (restricted) LSC in powers of y, to get (after some simplification) γ = 0, β α = 0, α + α/x = 0. There are no further constraints, so α = c 1 x 1 and β = c 1 x 2. Consequently, every generator of the restricted form is a multiple of X 1 = x 1 x x 2 y y. Exercise Use this symmetry generator to construct canonical coordinates and solve the ODE.

23 Beyond point symmetries Beyond point symmetries A nice version of the prolongation formula (given ξ and η) is η (k) = D k x Q + y (k+1) ξ, where Q = η y ξ is the characteristic; this is proved recursively. It leads to a very useful splitting of the infinitesimal generator: X = ξd x + Q y + (D x Q) y + (D 2 x Q) y +. The red term generates trivial Lie symmetries (for arbitrary ξ), because its restriction to solutions maps points along solution curves. So X is equivalent (up to a trivial symmetry generator) to X = Q y + (D x Q) y + (D 2 x Q) y +.

24 Beyond point symmetries X X (ˆx, ŷ) y (x, y) x ξd x ε Non-invariant solution curves: orbits cross these transversely For first-order ODEs, the geometrical picture is shown above: the generator X splits into a trivial tangential component, ξd x, and a nontrivial vertical component, X. Here, vertical means that there is no component in the direction of the independent variable x. The same picture holds (in the appropriate jet space) for higher order ODEs. The LSCs with X and X give the same characteristics.

25 Beyond point symmetries So far, we have restricted attention to point symmetries, which act as local diffeomorphisms on the (x, y)-plane. More symmetries are available, which act as local diffeomorphisms on various jet spaces. Given an ODE of order p, the types of symmetries are as follows. Point symmetries have a characteristic that depends only on x, y and y, and is linear in y, so Q = η(x, y) ξ(x, y)y. Contact symmetries have a characteristic Q(x, y, y ), with ξ = Q,y, η = Q y Q,y, η (1) = Q,x + y Q,y. Dynamical symmetries have a characteristic that depends on some or all of x, y, y,... y (p 1). It is usual to set ξ = 0, so that ˆx = x.

26 Beyond point symmetries To find any of these characteristics, first pick a suitable restriction on Q. Then solve the LSC in the form ( ) X y (p) ω(x, y,..., y (p 1) ) = 0 on solutions, using y (p) = ω(x, y,..., y (p 1) ) to eliminate y (p+k) for all k 0. Example Every characteristic of Lie contact symmetries for is a linear combination of y = x(x 1)y 3 2xy 2 + y Q 1 = 1, Q 2 = x, (point symmetry characteristics) Q 3 = e y, Q 4 = (xy y x)e y, Q 5 = e y xy +y +x. These are obtained by splitting the LSC in powers of y.

27 Beyond point symmetries Just as for first-order ODEs, there are infinitely many dynamical symmetry characteristics Q(x, y,..., y (p 1) ). To find any of them by splitting the LSC, one must make a further restriction. The most common restriction is to look for characteristics that are independent of y (p 1). All nontrivial symmetries, dynamical or otherwise, act by changing the arbitrary constants in the general solution of the ODE. (One can regard them as local point transformations on the p-dimensional space of first integrals.) The same idea applies to O Es.

28 The LSC for a given O E The LSC for a given O E The LSC for a given O E is written using the forward shift operator, S n, which replaces n by n + 1 in all functions of n: S n {f (n)} = f (n + 1), S n u = u 1, S n u i = u i+1. The last two equations can be combined by using u 0 as another notation for u; we will use this alternative whenever it is helpful. For simplicity, we will assume that u is defined for all n Z, though what follows can be adapted to any domain. For later use, we also introduce the identity operator I, which maps every function of n to itself.

29 The LSC for a given O E For scalar O Es, Lie symmetries are of the form ˆn = n, û = u + εq + O(ε 2 ). In particular, Lie point symmetries have û (and thus, Q) dependent on n and u only. The prolongation formula is obtained by replacing n by n + k: û k = u k + ε S k nq + O(ε 2 ). Therefore, the infinitesimal generator (prolonged to all orders) is X = k Z ( S k nq ) u k.

30 The LSC for a given O E The symmetry condition for is u p = ω(n, u,..., u p 1 ) (9) û p = ω(ˆn, û,..., û p 1 ) when (9) holds. Just as for ODEs, the LSC for Lie symmetries is obtained by expanding the symmetry condition to first order in ε: X { u p ω(n, u,..., u p 1 ) } = 0 when (9) holds. Using ω as shorthand for ω(n, u,..., u p 1 ), the LSC amounts to S p nq ω,1 Q ω,2 S n Q ω,p S p 1 n Q = 0 when (9) holds, where the subscript, i denotes the partial derivative with respect to the i th continuous argument. For example, ω,p = ω/ u (p 1).

31 The LSC for a given O E For most nonlinear O Es of order p > 1, one can find all Lie point symmetries from the LSC, without having to restrict Q. To see how, let us examine the LSC for a given second-order O E, u 2 = ω(n, u, u 1 ), (10) which is a functional equation (not an ODE or PDE) for Q: Q(n + 2, ω) ω,2 Q(n + 1, u 1 ) ω,1 Q(n, u) = 0. (11) Although functional equations are generally hard to solve, Lie symmetries are diffeomorphisms. Consequently, Q is a locally smooth function of its continuous argument and so the LSC can be solved by the method of differential elimination. Before discussing this method in general, we begin with a simple example.

32 The LSC for a given O E Example The LSC for Lie point symmetries of u 2 = (u 1 ) 2 /u is Q ( n + 2, (u 1 ) 2 /u ) 2(u 1 /u) Q(n + 1, u 1 ) + (u 1 /u) 2 Q(n, u) = 0. In each term, Q depends on only one continuous variable (and one discrete variable). We could differentiate the LSC in order to eliminate that term. In particular, { u2 (u 1 ) 2 u u 2u 1 u 1 } Q ( n + 2, (u 1 ) 2 /u ) = 0. By applying the differential operator in braces to the LSC, we eliminate the first term in the LSC and obtain Q (n + 1, u 1 ) 1 Q(n + 1, u 1 ) Q (n, u) + 1 Q(n, u) = 0. (12) u 1 u This has more terms than the LSC does, but Q and Q only take two pairs of arguments, instead of three.

33 The LSC for a given O E Continue the process, by differentiating (12) with respect to u, in order to eliminate Q (n + 1, u 1 ). This gives [ Q (n, u) + 1 u Q(n, u) ] = 0. (13) Having found a (parametrized) differential equation for Q, we can start to solve the system. Integrate (13) once, to obtain Q (n, u) 1 Q(n, u) = α(n), (14) u where α(n) is to be determined. Now substitute (14) into the intermediate equation (12), which gives the O E α 1 α = 0. Consequently α = c 1, so the general solution of (12) is Q(n, u) = c 1 u ln u + β(n) u. (15)

34 The LSC for a given O E Finally, substitute (15) into the LSC to obtain (after cancellation) β 2 2β 1 + β = 0. The general solution of this linear O E is β(n) = c 2 n + c 3, and hence the general solution of the LSC is Q(n, u) = c 1 u ln u + c 2 nu + c 3 u. In other words, every characteristic is a linear combination of Q 1 = u ln u, Q 2 = nu, Q 3 = u.

35 The LSC for a given O E This example illustrates the general technique for obtaining Lie point of any difference equation of order p 2: 1 Write down the LSC. 2 Apply appropriate differential operators to reduce the number of unknown functions. 3 Having reached a differential equation, back-substitute and solve the resulting linear difference equations. 4 Iterate, if necessary. Exactly the same approach may be used to obtain dynamical symmetries, with Q dependent on n, u, u 1,..., u p 1, provided that Q is restricted is some way. For first-order O Es, it is necessary to assume that Q has a particular dependence on u. A common assumption is Q(n, u) = α(n) u 2 + β(n) u + γ(n). (16)

36 The LSC for a given O E Example The LSC for the O E u 1 = u/(1 + nu) is Q ( n + 1, u/(1 + nu) ) Q(n, u)/(1 + nu) 2 = 0. With the restriction (16), the LSC reduces to u 2 α 1 + u(1 + nu)β 1 + (1 + nu) 2 γ 1 u 2 α u β γ = 0, which we split into a system of O Es by comparing powers of u: u 2 terms : α 1 + n β 1 + n 2 γ 1 α = 0, u terms : β 1 + 2n γ 1 β = 0, other terms : γ 1 γ = 0.

37 The LSC for a given O E Systems in this triangular form are solved by starting at the bottom and working up. In this case, each of the scalar equations is easy to solve, so there is no need to restrict α, β or γ. We obtain γ = c 1, β = c 2 c 1 n(n 1), α = c 3 c 2 n(n 1)/2 + c 1 n 2 (n 1) 2 /4. Therefore every characteristic that is quadratic in u is of the form ( n(n 1) 2 ( ) n(n 1) Q(n, u) = c 1 1 u) + c 2 u 1 u + c 3 u Note For systems not in triangular form, one may need to combine the O Es to solve them. If necessary, restrict α, β and γ.

38 The LSC for a given O E Potential problems with differential elimination 1. Generally, each elimination increases the number of terms, a phenomenon known as expression swell. This makes it difficult to find dynamical symmetries of high-order O Es. 2. In view of the above, it is wise to use computer algebra to keep track of the calculations. However, beware of computer algebra solvers for differential and difference equations. They may overlook some solutions! 3. It may be difficult (or even impossible) to solve some of the linear difference equations in Step 3. However, this is unusual for nonlinear O Es.

39 Inherited symmetries Inherited symmetries For ODEs and O Es, the set of Lie symmetry characteristics that satisfy a particular restriction is a vector space. However, for the set of all Lie point symmetries, there is an extra structure. Theorem If X i and X j generate Lie symmetries (of any type), so does [X i, X j ] = X i X j X j X i. Corollary The space L of all Lie point symmetry generators is a Lie algebra under the commutator [, ]. Note The space of all dynamical symmetry generators is also a Lie algebra, but it is easier to find the set of all Lie point symmetries (when p 2) and to construct their canonical coordinates!

40 Inherited symmetries An R-dimensional Lie algebra L is solvable if the derived series, L (0) = L, L (i+1) = [L (i), L (i) ], stops at {0}. If L is solvable, one can use the derived series to construct a chain of subalgebras, {0} = L 0 L 1 L R = L, such that each L k is a k-dimensional ideal of L k+1. Given such a chain, a canonical basis for L satisfies L k = Span(X 1,..., X k ), for k = 1,..., R. Theorem Given an ODE or O E of order p 2, suppose that L is solvable and that X 1,..., X R is a canonical basis. When the ODE or O E is reduced using canonical coordinates corresponding to X 1, the reduced equation inherits the symmetries generated by X 2,..., X R.

41 Inherited symmetries Example has a canonical basis The Lie algebra L of point symmetry generators for u 3 = u 2 + (u 2 u 1 ) 2 u 1 u X 1 = u, X 2 = u u. Reduce the order by using X 1, whose canonical coordinate is s = u: v 2 = v 2 1 /v v = u 1 u. The reduced O E inherits the symmetries generated by X 2, as X 2 v = u v u + u v 1 = u 1 u = v. u 1 Therefore X 2 reappears as the scaling symmetry generator X 2 = v v, which enables one to carry out a second reduction of order.

42 Inherited symmetries By contrast, let us use X 2 first, with s = ln u and ṽ = u 1 /u for simplicity. (For the same reason as for ODEs, we can take ṽ to be a locally invertible function of s 1 s.) This reduction of order works (as it is guaranteed to do), yielding ṽ 2 = 1 + ṽ (ṽ 1 1) 2 ṽ 1 (ṽ 1). (17) However, the reduced O E (17) has no Lie point symmetries. What has happened to the symmetries generated by X 1? As before, apply X 1 (prolonged) to ṽ, to get X 1 ṽ = ṽ u + ṽ = 1 u 1 u u 1 u 2 = 1 ṽ. u This is not a function of ṽ only, so it cannot be a generator of Lie point symmetries for the reduced O E (17).

43 Summary: the main results in Lecture 2 Summary: the main results in Lecture 2 Given an ODE or O E of order p 2, canonical coordinates reduce the order by one. If the reduced equation can be solved, the original equation is solved by one more integration or summation. The infinitesimal generator X determines the local behaviour of Lie symmetries. A function is invariant if and only if it is in the kernel of X. The prolonged infinitesimal generator yields the LSC. For an ODE, this is a differential equation; for an O E, the LSC is a functional equation. By restricting the characteristic, Q, one can seek all Lie symmetries that satisfy the restriction. Point symmetries act on the space of independent and dependent variables; dynamical symmetries act on the space of first integrals (or arbitrary constants).

44 Summary: the main results in Lecture 2 Conclusion: There are substantial differences between O E and ODE symmetry methods!

Applications of Symmetries and Conservation Laws to the Study of Nonlinear Elasticity Equations

Applications of Symmetries and Conservation Laws to the Study of Nonlinear Elasticity Equations Applications of Symmetries and Conservation Laws to the Study of Nonlinear Elasticity Equations A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment of the Requirements

More information

First Integrals/Invariants & Symmetries for Autonomous Difference Equations

First Integrals/Invariants & Symmetries for Autonomous Difference Equations Proceedings of Institute of Mathematics of NAS of Ukraine 2004, Vol. 50, Part 3, 1253 1260 First Integrals/Invariants & Symmetries for Autonomous Difference Equations Leon ARRIOLA Department of Mathematical

More information

B.7 Lie Groups and Differential Equations

B.7 Lie Groups and Differential Equations 96 B.7. LIE GROUPS AND DIFFERENTIAL EQUATIONS B.7 Lie Groups and Differential Equations Peter J. Olver in Minneapolis, MN (U.S.A.) mailto:olver@ima.umn.edu The applications of Lie groups to solve differential

More information

Algorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations

Algorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations dmitry.lyakhov@kaust.edu.sa Symbolic Computations May 4, 2018 1 / 25 Algorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations Dmitry A. Lyakhov 1 1 Computational Sciences

More information

Classification of cubics up to affine transformations

Classification of cubics up to affine transformations Classification of cubics up to affine transformations Mehdi Nadjafikah and Ahmad-Reza Forough Abstract. Classification of cubics (that is, third order planar curves in the R 2 ) up to certain transformations

More information

Symmetry Methods for Differential Equations and Conservation Laws. Peter J. Olver University of Minnesota

Symmetry Methods for Differential Equations and Conservation Laws. Peter J. Olver University of Minnesota Symmetry Methods for Differential Equations and Conservation Laws Peter J. Olver University of Minnesota http://www.math.umn.edu/ olver Santiago, November, 2010 Symmetry Groups of Differential Equations

More information

Linearization of Two Dimensional Complex-Linearizable Systems of Second Order Ordinary Differential Equations

Linearization of Two Dimensional Complex-Linearizable Systems of Second Order Ordinary Differential Equations Applied Mathematical Sciences, Vol. 9, 2015, no. 58, 2889-2900 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.4121002 Linearization of Two Dimensional Complex-Linearizable Systems of

More information

2. Intersection Multiplicities

2. Intersection Multiplicities 2. Intersection Multiplicities 11 2. Intersection Multiplicities Let us start our study of curves by introducing the concept of intersection multiplicity, which will be central throughout these notes.

More information

z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.

z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables. Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These

More information

Fundamentals of Linear Algebra. Marcel B. Finan Arkansas Tech University c All Rights Reserved

Fundamentals of Linear Algebra. Marcel B. Finan Arkansas Tech University c All Rights Reserved Fundamentals of Linear Algebra Marcel B. Finan Arkansas Tech University c All Rights Reserved 2 PREFACE Linear algebra has evolved as a branch of mathematics with wide range of applications to the natural

More information

A Lie-Group Approach for Nonlinear Dynamic Systems Described by Implicit Ordinary Differential Equations

A Lie-Group Approach for Nonlinear Dynamic Systems Described by Implicit Ordinary Differential Equations A Lie-Group Approach for Nonlinear Dynamic Systems Described by Implicit Ordinary Differential Equations Kurt Schlacher, Andreas Kugi and Kurt Zehetleitner kurt.schlacher@jku.at kurt.zehetleitner@jku.at,

More information

1 The functional equation for ζ

1 The functional equation for ζ 18.785: Analytic Number Theory, MIT, spring 27 (K.S. Kedlaya) The functional equation for the Riemann zeta function In this unit, we establish the functional equation property for the Riemann zeta function,

More information

Oblique derivative problems for elliptic and parabolic equations, Lecture II

Oblique derivative problems for elliptic and parabolic equations, Lecture II of the for elliptic and parabolic equations, Lecture II Iowa State University July 22, 2011 of the 1 2 of the of the As a preliminary step in our further we now look at a special situation for elliptic.

More information

MA22S3 Summary Sheet: Ordinary Differential Equations

MA22S3 Summary Sheet: Ordinary Differential Equations MA22S3 Summary Sheet: Ordinary Differential Equations December 14, 2017 Kreyszig s textbook is a suitable guide for this part of the module. Contents 1 Terminology 1 2 First order separable 2 2.1 Separable

More information

Lie Theory of Differential Equations and Computer Algebra

Lie Theory of Differential Equations and Computer Algebra Seminar Sophus Lie 1 (1991) 83 91 Lie Theory of Differential Equations and Computer Algebra Günter Czichowski Introduction The aim of this contribution is to show the possibilities for solving ordinary

More information

On the Linearization of Second-Order Dif ferential and Dif ference Equations

On the Linearization of Second-Order Dif ferential and Dif ference Equations Symmetry, Integrability and Geometry: Methods and Applications Vol. (006), Paper 065, 15 pages On the Linearization of Second-Order Dif ferential and Dif ference Equations Vladimir DORODNITSYN Keldysh

More information

7 Planar systems of linear ODE

7 Planar systems of linear ODE 7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

More information

Conservation laws for the geodesic equations of the canonical connection on Lie groups in dimensions two and three

Conservation laws for the geodesic equations of the canonical connection on Lie groups in dimensions two and three Appl Math Inf Sci 7 No 1 311-318 (013) 311 Applied Mathematics & Information Sciences An International Journal Conservation laws for the geodesic equations of the canonical connection on Lie groups in

More information

Symmetry Preserving Numerical Methods via Moving Frames

Symmetry Preserving Numerical Methods via Moving Frames Symmetry Preserving Numerical Methods via Moving Frames Peter J. Olver University of Minnesota http://www.math.umn.edu/ olver = Pilwon Kim, Martin Welk Cambridge, June, 2007 Symmetry Preserving Numerical

More information

Linear Algebra. Min Yan

Linear Algebra. Min Yan Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................

More information

Relevant sections from AMATH 351 Course Notes (Wainwright): Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls):

Relevant sections from AMATH 351 Course Notes (Wainwright): Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls): Lecture 5 Series solutions to DEs Relevant sections from AMATH 35 Course Notes (Wainwright):.4. Relevant sections from AMATH 35 Course Notes (Poulin and Ingalls): 2.-2.3 As mentioned earlier in this course,

More information

Preliminaries Lectures. Dr. Abdulla Eid. Department of Mathematics MATHS 101: Calculus I

Preliminaries Lectures. Dr. Abdulla Eid. Department of Mathematics   MATHS 101: Calculus I Preliminaries 2 1 2 Lectures Department of Mathematics http://www.abdullaeid.net/maths101 MATHS 101: Calculus I (University of Bahrain) Prelim 1 / 35 Pre Calculus MATHS 101: Calculus MATHS 101 is all about

More information

Symmetries and reduction techniques for dissipative models

Symmetries and reduction techniques for dissipative models Symmetries and reduction techniques for dissipative models M. Ruggieri and A. Valenti Dipartimento di Matematica e Informatica Università di Catania viale A. Doria 6, 95125 Catania, Italy Fourth Workshop

More information

2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1

2tdt 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 information

Chapter 3. Second Order Linear PDEs

Chapter 3. Second Order Linear PDEs Chapter 3. Second Order Linear PDEs 3.1 Introduction The general class of second order linear PDEs are of the form: ax, y)u xx + bx, y)u xy + cx, y)u yy + dx, y)u x + ex, y)u y + f x, y)u = gx, y). 3.1)

More information

On universality of critical behaviour in Hamiltonian PDEs

On universality of critical behaviour in Hamiltonian PDEs Riemann - Hilbert Problems, Integrability and Asymptotics Trieste, September 23, 2005 On universality of critical behaviour in Hamiltonian PDEs Boris DUBROVIN SISSA (Trieste) 1 Main subject: Hamiltonian

More information

An Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations

An Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical

More information

Chapter 3 Second Order Linear Equations

Chapter 3 Second Order Linear Equations Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,

More information

Vector Spaces, Orthogonality, and Linear Least Squares

Vector Spaces, Orthogonality, and Linear Least Squares Week Vector Spaces, Orthogonality, and Linear Least Squares. Opening Remarks.. Visualizing Planes, Lines, and Solutions Consider the following system of linear equations from the opener for Week 9: χ χ

More information

CLASSIFICATIONS OF THE FLOWS OF LINEAR ODE

CLASSIFICATIONS OF THE FLOWS OF LINEAR ODE CLASSIFICATIONS OF THE FLOWS OF LINEAR ODE PETER ROBICHEAUX Abstract. The goal of this paper is to examine characterizations of linear differential equations. We define the flow of an equation and examine

More information

MATH 304 Linear Algebra Lecture 19: Least squares problems (continued). Norms and inner products.

MATH 304 Linear Algebra Lecture 19: Least squares problems (continued). Norms and inner products. MATH 304 Linear Algebra Lecture 19: Least squares problems (continued). Norms and inner products. Orthogonal projection Theorem 1 Let V be a subspace of R n. Then any vector x R n is uniquely represented

More information

Multiple 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 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

arxiv: v1 [math-ph] 7 Dec 2017

arxiv: v1 [math-ph] 7 Dec 2017 Lie group classification of first-order delay ordinary differential equations arxiv:1712.02581v1 [math-ph] 7 Dec 2017 Vladimir A. Dorodnitsyn, Roman Kozlov, Sergey V. Meleshko and Pavel Winternitz Keldysh

More information

THEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)

THEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.) 4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M

More information

Multiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015

Multiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015 Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction

More information

Math 302 Outcome Statements Winter 2013

Math 302 Outcome Statements Winter 2013 Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a

More information

Chapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey

Chapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix

More information

SOME HYNDON S GENERALIZATIONS TO STARRETT S METHOD OF SOLVING FIRST ORDER ODES BY LIE GROUP SYMMETRY. Z.M. Mwanzia, K.C. Sogomo

SOME HYNDON S GENERALIZATIONS TO STARRETT S METHOD OF SOLVING FIRST ORDER ODES BY LIE GROUP SYMMETRY. Z.M. Mwanzia, K.C. Sogomo SOME HYNDON S GENERALIZATIONS TO STARRETT S METHOD OF SOLVING FIRST ORDER ODES BY LIE GROUP SYMMETRY ZM Mwanzia, KC Sogomo Zablon M Mwanzia, Department of Mathematics, PO Box 536, Egerton zmusyoka@egertonacke

More information

Vectors. January 13, 2013

Vectors. January 13, 2013 Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,

More information

Elementary realization of BRST symmetry and gauge fixing

Elementary realization of BRST symmetry and gauge fixing Elementary realization of BRST symmetry and gauge fixing Martin Rocek, notes by Marcelo Disconzi Abstract This are notes from a talk given at Stony Brook University by Professor PhD Martin Rocek. I tried

More information

Department of Mathematics Luleå University of Technology, S Luleå, Sweden. Abstract

Department of Mathematics Luleå University of Technology, S Luleå, Sweden. Abstract Nonlinear Mathematical Physics 1997, V.4, N 4, 10 7. Transformation Properties of ẍ + f 1 t)ẋ + f t)x + f t)x n = 0 Norbert EULER Department of Mathematics Luleå University of Technology, S-971 87 Luleå,

More information

UNDERGROUND LECTURE NOTES 1: Optimality Conditions for Constrained Optimization Problems

UNDERGROUND LECTURE NOTES 1: Optimality Conditions for Constrained Optimization Problems UNDERGROUND LECTURE NOTES 1: Optimality Conditions for Constrained Optimization Problems Robert M. Freund February 2016 c 2016 Massachusetts Institute of Technology. All rights reserved. 1 1 Introduction

More information

LECTURE 7, WEDNESDAY

LECTURE 7, WEDNESDAY LECTURE 7, WEDNESDAY 25.02.04 FRANZ LEMMERMEYER 1. Singular Weierstrass Curves Consider cubic curves in Weierstraß form (1) E : y 2 + a 1 xy + a 3 y = x 3 + a 2 x 2 + a 4 x + a 6, the coefficients a i

More information

Linear Algebra 1 Exam 2 Solutions 7/14/3

Linear Algebra 1 Exam 2 Solutions 7/14/3 Linear Algebra 1 Exam Solutions 7/14/3 Question 1 The line L has the symmetric equation: x 1 = y + 3 The line M has the parametric equation: = z 4. [x, y, z] = [ 4, 10, 5] + s[10, 7, ]. The line N is perpendicular

More information

A symmetry-based method for constructing nonlocally related partial differential equation systems

A symmetry-based method for constructing nonlocally related partial differential equation systems A symmetry-based method for constructing nonlocally related partial differential equation systems George W. Bluman and Zhengzheng Yang Citation: Journal of Mathematical Physics 54, 093504 (2013); doi:

More information

Vector Space Concepts

Vector Space Concepts Vector Space Concepts ECE 174 Introduction to Linear & Nonlinear Optimization Ken Kreutz-Delgado ECE Department, UC San Diego Ken Kreutz-Delgado (UC San Diego) ECE 174 Fall 2016 1 / 25 Vector Space Theory

More information

LECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS. 1. Lie groups

LECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS. 1. Lie groups LECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS 1. Lie groups A Lie group is a special smooth manifold on which there is a group structure, and moreover, the two structures are compatible. Lie groups are

More information

REVIEW OF DIFFERENTIAL CALCULUS

REVIEW OF DIFFERENTIAL CALCULUS REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be

More information

The Geometry of Euler s equation. Introduction

The Geometry of Euler s equation. Introduction The Geometry of Euler s equation Introduction Part 1 Mechanical systems with constraints, symmetries flexible joint fixed length In principle can be dealt with by applying F=ma, but this can become complicated

More information

Introduction to Algebraic and Geometric Topology Week 14

Introduction to Algebraic and Geometric Topology Week 14 Introduction to Algebraic and Geometric Topology Week 14 Domingo Toledo University of Utah Fall 2016 Computations in coordinates I Recall smooth surface S = {f (x, y, z) =0} R 3, I rf 6= 0 on S, I Chart

More information

Lecture Notes in Mathematics. Arkansas Tech University Department of Mathematics. The Basics of Linear Algebra

Lecture Notes in Mathematics. Arkansas Tech University Department of Mathematics. The Basics of Linear Algebra Lecture Notes in Mathematics Arkansas Tech University Department of Mathematics The Basics of Linear Algebra Marcel B. Finan c All Rights Reserved Last Updated November 30, 2015 2 Preface Linear algebra

More information

Definition (T -invariant subspace) Example. Example

Definition (T -invariant subspace) Example. Example Eigenvalues, Eigenvectors, Similarity, and Diagonalization We now turn our attention to linear transformations of the form T : V V. To better understand the effect of T on the vector space V, we begin

More information

Numerical Approximation of Phase Field Models

Numerical Approximation of Phase Field Models Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School

More information

Notes on SU(3) and the Quark Model

Notes on SU(3) and the Quark Model Notes on SU() and the Quark Model Contents. SU() and the Quark Model. Raising and Lowering Operators: The Weight Diagram 4.. Triangular Weight Diagrams (I) 6.. Triangular Weight Diagrams (II) 8.. Hexagonal

More information

LS.1 Review of Linear Algebra

LS.1 Review of Linear Algebra LS. LINEAR SYSTEMS LS.1 Review of Linear Algebra In these notes, we will investigate a way of handling a linear system of ODE s directly, instead of using elimination to reduce it to a single higher-order

More information

Regularity for Poisson Equation

Regularity for Poisson Equation Regularity for Poisson Equation OcMountain Daylight Time. 4, 20 Intuitively, the solution u to the Poisson equation u= f () should have better regularity than the right hand side f. In particular one expects

More information

Smooth Structure. lies on the boundary, then it is determined up to the identifications it 1 2

Smooth Structure. lies on the boundary, then it is determined up to the identifications it 1 2 132 3. Smooth Structure lies on the boundary, then it is determined up to the identifications 1 2 + it 1 2 + it on the vertical boundary and z 1/z on the circular part. Notice that since z z + 1 and z

More information

Construction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem

Construction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem Proceedings of 4th Workshop Group Analysis of Differential Equations & Integrability 2009, 1 23 Construction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem George W. BLUMAN,

More information

Lie Symmetries Analysis for SIR Model of Epidemiology

Lie Symmetries Analysis for SIR Model of Epidemiology Applied Mathematical Sciences, Vol. 7, 2013, no. 92, 4595-4604 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2013.36348 Lie Symmetries Analysis for SIR Model of Epidemiology A. Ouhadan 1,

More information

APPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.

APPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2. APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product

More information

arxiv: v1 [math.ap] 25 Aug 2009

arxiv: v1 [math.ap] 25 Aug 2009 Lie group analysis of Poisson s equation and optimal system of subalgebras for Lie algebra of 3 dimensional rigid motions arxiv:0908.3619v1 [math.ap] 25 Aug 2009 Abstract M. Nadjafikhah a, a School of

More information

Inductive and Recursive Moving Frames for Lie Pseudo-Groups

Inductive and Recursive Moving Frames for Lie Pseudo-Groups Inductive and Recursive Moving Frames for Lie Pseudo-Groups Francis Valiquette (Joint work with Peter) Symmetries of Differential Equations: Frames, Invariants and Applications Department of Mathematics

More information

Definitions and Properties of R N

Definitions and Properties of R N Definitions and Properties of R N R N as a set As a set R n is simply the set of all ordered n-tuples (x 1,, x N ), called vectors. We usually denote the vector (x 1,, x N ), (y 1,, y N ), by x, y, or

More information

New methods of reduction for ordinary differential equations

New methods of reduction for ordinary differential equations IMA Journal of Applied Mathematics (2001) 66, 111 125 New methods of reduction for ordinary differential equations C. MURIEL AND J. L. ROMERO Departamento de Matemáticas, Universidad de Cádiz, PO Box 40,

More information

Special Conformal Invariance

Special Conformal Invariance Chapter 6 Special Conformal Invariance Conformal transformation on the d-dimensional flat space-time manifold M is an invertible mapping of the space-time coordinate x x x the metric tensor invariant up

More information

CORRECTIONS TO SECOND PRINTING OF Olver, P.J., Equivalence, Invariants, and Symmetry, Cambridge University Press, Cambridge, 1995.

CORRECTIONS TO SECOND PRINTING OF Olver, P.J., Equivalence, Invariants, and Symmetry, Cambridge University Press, Cambridge, 1995. CORRECTIONS TO SECOND PRINTING OF Olver, PJ, Equivalence, Invariants, and Symmetry, Cambridge University Press, Cambridge, 1995 On back cover, line 17 18, change prospective geometry projective geometry

More information

Coordinate Systems and Canonical Forms

Coordinate Systems and Canonical Forms Appendix D Coordinate Systems and Canonical Forms D.1. Local Coordinates Let O be an open set in R n. We say that an n-tuple of smooth realvalued functions defined in O, (φ 1,...,φ n ), forms a local coordinate

More information

Renormalization for Lorenz maps

Renormalization for Lorenz maps Renormalization for Lorenz maps Denis Gaidashev, Matematiska Institutionen, Uppsala Universitet Tieste, June 5, 2012 D. Gaidashev, Uppsala Universitet () Renormalization for Lorenz maps Tieste, June 5,

More information

Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems

Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September

More information

Group Invariant Solutions of Complex Modified Korteweg-de Vries Equation

Group Invariant Solutions of Complex Modified Korteweg-de Vries Equation International Mathematical Forum, 4, 2009, no. 28, 1383-1388 Group Invariant Solutions of Complex Modified Korteweg-de Vries Equation Emanullah Hızel 1 Department of Mathematics, Faculty of Science and

More information

MATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES

MATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES MATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES J. WONG (FALL 2017) What did we cover this week? Basic definitions: DEs, linear operators, homogeneous (linear) ODEs. Solution techniques for some classes

More information

Symmetry classification of KdV-type nonlinear evolution equations

Symmetry classification of KdV-type nonlinear evolution equations arxiv:nlin/0201063v2 [nlin.si] 16 Sep 2003 Symmetry classification of KdV-type nonlinear evolution equations F. Güngör Department of Mathematics, Faculty of Science and Letters, Istanbul Technical University,

More information

CALCULUS ON MANIFOLDS. 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M =

CALCULUS ON MANIFOLDS. 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M = CALCULUS ON MANIFOLDS 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M = a M T am, called the tangent bundle, is itself a smooth manifold, dim T M = 2n. Example 1.

More information

SOLUTIONS FOR PROBLEMS 1-30

SOLUTIONS FOR PROBLEMS 1-30 . Answer: 5 Evaluate x x + 9 for x SOLUTIONS FOR PROBLEMS - 0 When substituting x in x be sure to do the exponent before the multiplication by to get (). + 9 5 + When multiplying ( ) so that ( 7) ( ).

More information

Series Solutions of ODEs. Special Functions

Series Solutions of ODEs. Special Functions C05.tex 6/4/0 3: 5 Page 65 Chap. 5 Series Solutions of ODEs. Special Functions We continue our studies of ODEs with Legendre s, Bessel s, and the hypergeometric equations. These ODEs have variable coefficients

More information

LECTURE 14: REGULAR SINGULAR POINTS, EULER EQUATIONS

LECTURE 14: REGULAR SINGULAR POINTS, EULER EQUATIONS LECTURE 14: REGULAR SINGULAR POINTS, EULER EQUATIONS 1. Regular Singular Points During the past few lectures, we have been focusing on second order linear ODEs of the form y + p(x)y + q(x)y = g(x). Particularly,

More information

Introduction to String Theory Prof. Dr. Lüst

Introduction to String Theory Prof. Dr. Lüst Introduction to String Theory Prof. Dr. Lüst Summer 2006 Assignment # 7 Due: July 3, 2006 NOTE: Assignments #6 and #7 have been posted at the same time, so please check the due dates and make sure that

More information

Last Update: April 7, 201 0

Last Update: April 7, 201 0 M ath E S W inter Last Update: April 7, Introduction to Partial Differential Equations Disclaimer: his lecture note tries to provide an alternative approach to the material in Sections.. 5 in the textbook.

More information

Symmetry Reduction of Two-Dimensional Damped Kuramoto Sivashinsky Equation

Symmetry Reduction of Two-Dimensional Damped Kuramoto Sivashinsky Equation Commun. Theor. Phys. 56 (2011 211 217 Vol. 56 No. 2 August 15 2011 Symmetry Reduction of Two-Dimensional Damped Kuramoto Sivashinsky Equation Mehdi Nadjafikhah and Fatemeh Ahangari School of Mathematics

More information

Generators for Continuous Coordinate Transformations

Generators for Continuous Coordinate Transformations Page 636 Lecture 37: Coordinate Transformations: Continuous Passive Coordinate Transformations Active Coordinate Transformations Date Revised: 2009/01/28 Date Given: 2009/01/26 Generators for Continuous

More information

Group Actions and Cohomology in the Calculus of Variations

Group Actions and Cohomology in the Calculus of Variations Group Actions and Cohomology in the Calculus of Variations JUHA POHJANPELTO Oregon State and Aalto Universities Focused Research Workshop on Exterior Differential Systems and Lie Theory Fields Institute,

More information

Symbolic-Algebraic Methods for Linear Partial Differential Operators (LPDOs) RISC Research Institute for Symbolic Computation

Symbolic-Algebraic Methods for Linear Partial Differential Operators (LPDOs) RISC Research Institute for Symbolic Computation Symbolic-Algebraic Methods for Linear Partial Differential Operators (LPDOs) Franz Winkler joint with Ekaterina Shemyakova RISC Research Institute for Symbolic Computation Johannes Kepler Universität.

More information

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.

More information

8 A pseudo-spectral solution to the Stokes Problem

8 A pseudo-spectral solution to the Stokes Problem 8 A pseudo-spectral solution to the Stokes Problem 8.1 The Method 8.1.1 Generalities We are interested in setting up a pseudo-spectral method for the following Stokes Problem u σu p = f in Ω u = 0 in Ω,

More information

Lecture 2: The Simplex method

Lecture 2: The Simplex method Lecture 2 1 Linear and Combinatorial Optimization Lecture 2: The Simplex method Basic solution. The Simplex method (standardform, b>0). 1. Repetition of basic solution. 2. One step in the Simplex algorithm.

More information

Élie Cartan s Theory of Moving Frames

Élie Cartan s Theory of Moving Frames Élie Cartan s Theory of Moving Frames Orn Arnaldsson Department of Mathematics University of Minnesota, Minneapolis Special Topics Seminar, Spring 2014 The Story of Symmetry Felix Klein and Sophus Lie

More information

Lecture Notes 6: Dynamic Equations Part A: First-Order Difference Equations in One Variable

Lecture Notes 6: Dynamic Equations Part A: First-Order Difference Equations in One Variable University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 54 Lecture Notes 6: Dynamic Equations Part A: First-Order Difference Equations in One Variable Peter J. Hammond latest revision 2017

More information

JUST THE MATHS UNIT NUMBER 1.5. ALGEBRA 5 (Manipulation of algebraic expressions) A.J.Hobson

JUST THE MATHS UNIT NUMBER 1.5. ALGEBRA 5 (Manipulation of algebraic expressions) A.J.Hobson JUST THE MATHS UNIT NUMBER 1.5 ALGEBRA 5 (Manipulation of algebraic expressions) by A.J.Hobson 1.5.1 Simplification of expressions 1.5.2 Factorisation 1.5.3 Completing the square in a quadratic expression

More information

We denote the space of distributions on Ω by D ( Ω) 2.

We denote the space of distributions on Ω by D ( Ω) 2. Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below. 54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and

More information

EXTENSION OF HÖLDER S THEOREM IN Diff 1+ɛ

EXTENSION OF HÖLDER S THEOREM IN Diff 1+ɛ EXTENSION OF HÖLDER S THEOREM IN Diff 1+ɛ + (I) Azer Akhmedov Abstract: We prove that if Γ is subgroup of Diff 1+ɛ + (I) and N is a natural number such that every non-identity element of Γ has at most

More information

Reminder on basic differential geometry

Reminder on basic differential geometry Reminder on basic differential geometry for the mastermath course of 2013 Charts Manifolds will be denoted by M, N etc. One should think of a manifold as made out of points (while the elements of a vector

More information

On the Stability of the Best Reply Map for Noncooperative Differential Games

On the Stability of the Best Reply Map for Noncooperative Differential Games On the Stability of the Best Reply Map for Noncooperative Differential Games Alberto Bressan and Zipeng Wang Department of Mathematics, Penn State University, University Park, PA, 68, USA DPMMS, University

More information

u n 2 4 u n 36 u n 1, n 1.

u n 2 4 u n 36 u n 1, n 1. Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then

More information

Some Notes on Linear Algebra

Some Notes on Linear Algebra Some Notes on Linear Algebra prepared for a first course in differential equations Thomas L Scofield Department of Mathematics and Statistics Calvin College 1998 1 The purpose of these notes is to present

More information

7.1 Another way to find scalings: breakdown of ordering

7.1 Another way to find scalings: breakdown of ordering 7 More matching! Last lecture we looked at matched asymptotic expansions in the situation where we found all the possible underlying scalings first, located where to put the boundary later from the direction

More information

Symmetry Reduction of Chazy Equation

Symmetry Reduction of Chazy Equation Applied Mathematical Sciences, Vol 8, 2014, no 70, 3449-3459 HIKARI Ltd, wwwm-hikaricom http://dxdoiorg/1012988/ams201443208 Symmetry Reduction of Chazy Equation Figen AÇIL KİRAZ Department of Mathematics,

More information

1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0

1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0 Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =

More information

21 Linear State-Space Representations

21 Linear State-Space Representations ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may

More information