Symmetry Methods for Differential and Difference Equations. Peter Hydon
|
|
- Sheryl Griffith
- 5 years ago
- Views:
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 A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment of the Requirements
More informationFirst 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 informationB.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 informationAlgorithmic 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 informationClassification 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 informationSymmetry 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 informationLinearization 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 information2. 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 informationz 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 informationFundamentals 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 informationA 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 information1 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 informationOblique 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 informationMA22S3 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 informationLie 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 informationOn 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 information7 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 informationConservation 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 informationSymmetry 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 informationLinear 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 informationRelevant 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 informationPreliminaries 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 informationSymmetries 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 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 informationChapter 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 informationOn 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 informationAn 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 informationChapter 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 informationVector 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 informationCLASSIFICATIONS 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 informationMATH 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 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 informationarxiv: 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 informationTHEODORE 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 informationMultiple 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 informationMath 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 informationChapter 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 informationSOME 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 informationVectors. 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 informationElementary 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 informationDepartment 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 informationUNDERGROUND 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 informationLECTURE 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 informationLinear 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 informationA 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 informationVector 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 informationLECTURE 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 informationREVIEW 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 informationThe 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 informationIntroduction 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 informationLecture 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 informationDefinition (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 informationNumerical 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 informationNotes 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 informationLS.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 informationRegularity 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 informationSmooth 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 informationConstruction 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 informationLie 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 informationAPPENDIX 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 informationarxiv: 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 informationInductive 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 informationDefinitions 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 informationNew 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 informationSpecial 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 informationCORRECTIONS 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 informationCoordinate 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 informationRenormalization 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 informationLecture 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 informationGroup 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 informationMATH 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 informationSymmetry 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 informationCALCULUS 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 informationSOLUTIONS 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 informationSeries 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 informationLECTURE 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 informationIntroduction 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 informationLast 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 informationSymmetry 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 informationGenerators 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 informationGroup 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 informationSymbolic-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 informationBIFURCATION 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 information8 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 informationLecture 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 Orn Arnaldsson Department of Mathematics University of Minnesota, Minneapolis Special Topics Seminar, Spring 2014 The Story of Symmetry Felix Klein and Sophus Lie
More informationLecture 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 informationJUST 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 informationWe 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 informationStability 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 information154 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 informationEXTENSION 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 informationReminder 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 informationOn 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 informationu 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 informationSome 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 information7.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 informationSymmetry 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 information1. 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 information21 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