On the Number of Limit Cycles for Discontinuous Generalized Liénard Polynomial Differential Systems

Size: px
Start display at page:

Download "On the Number of Limit Cycles for Discontinuous Generalized Liénard Polynomial Differential Systems"

Transcription

1 Internationa Journa of Bifurcation and Chaos Vo. 25 No pages c Word Scientific Pubishing Company DOI: /S X On the Number of Limit Cyces for Discontinuous Generaized Liénard Poynomia Differentia Systems Int. J. Bifurcation Chaos Downoaded from by CITY UNIVERSITY OF HONG KONG on 09/29/15. For persona use ony. Fangfang Jiang Department of Mathematics Tongi University Shanghai P. R. China iangfangfang87@126.com Junping Shi Department of Mathematics Coege of Wiiam and Mary Wiiamsburg VA USA shi@math.wm.edu Jitao Sun Department of Mathematics Tongi University Shanghai P. R. China sunt@sh163.net Received November ; Revised March In this paper we investigate the number of imit cyces for a cass of discontinuous panar differentia systems with mutipe sectors separated by many rays originating from the origin. In each sector it is a smooth generaized Liénard poynomia differentia system x = y g 1 xf 1 xy and y = x g 2 xf 2 xy wheref i x andg i x fori =1 2arepoynomias of variabe x with any given degree. By the averaging theory of first-order for discontinuous differentia systems we provide the criteria on the maximum number of medium ampitude imit cyces for the discontinuous generaized Liénard poynomia differentia systems. The upper bound for the number of medium ampitude imit cyces can be attained by specific exampes. Keywords: Discontinuous panar system; Liénard poynomia system; averaging theory; number of imit cyces. 1. Introduction One of main topics of the quaitative theory of panar differentia systems is the number and reative position of imit cyces. This probem restricted to continuous panar poynomia differentia systems is the we known Hibert s 16th probem see for exampe Li 2003]. Up to now there have been many achievements concerning the existence uniqueness and the number of imit cyces see for exampe García et a. 201; Justino & Jorge 2012; Li & Libre 2012; Libre 2010; Libre & Mereu ; Libre et a. 2015; Libre & Vas a 2013b; Loyd & Lynch 1988; Martins & Mereu 201; Shen & Han 2013; Sun 1992; Xiong & Zhong 2013] and references therein. A imit cyce bifurcating from a singe degenerate singuar point is caed a sma ampitude imit cyce and the one bifurcating from periodic orbits of a inear center is caed a medium ampitude imit cyce. In Libre & Vas 2012] the authors studied the number Author for correspondence

2 Int. J. Bifurcation Chaos Downoaded from by CITY UNIVERSITY OF HONG KONG on 09/29/15. For persona use ony. F. Jiang et a. of medium ampitude imit cyces for a cass of poynomia differentia systems of the form x = y g 1 x f 1 xy 1 y = x g 2 x f 2 xy where f i x and g i x i = 1 2 are poynomias of a given degree. For system 1 when f 1 x = g 1 x = 0 it is the generaized Liénard poynomia differentia system. Note that the cassica generaized Liénard poynomia differentia equation is x fxx gx = 0 or equivaenty x = y F x y = gx 2 where F x = x 0 fsds fx andgx arepoynomias in the variabe x. For system 2 when poynomias F x and gx have degrees n and m respectivey the generaization to discontinuous generaized Liénard poynomia differentia systems was studied in Libre & Mereu 2013]. To the best of our knowedge there are very few resuts concerning the number of imit cyces which bifurcate from a continuum of periodic orbits of differentia systems perturbed by discontinuous Liénard poynomia differentia systems. So we restrict our attention to a cass of discontinuous and piecewise smooth panar poynomia differentia systems with mutipe sectors separated by rays starting from the origin. In this paper we investigate the maximum number of medium ampitude imit cyces for a discontinuous generaized Liénard poynomia differentia system with even number sectors separated by /2 straight ines passing through the origin 0 0 where the discontinuous set is of the form /2 1 x y :y =tan α 2kπ } x. k=0 We first consider the case where two different smooth generaized Liénard poynomia differentia systems are defined aternativey in every two neighboring sectors on the pane. By appying the averaging theory of first-order for discontinuous differentia systems we obtain a ower bound for the maximum number of imit cyces and we show that the upper bound for the number of medium ampitude imit cyces can be attained by some constructed exampes see Exampe.1 for an exampe and Matab numerica simuation. The obtained imit cyces bifurcate from the periodic orbits of the inear center x = y y = x. Moreover we give the discussion on the resut if there are different vector fieds in different sectors respectivey. The paper is organized as foows. In the next section we present some reevant preiminaries for the discontinuous generaized Liénard poynomia differentia system. In Sec. 3 we first present the averaging theory of first-order for discontinuous differentia systems and then we prove our main resuts on the maximum number of medium ampitude imit cyces. In Sec. we give two exampes to iustrate the proved resuts. Concusion is presented in Sec Preiminaries Consider the foowing discontinuous generaized Liénard poynomia differentia system x = y εg 11 xf 11 xy y = x εg 12 xf 12 xy hx y > 0 x = y εg 21 xf 21 xy y = x εg 22 xf 22 xy hx y < 0 3 where ε is a sufficienty sma parameter and g i x and f i x fori =1 2 are poynomias of variabe x with degrees m and n respectivey. Let a function h : R 2 R be given by 2 1 hx y = k=0 y tan α 2kπ x with >1 being an even integer and α R. Then h 1 0 = 2 1 k=0 x y :y =tan α 2kπ } x 5 is a discontinuous set consisting of /2 straight ines passing through the origin O0 0 which divides the pane R 2 into sectors with ange 2π/. When ε = 0 the system 3 becomes a panar inear center as foows x = y y = x. 6 Note that the system 6 has a continuum of periodic orbits surrounding the origin and the orbits are countercockwise oriented in the phase pane. In this paper by using the averaging theory of first-order for discontinuous differentia systems we investigate the maximum number of medium ampitude imit cyces for the discontinuous generaized

