Reaction-Diusion Systems with. 1-Homogeneous Non-linearity. Matthias Buger. Mathematisches Institut der Justus-Liebig-Universitat Gieen

Size: px
Start display at page:

Download "Reaction-Diusion Systems with. 1-Homogeneous Non-linearity. Matthias Buger. Mathematisches Institut der Justus-Liebig-Universitat Gieen"

Transcription

1 Reaction-Dision Systems ith 1-Homogeneos Non-linearity Matthias Bger Mathematisches Institt der Jsts-Liebig-Uniersitat Gieen Arndtstrae 2, D Gieen, Germany Abstract We describe the dynamics of a system of to reaction-dision eqations ith 1-homogeneos non-linearity We sho that either an order-presering property holds and can be sed in order to determine the limiting behaior in some (inariant) sets or the long time behaior of all soltions can be described by looking at one scalar reaction-dision eqation only AMS No 35K57 Key ords: Reaction-dision eqations, order-presering, redction to one eqation, oscillation nmber 1 Introdction The dynamics of soltions of one reaction-dision eqation d dt y = y xx + g(t; x; y) (1) ith Dirichlet or other bondary conditions on an interal (0; L), L > 0, has been examined by many athors [10], [11], [12], [15], [18] Far less is knon for a system of to reaction-dision eqations d dt = xx xx + f t; x; ; x 2 (0; L); t > 0: (2) In this paper e examine system (2) ith a 1-homogeneos non-linearity f(; ) = h 1 (; ), hich means that h 1 (; ) = h 1 (; ) for all real Frthermore, e assme that h 1 is C 1 in IR 2 n f(0; 0)g It trns ot that there are to dierent cases, depending on hether the set E := ( 0 2 IR :? sin 0 cos 0 1 h 1 cos 0 sin 0 = 0 )

2 2 MBger is empty or not Case 1 E 6= ; Then E has innitely many elements and e can take ; 2 E, <, sch that 62 E for all < < We rite the initial data ( 0 ; 0 ) in polar coordinates 0 0 = r 0 cos ' 0 sin ' 0 : (3) If r 0 is non-negatie and ' 0 has its ales in the interal (; ), then e sho that the soltion (; ) starting at ( 0 ; 0 ) is attened ot and tends to a planar soltion r 0 cos ' 0 sin ' 0 here ' 0 is a constant and eqals either or We call a state ( 0 ; 0 ) planar, if it can be ritten in the form (4) ith a constant, ie x-independent angle ' 0 We note that ( 0 ; 0 ) is planar if and only if 0 and 0 are linearly dependent A soltion (; ) is called planar if (; )(t) is planar for all t This means that the angle fnction ' might depend on t bt not on the space ariable x It is easy to obsere that (; )(t) is a planar soltion if and only if the initial ale is planar This notation is motiated by the fact that for a planar soltion (; ) the cre (4) t : [0; 1] 3 x 7 (x; (t; x); (t; x)) 2 [0; 1] IR 2 lies in a plane P t (hich depends on t), as one can see in the folloing gre,-axis 6 P t 0 1 Figre 1: cre t for a planar soltion (; ) The proof of the fact that the anglar fnction '(t; ) associated ith (; )(t) conerges either to or to is based on the folloing order-presering property: If (? ;? )(t) is a planar soltion ith '? (0) '(0; ), then e get '? (t) '(t; ) for all t 0 Similarly, e get a fnction ' + Hence, the soltion (; ) lies beteen a loer and an pper planar soltion Then it only remains to sho that the angle fnctions ' (t) of the planar soltions both conerge either to or to

3 Reaction-Dision Systems ith 1-Homogeneos Non-linearity 3 6? 2 ẏ graph of (0; L) 3 x 7 0(x) 0 (x) 2 IR 2 2? '+ O Figre 2: a sitation in hich '?, ' + and ' conerge to We note that in the case E = 2 Z the non-linearity h 1 can be ritten in the form h 1 (; ) = 1(; ) : 2 (; ) Then the maximm priciple yields that soltions (; ) ith initial ale ( 0 ; 0 ), 0 > 0; 0 > 0, stay positie for all t > 0, and or reslt follos from ell knon order-presering properties [9], [14] We note that in this paper e se a dierent ordering relation hich has the adantage that it does also ork if E 6= 2 Z Therefore, or reslt in case 1 is a kind of generaliation of the ell knon orderpresering reslts for 1-homogeneos non-linearities Case 2 E = ; Another concept, hoeer, is needed in the case E = ; In this case, loosely speaking, the part of the non-linearity hich indces a rotation on IR 2 arond the origin does neer anish Hence, e cannot expect the existence of any stationary soltion besides ero, bt e shold be able to obsere periodic (or more diclt) motion heneer ero is nstable The description of the dynamics in this case seems to be really complicated, bt e sho that a major part can be done by looking at planar soltions hich are obtained from a scalar (time-periodic) reaction-dision eqation instead of a system of to eqations This method can be looked at as a generaliation of the approach sed in [3], [4], [5] for the model system d = xx +? + F (5) dt xx As a main reslt e obtain that each bonded soltion of (2) ith 1-homogeneos non-linearity h 1 tends either to a stationary or periodic (planar) soltion, ie e get a Poincare-Bendixson reslt Remarks on both cases We note that the fact that the dision rates coincide in both eqations is a crcial '? ū

4 4 MBger point Nmerical stdies sho that the dynamics can be totally dierent from the cases mentioned aboe for dierent dision rates, in general The restrictions on the non-linearity, hoeer, can be eakened It trns ot that all of or reslts hold for a non-linearity f of the form f(t; x; ; ) = h 1 (; ) + F (t; x; ; ) here F : [0; 1) [0; L] IR 2 IR is C 1 We ill ork ith the more general system d = xx + h dt 1 (; ) + F (t; x; ; ) (6) xx for the rest of this paper 2 Notation and preliminaries Let IR 2 := IR 2 n f(0; 0)g For eery (; ) 2 IR 2 the set f(; ); (?; )g is a base of IR 2 and, ths, h 1 (; ) 2 IR 2 can be ritten in the form h 1 = f 0 +? here f 0 (; ); g 0 (; ) 2 IR are niqely determined Since h 1 2 C 1 (IR 2 ; IR 2 ) and (7) implies that f 0 g 0 = = ? h 1 h 1 g 0 ; ; (7) e get f 0 ; g 0 2 C 1 (IR 2 ; IR) Since for all 2 IR = h 1 f 0 = + f 0? g 0 = h 1? + g 0 ; e obtain that f 0 = f 0 ; g 0 = g 0 : Therefore, it is scient to dene f 0 ; g 0 on S 1 here S 1 is dened by S 1 := f(x; y) 2 IR 2 : x 2 + y 2 = 1g :

5 Reaction-Dision Systems ith 1-Homogeneos Non-linearity 5 In particlar, this implies that f 0 ; g 0 and their deriaties are bonded (on S 1 and, hence, on IR 2 ) For 0 2 IR e introdce the soltion (; 0 ) : IR IR of the ODE d dt = g 0(cos ; sin ) ; (0) = 0 : (8) Since g 0 is C 1 on S 1 (and has, ths, bonded deriaties), is ell dened Lemma 1 The fnction (; 0 ) is either constant or strictly monotone Proof If g 0 (cos 0 ; sin 0 ) = 0, then (t; 0 ) = 0 for all t We assme that g 0 (cos 0 ; sin 0 ) 6= 0, log e only deal ith the case g 0 (cos 0 ; sin 0 ) > 0 If as not strictly increasing, then there old be some t 0 2 IR ith d dt (t 0; 0 ) = 0 hich old imply that g 0 (cos (t 0 ; 0 ); sin (t 0 ; 0 )) = 0 and, ths, old be constant, contradicting the assmption Lemma 2 If g 0 (cos ; sin ) = 0 for some 2 [0; 2), ie 2 E, then (; 0 ) is bonded for all 0 ; if there is no sch, ie E = ;, then (; 0 ) is nbonded Proof 1 If there is as aboe, then eery non-constant fnction (; 0 ) satises (t; 0 ) 6= mod 2 for all t Since (; 0 ) is continos, (IR; 0 ) is contained in an interal of length 2; in particlar, (; 0 ) is bonded 2 If (; 0 ) is bonded, then the limit := lim t1 (t; 0 ) 2 IR exists by Lemma 1 Ths, there is a seqence t n % 1 sch that d dt (t n; 0 ) 0 (n 1) This implies that g 0 (cos ; sin ) = 0 : 3 Main reslts 31 The case E = ; An element ( 0 ; 0 ) of C 1 ([0; L]; IR 2 ) is called planar, if 0 and 0 are linearly dependent oer IR, ie if there are r 0 2 C 1 ([0; L]; IR) and 0 2 IR sch that 0 0 = cos 0 sin 0 The set of all planar elements of C 1 ([0; L]; IR 2 ) is denoted by P r 0 :

