Invasive Spread of an Advantageous Mutation under Preemptive Competition

Size: px
Start display at page:

Download "Invasive Spread of an Advantageous Mutation under Preemptive Competition"

Transcription

1 Invasive Spread of an Advanageous Muaion under Preempive Compeiion Lauren O Malley, Joseph Yasi, James Basham, Balázs Kozma, and G. Korniss (Rensselaer) Andy Allsad and Thomas Caraco (SUNY Albany) Zolán Rácz (Eövös Univ.) Suppored by NSF, Research Corporaion

2 Moivaion invasive spread of a selecively favored muaion (advanageous allelle) R.A. Fisher, 37, Kolmogorov, Perovsky & Piscounov, 37 propagaing fron, iniially exiss and separaes he wo spaial regions occupied separaely by he wo alleles ( domain-wall moion, propagaion ino unsable saes) v + ( ) (FKPP) velociy selecion /marginal sabiliy: Aronson and Weinberger 78, Dee and Langer, 83, van Saarloos 87)

3 Model feaures Through muaions, an invasive allele appears in a habia originally dominaed by a common residen allele muaion is a rare sochasic process residens and invaders compee for common limiing resources hrough clonal propagaion (plans) compeiion is pre-empive (invader allele has an individual-level advanage, bu canno displace residens already presen, (Amarasekare, 003, Shurin e al., 004)

4 Laice Model 0: empy laice sie (available resource) allele : ( residen ) allele : ( invader ) n( x) 0 if allele is presen a sie x oherwise n( x) 0 if allele is presen a sie x oherwise common limiing resources excluded volume consrain A laice sie represens he minimum level of locally available resources required o susain an individual organism.

5 0: empy laice sie (available resource) : allele ( residen ) : allele ( invader ) local ransiion raes for an arbirary sie x: η( x η( x ϕ µ µ ) ) local spread of residens local spread of invaders forward-backward recurren muaion deah η ( x) i n 4 i x' nn( x) ( x') ϕ << µ < <

6 MF Equilibrium (, ) phase diagram 0 ϕ µ µ µ 0 µ ( ) ( ) ( ) ( ) φ µ φ µ + + ( ) ( ) µ µ ( µ) c

7 Saionary-Sae MC Simulaions (single-allele clonal plans: dispersal-limied exincion, Oborny e al., 005) (.65µ ) c * β ( β c )

8 Invasion ime (lifeime) ( ) i n (, ) i, i x L x ime-dependen global densiies ( ) τ / (firs-passage ime o a suiably chosen cu-off densiy) measable or quasi-equilibrium densiy of he residens

9 Single-cluser invasion L << R o! sochasic L L sysem wih p.b.c.

10 Single-cluser invasion L << R o ime 0 τ P no ( ) exp[ ( g ) / n ] for for > g g ( L I n ) g ~ L / v growh ime τ n + g n average ime beween nucleaion evens average lifeime Rikvold e al., 94 Richards e al., 95 Ramos e al., 99 Machado e al., 05 GK and Caraco, 05

11 Muli-cluser invasion L >> R o! self-averaging (near-deerminisic global densiies)

12 Muli-cluser invasion L >> R o 0 ime τ KJMA/Avrami s law (for homogeneous nucleaion, d): () e ln()( / τ ) 3 I( vτ ) τ ~ τ ~ τ ~ ( Iv /3 ) (average) lifeime R o ~ ( v / I / 3 ) average disance beween clusers

13 MC Resuls - Single-cluser invasion L µ 0.0 ϕ I n ~ ~ n ~ ϕ

14 MC Resuls - Muli-cluser invasion τ ~ ϕ 0.30 from Avrami s law: τ ~ I /3 ~ ϕ /3 From simulaions: τ ~ ϕ 0.30

15 Summary: Finie-size effecs slope slope 0.30

16 Surface/inerface properies