3 The Number of Limit Cyces for Discontinuous Systems Int. J. Bifurcation Chaos Downoaded from by CITY UNIVERSITY OF HONG KONG on 09/29/15. For persona use ony. Liénard poynomia differentia system 3 with ε > 0 sufficienty sma. 3. Main Resuts First we reca the foowing resut concerning the existence of imit cyces from the averaging theory of first-order for discontinuous differentia systems. Lemma 3.1 Libre & Mereu 201]. Consider the foowing discontinuous differentia equation x t =εht xε 2 Gt x ε 7 where ε is a sufficienty sma rea number functions Ht xgt x ε are given by Ht x =H 1 t xsgnht xh 2 t x Gt x ε =G 1 t x εsgnht xg 2 t x ε functions H 1 H 2 : R D R n G 1 G 2 : R D R R n and h : R D R are continuous T- periodic in the variabe t and D R n is an open subset. Let h be a C 1 function having 0 as a reguar vaue. Denote by M = h 1 0 N = 0} D M and K = N\M. Moreover define the averaging function d : D R n of 7 to be Assume that dx = T 0 Ht xdt. i functions H 1 H 2 G 1 G 2 and h are ocay Lipschitz with respect to the variabe x; ii for z 0 K with dz 0 = 0 there exists a neighborhood V of z 0 such that dz 0for every z V \z 0 } and d B d V z 0 0 where d B d V z 0 denotes the Brouwer degree of the function d with respect to the bounded open subset V D and the fixed point z 0 ; iii if h t t 0z 0 =0for some t 0 z 0 M then x h H 1 ] 2 x h H 2 t 0 z 0 > 0 where A B] denotes the inner product of the vectors A and B x h denotes the gradient of ht x with respect to the variabe x. Then for 0 < ε sufficienty sma there exists a T -periodic soution x ε of 7 such that x0ε z 0 as ε 0. By taking the ange θ as the new independent variabe and using Tayor series expansion the sys- Our main resut of this paper is as foows. tem 3 becomes m dρ dθ = ε a 1i ρ i cos i1 θ c 1i ρ i1 cos i1 θ sin θ b 1i ρ i cos i θ sin θ Oε 2 Theorem 3.2. Consider the system 3 with ε > 0 sufficienty sma. Assume that the poynomias g i x and f i x i=1 2 have degrees m 1 and n 1 respectivey then the maximum number of medium ampitude imit cyces bifurcating from the periodic orbits of system 6 is λ =maxm n 1}. Proof. Let x = ρ cos θ y = ρ sin θ with ρ>0 then the system 3 can be transformed into dρ dt = εcos θg 11ρ cos θf 11 ρ cos θρ sin θ] sinθg 12 ρ cos θf 12 ρ cos θρ sin θ]} dθ dt =1ε ρ cos θg 12ρ cos θf 12 ρ cos θρ sin θ] sin θg 11 ρ cos θf 11 ρ cos θρ sin θ]} for θ α kπ α k2π k = and dρ dt = εcos θg 21ρ cos θf 21 ρ cos θρ sin θ] sinθg 22 ρ cos θf 22 ρ cos θρ sin θ]} dθ dt =1ε ρ cos θg 22ρ cos θf 22 ρ cos θρ sin θ] sin θg 21 ρ cos θf 21 ρ cos θρ sin θ]} for θ α k2π α k1π k = Denote by g 11 x = g 21 x = f 11 x = f 21 x = a 1i x i g 12 x = a 2i x i g 22 x = c 1i x i f 12 x = c 2i x i f 22 x = b 1i x i b 2i x i d 1i x i d 2i x i. d 1i ρ i1 cos i θ sin 2 θ

4 F. Jiang et a. for θ α kπ α k2π k = and m dρ dθ = ε a 2i ρ i cos i1 θ c 2i ρ i1 cos i1 θ sin θ b 2i ρ i cos i θ sin θ d 2i ρ i1 cos i θ sin 2 θ Oε 2 9 Int. J. Bifurcation Chaos Downoaded from by CITY UNIVERSITY OF HONG KONG on 09/29/15. For persona use ony. for θ α k2π h 1 0 = h θ ρ = θ α k1π 2 1 k=0 2 1 p=0 k = From and 5 then θ ρ :sinθ =cosθtan α 2kπ } ρ 2 cos θ sinθtan α 2pπ ] 2 1 k p sin θ cos θ tan α 2kπ ]. By some tedious but simpe computations one has that h θ θ ρ 0foranyθ ρ h 1 0. On the other hand for the purpose of the averaging function we need to use the foowing formuas Gradshteyn & Ryshik 199] cos 2m θdθ = 1 2 2m cos 2m1 θdθ = 1 2 2m 2m θ 1 m 2 2m 1 m 1 2m sin2m 2θ 10 2m 2 2m 1 sin2m 2 1θ. 11 2m 2 1 Then the averaging function dρ of 8 and 9 is as foows where d 1 ρ = d 2 ρ = d 3 ρ = d ρ = 2 2 k=0 2 2 k=0 2 2 k=0 2 2 k=0 α k2π α kπ α k2π α kπ α k2π α kπ α k2π α kπ dρ =d 1 ρd 2 ρd 3 ρd ρ Substituting 10 and 11 into the above equaities then k1π α a 1i ρ i cos i1 θdθ a 2i ρ i cos i1 θdθ α k2π α k1π c 1i ρ i1 cos i1 θ sin θdθ c 2i ρ i1 cos i1 θ sin θdθ α k2π k1π α b 1i ρ i cos i θ sin θdθ b 2i ρ i cos i θ sin θdθ α k2π α k1π d 1i ρ i1 cos i θ sin 2 θdθ d 2i ρ i1 cos i θ sin 2 θdθ. α k2π d 1 ρ = 2 2 k=0 ρ i 2 i i i 1 a 1i A ik a 2i A ik m 1 a 12i1 a 22i1 ρ2i1 2 2i2 2i 2 i 1 2π

5 with A ik = A ik = sin i 2 1 α sin i 2 1 α The Number of Limit Cyces for Discontinuous Systems ] sin i 2 1 α kπ k 2π i 2 1 ] k 1π sin i 2 1 i 2 1 α ] ] k 2π Int. J. Bifurcation Chaos Downoaded from by CITY UNIVERSITY OF HONG KONG on 09/29/15. For persona use ony. with d 2 ρ = d 3 ρ = 2 2 k=0 2 2 k=0 B ik = B ik = C ik = C ik = c 1i ρ i1 B ik c 2i ρ i1 B ik b 1i ρ i C ik b 2i ρ i C ik cos i2 α cos i2 α cos i1 α cos i1 α k 2π i 2 k 1π k 1π cos i2 α kπ cos i2 α k 2π i 2 k 2π cos i1 α kπ i 1 cos i1 α i 1 n 1 2i d ρ = d 12i d 22i ρ 2i1 2 2i 1 i 2 2i2 2 2 k=0 2 2 k=0 ρ i1 2 i 1 ρ i1 2 i1 i i 1 i 2 i1 k 2π 2i 2 i 1 d 1i D ik d 2i D ik d 1i E ik d 2i E ik ] 2π with D ik = sin i 2 α k 2π i 2 ] sin i 2 α kπ ]

