Lipschitz stability for a coefficient inverse problem for the non-stationary transport equation via Carleman estimate

Size: px
Start display at page:

Download "Lipschitz stability for a coefficient inverse problem for the non-stationary transport equation via Carleman estimate"

Transcription

1 Lipschiz sabiliy for a coefficien inverse problem for he non-saionary ranspor equaion via Carleman esimae Michael V. Klibanov Sergey E. Pamyanykh Deparmen of Mahemaics Saisics, The Universiy of Norh Carolina a Charloe, Charloe, NC 83, U.S.A. s: mklibanv@ .uncc.edu, spamyan@uncc.edu Absrac The Lipschiz sabiliy esimae for a coefficien inverse problem for he non-saionary single-speed ranspor equaion wih he laeral boundary daa is obained. The mehod of Carleman esimaes is used. Uniqueness of he soluion follows.. Inroducion The ranspor equaion is used o model a variey of diffusion processes, such as diffusion of neurons in medium, scaering of ligh in he urbulen amosphere, propagaion of rays in a scaering medium, ec. (see, e.g., he book of Case Zweifel [6]). Coefficien inverse problems (CIPs) for he ranspor equaion are he problems of deermining of he absorpion coefficien, angular densiy of sources or scaering indicarix from an exra boundary daa. They find a variey of applicaions in opical omography, heory of nuclear reacors, ec. (see, e.g., he book of Anikonov, Kovanyuk Prokhorov [], [6]). This paper addresses he quesion of he Lipschiz sabiliy for a CIP for he non-saionary single-speed ranspor equaion wih he laeral boundary daa. In general, sabiliy esimaes for CIPs provide guidelines for he sabiliy of corresponding numerical mehods. Sabiliy, uniqueness exisence resuls references o such resuls for CIPs for he saionary ranspor equaion can be found, e. g., in [] in he book of Romanov [3]. Uniqueness exisence resuls for CIPs for he non-saionary ranspor equaion were obained in he works of Prilepko Ivankov [], [] []. Uniqueness exisence resuls in [] [] were obained for special forms of he unknown coefficien using he overdeerminaion a a poin. Also, uniqueness exisence resuls were obained for an inverse problem wih he final overdeerminaion, i.e. where complee laeral boundary daa is no presen bu boh iniial end condiions (a T) are given; see []. For some recen publicaions on inverse problems for he ranspor equaion see Tamasan [5] Sefanov [4]. A derivaion of he ranspor equaion for he non-saionary case can be found, for example, in [6]. The proof of he main resul of his paper is based on a Carleman esimae, obained by Klibanov Pamyanykh [6]. Tradiionally, Carleman esimaes have been used for proofs of sabiliy uniqueness resuls for non-sard Cauchy problems for PDEs. They were firs inroduced by Carleman in 939 [5], also see, e.g., books of Hörmer [7], Klibanov Timonov [7] Lavrenev, Romanov Shishaskii [9]. Bukhgeim Klibanov [4] have inroduced he ool of Carleman esimaes in he field of CIPs for proofs of global uniqueness sabiliy resuls; also, see Klibanov [], [3] [4], Klibanov Timonov [7], [8]. This mehod works for

2 non-overdeermined CIPs, as long as he iniial condiion is no vanishing he Carleman esimae holds for he corresponding differenial operaor (see Chaper in [7] for he definiion of non-overdeermined CIPs). Recenly, Klibanov Timonov have exended he original idea of [4] [] - [4] for he consrucion of numerical mehods for CIPs, including he case when he iniial condiion is he -funcion; see [7] for deails more references. Klibanov Malinsky [5] Kazemi Klibanov [] have proposed o use he Carleman esimaes for proofs of he Lipschiz sabiliy esimaes for hyperbolic equaions wih he laeral Cauchy daa; also see [7]. The mehod of [4], []-[4] [7] has generaed many publicaions, see, for example, Bellassoued [], [3], Imanuvilov Yamamoo [8], [9] [] he references cied herein. The Lipschiz sabiliy of he soluion of he non-saionary ranspor equaion wih he laeral daa was proved in [6]. In his paper he ideas of [] [5] are combined wih he ideas of [8], [9], [6]-[8]. In Secion he saemens of he resuls are given; in Secion 3, 4 5 he proofs of hese resuls are provided.. Saemens of resuls.. Saemens of resuls Le T R be posiive numbers. Denoe x R n : x R, R n :, H T,T, T,T, Z. Also, denoe C k H s C k H : D x, u x,, v C H, k The ranspor equaion in H has he form [6] u, u a x, u g x,,, u x,, d F x,,, (.) where is he uni vecor of paricle velociy, u x,,v C 3 H is he densiy of paricle flow, a x,v is he absorpion coefficien, F x,, v is he angular densiy of sources, g x,, v, is he scaering indicarix, u denoes he scalar produc of wo vecors u. Consider he following boundary condiion u p x,,, where! x,, v T, T" n, v#. (.) Here n,v is he scalar produc of he ouer uni normal vecor n o he surface he direcion v of he velociy. So, only incoming radiaion is given a he boundary in his case. Equaion (.) wih he boundary condiion (.) he iniial condiion a u x,, v f x,, $ x,v Z, form he classical forward problem for he ranspor equaion in any direcion of (posiive or negaive). Uniqueness, exisence sabiliy resuls for his problem are well known, see, e. g., Prilepko Ivankov []. Suppose now ha he absorpion coefficien a x,v is unknown, bu he following addiional (.3)

3 boundary condiion is given: u q x,,, where x,, v! T, T" n,v. The funcion q x,,v describes he ougoing radiaion a he boundary. Inroduce he funcion x,, x,, p x,,, q x,,, if if n, v#, n, v. (.4) Hence u x,,, $! x,, v T, T". (.5) Thus, we obain he following coefficien inverse problem for he non-saionary ranspor equaion: Inverse Problem: Given he iniial condiion (.3) he laeral daa (.5), deermine he coefficien a x,v of he equaion (.). For a posiive consan M, denoe D M s x C Z : s C Z M. Theorem. [Lipschiz sabiliy uniqueness] Le T R. Suppose ha derivaives k g exis in H k g C H r for k,,, where r is a posiive consan. Le f x,v r f x,v C Z r 3, where r 3 r. Suppose ha he coefficiens a,a D M correspond o he boundary daa x,, v x,,v, respecively, funcions k i,. Then he following Lipschiz sabiliy esimae holds a a L Z K L i L L for k,,, L ", where K K,T,r,r,r 3,M is he posiive consan depending on, T, r, r, r 3, M independen on he funcions a, a,,. In paricular, when, hen a x,v a x,v which implies ha he Inverse Problem has a mos one soluion. (.6) Below K K,T,r, r,r 3,M denoes differen posiive consans, depending on, T, r, r, r 3, M independen on funcions a, a,,, condiions of Theorem are assumed o be saisfied. The proof of Theorem is based on he Carleman esimae formulaed in Lemma. Le L u u, u u where u i u/ x i. Le x R n. Inroduce he funcion n iu i, i 3