17 Curren work: Wave propagaion PDE approach vs. Mone Carlo ( ) ( ) ( ) ( ) 4 4 µ δ µ δ + + ) ( ),, ( µ µ v Propagaion ino unsable saes: Velociy selecion /marginal sabiliy (Aronson and Weinberger 78, Dee and Langer, 83, van Saarloos 87)

18 MC simulaion for he fron velociy 0 planar Circular wavefron

19 MC simulaion for he fron roughening 0

20 Summary and Oulook KJMA/Avrami s law applicable o a wide range of sysems: ferromagneic (Rikvold e al. 94, Ramos 99) and ferroelecric (Ishibashi & Takagi 7, Duiker & Beale 90) swiching flame propagaion in slow combusion (Karunen 98) chemical reacions (Machado e al. 04) invasive allele spread and ecological invasion (GK & Caraco, 04) asympoic linear spreading velociy v * (,,µ) comparison of coninuum PDE and discree Mone Carlo approaches properies of criical cluser R c (,,µ) G. Korniss and T. Caraco, J. Theor. Biol. 33, 37 (005). J.A. Yasi, G. Korniss, and T. Caraco, cond-ma/ hp://

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

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi

More information

Spatial Ecology: Lecture 4, Integrodifference equations. II Southern-Summer school on Mathematical Biology

Spatial Ecology: Lecture 4, Integrodifference equations. II Southern-Summer school on Mathematical Biology Spaial Ecology: Lecure 4, Inegrodifference equaions II Souhern-Summer school on Mahemaical Biology Inegrodifference equaions Diffusion models assume growh and dispersal occur a he same ime. When reproducion

More information

From Particles to Rigid Bodies

From Particles to Rigid Bodies Rigid Body Dynamics From Paricles o Rigid Bodies Paricles No roaions Linear velociy v only Rigid bodies Body roaions Linear velociy v Angular velociy ω Rigid Bodies Rigid bodies have boh a posiion and

More information

Modeling Economic Time Series with Stochastic Linear Difference Equations

Modeling Economic Time Series with Stochastic Linear Difference Equations A. Thiemer, SLDG.mcd, 6..6 FH-Kiel Universiy of Applied Sciences Prof. Dr. Andreas Thiemer e-mail: andreas.hiemer@fh-kiel.de Modeling Economic Time Series wih Sochasic Linear Difference Equaions Summary:

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Simulating models with heterogeneous agents

Simulating models with heterogeneous agents Simulaing models wih heerogeneous agens Wouer J. Den Haan London School of Economics c by Wouer J. Den Haan Individual agen Subjec o employmen shocks (ε i, {0, 1}) Incomplee markes only way o save is hrough

More information

2) Of the following questions, which ones are thermodynamic, rather than kinetic concepts?

2) Of the following questions, which ones are thermodynamic, rather than kinetic concepts? AP Chemisry Tes (Chaper 12) Muliple Choice (40%) 1) Which of he following is a kineic quaniy? A) Enhalpy B) Inernal Energy C) Gibb s free energy D) Enropy E) Rae of reacion 2) Of he following quesions,

More information

Limit theorems for discrete-time metapopulation models

Limit theorems for discrete-time metapopulation models MASCOS Insiu de Mahémaiques de Toulouse, June 2010 - Page 1 Limi heorems for discree-ime meapopulaion models Phil Polle Deparmen of Mahemaics The Universiy of Queensland hp://www.mahs.uq.edu.au/ pkp AUSTRALIAN

More information

Today in Physics 218: radiation reaction

Today in Physics 218: radiation reaction Today in Physics 18: radiaion reacion Radiaion reacion The Abraham-Lorenz formula; radiaion reacion force The pah of he elecron in oday s firs example (radial decay grealy exaggeraed) 6 March 004 Physics

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

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

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

Problem Set on Differential Equations

Problem Set on Differential Equations Problem Se on Differenial Equaions 1. Solve he following differenial equaions (a) x () = e x (), x () = 3/ 4. (b) x () = e x (), x (1) =. (c) xe () = + (1 x ()) e, x () =.. (An asse marke model). Le p()

More information

ME425/525: Advanced Topics in Building Science

ME425/525: Advanced Topics in Building Science ME425/525: Advanced Topics in Building Science Indoor environmenal qualiy for susainable buildings: Lecure 6 Dr. Ellio T. Gall, Ph.D. Lecure 6 Today s objecives o Error propagaion Apply o SS soluion (venilaion,

More information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

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

Front propagation in A+B\2A reaction under subdiffusion

Front propagation in A+B\2A reaction under subdiffusion Fron propagaion in A+B\2A reacion under subdiffusion D. Froemberg, H. Schmid-Marens, and I. M. Sokolov Insiu für Physik, Humbold-Universiä zu Berlin, Newonsrasse 15, 12489 Berlin, Germany F. Sagués Deparamen

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

More information

Analytic Model and Bilateral Approximation for Clocked Comparator

Analytic Model and Bilateral Approximation for Clocked Comparator Analyic Model and Bilaeral Approximaion for Clocked Comparaor M. Greians, E. Hermanis, G.Supols Insiue of, Riga, Lavia, e-mail: gais.supols@edi.lv Research is suppored by: 1) ESF projec Nr.1DP/1.1.1.2.0/09/APIA/VIAA/020,

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

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 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

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