6 6 MBger Theorem 1 Let (; ) 2 C 1 ([0; 1)[0; L]; IR 2 ) be a (classical) soltion of system (6) ith Dirichlet or Nemann bondary conditions If E = ;, then e get dist C 1f(; )(t; ); Pg 0 (t 1) : This means that in the case E = ; the -limit set of any (bonded) soltion of (6) is contained in P By Theorem 1, the crcial part of the dynamics of (6) takes place in the set P of planar elements The folloing theorems deal ith the dynamics on P Theorem 2 The set P is positiely inariant, ie eery soltion (; ) of (6) hich starts in P is a planar soltion, ie (; )(t; ) 2 P for all t 0 If ( 0 ; 0 ) = (cos 0 ; sin 0 )r 0 ith 0 2 IR, then (; ) has the form (t; ) = cos (t; 0) sin (t; 0 ) r(t; ) (9) here r : [0; 1) [0; L] IR satises ths scalar eqation " d dt r = r cos (t; xx + r f 0 ) 0 sin (t; 0 ) + F t; x; cos (t; 0 ) sin (t; 0 ) r # (10) ith initial ale r 0 We note that the assertion of Theorem 2 is also tre in the case E 6= ;, bt the dynamics on P is important only if e hae the reslt of Theorem 1 hich ensres that all soltions nally tend to planar ones Theorem 3 Let (; ) be a (bonded) planar soltion and choose 0 ; r 0 ; r as in Theorem 2 If E = ;, then (; 0 ) mod 2 is periodic and the period p is gien by p = inff > 0 : j(; 0)j = 2g : In particlar, p depends only on g 0 If, in addition, (6) is atonomos or F is periodic in t ith a period p F sch that p F =p is rational, then the -limit set associated ith (; ) consists of periodic planar elements only, ie if e take an initial ale ( 0 ; 0 ) 2 (; ), then (; ) is a periodic planar soltion of (6) 32 The case E 6= ; Gien ; 2 IR, <, e call I ; := f( 0 ; 0 ) 2 C 1 ([0; L]; IR 2 ) : ( 0 ; 0 ) = r 0 (cos ' 0 ; sin ' 0 ) ith r 0 ; ' 0 2 C([0; L]; IR) sch that r 0 > 0 and < ' 0 < in (0; L)g the anglar space beteen and If, in addition, ; 2 E and

7 Reaction-Dision Systems ith 1-Homogeneos Non-linearity 7 62 E for all < <, then e call I ; an anglar space associated ith E We note that E 6= ; implies that E contains innitely many elements We sho that an anglar space associated ith E is positiely inariant nder the semio generated by (6) and all soltions ith initial ale in this anglar space tend to planar soltions ith (constant) angle or Theorem 4 Let I ; be an anglar space associated ith E and (; ) a (bonded) soltion ith initial ale (; )(0) 2 I ; Then e get: (a) (; )(t; ) 2 I ; for all t 0 (b) There are continos fnctions ' : (0; 1) C([0; L]; IR) and r : (0; 1) C 1 ([0; L]; IR) sch that = r cos ' sin ' (c) If g 0 (cos ; sin ) > 0 for < <, then '(t; ) conerges to in C([0; L]; IR), in case g 0 (cos ; sin ) < 0 to 4 Interpretation of the reslts If e hae E = ;, then the rotation on the ales of (; ) drien by g 0 does neer stop formally this means that (; 0 ) is strictly monotone and nbonded In this case, or main reslt, Theorem 1, ensres that a major part of the dynamics of (6) can be described by looking only at the dynamics on P, becase e kno that all soltions of (6) tend to P for t 1 Then Theorems 2 and 3 sho that the dynamics on P can be described by soling rst the ODE (8) and then looking at the scalar reaction-dision eqation (10) hich has been stdied in [6], [7] and [8] This is obiosly mch easier than dealing ith the fll system (6); in particlar e can se oscillation nmber reslts in order to attack eqation (10) hich are, in general, not aailable for systems Neertheless, e note that the global attractor of (6) (if it exists) contains, in general, more elements than jst the nion of ( 0 ; 0 ) for all points ( 0 ; 0 ) We can een sho that in some cases the global attractor ill contain non-planar elements The reslts [5] on nstable directions for periodic soltions of the model system (5) cold, for example, be sed in order to proe the existence of non-planar heteroclinic soltions If e hae E 6= ; and 2 E, then (t; ) = for all t and the set f( 0 ; 0 ) := r 0 (cos ; sin ) : r 0 2 C 1 ([0; L]; IR)g is inariant nder the semio indced by (6) In this case Theorem 4 describes the dynamics of all soltions starting in an anglar space I ; beteen to conseqtie elements of E, ie in a space hich is bonded by inariant sets :

8 8 MBger Theorem 4 says that all soltions hich start in the anglar space I ; ill be attened ot and conerge to the bondary of this anglar space (see Figre 2) In the proof of Theorem 4 e ill see that the crcial point is that? '(t 0 ; ) + for some t 0 yields (t? t 0 ;? ) '(t; ) (t? t 0 ; + ) for all t t 0 This (eak) order-presering property together ith the fact that (; 0 ) conerges for any 0 is the main idea of the proof of Theorem 4 We note that otside the anglar spaces een non-planar eqilibria may exist We can, for example, take L = and h 1 (; ) = (; =2) hich has the non-planar eqilibrim (sin x; sin(2x)) (and E = 2 Z) 5 Constrction of an appropriate transformation Let (; ) 2 C 1 ([0; 1) [0; L]; IR 2 ) be a soltion of (6) ith Dirichlet or Nemann bondary conditions For ' 2 IR e consider cos '? sin ' M(') := : sin ' cos ' Gien 0 2 IR e introdce (t; ) := M(?(t; 0 )) (t; ) : An elementary comptation shos that (; ) satises d = xx +? " g dt xx 0 M((t; 0 ))? d # dt (t; 0) + " F t; x; M((t; 0 )) # +f 0 M((t; 0 )) : (11) 6 Oscillation nmber reslts The concept of oscillation nmbers [1], [13], [16] (sometimes also called lap nmbers or Matano nmbers) deals ith soltions y : [0; 1)[0; L] IR of a linear eqation d dt y = y xx + q(t; x)y (12) here q : [0; 1) [0; L] IR is bonded on [0; T ] [0; L] for all T > 0 We dene the nmber of sign changes (the oscillation nmber) by Z(t) := spfk 2 IN : 9 0 < x 1 < : : : < x k < L :

9 Reaction-Dision Systems ith 1-Homogeneos Non-linearity 9 y(t; x j )y(t; x j+1 ) < j k? 1g : It is remarkable that Z(t) is nite for all t > 0, no matter ho many eros the initial state y(0; ) has, as Angenent [1] proed The most important reslt abot oscillation nmbers is that Z : (0; 1) IN is a non-increasing fnction Precisely, the ale of Z(t) decreases heneer y(t; ) has a mltiple ero, ie if y(t 0 ; ) has a mltiple ero, then Z(t 1 ) > Z(t 2 ) for all t 1 < t 0 < t 2 These oscillation nmber reslts play an important role in the examination of a scalar reaction-dision eqation bt can in general not be applied to systems of reaction-dision eqations In order to se these reslts for the fnction dened in the last section, e hae to sho that satises some eqation of the form (12) Since soles " d dt = xx + g 0 M((t; 0 )) " + F t; x; M((t; 0 ))? d # dt (t; 0) + f 0 M((t; 0 )) # by (11), it only remains to sho that q 1 := " g 0 M((t; 0 ))? d dt (t; 0) # is bonded on [0; T ] [0; L] for all T > 0 If q 1 is bonded, then ill sole eqation (12) ith q := q 1 + F t; x; M((t; 0 )) + f 0 M((t; 0 )) and e ill be able to apply oscillation nmber reslts to We assme that q 1 is not bonded Then there is T > 0 and a seqence (t n ; x n ) in [0; T ] [0; L] sch that jq 1 (t n ; x n )j 1 (n 1) and (t n ; x n ) 6= 0 for all n 2 IN We may assme that (t n ; x n ) conerges to (t 0 ; x 0 ) 2 [0; T ] [0; 1] Since g 0 and, by denition of, also d are bonded, this implies that dt (t n ; x n ) 0 (n 1) : (t n ; x n ) Hence, ( n ; n ) dened by n := n 1 q 2 (t n ; x n ) + 2 (t n ; n ) (t n ; x n ) (t n ; x n ) fll n 0 (n 1) This means that j n j 1 (n 1) and log e may assme that n 1 (n 1) Since the map [0; 2] S 1 3 ('; (; )) 7 g 0 M(') 2 IR