6 Int. J. Bifurcation Chaos Downoaded from by CITY UNIVERSITY OF HONG KONG on 09/29/15. For persona use ony. F. Jiang et a. E ik = D ik = E ik = sin i 2 2 α sin i 2 α Therefore we obtain that dρ = m 1 sin i 2 2 α 2 2 ] sin i 2 2 α kπ k 2π i 2 2 ] sin i 2 k 1π a 12i1 a 22i1 ρ 2i1 2i 2 2 2i2 i 1 k=0 2 2 k=0 2 2 k=0 ρ i 2 i i i 1 ρ i1 2 i 1 ρ i1 2 i1 i i 1 i 2 i1 From the above equaity we observe that the function dρ is a poynomia of variabe ρ with degree maxm n 1} so the poynomia dρ has at most λ = maxm n 1} positive roots that is the maximum number of bifurcating medium ampitude imit cyces is λ. Moreover if there exists some ρ 0 > 0 such that dρ 0 =0andd ρ 0 0 then the Brouwer degree d B d V ρ 0 0forsome sma open neighborhood V of ρ 0. Therefore by Lemma 3.1 for 0 < ε sufficienty sma the discontinuous system 3 has at most λ medium ampitude imit cyces bifurcating from the periodic orbits of the inear center x = y y = x. This competes the proof of Theorem 3.2. Remark 3.3. It shoud be noted that for some given degree m 1andn 1 there exist discontinuous generaized Liénard poynomia differentia systems having exacty λ =maxm n 1} medium ampitude imit cyces see Exampe.1 for an exampe with m = n =1. When α =0and =2wehave A ik = A ik =0 D ik = D ik =0 α i 2 ] k 1π sin i 2 2 i 2 2 2π n a 1i A ik a 2i A ik ] k 2π α d 12i d 22i ρ 2i1 2 2i1 i k=0 2 2 d 1i D ik d 2i D ik d 1i E ik d 2i E ik. k=0 ] ] k 2π 2i i 2π c 1i B ik c 2i B ik ρ i1 b 1i C ik b 2i C ik ρ i and E ik = E ik =0. In this case the function h : R 2 R is of the form hx y = y and the discontinuous generaized Liénard poynomia differentia system with two zones separated by the straight ine y = 0} is as foows x = y εg 11 xf 11 xy y = x εg 12 xf 12 xy y > 0 x = y εg 21 xf 21 xy y = x εg 22 xf 22 xy y < 0. We have the foowing resut:. 12 Coroary 3.. Consider the system 12 with 0 < ε sufficienty sma then the maximum number of medium ampitude imit cyces is λ =maxm n 1}. Proof. In the same way as Theorem 3.2 then the averaging function dρ isoftheform

7 dρ = = π 0 m a 1i ρ i cos i1 θ c 1i ρ i1 cos i1 θ sin θ The Number of Limit Cyces for Discontinuous Systems b 1i ρ i cos i θ sin θ 2π m d 1i ρ i1 cos i θ sin 2 θ dθ a 2i ρ i cos i1 θ π b 2i ρ i cos i θ sin θ d 2i ρ i1 cos i θ sin 2 θ dθ m 1 n 1 2i1 π2i 1!! a 12i1 a 22i1 ρ 2i 2!! c 2i ρ i1 cos i1 θ sin θ c 12i1 c 22i1 ρ 2i2 2 2i 3 Int. J. Bifurcation Chaos Downoaded from by CITY UNIVERSITY OF HONG KONG on 09/29/15. For persona use ony. m b 12i b 22i ρ 2i 2 2i 1 2i1 π2i 1!! d 12i d 22i ρ 2i 2!!. n Therefore we obtain that for ε > 0 sufficienty sma the discontinuous system 12 has at most λ = maxm n 1} that bifurcate medium ampitude imit cyces and then the concusion hods. When there are different vector fieds in the different sectors respectivey we consider the foowing discontinuous generaized Liénard poynomia differentia system of the form x = y εg i1 xf i1 xy y = x εg i2 xf i2 xy x y S i 13 where S i denotes the ith sector separated by /2 straight ines passing through the origin is a positive even integer and the poynomias g i x and f i x fori = =1 2are g i1 x = a i x g i2 x = b i x f i1 x = c i x f i2 x = d i x 1 satisfying m 1n 1. Then we have the foowing resut. Theorem 3.5. Consider the system 13 with 0 < ε sufficienty sma then the maximum number of medium ampitude imit cyces bifurcating from the periodic orbits of 6 is λ =maxm n 1}. Proof. Let x = ρ cos θ y = ρ sin θ with ρ>0 then the system 13 is transformed into dρ dθ = ε a i ρ cos 1 θ c i ρ 1 cos 1 θ sin θ b i ρ cos θ sin θ d i ρ 1 cos θ sin 2 θ Oε 2 for θ α 2i 2π α 2iπ i =1 2...With the same proof as Theorem 3.2 it foows from Lemma 3.1 that the concusion hods. Remark 3.6. For system 13 with 1 if S i i = is the ith sector separated by rays starting from the origin then the pane R 2 is divided by the discontinuous rays into sectors. In this case we have the same concusion as Theorem 3.5. Remark 3.7. Note that the maximum number of medium ampitude imit cyces is independent of the number of sectors but it depends on the degree of the poynomias g i x andf i x.. Exampes Exampe.1. Consider the foowing discontinuous generaized Liénard poynomia differentia system with two zones separated by the straight ine