Ordinary Differential Equations

Ordinary Differential Equations Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described

More information

Random Walk on Circle Imagine a Markov process governing the random motion of a particle on a circular

Random Walk on Circle Imagine a Markov process governing the random motion of a particle on a circular Random Walk on Circle Imagine a Markov process governing he random moion of a paricle on a circular laice: 1 2 γ γ γ The paricle moves o he righ or lef wih probabiliy γ and says where i is wih probabiliy

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

arxiv: v1 [cond-mat.stat-mech] 27 Nov 2012

arxiv: v1 [cond-mat.stat-mech] 27 Nov 2012 Emergence of flucuaing raveling fron soluions in macroscopic heory of noisy invasion frons Baruch Meerson and Arkady Vilenkin Racah Insiue of Physics, Hebrew Universiy of Jerusalem, Jerusalem 9194, Israel

More information

Numerical Solution of a Nonlinear Integro-Differential Equation

Numerical Solution of a Nonlinear Integro-Differential Equation EPJ Web of Conferences 18, 17 (16) DOI: 1.151/ epjconf/ 161817 C Owned by he auhors, published by EDP Sciences, 16 Numerical Soluion of a Nonlinear Inegro-Differenial Equaion Ján Buša 1,a, Michal Hnaič,3,4,b,

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

AP Calculus BC Chapter 10 Part 1 AP Exam Problems

AP Calculus BC Chapter 10 Part 1 AP Exam Problems AP Calculus BC Chaper Par AP Eam Problems All problems are NO CALCULATOR unless oherwise indicaed Parameric Curves and Derivaives In he y plane, he graph of he parameric equaions = 5 + and y= for, is a

More information

Exciting traffic jams: Nonlinear phenomena behind traffic jam formation on highways

Exciting traffic jams: Nonlinear phenomena behind traffic jam formation on highways PHYSICAL REVIEW E 8, 462 29 Exciing raffic jams: Nonlinear phenomena behind raffic jam formaion on highways Gábor Orosz,, * R. Eddie Wilson, 2, Róber Szalai, 2, and Gábor Sépán 3, Deparmen of Mechanical

More information

arxiv: v1 [nlin.ps] 19 Aug 2016

arxiv: v1 [nlin.ps] 19 Aug 2016 Rock-paper-scissors played wihin compeing domains in predaor-prey games Darka Labavić 1, and Hildegard Meyer-Ormanns 1 1 Physics and Earh Sciences, Jacobs Universiy Bremen, P.O.Bo 750561, 28725 Bremen,

More information

The Maxwell Equations, the Lorentz Field and the Electromagnetic Nanofield with Regard to the Question of Relativity

The Maxwell Equations, the Lorentz Field and the Electromagnetic Nanofield with Regard to the Question of Relativity The Maxwell Equaions, he Lorenz Field and he Elecromagneic Nanofield wih Regard o he Quesion of Relaiviy Daniele Sasso * Absrac We discuss he Elecromagneic Theory in some main respecs and specifically

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

Calculus and Seismic Analysis November Ken Mark Bechtel Project Manager

Calculus and Seismic Analysis November Ken Mark Bechtel Project Manager Calculus and Seismic Analysis November 2007 Ken Mark Bechel Projec Manager BART Today 5 Lines 104 Miles Underground 37 Miles Aerial 23 Miles Surface 44 Miles 43 Saions 330,000 Daily Ridership 150,000 Daily

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

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4) Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion

More information

School and Workshop on Market Microstructure: Design, Efficiency and Statistical Regularities March 2011

School and Workshop on Market Microstructure: Design, Efficiency and Statistical Regularities March 2011 2229-12 School and Workshop on Marke Microsrucure: Design, Efficiency and Saisical Regulariies 21-25 March 2011 Some mahemaical properies of order book models Frederic ABERGEL Ecole Cenrale Paris Grande

More information