10 10 MBger is (niformly) continosly dierentiable, it is Lipschit-continos in (; ) 2 S 1 niformly for all ' 2 [0; 2] ith Lipschit-constant L g Hence, e get g (t 0 M((t n ; 0 )) n ; x n )? d (t n ; x n ) dt (t n; 0 ) = g 0 M((t n ; 0 )) L n g? 1 n 0 = L g q( n? 1) 2? 2 n : n n? g 0 M((t n ; 0 )) 1 0 We dene n 2 IR by ( n ; n ) = (cos n ; sin n ) Then e get tan n = (t n ; x n )=(t n ; 0 ) and n 0 (n 1) Frthermore, it follos that q ( n? 1) 2? n 2 = 2 n sin = O(tan n ) = O (t n; x n ) (n 1) 2 (t n ; x n ) and, ths, q 1 (t n ; x n ) is bonded for n 1 contradicting or assmption 7 Proof of Theorem 1 Step 1 Assme that the assertion is false and constrct a non-planar element in the -limit set nder this assmption By standard argments [12], [17] e kno that the soltions of (6) bild a (local) semio on some Sobole space; in particlar, the orbit? := f(; )(t) : t 0g of the bonded soltion (; ) is a precompact sbset of the space X := D( ), 2 (0; 1) (in the notation of Henry [12]) Taking < 1 sciently large and sing the Sobole embedding theorem, e obtain that? lies precompact in C 1 ([0; L]; IR 2 ) This implies that, if the assertion of or theorem does not hold, then there is ( 0 0; 0) 0 2 C 1 ([0; L]; IR 2 ) n P and a seqence (t n ), t n % 1, sch that jj(; )(t n ; )? ( 0 0; 0 0)jj C 1 0 (n 1) : This means that ( 0 ; ) is contained in the -limit set associated ith (; ) Let ( 0 ; 0 ) be the soltion of (6) ith initial ale ( 0 ; 0 0 ) We note that 0 (0 ; 0 ) is dened for all t 2 IR since ( 0 0; 0) 0 is an element of the -limit set Step 2 Denition of 0 We ant to take 0 2 IR sch that 0 (0; ) dened sing the formla 0 0 (t; ) := M(?(t; 0 )) 0 0 (t; )

11 Reaction-Dision Systems ith 1-Homogeneos Non-linearity 11 has a mltiple ero If 0 or 0 0 hae a mltiple ero, then e can take 0 0 = =2 or 0 = 0 If neither 0 0 nor 0 0 hae mltiple eros, then e proceed as follos: Using the reslt of the last section, e can apply oscillation nmber reslts to 0 for eery 0 Since ( 0 ; 0 ) are also dened for negatie t, 0 soles an eqation of the form (12) for all t 2 IR Then Angenent's reslt [1] ensres that the nmber of eros of 0 (t; ) is nite for all t; in particlar for t = 0 Applying this reslt for 0 = =2 and 0 = 0, e obtain that 0 as ell as hae only nitely many eros (hich are all simple by assmption) Hence, there is a continos fnction ' 0 0 : [0; L] IR hich satises tan ' 0 0 (x) = 0 0(x) 0 0(x) for all x 2 [0; L] ith 0 0 (x) 6= 0 This reslt is an easy conseqence of de l'hospital's formla, and it can also be fond in [2,p696,Lemma 2] We set 0 := maxf' 0 0(x) : x 2 [0; L]g and take x 0 2 [0; L] sch that ' 0 0(x 0 ) = 0 Then e get 0 (0; x 0 ) = 0 sing the denition of 0 It remains to 0 (0; x 0 ) = 0 Case 1 0 (0; x 0 ) = 0 Then 0 (0; x 0 ) = 0 implies that 0 (x 0 0) = 0 (x 0 0) = 0 and, ths, Hence, it follos that tan 0 (x 0 0(x 0 0 (0; x 0 (0; x 0 ) = tan('0 (x 0 0)? 0 ) = 0 : In particlar, 0 (0; x 0 ) = 0 Case 2 0 (0; x 0 ) 6= 0 Let U 0 be a neighborhood of x 0 in [0; L] in hich 0 (0; ) has no eros Since ' 0 0(x) 0 for all x 2 [0; L], e get 0 tan(' 0 0(x)? 0 ) = 0 (0; x) 0 (0; x) for all x 2 U 0 Since 0 (0; ) has no sign change in U 0, it follos that 0 (0; ) mst not hae a sign change in U 0, too Hence x 0 mst be a mltiple ero of 0 (0; ) Step 3 Applying oscillation nmber reslts to 0 As shon in Section 6, e can apply oscillation nmber reslts to 0 Let Z 0 : IR IN be the corresponding oscillation nmber We note that e can dene Z 0 on the hole real axis becase 0 is dened for all t 2 IR (and soles a linear eqation of

12 12 MBger the form (12) on each interal [?; 1), 0) Since 0 (0; ) has a mltiple ero by Step 2, e get Z 0 (t? ) > Z 0 (t + ) for all t? < 0 < t + We note that 0 (t; ) can only hae mltiple eros for contably many t 2 IR since Z 0 has its ales in IN, is non-increasing and decreases each time 0 (t; ) has a mltiple ero Ths, e can choose t? < 0 < t + in a ay that 0 (t ; ) hae no mltiple eros Step 4 Denition of (; ) By assmption, E is empty and, ths, (; 0 ) is nbonded by Lemma 2 Therefore, (; 0 ) has to be periodic; e denote its period by p Wlog e may assme that t n mod p 2 [0; p) has the same ale for all n; otherise e take a sbseqence (t nk ) sch that m k := t nk mod p conerges in [0; p], denote its limit by and set k := t nk? (m k? ) Then e get (; )( k ) ( 0 0; 0), 0 and e can se ( k ) instead of (t n ) Wlog e may assme that = 0; otherise replace (; ) by (^; ^) := (; )(?) hich does also sole an eqation of the form (6) Then e introdce (; ) by (t; ) := M(?(t; 0 )) (t; ) for t 0 We note that (t n ) as chosen in a ay that t n is an integer mltiple of p for all n Ths, e get (t n ; 0 ) = 0 and (t n ; ) = M(?(t n ; 0 )) { } =M (? 0 ) Frthermore, an analogos comptation shos that (t n + t ; ) 0 0 (t n ; ) M(? 0 ) = (t ; ) in C 1 ([0; L]; IR) (0; ) (0; ) : (13) Since 0 (t + ; ) and 0 (t? ; ) hae no mltiple eros, it follos that the nmber of eros of (t n + t ; ) and 0 (t ; ) coincide from some index n on Since Z 0 (t? ) > Z 0 (t + ), the nmber Z(t) of eros of (t; ) accmlates at to dierent ales Z 0 (t ) Since oscillation nmber reslts hold for by Section 6, Z has to be non-increasing hich is a contradiction Remark We note that the assmption E = ; is sed in Step 4 only It ensres that (; 0 ) is nbonded and, ths, periodic Otherise (13) does not hold and e cannot conclde that ( 0 ; 0 )(0) is a limit point of some fnction (; ) dened as aboe