8 F. Jiang et a. y =0}: 2 x = y ε π x 3 2 xy y = x εx 1 y > 0 2 x = y ε π x 1y y = x ε x 2π y xy y < y From the system 15 we observe that g 11 x = 2 π x g 12x =x x Int. J. Bifurcation Chaos Downoaded from by CITY UNIVERSITY OF HONG KONG on 09/29/15. For persona use ony. g 21 x = 2 π x 1 g 22x =x f 11 x = 3 2 x f 12x =0 f 21 x =1 f 22 x = 2 π x so the functions g i x f i x fori =1 2have degrees m =1andn = 1 respectivey. By Lemma 3.1 the averaging function dρ of 15 is as foows π 2 dρ = cos θ 0 π ρ cos θ 3 2 ρ2 sin θ cos θ ] sinθρcos θ 1 dθ 2π π 2 cos θ ρ cos θ 1ρsin θ π sinθ ρ cos θ 2 π ρ sin θ ] ρ 2 sin θ cos θ dθ = ρ 2 3ρ 2= ρ 2ρ 1. It is easy to see that the zeros of dρ areρ 1 =2 and ρ 2 = 1 satisfying d 1 > 0andd 2 < 0. Therefore by Theorem 3.2 we obtain that for ε > 0 sufficienty sma the discontinuous system 15 has two medium ampitude imit cyces. One is a stabe imit cyce bifurcating from the periodic orbit of radius ρ 1 = 2 of the inear center x = y y = x and the other is an unstabe imit cyce bifurcating from the periodic orbit of radius ρ 2 =1.ByMatab we see that the discontinuous Fig. 1. Two medium ampitude imit cyces of 15 with ε =0.01. system 15 has indeed two medium ampitude imit cyces see Fig. 1 the bigger one is a stabe imit cyce and the smaer one is an unstabe imit cyce. Exampe.2. Let = and α = π and we consider the foowing discontinuous generaized Liénard poynomia differentia system x = y ε2x 2 xy y = x ε 1x x 2 8 π x2 y hx y > 0 x = y ε 13x 2 x π 2 y x2 y y = x εx 2 xy hx y < It foows that the scaar function h : R 2 R is defined as hx y =y xy x andthediscontinuity set is of the form h 1 0 = x y :y = x} x y :y = x}. Obviousy g 11 x =2x 2 g 12 x =1x x 2 g 21 x =13x 2 f 11 x =x f 12 x = 8 π x2 f 21 x =x 2 1 x π 2 g 22x =x 2 f 22 x =x then the poynomias g i x f i x fori =1 2 have degrees m =2andn = 2 respectivey

9 By Lemma 3.1 the averaging function dρ of 16 is as foows dρ = 3π π 5π 3π 7π 5π 9π 7π The Number of Limit Cyces for Discontinuous Systems cos θg 11 ρ cos θf 11 ρ cos θρ sin θsinθg 12 ρ cos θf 12 ρ cos θρ sin θ]dθ cos θg 21 ρ cos θf 21 ρ cos θρ sin θsinθg 22 ρ cos θf 22 ρ cos θρ sin θ]dθ cos θg 11 ρ cos θf 11 ρ cos θρ sin θsinθg 12 ρ cos θf 12 ρ cos θρ sin θ]dθ cos θg 21 ρ cos θf 21 ρ cos θρ sin θsinθg 22 ρ cos θf 22 ρ cos θρ sin θ]dθ Int. J. Bifurcation Chaos Downoaded from by CITY UNIVERSITY OF HONG KONG on 09/29/15. For persona use ony. = ρ 3 3ρ 2 2ρ = ρρ 2ρ 1. Therefore for ε > 0 sufficienty sma the discontinuous system 16 has two medium ampitude imit cyces. One is bifurcating from the periodic orbit of radius 2 centered at the origin and it is stabe. The other is bifurcating from the periodic orbit of radius 1 and it is unstabe. Remark.3. Note that for smooth generaized Liénard poynomia differentia system x = y εg 1 x f 1 xy y = x εg 2 x f 2 xy where f i x and g i x i = 1 2 are poynomias with degrees n and m respectivey. Then the maximum number of medium ampitude imit cyces by using the averaging theory of first-order is max m 1 n 2 ]}. So roughy speaking Theorem 3.2 in this paper shows that the discontinuous generaized Liénard poynomia differentia systems can have at east one more imit cyce than the corresponding smooth ones see for instance Exampes.1 and Concusion In this paper we have investigated the maximum number of medium ampitude imit cyces for a cass of discontinuous generaized Liénard poynomia differentia systems with sectors separated by /2 straight ines passing through the origin of coordinates. The discontinuous system consists of two different smooth generaized Liénard poynomia differentia systems of the form x = y g 1 x f 1 xy y = x g 2 x f 2 xy ocated aternativey in every two neighboring sectors where the poynomias f i x andg i x fori =1 2 have degrees n and m respectivey. By using the averaging theory of first-order for discontinuous differentia systems we have obtained that for ε > 0 sufficienty sma the discontinuous generaized Liénard poynomia differentia system has at most maxm n 1} medium ampitude imit cyces and the upper bound maxm n 1} of the number of medium ampitude imit cyces can be attained for some vaues m and n. Moreover we have aso obtained the same resut if there are different vector fieds in different sectors respectivey. Acknowedgments The authors are gratefu to the referee for his/her constructive comments. This work was revised when the first author visited Coege of Wiiam and Mary in Spring She expresses her gratitude to Department of Mathematics Coege of Wiiam and Mary for warm hospitaity and is aso very gratefu to Professor Shi for vauabe guidance. This work is supported by the NSF of China under grant References García B. Libre J. & Pérez de Río J. S. 201] Limit cyces of generaized Liénard poynomia differentia systems via averaging theory Chaos Soit. Fract Gradshteyn I. S. & Ryshik I. M. 199] Tabe of Integras Series and Products 5th edition ed. Jeffrey A. Academic Press NY. Justino A. R. & Jorge L. L. 2012] On the maximum number of imit cyces of a cass of generaized Liénard differentia systems Int. J. Bifurcation and Chaos

