arxiv: v1 [gr-qc] 16 Jun 2009

Size: px
Start display at page:

Download "arxiv: v1 [gr-qc] 16 Jun 2009"

Transcription

1 Lae-ime ails of wave maps copled o graviy Pior Bizoń, Tadesz Chmaj,,3 ndrzej Rosworowski, and Sanis law Zaj ac M. Smolchowski Insie of Physics, Jagiellonian Universiy, Kraków, Poland H. Niewodniczanski Insie of Nclear Physics, Polish cademy of Sciences, Kraków, Poland 3 Cracow Universiy of Technology, Kraków, Poland Daed: November, 8) We consider he lae-ime asympoic behavior for solions of Einsein s eqaions wih he wave map maer. Solions saring from small compacly sppored l-eqivarian iniial daa wih l are shown o decay as l+) a fre imelike infiniy and as l+) a fre nll infiniy. I. INTRODUCTION arxiv:96.99v gr-qc] 6 Jn 9 In his paper we conine or invesigaions, iniiaed in ], of he precise qaniaive descripion of he laeime asympoic behavior of self-graviaing massless fields. In ] we considered he simples case of a spherically symmeric massless scalar field. Using nonlinear perrbaion analysis we showed ha solions saring from small iniial daa decay as 3 a imelike infiniy and as a nll infiniy. We also derived a simple analyic formla for he amplide of he lae-ime ail in erms of iniial daa. Here we sdy he analogos problem for wave maps which are a naral geomeric generalizaion of he wave eqaion for he massless scalar field. This generalizaion seems ineresing becase in he so called eqivarian case he homoopy index l of he map plays he role similar o he mlipole index for spherical harmonics. However, in conras o he decomposiion of a scalar field ino spherical harmonics which makes sense only a he linearized level, i is consisen o sdy nonlinear evolion for he wave map wihin a fixed eqivariance class. In his sense l-eqivarian self-graviaing wave maps can serve as a poor man s oy-model of non-spherical collapse. The l = case redces o he spherically symmeric massless scalar field analyzed in ] so hereafer we assme ha l. We noe aside ha he l = case has been exensively sdied in he pas focsing on he criical behavior a he hreshold of black hole formaion, 3, 4, 5, 6], however, o or knowledge, he lae-ime behavior of wave maps copled o graviy has no been analyzed before. Using he same hird-order perrbaion mehod as in ] we show here ha for small compacly sppored iniial daa he lae-ime ail of he self-graviaing l-eqivarian wave map decays as l+) a fre imelike infiniy and as l+) a fre nll infiniy. We also compe he amplide of he ail in erms of iniial daa. These analyic resls are verified by he nmerical inegraion of he Einsein-wave map eqaions. II. SETUP Le U : M N be a map from a spaceime M,g ab ) ino a Riemannian manifold N,G B ). pair U,g ab ) is said o be a wave map copled o graviy if i is a criical poin of he acion fncional R S = 6πG λ ) gab a U b U B G B dv, ) M Rishescalarcrvareofhemericg ab, GisNewon sconsan,λishewavemapcoplingconsan,anddv is hevolmeelemenonm,g ab ). Thefieldeqaionsderivedfrom)areheEinseineqaionsR ab g abr = 8πGT ab wih he sress-energy ensor T ab = λ a U b U B ) g abg cd c U d U B ) G B, ) and he wave map eqaion g U +Γ BC U) au B b U C g ab =, 3) Γ BC are he Chrisoffel symbols of he arge meric G B and g is he wave operaor associaed wih he meric g ab. s a arge manifold we ake he hree-sphere wih he rond meric in polar coordinaes U = F,Ω) G B du du B = df +sin F dω. 4) For he for dimensional spaceime M we assme spherical symmery and se he following ansaz for he meric g ab dx a dx b = e α,r) e β,r) d +dr ) +r dω. 5)