13 Reaction-Dision Systems ith 1-Homogeneos Non-linearity 13 8 Proof of Theorems 2 and 3 Sppose that ( 0 ; 0 ), 0 and r 0 are gien as in Theorem 2 First dene r by (10) and then (; ) by (9) An elementary comptation shos that (; ) is a soltion of (6) ith initial ale ( 0 ; 0 ) Using the niqeness of the soltions of (6) (see, for example, [12]), the assertion of Theorem 2 follos By Lemma 2, (; 0 ) is nbonded In particlar, (; 0 ) is not constant and Lemma 1 implies that it is strictly monotone Wlog assme that (; 0 ) is strictly increasing Since (t; 0 ) is not bonded for t 1, there is a niqely determined p > 0 sch that (p ; 0 ) = : Ths, (; 0 ) is periodic ith period p Frthermore, there is t 0 2 [ 0 ; 0 + p ) sch that (t 0 ; 0 ) = 2k ith some integer k Then it follos that (t 0 + p j; 0 ) = 2(k + j) for all integers j In particlar, e get (t 0? p k; 0 ) = 0 and, ths, (p ; 0) = (p ; (t 0? p k; 0 )) = (t 0? p k + p ; 0 ) = (t 0? p (k? 1); 0 ) = 2 hich shos that p is gien by the expression mentioned in Theorem 3 We assme that F is t-periodic ith period p F here p =p F is rational This means that there is p > 0 ach that p=p as ell as p=p F are positie integers (in the atonomos case, one can jst take p = p ) In particlar, eqation (10) depends periodically on t (ith period p) Let ( 0 ; 0 ) be an element of (; ) and (; ) be the corresponding soltion of (6) Then (; ) is planar by Theorem 1 By Theorem 2, there is 0 2 IR and a fnction r sch that (9) and (10) hold Since ( 0 ; 0 ) = r(0)(cos 0 ; sin 0 ) is an element of (; ), it is a chain-recrrent point Then an elementary analysis like it as done in [5,Theorem 2] shos that r(0) is a also chain-reccrent point for eqation (10) Then [8] ensres that r is a p-periodic soltion of (10) and, since (; 0 ) is p-periodic, (; ) is a p-periodic soltion of (6) 9 Proof of Theorem 4 1 We ill proe the assertion only in the case =?=2 The general case can easily be redced to this sitation setting 0 0 (t) := M(?? =2) (t) : Then ( 0 ; 0 ) soles an eqation similar to (6) replacing F; f 0 ; g 0 by f 0 0 := f 0 M( + =2) ;

14 14 MBger F 0 t; x; g 0 0 := g 0 M( + =2) := F t; x; M( + =2) ; : Hence, e hae 0 =?=2, 0 =?? =2 and the assertion follos sing ' = ' =2 2 We note that =?=2 (Step 1) implies =2 2 E sing the fact that g 0 (0; 1) = g 0 (0;?1) (0-homogenicity) This gies 2 (?=2; =2] Hence, (; )(0) 2 I ; implies 0 = r 0 {} 0 cos ' {} 0 0 : 2(?=2;=2) Frthermore, e get g 0 (0;?1) = g 0 (cos ; sin ) = 0 Since the restriction of g 0 on S 1 is continosly dierentiable, then fnction ^g 0 gien by ^g 0 : S 1 3 (cos ; sin ) 7 g 0(cos ; sin ) cos sin 2 IR is ell dened and continos Setting ^g 0 (; ) := g 0(;) for all (; ) 2 IR 2, ^g 0 is 0-homogeneos, continos and bonded; in particlar, it is contained in L 1 (IR 2 ) Since satises d dt = xx + (F + f 0 )? g 0 = xx + (F + f 0? ^g 0 ) and (0; ) = 0 0, maximm principle argments imply that, for all t > 0, (t; x) > 0 for x 2 (0; L) and x (t; 0); x (t; 1) 6= 0 Ths, there is a continos fnction q : (0; 1) [0; L] IR sch that q(t; x) = (t; x) ; x 2 (0; L) : (t; x) (Note that q(t; 0); q(t; 1) can be dened sing de'l Hospital's formla) We dene ( ) ' : (0; 1) 3 t 7 [0; L] 3 x 7 arctan q(t; x) 2 (?=2; =2) 2 C([0; L]; IR) and r by r(t; ) := (t; )= cos '(t; ) Then, (b) follos We note that the positiity of and the fact that '(t; ) has its ales in (?=2; =2) implies that r(t; ) has only positie ales in (0; L) for all t > 0 3 We take t 0 > 0 and set? := min [0;L] '(t 0 ; ), + := max [0;L] '(t 0 ; ) Note that ' is continos by Step 2 Let = be dened as in Section 5 here the constant 0 shold hae the ale (?t 0 ; ) Then e get (t 0 ; ) = r(t 0 ; ) sin('(t 0 ; )? ) : Using the fact that r(t 0 ; ) is positie by Step 2 and '(t 0 ; )?? 2 [0; ), it follos that? (t 0 ; ) 0 We take x? 2 [0; L] sch that '(t 0 ; x? ) =? Then? (t 0 ; )