10 Int. J. Bifurcation Chaos Downoaded from by CITY UNIVERSITY OF HONG KONG on 09/29/15. For persona use ony. F. Jiang et a. Li J. B. 2003] Hibert s 16th probem and bifurcations of panar poynomia vector fieds Int. J. Bifurcation and Chaos Li C. Z. & Libre J. 2012] Uniqueness of imit cyces for Liénard differentia equations of degree four J. Diff. Eqs Libre J. 2010] The number of imit cyces in a Z2- equivariant Liénard system Math.Proc.Camb.Phi. Soc Libre J. & Vas C. 2012] On the number of imit cyces of a cass of poynomia differentia systems Proc. Roy. Soc. A Libre J. & Mereu A. C. 2013] Limit cyces for discontinuous generaized Liénard poynomia differentia equations Eectron.J.Diff.Eqs Libre J. & Vas C. 2013a] On the number of imit cyces for a generaization of Liénard poynomia differentia systems Int. J. Bifurcation and Chaos Libre J. & Vas C. 2013b] Limit cyces for a generaization of poynomia Liénard differentia systems Chaos Soit. Fract Libre J. & Mereu A. C. 201] Limit cyces for discontinuous quadratic differentia systems with two zones J. Math. Ana. App Libre J. Novaes D. D. & Teixeira M. A. 2015] On the birth of imit cyces for non-smooth dynamica systems Bu. Sci. Math Loyd N. D. & Lynch S. 1988] Sma-ampitude imit cyces of certain Liénard systems Proc. Roy. Soc. Lond. A Martins R. M. & Mereu A. C. 201] Limit cyces in discontinuous cassica Liénard equations Nonin. Ana: Rea Word App Shen J. H. & Han M. A. 2013] Bifurcations of canard imit cyces in severa singuary perturbed generaized poynomia Liénard systems Discr. Contin. Dyn. Syst.Ser.A Sun J. T. 1992] Boundedness of soutions of Liénard systems and existence of imit cyces Gaoxiao Yingyong Shuxue Xuebao Xiong Y. Q. & Zhong H. 2013] The number of imit cyces in a Z2-equivariant Liénard system Int. J. Bifurcation and Chaos

Theory of Generalized k-difference Operator and Its Application in Number Theory

Theory of Generalized k-difference Operator and Its Application in Number Theory Internationa Journa of Mathematica Anaysis Vo. 9, 2015, no. 19, 955-964 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ijma.2015.5389 Theory of Generaized -Difference Operator and Its Appication

More information

Week 6 Lectures, Math 6451, Tanveer

Week 6 Lectures, Math 6451, Tanveer Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n

More information

221B Lecture Notes Notes on Spherical Bessel Functions

221B Lecture Notes Notes on Spherical Bessel Functions Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Solution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima

Solution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution

More information

An approximate method for solving the inverse scattering problem with fixed-energy data

An approximate method for solving the inverse scattering problem with fixed-energy data J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999

More information

Physics 505 Fall Homework Assignment #4 Solutions

Physics 505 Fall Homework Assignment #4 Solutions Physics 505 Fa 2005 Homework Assignment #4 Soutions Textbook probems: Ch. 3: 3.4, 3.6, 3.9, 3.0 3.4 The surface of a hoow conducting sphere of inner radius a is divided into an even number of equa segments

More information

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,

More information

LECTURE 10. The world of pendula

LECTURE 10. The world of pendula LECTURE 0 The word of pendua For the next few ectures we are going to ook at the word of the pane penduum (Figure 0.). In a previous probem set we showed that we coud use the Euer- Lagrange method to derive

More information

Discrete Bernoulli s Formula and its Applications Arising from Generalized Difference Operator

Discrete Bernoulli s Formula and its Applications Arising from Generalized Difference Operator Int. Journa of Math. Anaysis, Vo. 7, 2013, no. 5, 229-240 Discrete Bernoui s Formua and its Appications Arising from Generaized Difference Operator G. Britto Antony Xavier 1 Department of Mathematics,

More information

Formulas for Angular-Momentum Barrier Factors Version II

Formulas for Angular-Momentum Barrier Factors Version II BNL PREPRINT BNL-QGS-06-101 brfactor1.tex Formuas for Anguar-Momentum Barrier Factors Version II S. U. Chung Physics Department, Brookhaven Nationa Laboratory, Upton, NY 11973 March 19, 2015 abstract A

More information

Establishment of Weak Conditions for Darboux- Goursat-Beudon Theorem

Establishment of Weak Conditions for Darboux- Goursat-Beudon Theorem Georgia Southern University Digita Commons@Georgia Southern Mathematica Sciences Facuty Pubications Department of Mathematica Sciences 2009 Estabishment of Weak Conditions for Darboux- Goursat-Beudon Theorem

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,

More information

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE 3 th Word Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 38 DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE Bo JIN SUMMARY The dynamic responses

More information

Lecture Notes 4: Fourier Series and PDE s

Lecture Notes 4: Fourier Series and PDE s Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for

More information

A Two-Parameter Trigonometric Series

A Two-Parameter Trigonometric Series A Two-Parameter Trigonometric Series Xiang-Qian Chang Xiang.Chang@bos.mcphs.edu Massachusetts Coege of Pharmacy and Heath Sciences Boston MA 25 Let us consider the foowing dua questions. A. Givenafunction

More information

On Integrals Involving Universal Associated Legendre Polynomials and Powers of the Factor (1 x 2 ) and Their Byproducts

On Integrals Involving Universal Associated Legendre Polynomials and Powers of the Factor (1 x 2 ) and Their Byproducts Commun. Theor. Phys. 66 (216) 369 373 Vo. 66, No. 4, October 1, 216 On Integras Invoving Universa Associated Legendre Poynomias and Powers of the Factor (1 x 2 ) and Their Byproducts Dong-Sheng Sun ( 孙东升

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

More information

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS TONY ALLEN, EMILY GEBHARDT, AND ADAM KLUBALL 3 ADVISOR: DR. TIFFANY KOLBA 4 Abstract. The phenomenon of noise-induced stabiization occurs

More information

Mat 1501 lecture notes, penultimate installment

Mat 1501 lecture notes, penultimate installment Mat 1501 ecture notes, penutimate instament 1. bounded variation: functions of a singe variabe optiona) I beieve that we wi not actuay use the materia in this section the point is mainy to motivate the

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

830. Nonlinear dynamic characteristics of SMA simply supported beam in axial stochastic excitation

830. Nonlinear dynamic characteristics of SMA simply supported beam in axial stochastic excitation 8. Noninear dynamic characteristics of SMA simpy supported beam in axia stochastic excitation Zhi-Wen Zhu 1, Wen-Ya Xie, Jia Xu 1 Schoo of Mechanica Engineering, Tianjin University, 9 Weijin Road, Tianjin

More information

Strauss PDEs 2e: Section Exercise 1 Page 1 of 7