4 x, x x, cons,. Le c cons. Denoe G c x x, : x x c x R. (.7) Obviously, G c G c if c c. (.7_) Inroduce he Carleman Weigh Funcion (CWF) as C x, exp x, ". Lemma. Choose he number such ha, T R/. Also, choose he consan c such ha G c x T,T. Then here exis posiive consans G c x B B G c x, depending only on he domain G c x, such ha he following poinwise Carleman esimae holds in G c x for all funcions u x,, C G c x C for all G c x : where he vecor funcion U,V saisfies he esimae L u C u C U V, (.8) U,V B u C. (.9) The proof of his lemma can be found in [6]. Also, we will use he following Lipschiz sabiliy resul, proved in [6] Theorem. Suppose ha he funcion u C! T, T" C saisfies he condiions (.) (.4). Le funcions a x,, g x,,, be bounded, i.e. a x,, r 5 $ x,,v H g x,,, r 6 $ x,, v, H, where r 5 r 6 are posiive consans. Le funcions x,,v L, F x,,v L H le T R. Then he following Lipschiz sabiliy esimae holds: u L H K L F L H ", where K K,T,r 5,r 6 is he posiive consan independen on funcions u, F... Preliminaries Before proceeding wih he proof of he Theorem, we inroduce some new funcions formulae necessary lemmaa. Le funcions u u be soluions of equaion (.) wih he iniial condiion (.3) he laeral daa (.5) for a x,v a x,v, x,,v x,, v a x,v a x,v, x,,v x,, v, respecively. Denoe u u u, 4

5 a a a, (.). From relaions (.), (.3), (.5) (.), noicing ha a u a u a u u, u a x, u g x,,, u x,, d au, au, we obain (.) u x,,v, $ x,v Z, (.) u x,,, $ x,, v! T, T". (.3) Applying he Theorem o he equaion (.) wih laeral daa (.3), we obain he following esimae for he funcion u Denoe u L H K L a L Z. u. Differeniaing (.) (.3) wih respec o, we obain, a x, g u g d au (.4) (.5) x,,, $ x,,v! T, T" Seing in (.), we obain x,,v au x,,v a x, v f x,v, where x, v Z. Differeniaing (.5) (.6) wih respec o denoing w, we obain w, w a x, w. g u g gw d au, (.6) (.7) (.8) w x,,, $ x,,v! T,T" We will need he following lemma Lemma. Le funcions a x,, a x, D M. The following Lipschiz sabiliy esimaes. (.9) hold: L H K a L Z L L ", (.) w L H K a L Z L L L ". (.) These esimaes are similar o he Lipschiz sabiliy esimae ha was obained in [6], bu do no follow direcly from he resul of [6] due o he presence of he funcion u in (.5) (.8). 5

6 The following lemma provides an esimae from he above for an inegral conaining he CWF. Lemma 3. For all funcions s C G c x for all, he following esimae holds G c x s x, d C x, dxd G c x s C x, dxd. See Secion 3. in [7] for he proof. Lemma 4. Le T R. Then for any c, R here exiss a G c T,T for all R,T,c,. Proof. By he definiion of he domain G c G c T,T max x, T c., R,T,c, such ha i.e. when which leads o he following inequaliy R T c, Since c,r R T hen 3. Proof of Lemma R c T., we can choose. Denoe G c G c R, T, such ha Choose so ha Choose a small number for arbirary c cons. Since T R, we can choose a small number T R 3 x 3 R,T, / ( Lemma 4 ) le, for he sake of definieness, G / 3 G /, /, such ha R,T, / T,T. T, T ". x 3 Consider he domains G / 3 G / G / G /. (See (.7_) Fig. for a schemaic represenaion in he - D case),.. (3.) (3.) (3.3) 6

7 Fig.. Ses G / 3 G / G / G /. Also, consider he cu-off funcion x, C T, T, such ha x, beween in in G /, T,T \G /, in G / \G /. The equaions (.5) (.8) imply ha, K u d d a, (3.4) w, w K w u d d w d a (3.5) Le x,, x,, x,. Then 7

8 n n i i i i i i. i i i Derivaives, i, i,..., n equal o zero in G / in T,T \G / are bounded in G / \G /. So, using he inequaliy (3.4), we obain n i i i n K u d d a. (3.6) Similarly, for w x,, w x,, x,, we obain from (3.5) w n iw i i K w u d d w d a w (3.7) Denoe u u x,, x,. Then (3.6) (3.7) become n i i i K u d d a (3.8) w n iw i i K w u d d w d a w (3.9) Muliplying (3.8) (3.9) by he CWF squaring boh sides, we obain 8

9 n i i C K u d d a C i K " C, w n iw i i C K w u d d w d a w C. The Carleman esimae (.8) leads o C U V (3.) K u d d w C U a V C (3.) K w u d d w d a w C where x,, H /, H / G / funcions U, V U, V are he funcions U, V from he Carleman esimae (.8)-(.9) for he case, when he funcion u is replaced by he funcions w, respecively. Inegraing over H / applying he Gauss formula, we obain C dh K u d d a C dh H / H / K C dh U,V ds (3.) H / M / Similarly, we obain for w w C dh (3.3) H / 9

10 K w u d d w d a C dh H / w C dh U,V ds. H / where dh dxd vd, M / G / Noicing ha for any funcion s x,,v M / G / denoes he boundary of he domain G /. C H s d C dh A s C dh, H / H / where A is he area of he uni sphere, we remove he inner inegrals over in (3.) (3.3). So, (3.) (3.3) become C dh K u a C dh H / H / K C dh U,V ds H / M / w C dh H / K w u a C dh H / w C dh U,V ds. H / Choose such ha K/ /. Then for all M / C dh we have H / K a C dh u C dh C dh U,V ds H / H / H / M /

11 w C dh H / K a C dh u C dh C dh w C dh H / H / H / H / U,V ds. M / Using (.9), we obain C dh (3.4) H / K a C dh u C dh C ds H / H / C dh K M H / / w C dh (3.5) H / K a C dh u C dh C dh w C dh H / H / H / H / K M w C ds. The boundary M / of he domain G / consiss of wo pars M / M / M /, where M / / x,,v : x R G / M / x,,v : x / G /. Since x,,v x,,v w x,,v x,,v, for x,,v M /,