Simplest nonequilibrium phase transition into an absorbing state

Simplest nonequilibrium phase transition into an absorbing state Simples nonequilibrium phase ransiion ino an absorbing sae A. C. Barao, 1 Juan A. Bonachela, 2 C. E. Fiore, 3 H. Hinrichsen, 1 and Miguel A. Muñoz 2 1 Universiä Würzburg, Fakulä für Physik und Asronomie,

More information

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS Andrei Tokmakoff, MIT Deparmen of Chemisry, 2/22/2007 2-17 2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS The mahemaical formulaion of he dynamics of a quanum sysem is no unique. So far we have described

More information

PHYSICS 220 Lecture 02 Motion, Forces, and Newton s Laws Textbook Sections

PHYSICS 220 Lecture 02 Motion, Forces, and Newton s Laws Textbook Sections PHYSICS 220 Lecure 02 Moion, Forces, and Newon s Laws Texbook Secions 2.2-2.4 Lecure 2 Purdue Universiy, Physics 220 1 Overview Las Lecure Unis Scienific Noaion Significan Figures Moion Displacemen: Δx

More information

ENV 6015 Solution to Mixing Problem Set

ENV 6015 Solution to Mixing Problem Set EN 65 Soluion o ixing Problem Se. A slug of dye ( ) is injeced ino a single ank wih coninuous mixing. The flow in and ou of he ank is.5 gpm. The ank volume is 5 gallons. When will he dye concenraion equal

More information

Variable acceleration, Mixed Exercise 11

Variable acceleration, Mixed Exercise 11 Variable acceleraion, Mixed Exercise 11 1 a v 1 P is a res when v 0. 0 1 b s 0 0 v d (1 ) 1 0 1 0 7. The disance ravelled by P is 7. m. 1 a v 6+ a d v 6 + When, a 6+ 0 The acceleraion of P when is 0 m

More information

Introduction to choice over time

Introduction to choice over time Microeconomic Theory -- Choice over ime Inroducion o choice over ime Individual choice Income and subsiuion effecs 7 Walrasian equilibrium ineres rae 9 pages John Riley Ocober 9, 08 Microeconomic Theory

More information

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II Module 3: The Damped Oscillaor-II Lecure 3: The Damped Oscillaor-II 3. Over-damped Oscillaions. This refers o he siuaion where β > ω (3.) The wo roos are and α = β + α 2 = β β 2 ω 2 = (3.2) β 2 ω 2 = 2

More information

Keywords: competition models; density-dependence; ecology; population dynamics; predation models; stochastic models UNESCO EOLSS

Keywords: competition models; density-dependence; ecology; population dynamics; predation models; stochastic models UNESCO EOLSS POPULATIO MODELS Michael B. Bonsall Deparmen of Zoology, Universiy of Oxford, Oxford, UK Keywords: compeiion models; densiy-dependence; ecology; populaion dynamics; predaion models; sochasic models Conens.

More information

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17 EES 16A Designing Informaion Devices and Sysems I Spring 019 Lecure Noes Noe 17 17.1 apaciive ouchscreen In he las noe, we saw ha a capacior consiss of wo pieces on conducive maerial separaed by a nonconducive

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

MEI STRUCTURED MATHEMATICS 4758

MEI STRUCTURED MATHEMATICS 4758 OXFORD CAMBRIDGE AND RSA EXAMINATIONS Advanced Subsidiary General Cerificae of Educaion Advanced General Cerificae of Educaion MEI STRUCTURED MATHEMATICS 4758 Differenial Equaions Thursday 5 JUNE 006 Afernoon

More information

RL Lecture 7: Eligibility Traces. R. S. Sutton and A. G. Barto: Reinforcement Learning: An Introduction 1

RL Lecture 7: Eligibility Traces. R. S. Sutton and A. G. Barto: Reinforcement Learning: An Introduction 1 RL Lecure 7: Eligibiliy Traces R. S. Suon and A. G. Baro: Reinforcemen Learning: An Inroducion 1 N-sep TD Predicion Idea: Look farher ino he fuure when you do TD backup (1, 2, 3,, n seps) R. S. Suon and

More information

h[n] is the impulse response of the discrete-time system:

h[n] is the impulse response of the discrete-time system: Definiion Examples Properies Memory Inveribiliy Causaliy Sabiliy Time Invariance Lineariy Sysems Fundamenals Overview Definiion of a Sysem x() h() y() x[n] h[n] Sysem: a process in which inpu signals are

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

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

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