2 In addiion we assme ha he map U is spherically l-eqivarian, ha is F = F,r), Ω = χ l ω), 6) χ l ω) is a homogeneos harmonic polynomial of degree l. For his ansaz, he energy-momenm ensor ) does no depend on angles and hs can be consisenly copled o he spherically symmeric Einsein eqaions. We noe in passing ha a very similar idea of inrodcing he anglar momenm ino spherical collapse was p forward by Olabarriea e al. 7] in he conex of criical phenomena. In erms of he mass fncion m,r) = r e α ) he Einsein eqaions ake he following form hereafer primes and dos denoe parial derivaives wih respec o r and, respecively) m = r e α F +e β F ) + ll+) sin F, 7) ṁ = r e α F F, 8) β = m r eα ll+)e αsin F, 9) r = 4πGλ is a dimensionless parameer. The wave map eqaion 3) akes he form e β ) F r r e β F ) +e β+α ll+) sinf r =. ) For l = he above eqaions redce o he Einsein-massless scalar field eqaions analyzed by s in ]. For = no graviy) eqaions 7)-9) are rivially solved by m = and β =, while eqaion ) redces o he fla space wave map eqaion. III. ITERTIVE SCHEME We assme ha iniial daa are small, smooh, and compacly sppored he las assmpion can be replaced by a siable fall-off condiion) F,r) = εgr), F,r) = εhr). ) We make he following perrbaive expansion m,r) = m,r)+εm,r)+ε m,r)+..., ) β,r) = β,r)+εβ,r)+ε β,r)+..., 3) F,r) = F,r)+εF,r)+ε F,r)+ε 3 F 3,r)+... 4) Sbsiing his expansion ino he field eqaions and groping erms wih he same power of ε we ge he ieraive scheme which can be solved recrsively. We consider perrbaions abo Minkowski spaceime, so m = β = F =. he firs order he meric fncions m = β = his follows from reglariy a r = ), while F saisfies he fla space radial wave eqaion for he l-h spherical harmonic l) F =, l) = r r r + ll+) r, 5) wih iniial daa F,r) = gr), F,r) = hr). The general every reglar solion of eqaion 5) is given by a sperposiion of ogoing and ingoing waves F,r) = F re,r)+f adv,r), 6) F re,r) = r l l k)! k!l k)! k= a k) ), Fadv v ) l k,r) = r l k+ l k)! a k) v) ), 7) k!l k)! v ) l k k= and =, v = +r are he rearded and advanced imes, respecively he sperscrip in rond brackes denoes he k-h derivaive). Noe ha for compacly sppored iniial daa he generaing fncion ax) can be chosen o have compac sppor as well his condiion deermines ax) niqely).

3 he second order l) F =, hence F = becase i has zero iniial daa), while he meric fncions saisfy he following eqaions m = r F +F + ll+) ) r F, 8) ṁ = r F F, 9) β = m r ll+) F r ) We emporarily pospone he discssion of his sysem and proceed now o he hird order, we have l) F 3 = β F + β F +β F ll+)m F r 3 + ll+)f3 3r. ) To solve his eqaion we se he Dhamel formla for he solion of he inhomogeneos wave eqaion l) F = N,r) wih zero iniial daa F,r) = r dτ +r τ τ ρp l µ)nτ,ρ)dρ, ) P l µ) are Legendre polynomials of degree l and µ = r +ρ τ) )/rρ noe ha µ wihin he inegraion range). pplying his formla o eqaion ), sing nll coordinaes η = τ ρ and ξ = τ +ρ, and he abbreviaion Km,β,F) = β F + β F +β F ll+)/r )mf/r F 3 /3), we obain F 3,r) = 8r +r dξ ξ η)p l µ)km ξ,η),β ξ,η),f ξ,η)), 3) ξ now µ = r +ξ ) η))/rξ η). If he iniial daa ) vanish oside a ball of radis R, hen for > r+r we may drop he advanced par of F,r) and inerchange he order of inegraion in 3) o ge F 3,r) = 8r +r ξ η)p l µ)km ξ,η),β ξ,η),f re ξ,η))dξ. 4) In order o deermine he lae-ime behavior of F 3,r) we need o know he behavior of he sorce erm K along he ligh cone for large vales of r he inersecion of he inegraion range in 4) wih he sppor of F re,r)). Having ha, we shall expand he fncion K in 4) in he inverse powers of ρ = ξ η)/ and calclae he inegrals sing he following ideniy see he appendix in 8] for he derivaion) +r dξ P l µ) ξ η) n = )ln )l l+)!! r l+ η) n l η) r ] n F l+ n, l+3 n l+3/ 3 ) ) r, 5) η n ) l = n )n 3) n l ). Now, we rern o he analysis of he second-order eqaions 8)-). Sbsiing he ogoing solion 7) ino 8) and inegraing, we ge r m,r) >R = ) a l+) ρ) k l+ n k ρ a l++n k) ρ)a l n) ρ) ρ k l+n) n l+k n) k n) l +l+k n)n+) ) k k k n)!n! ] dρ, 6) we sed ha m,r = ) =, which follows from reglariy of iniial daa a he origin and 9). Here and in he following we se repeaedly he fac ha ax) = for x > R, R being he radis of a ball on which he iniial daa ) are sppored. To describe he behavior of m,r) along he lighcone i is convenien o se he nll coordinae = insead of, and rewrie 6) as m,r) r+>r = ) a l+) l +l+ s) ds r a l) ) ) ll+)l +l+) 4 a l ) )a l) ) ) r +O r 3. 7)