12 x,,v w x,, v, for x,, v M /, hen C ds C ds w C ds C ds. M / M / M / M / Esimae boh sides of he inequaliy (3.4). Noe ha since in H / H / 3 H /, hen C dh C dh e / 3 dh. (3.6) Also, since Hence, H / H / 3 $ x, in G /, hen C e /, H / 3 x, H /. C dh e / H / dh. H / Therefore (3.4) (3.6) lead o e / 3 dh (3.7) H / 3 K a C dh u C dh e / dh C ds. H / H / H / M / Similarly, from (3.5) we obain e / 3 w dh (3.8) H / 3 K a C dh u C dh C dh e / w dh C ds. H / H / H / H / M / Le m sup x. Then (3.7) (3.8) yield G / e / 3 L H / 3 (3.9)

13 K e / L H / e m a L H / u L H / L M / e / 3 w L H / 3 (3.) K e / w L H / e m a L H / u L H / L H / /. L M Since u x,,v u x,,v x,,v x,, v $ x,, v H, hen (3.9) (3.) become e / 3 L H / 3 K e / L H / e m a L H / u L H / L M / e / 3 w L H / 3 K e / w L H / e m a L H / u L H / L H / /. L M Dividing hese inequaliies by exp / 3 ", we obain L H / 3 (3.) K e 5 L H / e m a L H / u L H / /, L M w L H / 3 K e 5 w L H / (3.) m K e a L H / u L H / L H / /. L M An inconvenience of he domain H / 3 for our goal is ha alhough he domain H / 3, bu H / 3. Thus, we now shif his domain. Choose an x such ha x 3 / consider he domain G / x, which is obained by a shif of he domain G /. Clearly one can choose R,T, / so small ha in addiion o (3.)-(3.3) Then G / x T, T G / 3 x T, T ". 3

14 G / 3 G / 3 x x x x 3 3 Consider now he ball B, / 3 : x : x / 3. By (3.) B, / 3, /. We prove now ha B G / 3 x. Le x of he ball B. Then Since Hence, by (3.4) B x x x x 3 x 3, /, hen 3 / 3. Hence, G / 3 Hence, here exiss a number E x x (3.3) (3.4), since B be an arbirary poin x. Therefore, using (3.3) (3.4), we obain ha G / 3 G / 3,T such ha he layer x, : x, The schemaic represenaion of he domains G /, G /. x. G / G / x. (3.5) x E in -D case is provided on Fig. 4

15 Fig.. G / Solid line, G / x Dashed line, E Shaded area. Since he Carleman esimae (.8)-(.9) is valid for he domain G / x, we can obain esimaes similar o (3.) (3.) L H / 3 x (3.6) K e 5 L H / x e m a L H / x u L H / x L M x / w L H / 3 x K e 5 w L H / x (3.7) m K e a L H / x u L H / x L H / x / x. L M where H / x G / x M / x G / x x, : x R. Consider now he layer E defined by (3.5) (see Fig.). Esimaes (3.), (3.6) (3.), (3.7) lead o he following esimaes in E : L E K e 5 L H e m w L E a L Z u L H L (3.8) (3.9) K e 5 w L H e m a L Z u L H L H L. Since for any funcion s x,,v C H here exiss, such ha hen (3.8) (3.9) lead o s x,,v dxd v s L E, x,, v dxd v N, (3.3) w x,,v dxd v N, 5

16 where N K e 5 L H e m a L Z u L H L (3.3) N K e 5 w L H e m a L Z u L H L H L. Le S, T, H, T, S T,, H T,. Denoe Y x,, n i i, (3.3) i x,, x,, S x,,. Esimae he L H norm of he funcion. Muliplying (3.3) by Z,, where, T, we obain inegraing over Consider he vecor funcion B so (3.33) becomes dxd vd n i v i idxd vd,,..., n. Then n i i i B, Ydxd vd. (3.33) x,,v dxd v x,,v dxd v B,n dsd vd K dxd vd Y dxd vd. Here B,n denoes he scalar produc of vecors B n, where n is he ouward normal vecor on. 6

17 Noicing ha B, where using he Cauchy-Schwarz inequaliy, we obain x,,v dxd v x,, v dxd v dsd vd (3.34) K dxd vd Y dxd vd, Esimae Y using (3.4) (3.3) Y K u d d a. (3.35) Esimaes (3.34) (3.35) lead o x,,v dxd v x,,v dxd v dsd vd K dxd vd u dxd vd a dxd vd. Using he Gronwall s inequaliy, we obain x,, v dxd v (3.36) K x,,v dxd v dsd vd u dxd vd a dxd vd. Subsiuing (3.3) (3.3) in he righ-h side of (3.36), we ge x,,v dxd v K N dsd vd u dxd vd a dxd vd K e 5 L H e m a L Z u L H L 7

18 K dsd vd u dxd vd a dxd vd K e 5 L H e m a L Z u L H L. Thus, L H K e 5 L H e m a L Z One can obain similar esimae for L H. Summing up ha esimae wih (3.37), we obain L H K e 5 e L H m a L Z u L H u L H L. L. (3.37) To remove he erm wih u from he laer formula we apply he esimae (.4). Hence L H K e 5 L H Consider, such ha e m a L Z Ke 5. L L. Then Choosing max,, we obain m L H K e a L Z 5 ln K. e m L L, (3.38) which implies he desired esimae (.). Applying he procedure, similar o (3.3)-(3.38), o he equaions depending on w, using he esimae (3.38), one can similarly obain he esimae (.). 4. Proof of he Theorem. This secion consiss of hree subsecions. In he subsecion 4. geomery is defined he proof of he Theorem is sared. In he subsecion 4. he supplemenary fac is proved. In he subsecion 4.3 he proof of he Theorem is finished. 4.. Beginning of he Proof of Theorem. The proof of he heorem is based on he Carleman esimae (.8)-(.9). The values of he 8

19 parameers, ha are used in he proof of his heorem are independen on he values of hese parameers used in he proof of he Lemma. Consider he problem (.)-(.3) in H. Also, consider he relaions, (.5)-(.7) (.8)-(.9). A equaion (.) becomes u x,,v au x,, v, (4.) Since u x,, v f x,v hen (4.) leads o Since f x,v r, a x,v K u x,, v. (4.) u x,,v u x,, v u x,, v d, we have u x,,v u x,,v u x,, v d. (4.3) $ $ Choose a poin x R n, R x R. Choose he number, such ha T R/ Denoe he domains P c G c x G c x, c, where he domains G c x are defined by (.7). Choose he consan c such ha x x c G c T, x R n : x R. Hence, x T. Define he domain b,b choose consans b such ha b P c 3 P c. (See fig. 3 for a schemaic represenaion in he - D case). 9