Strauss PDEs 2e: Section Exercise 1 Page 1 of 7 Strauss PDEs 2e: Section 4.3 - Exercise 1 Page 1 of 7 Exercise 1 Find the eigenvaues graphicay for the boundary conditions X(0) = 0, X () + ax() = 0. Assume that a 0. Soution The aim here is to determine

More information

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed

More information

An explicit Jordan Decomposition of Companion matrices

An explicit Jordan Decomposition of Companion matrices An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057

More information

2M2. Fourier Series Prof Bill Lionheart

2M2. Fourier Series Prof Bill Lionheart M. Fourier Series Prof Bi Lionheart 1. The Fourier series of the periodic function f(x) with period has the form f(x) = a 0 + ( a n cos πnx + b n sin πnx ). Here the rea numbers a n, b n are caed the Fourier

More information

Gaussian Curvature in a p-orbital, Hydrogen-like Atoms

Gaussian Curvature in a p-orbital, Hydrogen-like Atoms Advanced Studies in Theoretica Physics Vo. 9, 015, no. 6, 81-85 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/astp.015.5115 Gaussian Curvature in a p-orbita, Hydrogen-ike Atoms Sandro-Jose Berrio-Guzman

More information

Lecture 11. Fourier transform

Lecture 11. Fourier transform Lecture. Fourier transform Definition and main resuts Let f L 2 (R). The Fourier transform of a function f is a function f(α) = f(x)t iαx dx () The normaized Fourier transform of f is a function R ˆf =

More information

Course 2BA1, Section 11: Periodic Functions and Fourier Series

Course 2BA1, Section 11: Periodic Functions and Fourier Series Course BA, 8 9 Section : Periodic Functions and Fourier Series David R. Wikins Copyright c David R. Wikins 9 Contents Periodic Functions and Fourier Series 74. Fourier Series of Even and Odd Functions...........

More information

The EM Algorithm applied to determining new limit points of Mahler measures

The EM Algorithm applied to determining new limit points of Mahler measures Contro and Cybernetics vo. 39 (2010) No. 4 The EM Agorithm appied to determining new imit points of Maher measures by Souad E Otmani, Georges Rhin and Jean-Marc Sac-Épée Université Pau Veraine-Metz, LMAM,

More information

Research Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation

Research Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation Appied Mathematics and Stochastic Anaysis Voume 007, Artice ID 74191, 8 pages doi:10.1155/007/74191 Research Artice On the Lower Bound for the Number of Rea Roots of a Random Agebraic Equation Takashi

More information

On Some Basic Properties of Geometric Real Sequences

On Some Basic Properties of Geometric Real Sequences On Some Basic Properties of eometric Rea Sequences Khirod Boruah Research Schoar, Department of Mathematics, Rajiv andhi University Rono His, Doimukh-791112, Arunacha Pradesh, India Abstract Objective

More information

HYDROGEN ATOM SELECTION RULES TRANSITION RATES

HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS

More information

Assignment 7 Due Tuessday, March 29, 2016

Assignment 7 Due Tuessday, March 29, 2016 Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

Math 124B January 31, 2012

Math 124B January 31, 2012 Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

BASIC NOTIONS AND RESULTS IN TOPOLOGY. 1. Metric spaces. Sets with finite diameter are called bounded sets. For x X and r > 0 the set

BASIC NOTIONS AND RESULTS IN TOPOLOGY. 1. Metric spaces. Sets with finite diameter are called bounded sets. For x X and r > 0 the set BASIC NOTIONS AND RESULTS IN TOPOLOGY 1. Metric spaces A metric on a set X is a map d : X X R + with the properties: d(x, y) 0 and d(x, y) = 0 x = y, d(x, y) = d(y, x), d(x, y) d(x, z) + d(z, y), for a

More information

1D Heat Propagation Problems

1D Heat Propagation Problems Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2

More information

Two Kinds of Parabolic Equation algorithms in the Computational Electromagnetics

Two Kinds of Parabolic Equation algorithms in the Computational Electromagnetics Avaiabe onine at www.sciencedirect.com Procedia Engineering 9 (0) 45 49 0 Internationa Workshop on Information and Eectronics Engineering (IWIEE) Two Kinds of Paraboic Equation agorithms in the Computationa

More information

Approximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation

Approximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische

More information

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest

More information

General Decay of Solutions in a Viscoelastic Equation with Nonlinear Localized Damping

General Decay of Solutions in a Viscoelastic Equation with Nonlinear Localized Damping Journa of Mathematica Research with Appications Jan.,, Vo. 3, No., pp. 53 6 DOI:.377/j.issn:95-65...7 Http://jmre.dut.edu.cn Genera Decay of Soutions in a Viscoeastic Equation with Noninear Locaized Damping

More information

Tracking Control of Multiple Mobile Robots

Tracking Control of Multiple Mobile Robots Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart

More information

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18 Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms

More information

David Eigen. MA112 Final Paper. May 10, 2002