4 Nex, sing he gage freedom o se β,r = ) = and inegraing eqaion ), we ge r β,r) >R = ρ r ρ l +l+) a l+) s)) dsdρ a l) ρ) ) The firs inegral can be inegraed by pars giving β,r) r+>r = ) a a l+) l+) s) ) s) ds+ r r s ) ds l +l+) ρ 3 + ll+)3l +3l+) a l ) ρ)a l) ) ] ρ) ρ 4 +O ρ 5 dρ.8) a l) s) ) r s )) 3 + ll+)3l +3l+) 4 ] a l ) s)a l) ) s) r s )) 4 ds+o r 5. 9) IV. TILS Now, we shall apply he mehod described above o compe he lae-ime asympoics of solions in he hird-order approximaion. Hereafer, i is convenien o define he following inegrals for non-negaive inegers m, n) In m ) = s ) m a s)) n) ds. 3) From 9) we have. l = β,r) r+>r = r I)+ )] r I ) 5I))+O r, 3) β,r) r+>r = I r )+ )] r 4I ) 5a )) )+O r, 3) β,r) r+>r = r I) 5 )] r a )) +O r. 33) Sbsiing 7) and 3-33) ino 4) we obain F 3,r) = 4 r ξ η +r dξ P µ) ξ η) I η)a η) d U η) d I η)a η) ) ) +O ξ η) )], 34) U η) = 4I η)a η)+i η) 5I η))a η). 35) Performing he inner inegral over ξ in 34) wih he help of he ideniy 5) we ge he asympoic behavior which is valid for large rearded imes )] r F 3,r) = r ) +O, 36) = 8 3 a s)) as)ds. 37) From 36) we obain he lae-ime ails in boh asympoic regimes: F 3,r) r 4 a fre imelike infiniy r = cons, ) and rf 3 )v =,) ) a fre nll infiniy v =, ).

5 5 B. l We give he deailed calclaion only for l =. In his case we have from 9) r I 3 )+ r I 3 ) 3I ))+ r I3 3 ) 39I )+3a )) )+O )] r 3 β,r) r+>r =, 38) β,r) r+>r = I r 3)+ r 4I 3) 3a )) )+ )] r 6I 3) 39I) 6a )a ))+O r 3,39) β r+>r,r) = r I3 3 ) r a )) 6 )] r a )a )+O r 3. 4) Sbsiing 7) and 38-4) ino 4) we obain and F 3,r) = 4 r + 4 ξ η) +r dξ P l µ) ξ η) d a η)) 3 + ) I3 η)a3) η) + 5I3 ξ η η)a η)+ d )) ) aη)+ a η)) 3 3 a 3) η) 8 d V η) ) U η) +O ξ η) 3 )], 4) U η) = 8I 3η)a η)+i 3η) I η))a 3) η), 4) V η) = 4I 3η)aη)+36I 3η)a η)+ 7I η)+8i 3η) ) a η)+ 78I η)+4i 3 3η)+6a η)) ) a 3) η). 43) Performing he inner inegral over ξ in 4) wih he help of he ideniy 5) we ge he asympoic behavior for large rearded imes = 8 5 F 3,r) = r r ) 3 +B +O ) ] as) a s)) 3 3 a 3) s) ds and B = 64 5 )], 44) a s)) 3 ds. 45) For he general l i is easy o see ha he firs nonzero conribion o he ail comes from he erm wih n = l+ in he ideniy 5) which gives he following asympoics r l )] F 3,r) = r ) l+ l +B l +O. 46) The formla 46) gives he firs erm in he asympoic series approximaion of he solion for lae rearded imes, ha is for small ε we have r ) l+ r l F,r) ε 3 F 3,r) = Oε 5 ). 47) We havenoaemped oderiveageneralformlaforhecoefficiens l andb l hecompaionofhese coefficiens for each given l is sraighforward b as l increases he algebra becomes edios since i involves high-order expansions of he meric fncions along he ligh cone. nyway, i follows from 46) ha he ail behaves as F 3,r) r l l+) a fre imelike infiniy and as rf 3 )v =,) l+) a fre nll infiniy. Remark. For l he ail 46) has wo pars qanified by he coefficiens l and B l, respecively. The l -par comes from he graviaional self-ineracion of he wave map and vanishes for =. The B l -par comes from he cbic nonlineariy of he wave map eqaion and is presen wiho graviy as well. The case l = is special in he sense ha he B l -par is absen in 36) since i is sbdominan decaying as 5 ) wih respec o he leading order erm. Remark. I is insrcive o compare he ail 46) wih he ail for a es linear massless field propagaing on a fixed saionary asympoically fla backgrond. ccording o he Price law 9,, ] he l-h mlipole of his linear ail φ l,r) r l+ / r ) l+ for r. This decay is by one power faser han ha in 46). Of corse, his difference is no very srprising as he ail sdied here and Price s ail correspond o differen physical siaions, however we poin i o as anoher example of he inapplicabiliy of linearized heory in he sdy of radiaive relaxaion processes see, 3] for oher examples). We shall discss his isse in more deail else 4].