20 Fig. 3. The shaded area schemaically represens he domain P c. Consider he domains P c 3 P c P c P c. Also, consider he funcion x, C T,T, such ha x, beween in in P c, T, T \P c, in P c \P c, le x, be a non-increasing funcion of in he domain P c \P c non-decreasing funcion of in he domain P c \P c holds for any funcion s x,,v C H any x,,v H, a, so ha he following inequaliy x, s x,, v d x, s x,,v d. An example of such funcion is consruced in Appendix A. Denoe x,, x,, x, w x,, w x,, x,. Following he proof of Lemma from (3.4) o (3.5), we obain he analogs o esimaes (3.4) (3.5) for he domains C dh (4.4)

21 K a C dh u C dh C dh K Bc C ds, w C dh (4.5) K a C dh u C dh C dh w C dh K Bc w C ds, where B c is he boundary of he domain. Represen he inegrals u C dh C dh as a sums of inegrals u C dh u C dh u C dh \ C dh C dh C dh, consider he inegrals over he domain firs. Since u x,, v u x,, v \ u x,,v d x,,v x,,v x,,v d, we obain, using (.) (.7), u x,,v u x,, v u x,, v d u x,, v d (4.6)

22 x,, v x,,v x,, v d a f x,,v d, (4.7) Since u x,,v u x,, v, x,, v x,, v, $ x,,v, hen, applying (4.6) (4.7) o he inegrals u C dh C dh, we obain u C dh (4.8) K u x,,v d C dh K x,,v d C dh C dh K a C dh x,,v d C dh K a C dh w x,, v d C dh. (4.9) Applying Lemma 3 o (4.8) (4.9), we obain u C dh K Qc C dh (4.)

23 C dh K a C dh w C dh. (4.) Also, applying he esimae (4.7) o he righ-h side of (4.) using Lemma 3, we obain K u C dh Qc a C dh w C dh. (4.) Applying he esimaes (4.), (4.) (4.) o (4.4) (4.5), choosing sufficienly large, we obain C dh o be (4.3) K a C dh u C dh C dh K Bc C ds, \ w C dh (4.4) K a C dh u C dh C dh w C dh \ \ K Bc w C ds, Since u x,,v u x,,v x,,v x,, v $ x,, v H, (4.3) (4.4) become C dh (4.5) K a C dh u C dh C dh K Bc C ds \ \ 3

24 w C dh (4.6) K a C dh u C dh C dh w C dh \ \ K Bc w C ds. 4.. Proof of he Inegral Inequaliy. Here we esimae he inegral from he above hrough he inegral Consider he funcion c x Qc a C dh b a C dh. x x c. Then for any funcion s x,,v Hence, s x,,v dh C, which is even wih respec o he variable, we have Z c x s x,,v dd vdx c x a C dh a C dh Z c x s x,,v dd vdx. a C dh (4.7) (4.8) \ a C dh b a C dd vdx c x a C dd vdx. \ Z Z b Noe ha, since a x,v is independen of, we have b a C dd vdx c x a C d d vdx Z Z b 4

25 $ a b C x, dd vdx a c x C x, d d vdx. (4.9) Z Since he funcion is decreasing when, we have c x C x, d e x x b c x e b Z b e d c x b e x x e b c x b e x x b b e b d. Since hen b e b d b e d, c x b e x x b b e b d c So, by (4.8) (4.9), we obain x b e x x b a C dh b e a C dh d K b b K C x, d. a C dh. (4.) \ Noe ha C x, e c x, P c \P c. (4.) From (4.7), (4.9) (4.), we obain a C dh e c a dh e c a dh \ \ 5

26 e c a d vdx Z c x c x d Thus, we have Ke c Since b P c 3 P c, hen So, (4.) implies ha \ e c a d vdx Z a C dh b d Ke c a C dh Ke c e c 3 C x, $ K e c a dh b b b a dh. x, K b b. a dh. a C dh. (4.) (4.3) \ Finally, by (4.) (4.3), we have a C dh 4.3. The End of he Proof of he Theorem. K b a C dh. (4.4) Consider now he esimaes (4.), (4.3) (4.5). By (4.) we have a x,v K x,, v, (4.5) (4.3) leads o x,,v x,, v w x,, v d. (4.6) Combining (4.5) (4.6), we obain a x,v x,, v w x,,v d. Muliplying he las inequaliy by he C x, inegraing over 3, we obain a x, v C dh (4.7) 3 6

27 x,,v C dh w x,,v d C dh. 3 3 Since 3, he esimaes (4.5) (4.6) lead o C dh 3 (4.8) K a C dh u C dh C dh K Bc C ds \ w C dh 3 (4.9) K a C dh u C dh C dh w C dh \ \ K Bc w C ds. Since x,,v x,,v, $ x,, v 3, hen, combining he esimaes (4.7) (4.8), we obain a C dh w x,,v d C dh (4.3) 3 3 K a C dh u C dh C dh K Bc C ds. \ By Lemma 3 w x,, v d C x, dh w x,,v C x, dh. 3 3 Hence, (4.3) leads o 7

28 3 a C dh Qc 3 w C dh (4.3) K a C dh u C dh C dh K Bc C ds. \ Summing up he esimaes (4.3) (4.9), noicing ha w x,,v w x,, v, $ x,,v 3, aking, we obain a C dh w C dh (4.3) 3 3 K a C dh u C dh C dh \ \ K C dh w C dh K Bc C ds K Bc w C ds. The boundary B c consiss of wo pars. Denoe B c x, : x R P c, B c x, : x x c P c. Then B c B c B c. Since x,,v x,,v w x,,v x,, v, if x,, v B c, x,,v w x,, v, if x,,v B c, hen C ds C ds w C ds C ds. B c B c B c B c Thus, (4.3) leads o 8

29 3 a x, v C dh K a C dh u C dh C dh w C dh \ \ \ Noicing ha b K Bc C ds K B c C ds. 3 applying (4.4) o he las inequaliy, we obain a x, v C dh b (4.33) K b a C dh \ u C dh C dh \ w C dh \ K Bc C ds K B c C ds. Taking K in (4.33), we obain a x, v C dh b (4.34) K u C dh C dh w C dh \ \ \ Le m sup x x K Bc. Then, since C ds K B c C ds. max C x, : x, \ e c, inequaliy (4.34) yields b a x, v C dh 9

30 K e c u L H L H w L H Le d inf x x b K e m. Then (4.35) becomes L e d a L Z L. (4.35) K e c u L H L H w L H K e m L L. Using he esimaes for L H w L H, given by Lemma he esimae (.4) for u L H, we obain e d a L Z K e c a L Z L L L Since d K e m L c, hen dividing (4.36) by e d aking we obain he desired esimae (.6). K exp d L. c ", o be so large ha (4.36) Appendix A. Here we consruc supplemenary funcion. Consider consans C i, i,...,6, ha will be chosen laer, denoe he surfaces in R n, corresponding o hese consans, S i x, : x Le C C. Consider he funcion C C i, i,..., 6. C, C C e exp C C C C C C, C C C, C C. This is a non-increasing funcion of he parameer C. Consider he funcion 3