IS 709/809: Computational Methods in IS Research. Queueing Theory Introduction

IS 709/809: Computational Methods in IS Research. Queueing Theory Introduction IS 709/809: Compuaional Mehods in IS Research Queueing Theory Inroducion Nirmalya Roy Deparmen of Informaion Sysems Universiy of Maryland Balimore Couny www.umbc.edu Inroducion: Saisics of hings Waiing

More information

Crossing the Bridge between Similar Games

Crossing the Bridge between Similar Games Crossing he Bridge beween Similar Games Jan-David Quesel, Marin Fränzle, and Werner Damm Universiy of Oldenburg, Deparmen of Compuing Science, Germany CMACS Seminar CMU, Pisburgh, PA, USA 2nd December

More information

Signals and Systems Profs. Byron Yu and Pulkit Grover Fall Midterm 1 Solutions

Signals and Systems Profs. Byron Yu and Pulkit Grover Fall Midterm 1 Solutions 8-90 Signals and Sysems Profs. Byron Yu and Pulki Grover Fall 07 Miderm Soluions Name: Andrew ID: Problem Score Max 0 8 4 6 5 0 6 0 7 8 9 0 6 Toal 00 Miderm Soluions. (0 poins) Deermine wheher he following

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

V AK (t) I T (t) I TRM. V AK( full area) (t) t t 1 Axial turn-on. Switching losses for Phase Control and Bi- Directionally Controlled Thyristors

V AK (t) I T (t) I TRM. V AK( full area) (t) t t 1 Axial turn-on. Switching losses for Phase Control and Bi- Directionally Controlled Thyristors Applicaion Noe Swiching losses for Phase Conrol and Bi- Direcionally Conrolled Thyrisors V AK () I T () Causing W on I TRM V AK( full area) () 1 Axial urn-on Plasma spread 2 Swiching losses for Phase Conrol

More information

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion

More information

Lecture Outline. Introduction Transmission Line Equations Transmission Line Wave Equations 8/10/2018. EE 4347 Applied Electromagnetics.

Lecture Outline. Introduction Transmission Line Equations Transmission Line Wave Equations 8/10/2018. EE 4347 Applied Electromagnetics. 8/10/018 Course Insrucor Dr. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: rcrumpf@uep.edu EE 4347 Applied Elecromagneics Topic 4a Transmission Line Equaions Transmission These Line noes

More information

Chapter 7 Response of First-order RL and RC Circuits

Chapter 7 Response of First-order RL and RC Circuits Chaper 7 Response of Firs-order RL and RC Circuis 7.- The Naural Response of RL and RC Circuis 7.3 The Sep Response of RL and RC Circuis 7.4 A General Soluion for Sep and Naural Responses 7.5 Sequenial

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

Online Appendix to Solution Methods for Models with Rare Disasters

Online Appendix to Solution Methods for Models with Rare Disasters Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,

More information

arxiv: v1 [physics.optics] 9 Apr 2018

arxiv: v1 [physics.optics] 9 Apr 2018 Enhancemen of ransmission qualiy in solion-based opical waveguide sysems by frequency dependen linear gain-loss and he Raman self-frequency shif Avner Peleg, Debananda Chakrabory 2 Deparmen of Exac Sciences,

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 3 Signals & Sysems Prof. Mark Fowler Noe Se # Wha are Coninuous-Time Signals??? /6 Coninuous-Time Signal Coninuous Time (C-T) Signal: A C-T signal is defined on he coninuum of ime values. Tha is:

More information

arxiv:cond-mat/ May 2002

arxiv:cond-mat/ May 2002 -- uadrupolar Glass Sae in para-hydrogen and orho-deuerium under pressure. T.I.Schelkacheva. arxiv:cond-ma/5538 6 May Insiue for High Pressure Physics, Russian Academy of Sciences, Troisk 49, Moscow Region,

More information

Reliability of Technical Systems

Reliability of Technical Systems eliabiliy of Technical Sysems Main Topics Inroducion, Key erms, framing he problem eliabiliy parameers: Failure ae, Failure Probabiliy, Availabiliy, ec. Some imporan reliabiliy disribuions Componen reliabiliy

More information

Cryptanalysis of Reduced NORX initially discussed at ASK 2016