6 6 V. NUMERICS In his secion we compare he above analyic predicions wih he resls of nmerical solions of Einsein-wave map eqaions 7-) for varios iniial daa. The deails of he nmerical mehod were given in ] for he case l =. The only difference for higher l is he bondary condiion F,r) r l for small r which garanees reglariy a he origin. The iniial daa were generaed by he gassian for differen vales of ε. For hese iniial daa he formla 37) gives for l = and he formla 45) gives for l = = 4 8 εax) = εexp x ) 48) = 64 3π 7.769, 49) 7 3π 388., B = 48 3π ) 45 In order o exrac he parameers of he ails a imelike infiniy we fi or nmerical daa wih he formla F,r) = γ exp B/+C/ ). 5) imelike infiniy nll infiniy. e-7 F,r) e-6 e-9 l= F,r)/r l e- e-3 e-6 l= l= e- l= e-9 e-5 l=3 e- l= FIG. : Lef panel: The log-log plo of F,r) for fixed r = 5. Fiing 5) we ge power-law exponens γ = 4.96 l = ), 6.9 l = ), 8.49 l = 3), in agreemen wih he analyic predicion 46). Righ panel: The log-log plo of F,r)/r l for fixed large advanced ime v = +r = as he fncion of rearded ime =. The analogos fi o 5) yields he exponens.36 l = ), 3.4 l = ), 4.95 l = 3), in accordance wih 46). In boh panels =. and ε =. l = ), ε =.7 l = ), ε =.3 l = 3). l = ) l = ) ε heory nmerics ε heory nmerics e-5 9.5e e e TBLE I: The comparison of analyic and nmerical amplides of he ails a imelike infiniy. Here =. and r = 5. The hird-order approximaion reads = ε 3 r for l = and = ε 3 r +B ) for l =. The resls and heir confronaion wih analyic predicions are smmarized in Table and Figres,, and 3. From his comparison we conclde ha he hird-order approximaion is excellen for sfficienly small iniial daa. For large daa approaching he black-hole hreshold he hird-order approximaion breaks down his is seen in Fig. as he deviaion from he scaling ε 3 and in Fig. 3 as he deviaion from he linear dependence of on.

7 7 l= l= nmerics 3-order approximaion nmerics 3-order approximaion ε ε FIG. : The log-log plo of he amplide of he ail a imelike infiniy as a fncion of he amplide of iniial daa black dos) for fixed =. and r = 5. The hird-order approximaion dashed line) is excellen for small daa, b i breaks down for large daa lying near he hreshold of black hole formaion. l= l=. nmerics 3-order approximaion.5 nmerics 3-order approximaion FIG. 3: The plo of he amplide of he ail a imelike infiniy as a fncion of he copling consan black dos) for fixed ε =.5 and r = 5. s increases we leave he small-daa regime and conseqenly he hird-order approximaion dashed line) deerioraes. I shold be emphasized ha we ge he same decay raes l+) a imelike infiniy) and l+) a nll infiniy) for all sbcriical evolions, regardless of wheher or hird-order formla reprodces accraely he amplide of he ail for small daa) or fails for large daa). cknowledgmens: We acknowledge sppor by he MNII grans: NN 7935 and 89/6.PRUE/7/7. ] P. Bizoń, T. Chmaj, and. Rosworowski, gr-qc/ ] P. Bizoń and. Wasserman, Phys. Rev. D6, 843 ). 3] S. Hsa e al., Phys. Rev. D6, 47 ). 4] Ch. Lechner e al., Phys. Rev. D65, 85 ). 5] P. Bizoń and. Wasserman, Class. Qan. Grav. 9, 339 ). 6] P. ichelbrg, P. Bizoń, and Z. Tabor, Class. Qan. Grav. 3, S99 6). 7] I. Olabarriea e al., Phys. Rev. D76, 44 7). 8] P. Bizoń, T. Chmaj, and. Rosworowski, Phys. Rev. D78, 444 8). 9] R.H. Price, Phys. Rev. D 5, 49 97). ] C. Gndlach, R. Price, and J. Pllin, Phys. Rev. D 49, ). ] E. Poisson, Phys. Rev. D 66, 448 ). ] P. Bizoń, T. Chmaj, and. Rosworowski, Phys. Rev. D76, 435 7). 3] P. Bizoń, T. Chmaj, and. Rosworowski, Class. Qanm Grav. 4, F55 7). 4] P. Bizoń and. Rosworowski, in preparaion

Scalar Conservation Laws

Scalar Conservation Laws MATH-459 Nmerical Mehods for Conservaion Laws by Prof. Jan S. Heshaven Solion se : Scalar Conservaion Laws Eercise. The inegral form of he scalar conservaion law + f ) = is given in Eq. below. ˆ 2, 2 )