31 x, x, x, R n T,T. Consider any x R n, such ha he line x x in R n T,T crosses boh surfaces S S. Le firs. Choose arbirary,, T",, such ha he poins x, x, are locaed beween he surfaces S S. Clearly, he poins x, x, correspond o differen level surfaces of he funcion x,, S 3 S 4, respecively, ha have corresponding consans C 3 C 4, such ha C C 4 C 3 C (see. Fig.4). Since Fig.4. Schemaic represenaion of level surfaces for -D case. C is a non-increasing funcion, we have x, x,. Thus, he funcion x, is non-increasing wih respec o, when. Le. Choose arbirary 3, 4 T,", 3 4, such ha he poins x, 3 x, 4 are locaed beween he surfaces S S. Clearly, he poins x, 3 x, 4 correspond o differen level surfaces of funcion x,, S 5 S 6, respecively, ha have corresponding consans C 5 C 6, such ha C C 6 C 5 C (see Fig. 4). Since he funcion C is a non-increasing funcion, we have x, 3 x, 4. Thus, he funcion x, is non-decreasing wih respec o, when. So, since he funcion x, is coninuously differeniable in R n T,T, we can ake i as he funcion. Acknowledgmen The work of M.V. Klibanov was suppored by, or in par by, he U.S. Army Research Laboraory U.S. Army Research Office under conrac/gran number W9NF His work was also 3

32 parially suppored by NATO under he gran number PDD(CP)-(PST.NR.CLG 9863). References. Anikonov D.S., Kovanyuk A.E. Prokhorov I.V. Transpor Equaion Tomography. VSP Publ., Urech, The Neherls,.. Bellassoued M. Global logarihmic sabiliy for inverse hyperbolic problem by arbirary boundary observaion. Inverse Problems, v., (4), pp Bellassoued M. Uniqueness sabiliy in deermining he speed of propagaion of second-order hyperbolic equaion wih variable coefficiens. Applicable Analysis, v. 83, (4), pp Bukhgeim A.L. Klibanov M.V. Uniqueness in he large of a class of mulidimensional inverse problems. Sovie Mah. Doklady, v. 7, (98), pp (Also available online a hp:// 5. Carleman T. Sur un problème d unicié pour les sysèmes d équaions aux dérivées parielles à deux variables indépendanes. Ark. Ma. Asr. Fys., 6B, No. 7, (939), pp Case K.M. Zweifel P.F. Linear Transpor Theory. Addison-Wesley, Reading, MA, Hörmer L. Linear Parial Differenial Operaors. Springer-Verlag, Berlin, Imanuvilov O.Y., Yamamoo M. Global Lipschiz sabiliy in an inverse hyperbolic problem by inerior observaions. Inverse Problems, v. 7, (), pp Imanuvilov O.Y., Yamamoo M. Global uniqueness sabiliy in deermining coefficiens of wave equaions. Commun. in Par. Diff. Eqs., v. 6, (), pp Imanuvilov O.Y., Yamamoo M. Deerminaion of a coefficien in an acousic equaion wih a single measuremen. Inverse Problems, v. 9, (3), pp Kazemi M. Klibanov M.V. Sabiliy esimaes for ill-posed Cauchy problem involving hyperbolic equaions inequaliies. Applicable Analysis, v. 5, (993), pp Klibanov M.V. Uniqueness in he large of some mulidimensional inverse problems, in Non-Classical Problems of Mahemaical Physics, Proc. Compuing Cener of he Siberian Branch of he Russian Academy of Science, Novosibrisk, 98, pp. -4. (In Russian). 3. Klibanov M.V. Inverse problems in he large Carleman bounds. Differenial Equaions, v., (984), pp (Also available online a hp:// 4. Klibanov M.V. Inverse problems Carleman esimaes. Inverse Problems, v. 8, (99) pp (Also available online a hp:// 5. Klibanov M.V. Malinsky J. Newon-Kanorovich mehod for 3-dimensional poenial inverse scaering problem sabiliy of he hyperbolic Cauchy problem wih ime dependen daa. Inverse Problems, v. 7, (99), pp Klibanov M.V. Pamyanykh S.E. Lipschiz sabiliy of a non-sard problem for he non-saionary ranspor equaion via Carleman esimae. Hp:// preprin number Klibanov M.V. Timonov A. Carleman Esimaes for Coeffcien Inverse Problems Numerical Applicaions. VSP Publ., Urech, The Neherls, Klibanov M.V. Timonov A. Global uniqueness for a 3d/d inverse conduciviy problem via he modified mehod of Carleman esimaes. J. Inverse Ill-Posed Problems, v. 3, (5), pp Lavrenev M.M., Romanov V.G. Shishaskii S.P. Ill-posed problems of mahemaical physics analysis. AMS, Providence, Rhode Isl,

33 . Prilepko A.I. Ivankov A.L. Inverse problems for non-saionary ranspor equaion. Sovie Mah. Doklady, v. 9, (984), No. 3, pp Prilepko A.I. Ivankov A.L. Inverse problems for he deerminaion of a coefficien he righ side of a non-saionary mulivelociy ranspor equaion wih overdeerminaion a a poin. Differenial equaions, v., (985), No., pp Prilepko A.I. Ivankov A.L. Inverse problems of finding a coefficien, he scaering indicarix, he righ side of a nonsaionary many-velociy ranspor equaion. Differenial equaions, v., (985), No. 5, pp Romanov V.G. Invesigaion Mehods for Inverse Problems. VSP Publ., Urech, The Neherls,. 4. Sefanov P. Inverse problems in ranspor heory. Inside ou: inverse problems applicaions, pp.-3, Mah. Sci. Res. Ins. Publ., 47, Cambridge Univ. Press, Cambridge, Tamasan A. An inverse boundary value problem in wo-dimensional ranspor. Inverse Problems, v. 8, (), No., pp

Lipschitz stability for a coefficient inverse problem for the non-stationary transport equation via Carleman estimate

Lipschitz stability for a coefficient inverse problem for the non-stationary transport equation via Carleman estimate Lipschiz sabiliy for a coefficien inverse problem for he non-saionary ranspor equaion via Carleman esimae Sergey E Pamyanykh Deparmen of Mahemaics Saisics, The Universiy of Norh Carolina a Charloe, Charloe,

More information

Global uniqueness for a coefficient inverse problem for the non-stationary transport equation via Carleman estimate

Global uniqueness for a coefficient inverse problem for the non-stationary transport equation via Carleman estimate Global uniqueness for a coefficien inverse problem for he non-saionary ranspor equaion via Carleman esimae Michael V. Klibanov and Sergey E. Pamyanykh Deparmen of Mahemaics and Saisics, The Universiy of

More information

Estimates of Initial Conditions of Parabolic Equations and Inequalities Via the Lateral Cauchy Data

Estimates of Initial Conditions of Parabolic Equations and Inequalities Via the Lateral Cauchy Data Esimaes of Iniial Condiions of Parabolic Equaions and Inequaliies Via he Laeral Cauchy Daa Michael V. Klibanov Deparmen of Mahemaics and Saisics Universiy of Norh Carolina a Charloe, Charloe, NC 28223,

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

t 2 B F x,t n dsdt t u x,t dxdt

t 2 B F x,t n dsdt t u x,t dxdt Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Convergence of the Neumann series in higher norms

Convergence of the Neumann series in higher norms Convergence of he Neumann series in higher norms Charles L. Epsein Deparmen of Mahemaics, Universiy of Pennsylvania Version 1.0 Augus 1, 003 Absrac Naural condiions on an operaor A are given so ha he Neumann

More information

Undetermined coefficients for local fractional differential equations

Undetermined coefficients for local fractional differential equations Available online a www.isr-publicaions.com/jmcs J. Mah. Compuer Sci. 16 (2016), 140 146 Research Aricle Undeermined coefficiens for local fracional differenial equaions Roshdi Khalil a,, Mohammed Al Horani

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

On Gronwall s Type Integral Inequalities with Singular Kernels

On Gronwall s Type Integral Inequalities with Singular Kernels Filoma 31:4 (217), 141 149 DOI 1.2298/FIL17441A Published by Faculy of Sciences and Mahemaics, Universiy of Niš, Serbia Available a: hp://www.pmf.ni.ac.rs/filoma On Gronwall s Type Inegral Inequaliies

More information

Heat kernel and Harnack inequality on Riemannian manifolds

Heat kernel and Harnack inequality on Riemannian manifolds Hea kernel and Harnack inequaliy on Riemannian manifolds Alexander Grigor yan UHK 11/02/2014 onens 1 Laplace operaor and hea kernel 1 2 Uniform Faber-Krahn inequaliy 3 3 Gaussian upper bounds 4 4 ean-value

More information

STABILITY OF NONLINEAR NEUTRAL DELAY DIFFERENTIAL EQUATIONS WITH VARIABLE DELAYS

STABILITY OF NONLINEAR NEUTRAL DELAY DIFFERENTIAL EQUATIONS WITH VARIABLE DELAYS Elecronic Journal of Differenial Equaions, Vol. 217 217, No. 118, pp. 1 14. ISSN: 172-6691. URL: hp://ejde.mah.xsae.edu or hp://ejde.mah.un.edu STABILITY OF NONLINEAR NEUTRAL DELAY DIFFERENTIAL EQUATIONS

More information

Homogenization of random Hamilton Jacobi Bellman Equations

Homogenization of random Hamilton Jacobi Bellman Equations Probabiliy, Geomery and Inegrable Sysems MSRI Publicaions Volume 55, 28 Homogenizaion of random Hamilon Jacobi Bellman Equaions S. R. SRINIVASA VARADHAN ABSTRACT. We consider nonlinear parabolic equaions

More information

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a

More information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,

More information

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL: Ann. Func. Anal. 2 2011, no. 2, 34 41 A nnals of F uncional A nalysis ISSN: 2008-8752 elecronic URL: www.emis.de/journals/afa/ CLASSIFICAION OF POSIIVE SOLUIONS OF NONLINEAR SYSEMS OF VOLERRA INEGRAL EQUAIONS

More information

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity ANNALES POLONICI MATHEMATICI LIV.2 99) L p -L q -Time decay esimae for soluion of he Cauchy problem for hyperbolic parial differenial equaions of linear hermoelasiciy by Jerzy Gawinecki Warszawa) Absrac.