Cryptanalysis of Reduced NORX initially discussed at ASK 2016 Crypanalysis of Reduced NORX iniiy discussed a ASK 2016 Nasour Bagheri (SRTTU, Iran), Tao Huang (NTU, Singaopre), Keing Jia (Tsinghua Univ., Sae Key Lab, China), Florian Mendel (TUGraz, Ausria), Yu Sasaki

More information

in Engineering Prof. Dr. Michael Havbro Faber ETH Zurich, Switzerland Swiss Federal Institute of Technology

in Engineering Prof. Dr. Michael Havbro Faber ETH Zurich, Switzerland Swiss Federal Institute of Technology Risk and Saey in Engineering Pro. Dr. Michael Havbro Faber ETH Zurich, Swizerland Conens o Today's Lecure Inroducion o ime varian reliabiliy analysis The Poisson process The ormal process Assessmen o he

More information

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

More information

The English physicist James Joule was perhaps the first to use the

The English physicist James Joule was perhaps the first to use the The sochasic pison problem G. Lin, C. H. Su, and G. E. Karniadakis* Division of Applied Mahemaics, Brown Universiy, 182 George Sree, Providence, RI 02912 Edied by Richard V. Kadison, Universiy of Pennsylvania,

More information

Online publication date: 04 August 2010 PLEASE SCROLL DOWN FOR ARTICLE

Online publication date: 04 August 2010 PLEASE SCROLL DOWN FOR ARTICLE This aricle was downloaded by: [Lanchier, Nicolas] On: 6 Augus 2010 Access deails: Access Deails: [subscripion number 925212420] Publisher Taylor & Francis Informa Ld Regisered in England and Wales Regisered

More information

1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population

1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population Scienific Comuing I Module : Poulaion Modelling Coninuous Models Michael Bader Par I ODE Models Lehrsuhl Informaik V Winer 7/ Discree vs. Coniuous Models d d = F,,...) ) =? discree model: coninuous model:

More information

arxiv: v1 [math.na] 23 Feb 2016

arxiv: v1 [math.na] 23 Feb 2016 EPJ Web of Conferences will be se by he publisher DOI: will be se by he publisher c Owned by he auhors, published by EDP Sciences, 16 arxiv:163.67v1 [mah.na] 3 Feb 16 Numerical Soluion of a Nonlinear Inegro-Differenial

More information

Lecture Notes 3: Quantitative Analysis in DSGE Models: New Keynesian Model

Lecture Notes 3: Quantitative Analysis in DSGE Models: New Keynesian Model Lecure Noes 3: Quaniaive Analysis in DSGE Models: New Keynesian Model Zhiwei Xu, Email: xuzhiwei@sju.edu.cn The moneary policy plays lile role in he basic moneary model wihou price sickiness. We now urn

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

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy Answer 4 of he following 5 quesions. 1. Consider a pure-exchange economy wih sochasic endowmens. The sae of he economy in period, 0,1,..., is he hisory of evens s ( s0, s1,..., s ). The iniial sae is given.