15 Reaction-Dision Systems ith 1-Homogeneos Non-linearity 15 has a mltiple ero at x = x? By Section 5, e can apply oscillation nmber reslts to? Hence,? (t; ) is either the ero fnction or strictly positie for all t > t 0 In particlar, e get Analogosly, e get min '(t; ) (t? t 0 ;? ) for all t > t 0 [0;L] max '(t; ) (t? t 0 ; + ) for all t > t 0 [0;L] 4 The soltions of (6) form a (local) semio on the space C 1 ([0; L]; IR 2 ) Hence, (; )(0) 2 I ; implies that there is > 0 sch that < '(t; ) < for all 0 t < Ths, e hae (; )(t) 2 I ; for t 2 [0; ) Take t 0 := =2 and dene as in Step 3 Then e get <? + < and, by denition of the fnction (; ), < (t;? ) (t; + ) < for all t 0 Ths, Step 3 implies that (; )(t) 2 I ; for all t =2 and, sing the reslt mentioned aboe, for all t 0 This proes (a) 5 Take t 0 > 0 and dene as in Step 3 Then e get (; )(t 0 ) 2 I ; by Step 4 and, ths, <? + < : If g 0 (cos 0 ; sin 0 ) > 0 for < 0 <, then Lemmata 1 and 2 yield (t; ) (t 1) Since (t? t 0 ;? ) min [0;L] '(t; ) max '(t; ) (t? t 0 ; + ) for all t > t 0 [0;L] by Step 3, it follos that '(t; ) conerges to in C([0; L]; IR) Analogosly, e get conergence to in the case g 0 (cos 0 ; sin 0 ) < 0 for < 0 < This proes (c) Acknoledgement The athor thanks Prof B Fiedler, Freie Uniersitat Berlin, for helpfl sggestions hich led to the examination of the problem dealt ith in this paper References [1] S Angenent, The ero set of a soltion of a parabolic eqation, J reine ange Math 390 (1988), 79{96 [2] M Bger, Torsion nmbers, a tool for the examination of symmetric reaction-dision systems related to oscillation nmbers, Disc Cont Dyn Sys 4 (1998), 691{708

16 16 MBger [3] M Bger, Periodic soltions of an atonomos reaction-dision eqation, preprint 1997 [4] M Bger, On the existence of stationary and periodic soltions for a class of atonomos reaction-dision systems, to appear in Ad Di Eq [5] M Bger, Conergence to periodic soltions in atonomos reactiondision systems on one space ariable, to appear in J Di Eq [6] P Brnosky, B Sanstede, Conergence in general periodic parabolic eqations in one space dimension, Nonlin Anal 18 (1992), 209{215 [7] X-Y Chen, H Matano, Conergence, asymptotic periodicity, and nitepoint blo-p in one-dimensional semilinear heat eqations, J Di Eq 78 (1989), 160{190 [8] X-Y Chen, P Polacik, Gradient-like strctre and Morse decompositions for time-periodic one-dimensional parabolic eqations, J Dynam Di Eq 7 (1995), 73{107 [9] E Dancer, P Hess, Stability of xed points for order-presering discrete-time dynamical systems, J Reine Ange Math 419 (1991), 125{139 [10] B Fiedler, J Mallet-Paret, A Poincare-Bendixson theorem for scalar reaction dision eqations, Arch Rat Mech Anal 107:4 (1989), 325{345 [11] JK Hale, G Ragel, Conergence in gradient-like systems ith applications to PDE, Z Ange Math Phys 43 (1992), 63{124 [12] D Henry, Geometric theory of semilinear parabolic eqations, Springer- Verlag, Berlin-Heidelberg 1981 [13] D Henry, Some innite dimensional Morse Smale systems dened by parabolic dierential eqations, J D E 59 (1985), 165{205 [14] SB Hs, HL Smith, P Waltman, Competitie exclsion and coexistence for competitie systems on ordered Banach spaces, Trans AMS 348 (1996), 4083{4094 [15] H Matano, Conergence of soltions of one-dimensional semilinear eqations, J Math Kyoto Uni 18 (1978), 221{227 [16] H Matano, Nonincrease of the lap nmber of a soltion for a one-dimensional semi-linear parabolic eqation, J Fac Sci Uni Tokyo, Sec IA 29 (1982), 401{441 [17] A Pay, Semigrops of linear operators and applications to partial dierential eqations, Springer-Verlag, Ne York-Berlin-Heidelberg-Tokio 1983

17 Reaction-Dision Systems ith 1-Homogeneos Non-linearity 17 [18] TI Zelenyak, Stabiliation of soltions of bondary ale problems for a second order parabolic eqation ith one space ariable, Di Eq 4 (1968), 17{22

Subcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany

Subcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2

More information

1 The space of linear transformations from R n to R m :

1 The space of linear transformations from R n to R m : Math 540 Spring 20 Notes #4 Higher deriaties, Taylor s theorem The space of linear transformations from R n to R m We hae discssed linear transformations mapping R n to R m We can add sch linear transformations

More information

Spanning Trees with Many Leaves in Graphs without Diamonds and Blossoms

Spanning Trees with Many Leaves in Graphs without Diamonds and Blossoms Spanning Trees ith Many Leaes in Graphs ithot Diamonds and Blossoms Pal Bonsma Florian Zickfeld Technische Uniersität Berlin, Fachbereich Mathematik Str. des 7. Jni 36, 0623 Berlin, Germany {bonsma,zickfeld}@math.t-berlin.de

More information

The Brauer Manin obstruction

The Brauer Manin obstruction The Braer Manin obstrction Martin Bright 17 April 2008 1 Definitions Let X be a smooth, geometrically irredcible ariety oer a field k. Recall that the defining property of an Azmaya algebra A is that,

More information

arxiv: v1 [math.co] 25 Sep 2016

arxiv: v1 [math.co] 25 Sep 2016 arxi:1609.077891 [math.co] 25 Sep 2016 Total domination polynomial of graphs from primary sbgraphs Saeid Alikhani and Nasrin Jafari September 27, 2016 Department of Mathematics, Yazd Uniersity, 89195-741,

More information

Concept of Stress at a Point

Concept of Stress at a Point Washkeic College of Engineering Section : STRONG FORMULATION Concept of Stress at a Point Consider a point ithin an arbitraril loaded deformable bod Define Normal Stress Shear Stress lim A Fn A lim A FS

More information

PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS. 1. Introduction

PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS. 1. Introduction PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS TANYA DEWLAND, JEROME WESTON, AND RACHEL WEYRENS Abstract. We will be determining qalitatie featres of a discrete dynamical system of homogeneos difference

More information

Reduction of over-determined systems of differential equations

Reduction of over-determined systems of differential equations Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical

More information

Graphs and Networks Lecture 5. PageRank. Lecturer: Daniel A. Spielman September 20, 2007

Graphs and Networks Lecture 5. PageRank. Lecturer: Daniel A. Spielman September 20, 2007 Graphs and Networks Lectre 5 PageRank Lectrer: Daniel A. Spielman September 20, 2007 5.1 Intro to PageRank PageRank, the algorithm reportedly sed by Google, assigns a nmerical rank to eery web page. More

More information

ON THE PERFORMANCE OF LOW

ON THE PERFORMANCE OF LOW Monografías Matemáticas García de Galdeano, 77 86 (6) ON THE PERFORMANCE OF LOW STORAGE ADDITIVE RUNGE-KUTTA METHODS Inmaclada Higeras and Teo Roldán Abstract. Gien a differential system that inoles terms

More information

J.A. BURNS AND B.B. KING redced order controllers sensors/actators. The kernels of these integral representations are called fnctional gains. In [4],

J.A. BURNS AND B.B. KING redced order controllers sensors/actators. The kernels of these integral representations are called fnctional gains. In [4], Jornal of Mathematical Systems, Estimation, Control Vol. 8, No. 2, 1998, pp. 1{12 c 1998 Birkhaser-Boston A Note on the Mathematical Modelling of Damped Second Order Systems John A. Brns y Belinda B. King

More information

Characterizations of probability distributions via bivariate regression of record values

Characterizations of probability distributions via bivariate regression of record values Metrika (2008) 68:51 64 DOI 10.1007/s00184-007-0142-7 Characterizations of probability distribtions via bivariate regression of record vales George P. Yanev M. Ahsanllah M. I. Beg Received: 4 October 2006

More information

NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS

NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS COLLEEN ACKERMANN AND WILL HARDESTY Abstract. The Cahn-Morral System has often been sed to model the dynamics of phase separation in mlti-component

More information

SUBJECT:ENGINEERING MATHEMATICS-I SUBJECT CODE :SMT1101 UNIT III FUNCTIONS OF SEVERAL VARIABLES. Jacobians

SUBJECT:ENGINEERING MATHEMATICS-I SUBJECT CODE :SMT1101 UNIT III FUNCTIONS OF SEVERAL VARIABLES. Jacobians SUBJECT:ENGINEERING MATHEMATICS-I SUBJECT CODE :SMT0 UNIT III FUNCTIONS OF SEVERAL VARIABLES Jacobians Changing ariable is something e come across er oten in Integration There are man reasons or changing

More information

THE CATEGORY OF CGWH SPACES

THE CATEGORY OF CGWH SPACES THE CATEGORY OF CGWH SPACES N. P. STRICKLAND It is ell-knon that the category U of compactly generated eak Hasdorff spaces is a conenient setting for homotopy theory and algebraic topology. In this paper

More information

Gradient Projection Anti-windup Scheme on Constrained Planar LTI Systems. Justin Teo and Jonathan P. How

Gradient Projection Anti-windup Scheme on Constrained Planar LTI Systems. Justin Teo and Jonathan P. How 1 Gradient Projection Anti-windp Scheme on Constrained Planar LTI Systems Jstin Teo and Jonathan P. How Technical Report ACL1 1 Aerospace Controls Laboratory Department of Aeronatics and Astronatics Massachsetts

More information

Chords in Graphs. Department of Mathematics Texas State University-San Marcos San Marcos, TX Haidong Wu

Chords in Graphs. Department of Mathematics Texas State University-San Marcos San Marcos, TX Haidong Wu AUSTRALASIAN JOURNAL OF COMBINATORICS Volme 32 (2005), Pages 117 124 Chords in Graphs Weizhen G Xingde Jia Department of Mathematics Texas State Uniersity-San Marcos San Marcos, TX 78666 Haidong W Department

More information

MAT389 Fall 2016, Problem Set 6

MAT389 Fall 2016, Problem Set 6 MAT389 Fall 016, Problem Set 6 Trigonometric and hperbolic fnctions 6.1 Show that e iz = cos z + i sin z for eer comple nmber z. Hint: start from the right-hand side and work or wa towards the left-hand

More information

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2 MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,

More information

Change of Variables. (f T) JT. f = U

Change of Variables. (f T) JT. f = U Change of Variables 4-5-8 The change of ariables formla for mltiple integrals is like -sbstittion for single-ariable integrals. I ll gie the general change of ariables formla first, and consider specific

More information

Restricted cycle factors and arc-decompositions of digraphs. J. Bang-Jensen and C. J. Casselgren

Restricted cycle factors and arc-decompositions of digraphs. J. Bang-Jensen and C. J. Casselgren Restricted cycle factors and arc-decompositions of digraphs J. Bang-Jensen and C. J. Casselgren REPORT No. 0, 0/04, spring ISSN 0-467X ISRN IML-R- -0-/4- -SE+spring Restricted cycle factors and arc-decompositions

More information

The Scalar Conservation Law

The Scalar Conservation Law The Scalar Conservation Law t + f() = 0 = conserved qantity, f() =fl d dt Z b a (t, ) d = Z b a t (t, ) d = Z b a f (t, ) d = f (t, a) f (t, b) = [inflow at a] [otflow at b] f((a)) f((b)) a b Alberto Bressan

More information

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation:

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation: Math 263 Assignment #3 Soltions 1. A fnction z f(x, ) is called harmonic if it satisfies Laplace s eqation: 2 + 2 z 2 0 Determine whether or not the following are harmonic. (a) z x 2 + 2. We se the one-variable

More information

We automate the bivariate change-of-variables technique for bivariate continuous random variables with

We automate the bivariate change-of-variables technique for bivariate continuous random variables with INFORMS Jornal on Compting Vol. 4, No., Winter 0, pp. 9 ISSN 09-9856 (print) ISSN 56-558 (online) http://dx.doi.org/0.87/ijoc.046 0 INFORMS Atomating Biariate Transformations Jeff X. Yang, John H. Drew,

More information

Part II. Martingale measres and their constrctions 1. The \First" and the \Second" fndamental theorems show clearly how \mar tingale measres" are impo

Part II. Martingale measres and their constrctions 1. The \First and the \Second fndamental theorems show clearly how \mar tingale measres are impo Albert N. Shiryaev (Stelov Mathematical Institte and Moscow State University) ESSENTIALS of the ARBITRAGE THEORY Part I. Basic notions and theorems of the \Arbitrage Theory" Part II. Martingale measres

More information

Krauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields.

Krauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields. Kraskopf, B, Lee,, & Osinga, H (28) odimension-one tangency bifrcations of global Poincaré maps of for-dimensional vector fields Early version, also known as pre-print Link to pblication record in Explore

More information

Complex Variables. For ECON 397 Macroeconometrics Steve Cunningham

Complex Variables. For ECON 397 Macroeconometrics Steve Cunningham Comple Variables For ECON 397 Macroeconometrics Steve Cnningham Open Disks or Neighborhoods Deinition. The set o all points which satis the ineqalit

More information

Upper Bounds on the Spanning Ratio of Constrained Theta-Graphs

Upper Bounds on the Spanning Ratio of Constrained Theta-Graphs Upper Bonds on the Spanning Ratio of Constrained Theta-Graphs Prosenjit Bose and André van Renssen School of Compter Science, Carleton University, Ottaa, Canada. jit@scs.carleton.ca, andre@cg.scs.carleton.ca

More information

B-469 Simplified Copositive and Lagrangian Relaxations for Linearly Constrained Quadratic Optimization Problems in Continuous and Binary Variables

B-469 Simplified Copositive and Lagrangian Relaxations for Linearly Constrained Quadratic Optimization Problems in Continuous and Binary Variables B-469 Simplified Copositive and Lagrangian Relaxations for Linearly Constrained Qadratic Optimization Problems in Continos and Binary Variables Naohiko Arima, Snyong Kim and Masakaz Kojima October 2012,

More information

Axiomatizing the Cyclic Interval Calculus

Axiomatizing the Cyclic Interval Calculus Axiomatizing the Cyclic Interal Calcls Jean-François Condotta CRIL-CNRS Uniersité d Artois 62300 Lens (France) condotta@cril.ni-artois.fr Gérard Ligozat LIMSI-CNRS Uniersité de Paris-Sd 91403 Orsay (France)

More information

The Minimal Estrada Index of Trees with Two Maximum Degree Vertices

The Minimal Estrada Index of Trees with Two Maximum Degree Vertices MATCH Commnications in Mathematical and in Compter Chemistry MATCH Commn. Math. Compt. Chem. 64 (2010) 799-810 ISSN 0340-6253 The Minimal Estrada Index of Trees with Two Maximm Degree Vertices Jing Li

More information

Distributed Weighted Vertex Cover via Maximal Matchings

Distributed Weighted Vertex Cover via Maximal Matchings Distribted Weighted Vertex Coer ia Maximal Matchings FABRIZIO GRANDONI Uniersità di Roma Tor Vergata JOCHEN KÖNEMANN Uniersity of Waterloo and ALESSANDRO PANCONESI Sapienza Uniersità di Roma In this paper

More information

The Linear Quadratic Regulator

The Linear Quadratic Regulator 10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.

More information

Direct linearization method for nonlinear PDE s and the related kernel RBFs

Direct linearization method for nonlinear PDE s and the related kernel RBFs Direct linearization method for nonlinear PDE s and the related kernel BFs W. Chen Department of Informatics, Uniersity of Oslo, P.O.Box 1080, Blindern, 0316 Oslo, Norway Email: wenc@ifi.io.no Abstract

More information

Predicting Popularity of Twitter Accounts through the Discovery of Link-Propagating Early Adopters

Predicting Popularity of Twitter Accounts through the Discovery of Link-Propagating Early Adopters Predicting Poplarity of Titter Acconts throgh the Discoery of Link-Propagating Early Adopters Daichi Imamori Gradate School of Informatics, Kyoto Uniersity Sakyo, Kyoto 606-850 Japan imamori@dl.soc.i.kyoto-.ac.jp

More information

2 Faculty of Mechanics and Mathematics, Moscow State University.

2 Faculty of Mechanics and Mathematics, Moscow State University. th World IMACS / MODSIM Congress, Cairns, Astralia 3-7 Jl 9 http://mssanz.org.a/modsim9 Nmerical eamination of competitie and predator behaior for the Lotka-Volterra eqations with diffsion based on the

More information

Setting The K Value And Polarization Mode Of The Delta Undulator

Setting The K Value And Polarization Mode Of The Delta Undulator LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions

More information

On Minima of Discrimination Functions

On Minima of Discrimination Functions On Minima of Discrimination Fnctions Janne V. Kjala Uniersity of Jyäskylä Ehtibar N. Dzhafaro Prde Uniersity Abstract A discrimination fnction (x; y) assigns a measre of discriminability to stimls pairs

More information

Non-Lecture I: Linear Programming. Th extremes of glory and of shame, Like east and west, become the same.

Non-Lecture I: Linear Programming. Th extremes of glory and of shame, Like east and west, become the same. The greatest flood has the soonest ebb; the sorest tempest the most sdden calm; the hottest love the coldest end; and from the deepest desire oftentimes enses the deadliest hate. Th extremes of glory and

More information

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1 y g j P(t) P(,y) r t0 i 4/4/006 Motion ( F.Robilliard) 1 Motion: We stdy in detail three cases of motion: 1. Motion in one dimension with constant acceleration niform linear motion.. Motion in two dimensions

More information

A scalar nonlocal bifurcation of solitary waves for coupled nonlinear Schrödinger systems

A scalar nonlocal bifurcation of solitary waves for coupled nonlinear Schrödinger systems INSTITUTE OF PHYSICS PUBLISHING Nonlinearity 5 (22) 265 292 NONLINEARITY PII: S95-775(2)349-4 A scalar nonlocal bifrcation of solitary waes for copled nonlinear Schrödinger systems Alan R Champneys and

More information

A Survey of the Implementation of Numerical Schemes for Linear Advection Equation

A Survey of the Implementation of Numerical Schemes for Linear Advection Equation Advances in Pre Mathematics, 4, 4, 467-479 Pblished Online Agst 4 in SciRes. http://www.scirp.org/jornal/apm http://dx.doi.org/.436/apm.4.485 A Srvey of the Implementation of Nmerical Schemes for Linear

More information

Nonlinear parametric optimization using cylindrical algebraic decomposition

Nonlinear parametric optimization using cylindrical algebraic decomposition Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic

More information

Complex Tire-Ground Interaction Simulation: Recent Developments Of An Advanced Shell Theory Based Tire Model

Complex Tire-Ground Interaction Simulation: Recent Developments Of An Advanced Shell Theory Based Tire Model . ozdog and W. W. Olson Complex Tire-Grond Interaction Simlation: ecent eelopments Of n danced Shell Theory ased Tire odel EFEECE: ozdog. and Olson W. W. Complex Tire-Grond Interaction Simlation: ecent

More information

MATH2715: Statistical Methods

MATH2715: Statistical Methods MATH275: Statistical Methods Exercises III (based on lectres 5-6, work week 4, hand in lectre Mon 23 Oct) ALL qestions cont towards the continos assessment for this modle. Q. If X has a niform distribtion

More information

3.3 Operations With Vectors, Linear Combinations

3.3 Operations With Vectors, Linear Combinations Operations With Vectors, Linear Combinations Performance Criteria: (d) Mltiply ectors by scalars and add ectors, algebraically Find linear combinations of ectors algebraically (e) Illstrate the parallelogram

More information

SECTION 6.7. The Dot Product. Preview Exercises. 754 Chapter 6 Additional Topics in Trigonometry. 7 w u 7 2 =?. 7 v 77w7

SECTION 6.7. The Dot Product. Preview Exercises. 754 Chapter 6 Additional Topics in Trigonometry. 7 w u 7 2 =?. 7 v 77w7 754 Chapter 6 Additional Topics in Trigonometry 115. Yo ant to fly yor small plane de north, bt there is a 75-kilometer ind bloing from est to east. a. Find the direction angle for here yo shold head the

More information

Essentials of optimal control theory in ECON 4140

Essentials of optimal control theory in ECON 4140 Essentials of optimal control theory in ECON 4140 Things yo need to know (and a detail yo need not care abot). A few words abot dynamic optimization in general. Dynamic optimization can be thoght of as

More information

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018 Lectre 3 The dot prodct Dan Nichols nichols@math.mass.ed MATH 33, Spring 018 Uniersity of Massachsetts Janary 30, 018 () Last time: 3D space Right-hand rle, the three coordinate planes 3D coordinate system:

More information

A Single Species in One Spatial Dimension

A Single Species in One Spatial Dimension Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,

More information

Minimizing Intra-Edge Crossings in Wiring Diagrams and Public Transportation Maps

Minimizing Intra-Edge Crossings in Wiring Diagrams and Public Transportation Maps Minimizing Intra-Edge Crossings in Wiring Diagrams and Pblic Transportation Maps Marc Benkert 1, Martin Nöllenbrg 1, Takeaki Uno 2, and Alexander Wolff 1 1 Department of Compter Science, Karlsrhe Uniersity,

More information

ECE Notes 4 Functions of a Complex Variable as Mappings. Fall 2017 David R. Jackson. Notes are adapted from D. R. Wilton, Dept.

ECE Notes 4 Functions of a Complex Variable as Mappings. Fall 2017 David R. Jackson. Notes are adapted from D. R. Wilton, Dept. ECE 638 Fall 017 Daid R. Jackson Notes 4 Fnctions of a Comple Variable as Mappings Notes are adapted from D. R. Wilton, Dept. of ECE 1 A Fnction of a Comple Variable as a Mapping A fnction of a comple

More information

STURM-LIOUVILLE PROBLEMS

STURM-LIOUVILLE PROBLEMS STURM-LIOUVILLE PROBLEMS ANTON ZETTL Mathematics Department, Northern Illinois University, DeKalb, Illinois 60115. Dedicated to the memory of John Barrett. ABSTRACT. Reglar and singlar Strm-Lioville problems

More information

EE2 Mathematics : Functions of Multiple Variables

EE2 Mathematics : Functions of Multiple Variables EE2 Mathematics : Fnctions of Mltiple Variables http://www2.imperial.ac.k/ nsjones These notes are not identical word-for-word with m lectres which will be gien on the blackboard. Some of these notes ma

More information

Existence of HCMU metrics on higher genus Riemann surfaces

Existence of HCMU metrics on higher genus Riemann surfaces Existence of HCMU metrics on higher gens Riemann srfaces October 4, 0 bstract We will generalize the reslt in [CCW05] and roe the existence of HCMU metrics on higher gens K-srfaces, i.e. Riemann srfaces

More information

CRITERIA FOR TOEPLITZ OPERATORS ON THE SPHERE. Jingbo Xia

CRITERIA FOR TOEPLITZ OPERATORS ON THE SPHERE. Jingbo Xia CRITERIA FOR TOEPLITZ OPERATORS ON THE SPHERE Jingbo Xia Abstract. Let H 2 (S) be the Hardy space on the nit sphere S in C n. We show that a set of inner fnctions Λ is sfficient for the prpose of determining

More information

Math 4A03: Practice problems on Multivariable Calculus

Math 4A03: Practice problems on Multivariable Calculus Mat 4A0: Practice problems on Mltiariable Calcls Problem Consider te mapping f, ) : R R defined by fx, y) e y + x, e x y) x, y) R a) Is it possible to express x, y) as a differentiable fnction of, ) near

More information

MAXIMUM AND ANTI-MAXIMUM PRINCIPLES FOR THE P-LAPLACIAN WITH A NONLINEAR BOUNDARY CONDITION. 1. Introduction. ν = λ u p 2 u.

MAXIMUM AND ANTI-MAXIMUM PRINCIPLES FOR THE P-LAPLACIAN WITH A NONLINEAR BOUNDARY CONDITION. 1. Introduction. ν = λ u p 2 u. 2005-Ojda International Conference on Nonlinear Analysis. Electronic Jornal of Differential Eqations, Conference 14, 2006, pp. 95 107. ISSN: 1072-6691. URL: http://ejde.math.txstate.ed or http://ejde.math.nt.ed

More information

Model (In-)Validation from a H and µ perspective

Model (In-)Validation from a H and µ perspective Model (In-)Validation from a H and µ perspectie Wolfgang Reinelt Department of Electrical Engineering Linköping Uniersity, S-581 83 Linköping, Seden WWW: http://.control.isy.li.se/~olle/ Email: olle@isy.li.se

More information

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching

More information

1. Tractable and Intractable Computational Problems So far in the course we have seen many problems that have polynomial-time solutions; that is, on

1. Tractable and Intractable Computational Problems So far in the course we have seen many problems that have polynomial-time solutions; that is, on . Tractable and Intractable Comptational Problems So far in the corse we have seen many problems that have polynomial-time soltions; that is, on a problem instance of size n, the rnning time T (n) = O(n

More information

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises LECTURE 1 State-Space Linear Sstems This lectre introdces state-space linear sstems, which are the main focs of this book. Contents 1. State-Space Linear Sstems 2. Block Diagrams 3. Exercises 1.1 State-Space

More information

Worst-case analysis of the LPT algorithm for single processor scheduling with time restrictions

Worst-case analysis of the LPT algorithm for single processor scheduling with time restrictions OR Spectrm 06 38:53 540 DOI 0.007/s009-06-043-5 REGULAR ARTICLE Worst-case analysis of the LPT algorithm for single processor schedling with time restrictions Oliver ran Fan Chng Ron Graham Received: Janary

More information

Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Porous Rotating Disk

Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Porous Rotating Disk Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Poros otating Disk Hemant Poonia * * Asstt. Prof., Deptt. of Math, Stat & Physics, CCSHAU, Hisar-54.. C. Chadhary etd. Professor, Deptt.

More information

ON THE SHAPES OF BILATERAL GAMMA DENSITIES

ON THE SHAPES OF BILATERAL GAMMA DENSITIES ON THE SHAPES OF BILATERAL GAMMA DENSITIES UWE KÜCHLER, STEFAN TAPPE Abstract. We investigate the for parameter family of bilateral Gamma distribtions. The goal of this paper is to provide a thorogh treatment

More information

A new integral transform on time scales and its applications

A new integral transform on time scales and its applications Agwa et al. Advances in Difference Eqations 202, 202:60 http://www.advancesindifferenceeqations.com/content/202//60 RESEARCH Open Access A new integral transform on time scales and its applications Hassan

More information

Applied Mathematics Letters

Applied Mathematics Letters Applied Mathematics Letters 23 (21 49 416 Contents lists available at ScienceDirect Applied Mathematics Letters jornal homepage: www.elsevier.com/locate/aml Exponential trichotomy and homoclinic bifrcation

More information

Elements of Coordinate System Transformations

Elements of Coordinate System Transformations B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and

More information

2 PTR W. MICHOR So we see that for xed (y; v) the transition fnctions are linear in (; w) 2 R m V. This describes the vector bndle strctre of the tang

2 PTR W. MICHOR So we see that for xed (y; v) the transition fnctions are linear in (; w) 2 R m V. This describes the vector bndle strctre of the tang TH JACOBI FLOW Peter W. Michor rwin Schrodinger Institt fr Mathematische Physik, Pastergasse 6/7, A-1090 Wien, Astria November 19, 1996 For Wlodek Tlczyjew, on the occasion of his 65th birthday. It is

More information

THE HOHENBERG-KOHN THEOREM FOR MARKOV SEMIGROUPS

THE HOHENBERG-KOHN THEOREM FOR MARKOV SEMIGROUPS THE HOHENBERG-KOHN THEOREM FOR MARKOV SEMIGROUPS OMAR HIJAB Abstract. At the basis of mch of comptational chemistry is density fnctional theory, as initiated by the Hohenberg-Kohn theorem. The theorem

More information

Intro to path analysis Richard Williams, University of Notre Dame, https://www3.nd.edu/~rwilliam/ Last revised April 6, 2015

Intro to path analysis Richard Williams, University of Notre Dame, https://www3.nd.edu/~rwilliam/ Last revised April 6, 2015 Intro to path analysis Richard Williams, Uniersity of Notre Dame, https://3.nd.ed/~rilliam/ Last reised April 6, 05 Sorces. This discssion dras heaily from Otis Ddley Dncan s Introdction to Strctral Eqation

More information

Discussion Papers Department of Economics University of Copenhagen

Discussion Papers Department of Economics University of Copenhagen Discssion Papers Department of Economics University of Copenhagen No. 10-06 Discssion of The Forward Search: Theory and Data Analysis, by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen,

More information

Figure 1 Probability density function of Wedge copula for c = (best fit to Nominal skew of DRAM case study).

Figure 1 Probability density function of Wedge copula for c = (best fit to Nominal skew of DRAM case study). Wedge Copla This docment explains the constrction and properties o a particlar geometrical copla sed to it dependency data rom the edram case stdy done at Portland State University. The probability density

More information

Sign-reductions, p-adic valuations, binomial coefficients modulo p k and triangular symmetries

Sign-reductions, p-adic valuations, binomial coefficients modulo p k and triangular symmetries Sign-redctions, p-adic valations, binomial coefficients modlo p k and trianglar symmetries Mihai Prnesc Abstract According to a classical reslt of E. Kmmer, the p-adic valation v p applied to a binomial

More information

Computational Geosciences 2 (1998) 1, 23-36

Computational Geosciences 2 (1998) 1, 23-36 A STUDY OF THE MODELLING ERROR IN TWO OPERATOR SPLITTING ALGORITHMS FOR POROUS MEDIA FLOW K. BRUSDAL, H. K. DAHLE, K. HVISTENDAHL KARLSEN, T. MANNSETH Comptational Geosciences 2 (998), 23-36 Abstract.

More information

J. Klein, R. Schwanzl, R. M. Vogt. It has been known for a long time that an (n? 1)-connected CW -complex

J. Klein, R. Schwanzl, R. M. Vogt. It has been known for a long time that an (n? 1)-connected CW -complex Co-Mltiplication and Sspension J. Klein, R. Schanzl, R. M. Vogt 1 Introdction It has been knon for a long time that an (n? 1)-connected CW -complex of dimension 2n? 1 has the homotopy type of a sspension.

More information

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE

More information

Spring, 2008 CIS 610. Advanced Geometric Methods in Computer Science Jean Gallier Homework 1, Corrected Version

Spring, 2008 CIS 610. Advanced Geometric Methods in Computer Science Jean Gallier Homework 1, Corrected Version Spring, 008 CIS 610 Adanced Geometric Methods in Compter Science Jean Gallier Homework 1, Corrected Version Febrary 18, 008; De March 5, 008 A problems are for practice only, and shold not be trned in.

More information

Decoder Error Probability of MRD Codes

Decoder Error Probability of MRD Codes Decoder Error Probability of MRD Codes Maximilien Gadolea Department of Electrical and Compter Engineering Lehigh University Bethlehem, PA 18015 USA E-mail: magc@lehighed Zhiyan Yan Department of Electrical

More information

Discussion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli

Discussion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli 1 Introdction Discssion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen Department of Economics, University of Copenhagen and CREATES,

More information

Change of Variables. f(x, y) da = (1) If the transformation T hasn t already been given, come up with the transformation to use.

Change of Variables. f(x, y) da = (1) If the transformation T hasn t already been given, come up with the transformation to use. MATH 2Q Spring 26 Daid Nichols Change of Variables Change of ariables in mltiple integrals is complicated, bt it can be broken down into steps as follows. The starting point is a doble integral in & y.

More information

Can a traveling wave connect two unstable states? The case of the nonlocal Fisher equation

Can a traveling wave connect two unstable states? The case of the nonlocal Fisher equation Can a traeling wae connect two nstable states? The case of the nonlocal Fisher eqation Grégoire Nadin, Benoît Perthame, Min Tang To cite this ersion: Grégoire Nadin, Benoît Perthame, Min Tang. Can a traeling

More information

more detailed way in comparison with two others. The system -- (b) seems to be nstable with respect to transformation to the - congration (a) [5]. The

more detailed way in comparison with two others. The system -- (b) seems to be nstable with respect to transformation to the - congration (a) [5]. The For varios string baryon models and Regge trajectories G.S. Sharov Tver state niversity, 17000, Sadovyj per., 35, Tver, Rssia Abstract Rotational motions for the ark-diark (-), linear (--) \threestring"

More information

DYNAMICAL LOWER BOUNDS FOR 1D DIRAC OPERATORS. 1. Introduction We consider discrete, resp. continuous, Dirac operators

DYNAMICAL LOWER BOUNDS FOR 1D DIRAC OPERATORS. 1. Introduction We consider discrete, resp. continuous, Dirac operators DYNAMICAL LOWER BOUNDS FOR D DIRAC OPERATORS ROBERTO A. PRADO AND CÉSAR R. DE OLIVEIRA Abstract. Qantm dynamical lower bonds for continos and discrete one-dimensional Dirac operators are established in

More information

Simpler Testing for Two-page Book Embedding of Partitioned Graphs

Simpler Testing for Two-page Book Embedding of Partitioned Graphs Simpler Testing for Two-page Book Embedding of Partitioned Graphs Seok-Hee Hong 1 Hiroshi Nagamochi 2 1 School of Information Technologies, Uniersity of Sydney, seokhee.hong@sydney.ed.a 2 Department of

More information

arxiv: v1 [physics.flu-dyn] 11 Mar 2011

arxiv: v1 [physics.flu-dyn] 11 Mar 2011 arxiv:1103.45v1 [physics.fl-dyn 11 Mar 011 Interaction of a magnetic dipole with a slowly moving electrically condcting plate Evgeny V. Votyakov Comptational Science Laboratory UCY-CompSci, Department

More information

Introdction Finite elds play an increasingly important role in modern digital commnication systems. Typical areas of applications are cryptographic sc

Introdction Finite elds play an increasingly important role in modern digital commnication systems. Typical areas of applications are cryptographic sc A New Architectre for a Parallel Finite Field Mltiplier with Low Complexity Based on Composite Fields Christof Paar y IEEE Transactions on Compters, Jly 996, vol 45, no 7, pp 856-86 Abstract In this paper

More information

arxiv: v2 [cs.si] 27 Apr 2017

arxiv: v2 [cs.si] 27 Apr 2017 Opinion Dynamics in Networks: Conergence, Stability and Lack of Explosion arxi:1607.038812 [cs.si] 27 Apr 2017 Tng Mai Georgia Institte of Technology maitng89@gatech.ed Vijay V. Vazirani Georgia Institte

More information

Geometry of Span (continued) The Plane Spanned by u and v

Geometry of Span (continued) The Plane Spanned by u and v Geometric Description of Span Geometr of Span (contined) 2 Geometr of Span (contined) 2 Span {} Span {, } 2 Span {} 2 Geometr of Span (contined) 2 b + 2 The Plane Spanned b and If a plane is spanned b

More information

The Cryptanalysis of a New Public-Key Cryptosystem based on Modular Knapsacks

The Cryptanalysis of a New Public-Key Cryptosystem based on Modular Knapsacks The Cryptanalysis of a New Pblic-Key Cryptosystem based on Modlar Knapsacks Yeow Meng Chee Antoine Jox National Compter Systems DMI-GRECC Center for Information Technology 45 re d Ulm 73 Science Park Drive,

More information

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they

More information

Chapter 2 Difficulties associated with corners

Chapter 2 Difficulties associated with corners Chapter Difficlties associated with corners This chapter is aimed at resolving the problems revealed in Chapter, which are cased b corners and/or discontinos bondar conditions. The first section introdces

More information

Study of the diffusion operator by the SPH method

Study of the diffusion operator by the SPH method IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait

More information

MEAN VALUE ESTIMATES OF z Ω(n) WHEN z 2.

MEAN VALUE ESTIMATES OF z Ω(n) WHEN z 2. MEAN VALUE ESTIMATES OF z Ωn WHEN z 2 KIM, SUNGJIN 1 Introdction Let n i m pei i be the prime factorization of n We denote Ωn by i m e i Then, for any fixed complex nmber z, we obtain a completely mltiplicative

More information

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer 6. Similarit Soltions Chapter 6 Momentm Transfer in an Eternal Laminar Bondar Laer Consider a laminar incompressible bondar laer with constant properties. Assme the flow is stead and two-dimensional aligned

More information

Online Stochastic Matching: New Algorithms and Bounds

Online Stochastic Matching: New Algorithms and Bounds Online Stochastic Matching: New Algorithms and Bonds Brian Brbach, Karthik A. Sankararaman, Araind Sriniasan, and Pan X Department of Compter Science, Uniersity of Maryland, College Park, MD 20742, USA

More information

Decoder Error Probability of MRD Codes

Decoder Error Probability of MRD Codes Decoder Error Probability of MRD Codes Maximilien Gadolea Department of Electrical and Compter Engineering Lehigh University Bethlehem, PA 18015 USA E-mail: magc@lehigh.ed Zhiyan Yan Department of Electrical

More information

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

More information

A Theory of Markovian Time Inconsistent Stochastic Control in Discrete Time

A Theory of Markovian Time Inconsistent Stochastic Control in Discrete Time A Theory of Markovian Time Inconsistent Stochastic Control in Discrete Time Tomas Björk Department of Finance, Stockholm School of Economics tomas.bjork@hhs.se Agatha Mrgoci Department of Economics Aarhs

More information

A Regulator for Continuous Sedimentation in Ideal Clarifier-Thickener Units

A Regulator for Continuous Sedimentation in Ideal Clarifier-Thickener Units A Reglator for Continos Sedimentation in Ideal Clarifier-Thickener Units STEFAN DIEHL Centre for Mathematical Sciences, Lnd University, P.O. Box, SE- Lnd, Sweden e-mail: diehl@maths.lth.se) Abstract. The

More information