More information

CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS

CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS SARAJEVO JOURNAL OF MATHEMATICS Vol.10 (22 (2014, 67 76 DOI: 10.5644/SJM.10.1.09 CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS ALMA OMERSPAHIĆ AND VAHIDIN HADŽIABDIĆ Absrac. This paper presens sufficien

More information

MATH 5720: Gradient Methods Hung Phan, UMass Lowell October 4, 2018

MATH 5720: Gradient Methods Hung Phan, UMass Lowell October 4, 2018 MATH 5720: Gradien Mehods Hung Phan, UMass Lowell Ocober 4, 208 Descen Direcion Mehods Consider he problem min { f(x) x R n}. The general descen direcions mehod is x k+ = x k + k d k where x k is he curren

More information

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,

More information

POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION

POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Novi Sad J. Mah. Vol. 32, No. 2, 2002, 95-108 95 POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Hajnalka Péics 1, János Karsai 2 Absrac. We consider he scalar nonauonomous neural delay differenial

More information

CHARACTERIZATION OF REARRANGEMENT INVARIANT SPACES WITH FIXED POINTS FOR THE HARDY LITTLEWOOD MAXIMAL OPERATOR

CHARACTERIZATION OF REARRANGEMENT INVARIANT SPACES WITH FIXED POINTS FOR THE HARDY LITTLEWOOD MAXIMAL OPERATOR Annales Academiæ Scieniarum Fennicæ Mahemaica Volumen 31, 2006, 39 46 CHARACTERIZATION OF REARRANGEMENT INVARIANT SPACES WITH FIXED POINTS FOR THE HARDY LITTLEWOOD MAXIMAL OPERATOR Joaquim Marín and Javier

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Differential Harnack Estimates for Parabolic Equations

Differential Harnack Estimates for Parabolic Equations Differenial Harnack Esimaes for Parabolic Equaions Xiaodong Cao and Zhou Zhang Absrac Le M,g be a soluion o he Ricci flow on a closed Riemannian manifold In his paper, we prove differenial Harnack inequaliies

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Positive continuous solution of a quadratic integral equation of fractional orders

Positive continuous solution of a quadratic integral equation of fractional orders Mah. Sci. Le., No., 9-7 (3) 9 Mahemaical Sciences Leers An Inernaional Journal @ 3 NSP Naural Sciences Publishing Cor. Posiive coninuous soluion of a quadraic inegral equaion of fracional orders A. M.

More information

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

Existence Theory of Second Order Random Differential Equations

Existence Theory of Second Order Random Differential Equations Global Journal of Mahemaical Sciences: Theory and Pracical. ISSN 974-32 Volume 4, Number 3 (22), pp. 33-3 Inernaional Research Publicaion House hp://www.irphouse.com Exisence Theory of Second Order Random

More information

A New Perturbative Approach in Nonlinear Singularity Analysis

A New Perturbative Approach in Nonlinear Singularity Analysis Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

TO our knowledge, most exciting results on the existence

TO our knowledge, most exciting results on the existence IAENG Inernaional Journal of Applied Mahemaics, 42:, IJAM_42 2 Exisence and Uniqueness of a Periodic Soluion for hird-order Delay Differenial Equaion wih wo Deviaing Argumens A. M. A. Abou-El-Ela, A. I.

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

Theory of! Partial Differential Equations-I!

Theory of! Partial Differential Equations-I! hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

An Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation

An Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation Commun Theor Phys Beijing, China 43 2005 pp 591 596 c Inernaional Academic Publishers Vol 43, No 4, April 15, 2005 An Invariance for 2+1-Eension of Burgers Equaion Formulae o Obain Soluions of KP Equaion

More information

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x . 1 Mah 211 Homework #3 February 2, 2001 2.4.3. y + (2/x)y = (cos x)/x 2 Answer: Compare y + (2/x) y = (cos x)/x 2 wih y = a(x)x + f(x)and noe ha a(x) = 2/x. Consequenly, an inegraing facor is found wih

More information

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11. 1 Mah 334 Tes 1 KEY Spring 21 Secion: 1 Insrucor: Sco Glasgow Daes: Ma 1 and 11. Do NOT wrie on his problem saemen bookle, excep for our indicaion of following he honor code jus below. No credi will be

More information

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,

More information

On Oscillation of a Generalized Logistic Equation with Several Delays

On Oscillation of a Generalized Logistic Equation with Several Delays Journal of Mahemaical Analysis and Applicaions 253, 389 45 (21) doi:1.16/jmaa.2.714, available online a hp://www.idealibrary.com on On Oscillaion of a Generalized Logisic Equaion wih Several Delays Leonid

More information

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations! hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

On a Fractional Stochastic Landau-Ginzburg Equation

On a Fractional Stochastic Landau-Ginzburg Equation Applied Mahemaical Sciences, Vol. 4, 1, no. 7, 317-35 On a Fracional Sochasic Landau-Ginzburg Equaion Nguyen Tien Dung Deparmen of Mahemaics, FPT Universiy 15B Pham Hung Sree, Hanoi, Vienam dungn@fp.edu.vn

More information

Lecture 10: The Poincaré Inequality in Euclidean space

Lecture 10: The Poincaré Inequality in Euclidean space Deparmens of Mahemaics Monana Sae Universiy Fall 215 Prof. Kevin Wildrick n inroducion o non-smooh analysis and geomery Lecure 1: The Poincaré Inequaliy in Euclidean space 1. Wha is he Poincaré inequaliy?

More information

Lecture 10: Wave equation, solution by spherical means

Lecture 10: Wave equation, solution by spherical means Lecure : Wave equaion, soluion by spherical means Physical modeling eample: Elasodynamics u (; ) displacemen vecor in elasic body occupying a domain U R n, U, The posiion of he maerial poin siing a U in

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

Boundedness and Exponential Asymptotic Stability in Dynamical Systems with Applications to Nonlinear Differential Equations with Unbounded Terms

Boundedness and Exponential Asymptotic Stability in Dynamical Systems with Applications to Nonlinear Differential Equations with Unbounded Terms Advances in Dynamical Sysems and Applicaions. ISSN 0973-531 Volume Number 1 007, pp. 107 11 Research India Publicaions hp://www.ripublicaion.com/adsa.hm Boundedness and Exponenial Asympoic Sabiliy in Dynamical

More information

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx. . Use Simpson s rule wih n 4 o esimae an () +. Soluion: Since we are using 4 seps, 4 Thus we have [ ( ) f() + 4f + f() + 4f 3 [ + 4 4 6 5 + + 4 4 3 + ] 5 [ + 6 6 5 + + 6 3 + ]. 5. Our funcion is f() +.

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3. Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies

More information

Y 0.4Y 0.45Y Y to a proper ARMA specification.

Y 0.4Y 0.45Y Y to a proper ARMA specification. HG Jan 04 ECON 50 Exercises II - 0 Feb 04 (wih answers Exercise. Read secion 8 in lecure noes 3 (LN3 on he common facor problem in ARMA-processes. Consider he following process Y 0.4Y 0.45Y 0.5 ( where

More information

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary

More information

Fractional Method of Characteristics for Fractional Partial Differential Equations

Fractional Method of Characteristics for Fractional Partial Differential Equations Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics

More information

arxiv:math/ v1 [math.nt] 3 Nov 2005

arxiv:math/ v1 [math.nt] 3 Nov 2005 arxiv:mah/0511092v1 [mah.nt] 3 Nov 2005 A NOTE ON S AND THE ZEROS OF THE RIEMANN ZETA-FUNCTION D. A. GOLDSTON AND S. M. GONEK Absrac. Le πs denoe he argumen of he Riemann zea-funcion a he poin 1 + i. Assuming

More information

Symmetry and Numerical Solutions for Systems of Non-linear Reaction Diffusion Equations

Symmetry and Numerical Solutions for Systems of Non-linear Reaction Diffusion Equations Symmery and Numerical Soluions for Sysems of Non-linear Reacion Diffusion Equaions Sanjeev Kumar* and Ravendra Singh Deparmen of Mahemaics, (Dr. B. R. Ambedkar niversiy, Agra), I. B. S. Khandari, Agra-8

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

More information

EXISTENCE OF NON-OSCILLATORY SOLUTIONS TO FIRST-ORDER NEUTRAL DIFFERENTIAL EQUATIONS

EXISTENCE OF NON-OSCILLATORY SOLUTIONS TO FIRST-ORDER NEUTRAL DIFFERENTIAL EQUATIONS Elecronic Journal of Differenial Equaions, Vol. 206 (206, No. 39, pp.. ISSN: 072-669. URL: hp://ejde.mah.xsae.edu or hp://ejde.mah.un.edu fp ejde.mah.xsae.edu EXISTENCE OF NON-OSCILLATORY SOLUTIONS TO

More information

Lab #2: Kinematics in 1-Dimension

Lab #2: Kinematics in 1-Dimension Reading Assignmen: Chaper 2, Secions 2-1 hrough 2-8 Lab #2: Kinemaics in 1-Dimension Inroducion: The sudy of moion is broken ino wo main areas of sudy kinemaics and dynamics. Kinemaics is he descripion

More information

THE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant).

THE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant). THE WAVE EQUATION 43. (S) Le u(x, ) be a soluion of he wave equaion u u xx = 0. Show ha Q43(a) (c) is a. Any ranslaion v(x, ) = u(x + x 0, + 0 ) of u(x, ) is also a soluion (where x 0, 0 are consans).

More information

Solution of Integro-Differential Equations by Using ELzaki Transform

Solution of Integro-Differential Equations by Using ELzaki Transform Global Journal of Mahemaical Sciences: Theory and Pracical. Volume, Number (), pp. - Inernaional Research Publicaion House hp://www.irphouse.com Soluion of Inegro-Differenial Equaions by Using ELzaki Transform

More information

BEng (Hons) Telecommunications. Examinations for / Semester 2

BEng (Hons) Telecommunications. Examinations for / Semester 2 BEng (Hons) Telecommunicaions Cohor: BTEL/14/FT Examinaions for 2015-2016 / Semeser 2 MODULE: ELECTROMAGNETIC THEORY MODULE CODE: ASE2103 Duraion: 2 ½ Hours Insrucions o Candidaes: 1. Answer ALL 4 (FOUR)

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

More information

Research Article Existence and Uniqueness of Periodic Solution for Nonlinear Second-Order Ordinary Differential Equations

Research Article Existence and Uniqueness of Periodic Solution for Nonlinear Second-Order Ordinary Differential Equations Hindawi Publishing Corporaion Boundary Value Problems Volume 11, Aricle ID 19156, 11 pages doi:1.1155/11/19156 Research Aricle Exisence and Uniqueness of Periodic Soluion for Nonlinear Second-Order Ordinary

More information

Some New Uniqueness Results of Solutions to Nonlinear Fractional Integro-Differential Equations

Some New Uniqueness Results of Solutions to Nonlinear Fractional Integro-Differential Equations Annals of Pure and Applied Mahemaics Vol. 6, No. 2, 28, 345-352 ISSN: 2279-87X (P), 2279-888(online) Published on 22 February 28 www.researchmahsci.org DOI: hp://dx.doi.org/.22457/apam.v6n2a Annals of

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Hamilton Jacobi equations

Hamilton Jacobi equations Hamilon Jacobi equaions Inoducion o PDE The rigorous suff from Evans, mosly. We discuss firs u + H( u = 0, (1 where H(p is convex, and superlinear a infiniy, H(p lim p p = + This by comes by inegraion

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS Xinping Guan ;1 Fenglei Li Cailian Chen Insiue of Elecrical Engineering, Yanshan Universiy, Qinhuangdao, 066004, China. Deparmen

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

POSITIVE AND MONOTONE SYSTEMS IN A PARTIALLY ORDERED SPACE

POSITIVE AND MONOTONE SYSTEMS IN A PARTIALLY ORDERED SPACE Urainian Mahemaical Journal, Vol. 55, No. 2, 2003 POSITIVE AND MONOTONE SYSTEMS IN A PARTIALLY ORDERED SPACE A. G. Mazo UDC 517.983.27 We invesigae properies of posiive and monoone differenial sysems wih

More information

Math 333 Problem Set #2 Solution 14 February 2003

Math 333 Problem Set #2 Solution 14 February 2003 Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial

More information

Endpoint Strichartz estimates

Endpoint Strichartz estimates Endpoin Sricharz esimaes Markus Keel and Terence Tao (Amer. J. Mah. 10 (1998) 955 980) Presener : Nobu Kishimoo (Kyoo Universiy) 013 Paricipaing School in Analysis of PDE 013/8/6 30, Jeju 1 Absrac of he

More information

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,

More information

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4. PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence

More information

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

Stochastic Model for Cancer Cell Growth through Single Forward Mutation Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

After the completion of this section the student. Theory of Linear Systems of ODEs. Autonomous Systems. Review Questions and Exercises

After the completion of this section the student. Theory of Linear Systems of ODEs. Autonomous Systems. Review Questions and Exercises Chaper V ODE V.5 Sysems of Ordinary Differenial Equaions 45 V.5 SYSTEMS OF FIRST ORDER LINEAR ODEs Objecives: Afer he compleion of his secion he suden - should recall he definiion of a sysem of linear

More information

On the probabilistic stability of the monomial functional equation

On the probabilistic stability of the monomial functional equation Available online a www.jnsa.com J. Nonlinear Sci. Appl. 6 (013), 51 59 Research Aricle On he probabilisic sabiliy of he monomial funcional equaion Claudia Zaharia Wes Universiy of Timişoara, Deparmen of

More information

arxiv: v1 [math.fa] 9 Dec 2018

arxiv: v1 [math.fa] 9 Dec 2018 AN INVERSE FUNCTION THEOREM CONVERSE arxiv:1812.03561v1 [mah.fa] 9 Dec 2018 JIMMIE LAWSON Absrac. We esablish he following converse of he well-known inverse funcion heorem. Le g : U V and f : V U be inverse

More information

A Necessary and Sufficient Condition for the Solutions of a Functional Differential Equation to Be Oscillatory or Tend to Zero

A Necessary and Sufficient Condition for the Solutions of a Functional Differential Equation to Be Oscillatory or Tend to Zero JOURNAL OF MAEMAICAL ANALYSIS AND APPLICAIONS 24, 7887 1997 ARICLE NO. AY965143 A Necessary and Sufficien Condiion for he Soluions of a Funcional Differenial Equaion o Be Oscillaory or end o Zero Piambar

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

CONTRIBUTION TO IMPULSIVE EQUATIONS

CONTRIBUTION TO IMPULSIVE EQUATIONS European Scienific Journal Sepember 214 /SPECIAL/ ediion Vol.3 ISSN: 1857 7881 (Prin) e - ISSN 1857-7431 CONTRIBUTION TO IMPULSIVE EQUATIONS Berrabah Faima Zohra, MA Universiy of sidi bel abbes/ Algeria

More information

THE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2.

THE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2. THE BERNOULLI NUMBERS The Bernoulli numbers are defined here by he exponenial generaing funcion ( e The firs one is easy o compue: (2 and (3 B 0 lim 0 e lim, 0 e ( d B lim 0 d e +e e lim 0 (e 2 lim 0 2(e

More information

LINEAR INVARIANCE AND INTEGRAL OPERATORS OF UNIVALENT FUNCTIONS

LINEAR INVARIANCE AND INTEGRAL OPERATORS OF UNIVALENT FUNCTIONS LINEAR INVARIANCE AND INTEGRAL OPERATORS OF UNIVALENT FUNCTIONS MICHAEL DORFF AND J. SZYNAL Absrac. Differen mehods have been used in sudying he univalence of he inegral ) α ) f) ) J α, f)z) = f ) d, α,

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

arxiv: v1 [math.ca] 15 Nov 2016

arxiv: v1 [math.ca] 15 Nov 2016 arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University THE MYSTERY OF STOCHASTIC MECHANICS Edward Nelson Deparmen of Mahemaics Princeon Universiy 1 Classical Hamilon-Jacobi heory N paricles of various masses on a Euclidean space. Incorporae he masses in he

More information

Testing for a Single Factor Model in the Multivariate State Space Framework

Testing for a Single Factor Model in the Multivariate State Space Framework esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics

More information