David Eigen. MA112 Final Paper. May 10, 2002 David Eigen MA112 Fina Paper May 1, 22 The Schrodinger equation describes the position of an eectron as a wave. The wave function Ψ(t, x is interpreted as a probabiity density for the position of the eectron.

More information

(This is a sample cover image for this issue. The actual cover is not yet available at this time.)

(This is a sample cover image for this issue. The actual cover is not yet available at this time.) (This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna

More information

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated

More information

arxiv: v1 [math.fa] 23 Aug 2018

arxiv: v1 [math.fa] 23 Aug 2018 An Exact Upper Bound on the L p Lebesgue Constant and The -Rényi Entropy Power Inequaity for Integer Vaued Random Variabes arxiv:808.0773v [math.fa] 3 Aug 08 Peng Xu, Mokshay Madiman, James Mebourne Abstract

More information

Path planning with PH G2 splines in R2

Path planning with PH G2 splines in R2 Path panning with PH G2 spines in R2 Laurent Gajny, Richard Béarée, Eric Nyiri, Oivier Gibaru To cite this version: Laurent Gajny, Richard Béarée, Eric Nyiri, Oivier Gibaru. Path panning with PH G2 spines

More information

The Group Structure on a Smooth Tropical Cubic

The Group Structure on a Smooth Tropical Cubic The Group Structure on a Smooth Tropica Cubic Ethan Lake Apri 20, 2015 Abstract Just as in in cassica agebraic geometry, it is possibe to define a group aw on a smooth tropica cubic curve. In this note,

More information

Module 22: Simple Harmonic Oscillation and Torque

Module 22: Simple Harmonic Oscillation and Torque Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque

More information

Legendre Polynomials - Lecture 8

Legendre Polynomials - Lecture 8 Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha

More information

Integrating Factor Methods as Exponential Integrators

Integrating Factor Methods as Exponential Integrators Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been

More information

ANALOG OF HEAT EQUATION FOR GAUSSIAN MEASURE OF A BALL IN HILBERT SPACE GYULA PAP

ANALOG OF HEAT EQUATION FOR GAUSSIAN MEASURE OF A BALL IN HILBERT SPACE GYULA PAP ANALOG OF HEAT EQUATION FOR GAUSSIAN MEASURE OF A BALL IN HILBERT SPACE GYULA PAP ABSTRACT. If µ is a Gaussian measure on a Hibert space with mean a and covariance operator T, and r is a} fixed positive

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

EXISTENCE OF SOLUTIONS FOR ONE-DIMENSIONAL WAVE EQUATIONS WITH NONLOCAL CONDITIONS

EXISTENCE OF SOLUTIONS FOR ONE-DIMENSIONAL WAVE EQUATIONS WITH NONLOCAL CONDITIONS Eectronic Journa of Differentia Equations, Vo. 21(21), No. 76, pp. 1 8. ISSN: 172-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (ogin: ftp) EXISTENCE OF SOLUTIONS

More information

In Coulomb gauge, the vector potential is then given by

In Coulomb gauge, the vector potential is then given by Physics 505 Fa 007 Homework Assignment #8 Soutions Textbook probems: Ch. 5: 5.13, 5.14, 5.15, 5.16 5.13 A sphere of raius a carries a uniform surface-charge istribution σ. The sphere is rotate about a

More information

EXPONENTIAL DECAY OF SOLUTIONS TO A VISCOELASTIC EQUATION WITH NONLINEAR LOCALIZED DAMPING

EXPONENTIAL DECAY OF SOLUTIONS TO A VISCOELASTIC EQUATION WITH NONLINEAR LOCALIZED DAMPING Eectronic Journa of Differentia Equations, Vo. 24(24), No. 88, pp. 1 1. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (ogin: ftp) EXPONENTIAL DECAY

More information

Math-Net.Ru All Russian mathematical portal

Math-Net.Ru All Russian mathematical portal Math-Net.Ru A Russian mathematica porta D. Zaora, On properties of root eements in the probem on sma motions of viscous reaxing fuid, Zh. Mat. Fiz. Ana. Geom., 217, Voume 13, Number 4, 42 413 DOI: https://doi.org/1.1547/mag13.4.42

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

ON CERTAIN SUMS INVOLVING THE LEGENDRE SYMBOL. Borislav Karaivanov Sigma Space Inc., Lanham, Maryland

ON CERTAIN SUMS INVOLVING THE LEGENDRE SYMBOL. Borislav Karaivanov Sigma Space Inc., Lanham, Maryland #A14 INTEGERS 16 (2016) ON CERTAIN SUMS INVOLVING THE LEGENDRE SYMBOL Borisav Karaivanov Sigma Sace Inc., Lanham, Maryand borisav.karaivanov@sigmasace.com Tzvetain S. Vassiev Deartment of Comuter Science

More information

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c)

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c) A Simpe Efficient Agorithm of 3-D Singe-Source Locaization with Uniform Cross Array Bing Xue a * Guangyou Fang b Yicai Ji c Key Laboratory of Eectromagnetic Radiation Sensing Technoogy, Institute of Eectronics,

More information

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers Supporting Information for Suppressing Kein tunneing in graphene using a one-dimensiona array of ocaized scatterers Jamie D Was, and Danie Hadad Department of Chemistry, University of Miami, Cora Gabes,

More information

Generalized Bell polynomials and the combinatorics of Poisson central moments

Generalized Bell polynomials and the combinatorics of Poisson central moments Generaized Be poynomias and the combinatorics of Poisson centra moments Nicoas Privaut Division of Mathematica Sciences Schoo of Physica and Mathematica Sciences Nanyang Technoogica University SPMS-MAS-05-43,

More information

Product Cosines of Angles between Subspaces

Product Cosines of Angles between Subspaces Product Cosines of Anges between Subspaces Jianming Miao and Adi Ben-Israe Juy, 993 Dedicated to Professor C.R. Rao on his 75th birthday Let Abstract cos{l,m} : i cos θ i, denote the product of the cosines

More information

Laplace - Fibonacci transform by the solution of second order generalized difference equation

Laplace - Fibonacci transform by the solution of second order generalized difference equation Nonauton. Dyn. Syst. 017; 4: 30 Research Artice Open Access Sandra Pineas*, G.B.A Xavier, S.U. Vasantha Kumar, and M. Meganathan Lapace - Fibonacci transform by the soution of second order generaized difference

More information

6 Wave Equation on an Interval: Separation of Variables

6 Wave Equation on an Interval: Separation of Variables 6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.

More information

Published in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics

Published in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of

More information

A proposed nonparametric mixture density estimation using B-spline functions

A proposed nonparametric mixture density estimation using B-spline functions A proposed nonparametric mixture density estimation using B-spine functions Atizez Hadrich a,b, Mourad Zribi a, Afif Masmoudi b a Laboratoire d Informatique Signa et Image de a Côte d Opae (LISIC-EA 4491),

More information

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006 expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

arxiv: v1 [hep-th] 10 Dec 2018