More information

4.2 Continuous-Time Systems and Processes Problem Definition Let the state variable representation of a linear system be

4.2 Continuous-Time Systems and Processes Problem Definition Let the state variable representation of a linear system be 4 COVARIANCE ROAGAION 41 Inrodcion Now ha we have compleed or review of linear sysems and random processes, we wan o eamine he performance of linear sysems ecied by random processes he sandard approach

More information

, u denotes uxt (,) and u. mean first partial derivatives of u with respect to x and t, respectively. Equation (1.1) can be simply written as

, u denotes uxt (,) and u. mean first partial derivatives of u with respect to x and t, respectively. Equation (1.1) can be simply written as Proceedings of he rd IMT-GT Regional Conference on Mahemaics Saisics and Applicaions Universii Sains Malaysia ANALYSIS ON () + () () = G( ( ) ()) Jessada Tanhanch School of Mahemaics Insie of Science Sranaree

More information

Srednicki Chapter 20

Srednicki Chapter 20 Srednicki Chaper QFT Problems & Solions. George Ocober 4, Srednicki.. Verify eqaion.7. Using eqaion.7,., and he fac ha m = in his limi, or ask is o evalae his inegral:! x x x dx dx dx x sx + x + x + x

More information

A Mathematical model to Solve Reaction Diffusion Equation using Differential Transformation Method

A Mathematical model to Solve Reaction Diffusion Equation using Differential Transformation Method Inernaional Jornal of Mahemaics Trends and Technology- Volme Isse- A Mahemaical model o Solve Reacion Diffsion Eqaion sing Differenial Transformaion Mehod Rahl Bhadaria # A.K. Singh * D.P Singh # #Deparmen

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

A Comparison Among Homotopy Perturbation Method And The Decomposition Method With The Variational Iteration Method For Dispersive Equation

A Comparison Among Homotopy Perturbation Method And The Decomposition Method With The Variational Iteration Method For Dispersive Equation Inernaional Jornal of Basic & Applied Sciences IJBAS-IJENS Vol:9 No: A Comparison Among Homoopy Perrbaion Mehod And The Decomposiion Mehod Wih The Variaional Ieraion Mehod For Dispersive Eqaion Hasan BULUT*

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

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

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

PH2130 Mathematical Methods Lab 3. z x

PH2130 Mathematical Methods Lab 3. z x PH130 Mahemaical Mehods Lab 3 This scrip shold keep yo bsy for he ne wo weeks. Yo shold aim o creae a idy and well-srcred Mahemaica Noebook. Do inclde plenifl annoaions o show ha yo know wha yo are doing,

More information

A Model for the Expansion of the Universe

A Model for the Expansion of the Universe Issue 2 April PROGRESS IN PHYSICS Volume 0 204 A Model for he Expansion of he Universe Nilon Penha Silva Deparmeno de Física Reired Professor, Universidade Federal de Minas Gerais, Belo Horizone, MG, Brazil

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

Computers and Mathematics with Applications

Computers and Mathematics with Applications Compers and Mahemaics wih Applicaions 59 (00) 80 809 Conens liss available a ScienceDirec Compers and Mahemaics wih Applicaions jornal homepage: www.elsevier.com/locae/camwa Solving fracional bondary vale

More information

A Direct Method for Solving Nonlinear PDEs and. New Exact Solutions for Some Examples

A Direct Method for Solving Nonlinear PDEs and. New Exact Solutions for Some Examples In. J. Conemp. Mah. Sciences, Vol. 6, 011, no. 46, 83-90 A Direc Mehod for Solving Nonlinear PDEs and New Eac Solions for Some Eamples Ameina S. Nseir Jordan Universiy of Science and Technology Deparmen

More information

An impact of noise on invariant manifolds in nonlinear dynamical systems

An impact of noise on invariant manifolds in nonlinear dynamical systems JOURNAL OF MATHEMATICAL PHYSICS 51, 4272 21 An impac of noise on invarian manifolds in nonlinear dynamical sysems X Sn, a Jinqiao Dan, and Xiaofan Li Deparmen of Applied Mahemaics, Illinois Insie of Technology,

More information

ON THE OSCILLATION OF THIRD ORDER FUNCTIONAL DIFFERENTIAL EQUATIONS. Cairo University, Orman, Giza 12221, Egypt

ON THE OSCILLATION OF THIRD ORDER FUNCTIONAL DIFFERENTIAL EQUATIONS. Cairo University, Orman, Giza 12221, Egypt a 1/α s)ds < Indian J. pre appl. Mah., 396): 491-507, December 2008 c Prined in India. ON THE OSCILLATION OF THIRD ORDER FUNCTIONAL DIFFERENTIAL EQUATIONS SAID R. GRACE 1, RAVI P. AGARWAL 2 AND MUSTAFA