More information

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive. LECTURE 3 Linear/Nonnegaive Marix Models x ( = Px ( A= m m marix, x= m vecor Linear sysems of difference equaions arise in several difference conexs: Linear approximaions (linearizaion Perurbaion analysis

More information

Speaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis

Speaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis Speaker Adapaion Techniques For Coninuous Speech Using Medium and Small Adapaion Daa Ses Consaninos Boulis Ouline of he Presenaion Inroducion o he speaker adapaion problem Maximum Likelihood Sochasic Transformaions

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

Fishing limits and the Logistic Equation. 1

Fishing limits and the Logistic Equation. 1 Fishing limis and he Logisic Equaion. 1 1. The Logisic Equaion. The logisic equaion is an equaion governing populaion growh for populaions in an environmen wih a limied amoun of resources (for insance,

More information

6.003 Homework 1. Problems. Due at the beginning of recitation on Wednesday, February 10, 2010.

6.003 Homework 1. Problems. Due at the beginning of recitation on Wednesday, February 10, 2010. 6.003 Homework Due a he beginning of reciaion on Wednesday, February 0, 200. Problems. Independen and Dependen Variables Assume ha he heigh of a waer wave is given by g(x v) where x is disance, v is velociy,

More information

Case Study 1: Spreading Infectious Disease in Rabbits

Case Study 1: Spreading Infectious Disease in Rabbits Case Sudy 1: Spreading Infecious Disease in Rabbis Huaiong Huang Deparmen of Mahemaics & Saisics, York Universiy, Torono, ON, Canada Sepember 28, 24 Absrac In his repor we sudy a basic model for he spreading

More information

Sequential Importance Resampling (SIR) Particle Filter

Sequential Importance Resampling (SIR) Particle Filter Paricle Filers++ Pieer Abbeel UC Berkeley EECS Many slides adaped from Thrun, Burgard and Fox, Probabilisic Roboics 1. Algorihm paricle_filer( S -1, u, z ): 2. Sequenial Imporance Resampling (SIR) Paricle

More information

University of Cyprus Biomedical Imaging and Applied Optics. Appendix. DC Circuits Capacitors and Inductors AC Circuits Operational Amplifiers

University of Cyprus Biomedical Imaging and Applied Optics. Appendix. DC Circuits Capacitors and Inductors AC Circuits Operational Amplifiers Universiy of Cyprus Biomedical Imaging and Applied Opics Appendix DC Circuis Capaciors and Inducors AC Circuis Operaional Amplifiers Circui Elemens An elecrical circui consiss of circui elemens such as

More information

Oscillation death in diffusively coupled oscillators by local repulsive link

Oscillation death in diffusively coupled oscillators by local repulsive link Oscillaion deah in diffusively coupled oscillaors by local repulsive link C.R. Hens, Olasunkanmi I.Olusola,, Pinaki Pal 3, Syamal K. Dana CSIR-Indian Insiue of Chemical Biology, Jadavpur, Kolkaa 73, India

More information

Stationary Distribution. Design and Analysis of Algorithms Andrei Bulatov

Stationary Distribution. Design and Analysis of Algorithms Andrei Bulatov Saionary Disribuion Design and Analysis of Algorihms Andrei Bulaov Algorihms Markov Chains 34-2 Classificaion of Saes k By P we denoe he (i,j)-enry of i, j Sae is accessible from sae if 0 for some k 0

More information

ψ(t) = V x (0)V x (t)

ψ(t) = V x (0)V x (t) .93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in

More information

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance Ch: Moion along a sraigh line Moion Posiion and Displacemen Average Velociy and Average Speed Insananeous Velociy and Speed Acceleraion Consan Acceleraion: A Special Case Anoher Look a Consan Acceleraion

More information

Single-loop System Reliability-based Topology Optimization Accounting for Statistical Dependence between Limit-states

Single-loop System Reliability-based Topology Optimization Accounting for Statistical Dependence between Limit-states 11 h US Naional Congress on Compuaional Mechanics Minneapolis, Minnesoa Single-loop Sysem Reliabiliy-based Topology Opimizaion Accouning for Saisical Dependence beween imi-saes Tam Nguyen, Norheasern Universiy

More information

CHE302 LECTURE VI DYNAMIC BEHAVIORS OF REPRESENTATIVE PROCESSES. Professor Dae Ryook Yang

CHE302 LECTURE VI DYNAMIC BEHAVIORS OF REPRESENTATIVE PROCESSES. Professor Dae Ryook Yang CHE30 LECTURE VI DYNAMIC BEHAVIORS OF REPRESENTATIVE PROCESSES Professor Dae Ryook Yang Fall 00 Dep. of Chemical and Biological Engineering orea Universiy CHE30 Process Dynamics and Conrol orea Universiy

More information

18 Biological models with discrete time

18 Biological models with discrete time 8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so

More information

( ) = Q 0. ( ) R = R dq. ( t) = I t

( ) = Q 0. ( ) R = R dq. ( t) = I t ircuis onceps The addiion of a simple capacior o a circui of resisors allows wo relaed phenomena o occur The observaion ha he ime-dependence of a complex waveform is alered by he circui is referred o as

More information

Mechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel,

Mechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel, Mechanical Faigue and Load-Induced Aging of Loudspeaker Suspension Wolfgang Klippel, Insiue of Acousics and Speech Communicaion Dresden Universiy of Technology presened a he ALMA Symposium 2012, Las Vegas

More information

THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES

THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES Kragujevac J. Sci. 3 () 7-4. UDC 53.5:536. 4 THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES Hazem A. Aia Dep. of Mahemaics, College of Science,King Saud Universiy

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information