arxiv: v1 [hep-th] 10 Dec 2018 Casimir energy of an open string with ange-dependent boundary condition A. Jahan 1 and I. Brevik 2 1 Research Institute for Astronomy and Astrophysics of Maragha (RIAAM, Maragha, Iran 2 Department of Energy

More information

AXISYMMETRIC SOLUTIONS OF A TWO-DIMENSIONAL NONLINEAR WAVE SYSTEM WITH A TWO-CONSTANT EQUATION OF STATE

AXISYMMETRIC SOLUTIONS OF A TWO-DIMENSIONAL NONLINEAR WAVE SYSTEM WITH A TWO-CONSTANT EQUATION OF STATE Eectronic Journa of Differentia Equations, Vo. 07 (07), No. 56, pp. 8. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu AXISYMMETRIC SOLUTIONS OF A TWO-DIMENSIONAL NONLINEAR

More information

King Fahd University of Petroleum & Minerals

King Fahd University of Petroleum & Minerals King Fahd University of Petroeum & Mineras DEPARTMENT OF MATHEMATICAL SCIENCES Technica Report Series TR 369 December 6 Genera decay of soutions of a viscoeastic equation Saim A. Messaoudi DHAHRAN 3161

More information

Global Optimality Principles for Polynomial Optimization Problems over Box or Bivalent Constraints by Separable Polynomial Approximations

Global Optimality Principles for Polynomial Optimization Problems over Box or Bivalent Constraints by Separable Polynomial Approximations Goba Optimaity Principes for Poynomia Optimization Probems over Box or Bivaent Constraints by Separabe Poynomia Approximations V. Jeyakumar, G. Li and S. Srisatkunarajah Revised Version II: December 23,

More information

Research Article Numerical Range of Two Operators in Semi-Inner Product Spaces

Research Article Numerical Range of Two Operators in Semi-Inner Product Spaces Abstract and Appied Anaysis Voume 01, Artice ID 846396, 13 pages doi:10.1155/01/846396 Research Artice Numerica Range of Two Operators in Semi-Inner Product Spaces N. K. Sahu, 1 C. Nahak, 1 and S. Nanda

More information

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,

More information

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary

More information

HILBERT? What is HILBERT? Matlab Implementation of Adaptive 2D BEM. Dirk Praetorius. Features of HILBERT

HILBERT? What is HILBERT? Matlab Implementation of Adaptive 2D BEM. Dirk Praetorius. Features of HILBERT Söerhaus-Workshop 2009 October 16, 2009 What is HILBERT? HILBERT Matab Impementation of Adaptive 2D BEM joint work with M. Aurada, M. Ebner, S. Ferraz-Leite, P. Godenits, M. Karkuik, M. Mayr Hibert Is

More information

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation

More information

Small generators of function fields

Small generators of function fields Journa de Théorie des Nombres de Bordeaux 00 (XXXX), 000 000 Sma generators of function fieds par Martin Widmer Résumé. Soit K/k une extension finie d un corps goba, donc K contient un éément primitif

More information

Summation of p-adic Functional Series in Integer Points

Summation of p-adic Functional Series in Integer Points Fiomat 31:5 (2017), 1339 1347 DOI 10.2298/FIL1705339D Pubished by Facuty of Sciences and Mathematics, University of Niš, Serbia Avaiabe at: http://www.pmf.ni.ac.rs/fiomat Summation of p-adic Functiona

More information

Dual Integral Equations and Singular Integral. Equations for Helmholtz Equation

Dual Integral Equations and Singular Integral. Equations for Helmholtz Equation Int.. Contemp. Math. Sciences, Vo. 4, 9, no. 34, 1695-1699 Dua Integra Equations and Singuar Integra Equations for Hemhotz Equation Naser A. Hoshan Department of Mathematics TafiaTechnica University P.O.

More information

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix VOL. NO. DO SCHOOLS MATTER FOR HIGH MATH ACHIEVEMENT? 43 Do Schoos Matter for High Math Achievement? Evidence from the American Mathematics Competitions Genn Eison and Ashey Swanson Onine Appendix Appendix

More information

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization

More information

An H 2 type Riemannian metric on the space of planar curves

An H 2 type Riemannian metric on the space of planar curves An H 2 type Riemannian metric on the space of panar curves Jayant hah Mathematics Department, Northeastern University, Boston MA emai: shah@neu.edu Abstract An H 2 type metric on the space of panar curves

More information

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n. Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

Some Applications on Generalized Hypergeometric and Confluent Hypergeometric Functions

Some Applications on Generalized Hypergeometric and Confluent Hypergeometric Functions Internationa Journa of Mathematica Anaysis and Appications 0; 5(): 4-34 http://www.aascit.org/journa/ijmaa ISSN: 375-397 Some Appications on Generaized Hypergeometric and Confuent Hypergeometric Functions

More information

Binomial Transform and Dold Sequences

Binomial Transform and Dold Sequences 1 2 3 47 6 23 11 Journa of Integer Sequences, Vo. 18 (2015), Artice 15.1.1 Binomia Transform and Dod Sequences Kaudiusz Wójcik Department of Mathematics and Computer Science Jagieonian University Lojasiewicza

More information

Some identities of Laguerre polynomials arising from differential equations

Some identities of Laguerre polynomials arising from differential equations Kim et a. Advances in Difference Equations 2016) 2016:159 DOI 10.1186/s13662-016-0896-1 R E S E A R C H Open Access Some identities of Laguerre poynomias arising from differentia equations Taekyun Kim

More information

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing

More information

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27 Physics 55 Fa 7 Homework Assignment #5 Soutions Textook proems: Ch. 3: 3.3, 3.7, 3.6, 3.7 3.3 Sove for the potentia in Proem 3., using the appropriate Green function otained in the text, and verify that

More information

Introduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

Introduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization

More information

Higher dimensional PDEs and multidimensional eigenvalue problems

Higher dimensional PDEs and multidimensional eigenvalue problems Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u

More information

Input-to-state stability for a class of Lurie systems

Input-to-state stability for a class of Lurie systems Automatica 38 (2002) 945 949 www.esevier.com/ocate/automatica Brief Paper Input-to-state stabiity for a cass of Lurie systems Murat Arcak a;, Andrew Tee b a Department of Eectrica, Computer and Systems

More information

Another Class of Admissible Perturbations of Special Expressions

Another Class of Admissible Perturbations of Special Expressions Int. Journa of Math. Anaysis, Vo. 8, 014, no. 1, 1-8 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ijma.014.31187 Another Cass of Admissibe Perturbations of Specia Expressions Jerico B. Bacani

More information

Research Article Building Infinitely Many Solutions for Some Model of Sublinear Multipoint Boundary Value Problems

Research Article Building Infinitely Many Solutions for Some Model of Sublinear Multipoint Boundary Value Problems Abstract and Appied Anaysis Voume 2015, Artice ID 732761, 4 pages http://dx.doi.org/10.1155/2015/732761 Research Artice Buiding Infinitey Many Soutions for Some Mode of Subinear Mutipoint Boundary Vaue

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information