More information

Anomalous transport regimes and asymptotic concentration distributions in the presence of advection and diffusion on a comb structure

Anomalous transport regimes and asymptotic concentration distributions in the presence of advection and diffusion on a comb structure Anomalos ranspor regimes and asympoic concenraion disribions in he presence of advecion and diffsion on a comb srcre Olga A. Dvoreskaya and Peer S. Kondraenko Nclear Safey Insie, Rssian Academy of Sciences,

More information

1 First Order Partial Differential Equations

1 First Order Partial Differential Equations Firs Order Parial Differenial Eqaions The profond sdy of nare is he mos ferile sorce of mahemaical discoveries. - Joseph Forier (768-830). Inrodcion We begin or sdy of parial differenial eqaions wih firs

More information

DESIGN OF TENSION MEMBERS

DESIGN OF TENSION MEMBERS CHAPTER Srcral Seel Design LRFD Mehod DESIGN OF TENSION MEMBERS Third Ediion A. J. Clark School of Engineering Deparmen of Civil and Environmenal Engineering Par II Srcral Seel Design and Analysis 4 FALL

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

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

Probing anomalous off-shell electroweak dipole moment of the top quark

Probing anomalous off-shell electroweak dipole moment of the top quark EPJ Web of Conferences 18, 42 (17) QFTHEP 17 DOI: 1.11/epjconf/171842 Probing anomalos off-shell elecroweak dipole momen of he op qark Edward Boos 1,2, Vyacheslav Bnichev 1, Lev Ddko 1, and Georgii Voronikov

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

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

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

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du.

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du. MATH 3B: MIDTERM REVIEW JOE HUGHES. Inegraion by Pars. Evaluae 3 e. Soluion: Firs make he subsiuion u =. Then =, hence 3 e = e = ue u Now inegrae by pars o ge ue u = ue u e u + C and subsiue he definiion

More information

Asymptotic Solution of the Anti-Plane Problem for a Two-Dimensional Lattice

Asymptotic Solution of the Anti-Plane Problem for a Two-Dimensional Lattice Asympoic Solion of he Ani-Plane Problem for a Two-Dimensional Laice N.I. Aleksandrova N.A. Chinakal Insie of Mining, Siberian Branch, Rssian Academy of Sciences, Krasnyi pr. 91, Novosibirsk, 6391 Rssia,

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

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

THE SINE INTEGRAL. x dt t

THE SINE INTEGRAL. x dt t THE SINE INTEGRAL As one learns in elemenary calculus, he limi of sin(/ as vanishes is uniy. Furhermore he funcion is even and has an infinie number of zeros locaed a ±n for n1,,3 Is plo looks like his-

More information

Non-Gaussian Test Models for Prediction and State Estimation with Model Errors

Non-Gaussian Test Models for Prediction and State Estimation with Model Errors Non-Gassian Tes Models for Predicion and Sae Esimaion wih Model Errors M. Branicki, N. Chen, and A.J. Majda Deparmen of Mahemaics and Cener for Amosphere Ocean Science, Coran Insie of Mahemaical Sciences

More information

Modelling Traffic Flow with Constant Speed using the Galerkin Finite Element Method

Modelling Traffic Flow with Constant Speed using the Galerkin Finite Element Method Proceedings of he World Congress on Engineering 29 Vol II WCE 29, Jly - 3, 29, London, U.K. Modelling Traffic Flow wih Consan Speed sing he Galerin Finie Elemen Mehod Wesley Celemans, Magd A. Wahab, Kr

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

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

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

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

first-order circuit Complete response can be regarded as the superposition of zero-input response and zero-state response.

first-order circuit Complete response can be regarded as the superposition of zero-input response and zero-state response. Experimen 4:he Sdies of ransiional processes of 1. Prpose firs-order circi a) Use he oscilloscope o observe he ransiional processes of firs-order circi. b) Use he oscilloscope o measre he ime consan of

More information

Dispersive Systems. 1) Schrödinger equation 2) Cubic Schrödinger 3) KdV 4) Discreterised hyperbolic equation 5) Discrete systems.

Dispersive Systems. 1) Schrödinger equation 2) Cubic Schrödinger 3) KdV 4) Discreterised hyperbolic equation 5) Discrete systems. Dispersive Sysems 1) Schrödinger eqaion ) Cbic Schrödinger 3) KdV 4) Discreerised hyperbolic eqaion 5) Discree sysems KdV + + ε =, = ( ) ( ) d d + = d d =, =. ( ) = ( ) DISCONTINUITY, prescribed cri Collision

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

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

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.

More information

Conservation Laws and Hamiltonian Symmetries of Whitham-Broer-Kaup Equations

Conservation Laws and Hamiltonian Symmetries of Whitham-Broer-Kaup Equations Indian Jornal of Science and Technology Vol 8( 78 84 Janary 05 ISSN (Prin : 0974-84 ISSN (Online : 0974-545 DOI : 0.7485/ijs/05/8i/47809 Conseraion Laws and Hamilonian Symmeries of Whiham-Broer-Kap Eqaions

More information

Integration Over Manifolds with Variable Coordinate Density

Integration Over Manifolds with Variable Coordinate Density Inegraion Over Manifolds wih Variable Coordinae Densiy Absrac Chrisopher A. Lafore clafore@gmail.com In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure

More information

Operators related to the Jacobi setting, for all admissible parameter values

Operators related to the Jacobi setting, for all admissible parameter values Operaors relaed o he Jacobi seing, for all admissible parameer values Peer Sjögren Universiy of Gohenburg Join work wih A. Nowak and T. Szarek Alba, June 2013 () 1 / 18 Le Pn α,β be he classical Jacobi

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

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

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species

More information

Average Number of Lattice Points in a Disk

Average Number of Lattice Points in a Disk Average Number of Laice Poins in a Disk Sujay Jayakar Rober S. Sricharz Absrac The difference beween he number of laice poins in a disk of radius /π and he area of he disk /4π is equal o he error in he

More information

Huazhong Tang 1 and Gerald Warnecke Introduction ANOTEON(2K + 1)-POINT CONSERVATIVE MONOTONE SCHEMES

Huazhong Tang 1 and Gerald Warnecke Introduction ANOTEON(2K + 1)-POINT CONSERVATIVE MONOTONE SCHEMES ESAIM: MAN Vol. 38, N o, 4, pp. 345 357 DOI:.5/man:46 ESAIM: Mahemaical Modelling and Nmerical Analysis ANOTEON(K + )-POINT CONSERVATIVE MONOTONE SCHEMES Hazhong Tang and Gerald Warnecke Absrac. Firs order

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

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

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

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that ODEs, Homework #4 Soluions. Check ha y ( = is a soluion of he second-order ODE ( cos sin y + y sin y sin = 0 and hen use his fac o find all soluions of he ODE. When y =, we have y = and also y = 0, so

More information

ON JENSEN S INEQUALITY FOR g-expectation

ON JENSEN S INEQUALITY FOR g-expectation Chin. Ann. Mah. 25B:3(2004),401 412. ON JENSEN S INEQUALITY FOR g-expectation JIANG Long CHEN Zengjing Absrac Briand e al. gave a conerexample showing ha given g, Jensen s ineqaliy for g-expecaion sally

More information

Hf(x) := 1 π p.v. f(y)

Hf(x) := 1 π p.v. f(y) NOTES ON THE NUMERICAL SOLUTION OF THE BENJAMIN EQUATION arxiv:131.2749v3 [mah.na] 22 May 214 VASSILIOS A. DOUGALIS, ANGEL DURAN, AND DIMITRIOS MITSOTAKIS Absrac. In his paper we consider he Benjamin eqaion,

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

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

Effects of Coordinate Curvature on Integration

Effects of Coordinate Curvature on Integration Effecs of Coordinae Curvaure on Inegraion Chrisopher A. Lafore clafore@gmail.com Absrac In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure of he manifold

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

Generating Solutions of Einstein s Equations

Generating Solutions of Einstein s Equations WDS'09 Proceedings of Conribued Papers, Par III, 13 19, 2009. ISBN 978-80-7378-103-3 MATFYZPRESS Generaing Soluions of Einsein s Equaions J. Hruška Charles Universiy, Faculy of Mahemaics and Physics, Prague,

More information

ENGI 9420 Engineering Analysis Assignment 2 Solutions

ENGI 9420 Engineering Analysis Assignment 2 Solutions ENGI 940 Engineering Analysis Assignmen Soluions 0 Fall [Second order ODEs, Laplace ransforms; Secions.0-.09]. Use Laplace ransforms o solve he iniial value problem [0] dy y, y( 0) 4 d + [This was Quesion

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

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

Chapter 8 The Complete Response of RL and RC Circuits

Chapter 8 The Complete Response of RL and RC Circuits Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior

More information

Late-time tails of self-gravitating waves

Late-time tails of self-gravitating waves Late-time tails of self-gravitating waves (non-rigorous quantitative analysis) Piotr Bizoń Jagiellonian University, Kraków Based on joint work with Tadek Chmaj and Andrzej Rostworowski Outline: Motivation

More information

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

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

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

Well-posedness of the generalized Proudman-Johnson equation without viscosity. Hisashi Okamoto RIMS, Kyoto University

Well-posedness of the generalized Proudman-Johnson equation without viscosity. Hisashi Okamoto RIMS, Kyoto University Well-posedness o he generalized Prodman-Johnson eqaion wiho viscosiy Hisashi Okamoo RIMS, Kyoo Universiy okamoo@krims.kyoo-.ac.jp Generalized Prodman-Johnson eqaion Proposed in by Zh and O. in order o

More information

EE650R: Reliability Physics of Nanoelectronic Devices Lecture 9:

EE650R: Reliability Physics of Nanoelectronic Devices Lecture 9: EE65R: Reliabiliy Physics of anoelecronic Devices Lecure 9: Feaures of Time-Dependen BTI Degradaion Dae: Sep. 9, 6 Classnoe Lufe Siddique Review Animesh Daa 9. Background/Review: BTI is observed when he

More information

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

KEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow KEY Mah 334 Miderm III Fall 28 secions and 3 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering

More information

Riemann Function and Methods of Group Analysis

Riemann Function and Methods of Group Analysis American Research Jornal of Mahemaics Original Aricle ISSN 378-74X Volme Isse 3 5 Riemann Fncion and Mehods of Grop Analsis Akimov Andre Chernov Igor Abdllina Rfina 3 4533 Serliamak Rssia Lenina sree 47A

More information

Earthquake, Volcano and Tsunami

Earthquake, Volcano and Tsunami A. Merapi Volcano Erpion Earhqake, Volcano and Tsnami Qesion Answer Marks A. Using Black s Principle he eqilibrim emperare can be obained Ths,.5 A. For ideal gas, pv e e RTe, hs.3 A.3 The relaive velociy

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

Electrical and current self-induction

Electrical and current self-induction Elecrical and curren self-inducion F. F. Mende hp://fmnauka.narod.ru/works.hml mende_fedor@mail.ru Absrac The aricle considers he self-inducance of reacive elemens. Elecrical self-inducion To he laws of

More information

Generalized Chebyshev polynomials

Generalized Chebyshev polynomials Generalized Chebyshev polynomials Clemene Cesarano Faculy of Engineering, Inernaional Telemaic Universiy UNINETTUNO Corso Viorio Emanuele II, 39 86 Roma, Ialy email: c.cesarano@unineunouniversiy.ne ABSTRACT

More information

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Applied Mathematics Letters. Oscillation results for fourth-order nonlinear dynamic equations

Applied Mathematics Letters. Oscillation results for fourth-order nonlinear dynamic equations Applied Mahemaics Leers 5 (0) 058 065 Conens liss available a SciVerse ScienceDirec Applied Mahemaics Leers jornal homepage: www.elsevier.com/locae/aml Oscillaion resls for forh-order nonlinear dynamic

More information

Exact solitary-wave Special Solutions for the Nonlinear Dispersive K(m,n) Equations by Means of the Homotopy Analysis Method

Exact solitary-wave Special Solutions for the Nonlinear Dispersive K(m,n) Equations by Means of the Homotopy Analysis Method Available a hp://pva.ed/aa Appl. Appl. Mah. ISSN: 93-9466 Special Isse No. (Ags ) pp. 8 93 Applicaions Applied Maheaics: An Inernaional Jornal (AAM) Eac soliary-wave Special Solions for he Nonlinear Dispersive

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

Optimal Stochastic Control in Continuous Time with Wiener Processes: General Results and Applications to Optimal Wildlife Management

Optimal Stochastic Control in Continuous Time with Wiener Processes: General Results and Applications to Optimal Wildlife Management Iranian Jornal of Operaions Research Vol. 8, No., 017, pp. 58-67 Downloaded from iors.ir a 1:47 +0330 on Wednesday November 14h 018 [ DOI: 10.95/iors.8..58 ] Opimal Sochasic Conrol in Coninos Time wih

More information

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation. Mah 36. Rumbos Spring 1 1 Soluions o Assignmen #6 1. Suppose he growh of a populaion is governed by he differenial equaion where k is a posiive consan. d d = k (a Explain why his model predics ha he populaion

More information

4 3 a b (C) (a 2b) (D) (2a 3b)

4 3 a b (C) (a 2b) (D) (2a 3b) * A balloon is moving verically pwards wih a velociy of 9 m/s. A sone is dropped from i and i reaches he grond in 10 sec. The heigh of he balloon when he sone was dropped is (ake g = 9.8 ms - ) (a) 100

More information

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick

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

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial

More information

Optimizing heat exchangers

Optimizing heat exchangers Opimizing hea echangers Jean-Luc Thiffeaul Deparmen of Mahemaics, Universiy of Wisconsin Madison, 48 Lincoln Dr., Madison, WI 5376, USA wih: Florence Marcoe, Charles R. Doering, William R. Young (Daed:

More information

Math 10B: Mock Mid II. April 13, 2016

Math 10B: Mock Mid II. April 13, 2016 Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

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

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

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering

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

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

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