ASTR415: Problem Set #5
|
|
- Sarah Alexander
- 6 years ago
- Views:
Transcription
1 ASTR45: Problem Se #5 Curran D. Muhlberger Universi of Marland (Daed: April 25, 27) Three ssems of coupled differenial equaions were sudied using inegraors based on Euler s mehod, a fourh-order Runge-Kua mehod, and he mehod. For even a basic ssem of wo smooh equaions, Euler s mehod performed ver poorl and was slow o converge, as he error scaled as O(h). The Runge-Kua mehod performed well for all hree ssems, converging quickl even for moderae ime seps and ehibiing an error ha scaled as O(h 4 ). The mehod, where applicable, also converged quickl o a soluion, alhough he error onl scaled as O(h 2 ). The leapfrog mehod proved o be he mos sable, conserving he energ of an orbi hroughou a long inegraion. Finall, a bisecion echnique was used when suding a Loka-Volerra Predaor-Pre model o deermine ha he minimum huning rae o achieve eincion was q =.24. I. A ONE-DIMENSIONAL SECOND-ORDER ODE Consider he second-order ordinar differenial equaion d 2 = () d2 wih he iniial condiions () =, ẋ =. The soluion is given b () = sin(). Noe ha his is a conservaive equaion, since he second derivaive of is independen of and ẋ. Therefore, he mehod is a valid inegraion scheme. This second-order equaion can be wrien as a ssem of he coupled firs-order equaions d d = d d = (2) wih iniial condiions () =, () =. Here, he eac soluion for is given b () = cos(). This ssem is in a form ha can be inegraed using Euler s mehod or a Runge-Kua mehod. Euler s mehod, he fourh-order Runge-Kua mehod (as described in []), and he mehod can all be used wih a consan sep size h. When solving his equaion numericall, h was chosen o be,.3,.,.3, and., and he equaion was inegraed from = o = 5, represening abou 2.5 periods of he sine wave. The resuls from Euler s mehod are shown in figure, he resuls from he fourh-order Runge-Kua mehod are shown in figure 2, and he resuls from he mehod are shown in figure 3. In all cases, he eac soluion () = sin() is ploed as a dashed curve. Euler s mehod converges o he soluion ver slowl and canno mainain a consan ampliude even wih small values of h. The Runge-Kua mehod does a much beer job a approimaing he soluion, being visuall indisinguishable from he eac soluion as earl as h =.. Furhermore, he ampliude of is soluions does no appear o grow over ime on his ime scale. Finall, he soluions also converge ver quickl, being visuall idenical o he eac soluion b h =. and mainaining a consan ampliude. However, for he ver large ime sep h =, he soluion drifs ou of phase wih he eac soluion fairl quickl. This behavior was no observed for he Runge-Kua soluion using he same ime sep. For each of hese mehods, we epec he error a an ime in he soluion o decrease as some power of he sep size h. To confirm his behavior, (5) sin(5) is ploed agains h on a logarihmic scale in figure 4. For all hree mehods, he slope is consan for h.. In paricular, he slope for Euler s mehod is.2, he slope for he fourh-order Runge-Kua mehod is 3.92, and he slope for he mehod is 2.. From his, we can conclude ha he error in he inegraion is h for Euler s mehod, h 4 for he Runge-Kua mehod, and h 2 for he mehod. These are he epeced dependencies for he error for each of he mehods. Euler s mehod is a firs-order mehod, meaning ha he error in each sep is O(h 2 ). Since he error a f is he accumulaion of he error from ( f )/h seps, his final error is O(h), as was observed above. Similarl, he Runge-Kua mehod used is fourh-order, wih an error in each erm of O(h 5 ). Summed over all seps, he oal error should be O(h 4 ). Finall, he error in each erm of he mehod is O(h 3 ) for a oal accumulaed error of O(h 2 ).
2 2 Euler Mehod (h =.) Euler Mehod (h =.3) Euler Mehod (h =.) Euler Mehod (h =.3) Euler Mehod (h =.) FIG. : Soluions of obeing equaion 2 using Euler s mehod for various ime seps h. The eac soluion is represened b he dashed curve.
3 3 Runge-Kua Mehod (h =.) Runge-Kua Mehod (h =.3) Runge-Kua Mehod (h =.) Runge-Kua Mehod (h =.3) Runge-Kua Mehod (h =.) FIG. 2: Soluions of obeing equaion 2 using he fourh-order Runge-Kua mehod for various ime seps h. The eac soluion is represened b he dashed curve.
4 4 Mehod (h =.) Mehod (h =.3) Mehod (h =.) Mehod (h =.3) Mehod (h =.) FIG. 3: Soluions of equaion using he mehod for various ime seps h. The eac soluion is represened b he dashed curve.
5 5 Error in Soluion a = 5 Euler Runge-Kua. Error e-4 e-6 e-8 e-.. h FIG. 4: The error of all hree mehods a = 5 for each ime sep h. II. A TWO-DIMENSIONAL ORBIT Consider he wo-dimensional orbi of a mass in he graviaional poenial given b Φ = ( ) /2 (3) The resuling acceleraion is given b ẍ = Φ. Wriing he componens eplicil, d 2 d 2 d 2 d 2 = dφ d = 2( ) 3/2 = dφ d = 2( ) 3/2 (4) This ssem of wo coupled second-order conservaive differenial equaions can be recas ino a ssem of he four coupled firs-order equaions d d = a da d = 2( ) 3/2 d d = b db d = 2( ) 3/2 (5) This equaion was solved using boh he fourh-order Runge-Kua mehod and he Leepfrog mehod. The iniial condiions were aken o be () =, ẋ() =, () =, and ẏ() =., and he soluions were inegraed ou o =. Sep sizes of,.5,.25, and. were used. The soluions using he Runge-Kua mehod are shown in figure 5, and he soluions from he mehod are ploed in figure 6. Neiher mehod achieves good convergence unil h =.25.
6 6 Runge-Kua Mehod (h =.) Runge-Kua Mehod (h =.5) Runge-Kua Mehod (h =.25) Runge-Kua Mehod (h =.) FIG. 5: Orbis from he Runge-Kua soluion of equaion 5 for various ime seps h. The energ per uni mass of he paricle is given b ε = ẋ2 + ẏ Φ(, ) (6) The velociies needed o calculae he energ are compued as a and b in he Runge-Kua soluion of equaion 5, bu mus be eraced from he mehod b resncing he offse velociies using half Euler seps. The energ of he orbi is ploed agains for boh solvers in figure 7. The energ for h = is meaningless for boh solvers, as neiher produced sable orbis for his ime sep. However, for all smaller values of h, he energ of he soluion is sead over long periods of ime while he energ of he Runge-Kua soluion seadil decreases wih ime. Thus, while he local variaions in he energ of he soluion ma be much larger han hose of he Runge-Kua soluion, he energ is more sable over long periods of ime, indicaing ha he mehod has auomaicall enforced energ conservaion for his ssem.
7 7 Mehod (h =.) Mehod (h =.5) Mehod (h =.25) Mehod (h =.) FIG. 6: Orbis from he soluion of equaion 4 for various ime seps h. III. THE LOTKA-VOLTERRA PREDATOR-PREY MODEL The populaions of wo species in a predaor-pre relaionship, such as rabbis and foes, can be described b he firs-order Loka-Volerra equaions ẋ = (A d) B ẏ = ( C e) + D (7) where and are he populaion densiies of he pre and predaor respecivel. A is he pre s reproducion rae, B is he pre s consumpion rae b he predaor, C is he predaor s deah rae, D is he predaor s populaion growh rae due o consumpion of he pre, and d and e are he huning raes of he pre and predaor respecivel. Here we consider an ecossem conaining rabbis and foes wih parameers A =, B =., C =.5, D =.3, and d = e =. The iniial condiions are aken o be () = 3 and () = 3. The fourh-order Runge-Kua inegraor was used o inegrae he soluion ou o = for ime seps,.5,.25, and.. Phase diagrams for each ime sep are shown in figure 8. Supposing ha boh species are huned a an equal rae q (so d = e = q), here should be some value of q such ha boh populaions become einc b = (where eincion is defined b, < 9 ). To find his value of q, a
8 8 Energ Conservaion (h =.) Energ Conservaion (h =.5).4.2 Runge-Kua Runge-Kua Energ -.4 Energ Energ Conservaion (h =.25) Energ Conservaion (h =.) Runge-Kua Runge-Kua Energ Energ FIG. 7: Energ Conservaion. bisecion-sle mehod was used. The value was firs brackeed b observing ha for q =, (), () > 9, and for q = 5, (), () < 9. The ssem was hen solved choosing q a he midpoin of his inerval. If boh species were found o be einc a =, he upper bound was se equal o he midpoin. Oherwise he lower bound was se equal o he midpoin. This procedure was repeaed unil he size of he inerval fell below 4, obaining 4 digis of accurac. Using a ime sep of h =., i was found ha q =.2398 is he minimum huning rae necessar o drive boh species o eincion b =. [] Press, William H., Saul A. Teukolsk, William T. Veerling, Brian P. Flanner. Numerical Recipes in FORTRAN 77: The Ar of Scienific Compuing (Volume of Forran Numerical Recipes). Second Ediion. Cambridge Universi Press, 2.
9 9 Runge-Kua Mehod (h =.) Runge-Kua Mehod (h =.5) Predaor Densi (Foes) Predaor Densi (Foes) Pre Densi (Rabbis) Pre Densi (Rabbis) Runge-Kua Mehod (h =.25) Runge-Kua Mehod (h =.) Predaor Densi (Foes) 2 5 Predaor Densi (Foes) Pre Densi (Rabbis) Pre Densi (Rabbis) FIG. 8: Phase diagrams for he Loka-Volerra Predaor-Pre model (equaion 7) inegraed using he Runge-Kua mehod wih various ime seps.
ME 391 Mechanical Engineering Analysis
Fall 04 ME 39 Mechanical Engineering Analsis Eam # Soluions Direcions: Open noes (including course web posings). No books, compuers, or phones. An calculaor is fair game. Problem Deermine he posiion of
More informationSimulation-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 informationProblem Set 7-7. dv V ln V = kt + C. 20. Assume that df/dt still equals = F RF. df dr = =
20. Assume ha df/d sill equals = F + 0.02RF. df dr df/ d F+ 0. 02RF = = 2 dr/ d R 0. 04RF 0. 01R 10 df 11. 2 R= 70 and F = 1 = = 0. 362K dr 31 21. 0 F (70, 30) (70, 1) R 100 Noe ha he slope a (70, 1) is
More informationSolutions 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 informationKEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering
More informationln 2 1 ln y x c y C x
Lecure 14 Appendi B: Some sample problems from Boas Here are some soluions o he sample problems assigned for Chaper 8 8: 6 Soluion: We wan o find he soluion o he following firs order equaion using separaion
More information4.6 One Dimensional Kinematics and Integration
4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of
More informationEXERCISES 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 informationMath Wednesday March 3, , 4.3: First order systems of Differential Equations Why you should expect existence and uniqueness for the IVP
Mah 2280 Wednesda March 3, 200 4., 4.3: Firs order ssems of Differenial Equaions Wh ou should epec eisence and uniqueness for he IVP Eample: Consider he iniial value problem relaed o page 4 of his eserda
More informationThe fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation
Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or
More information8. Basic RL and RC Circuits
8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics
More information1. 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 informationElementary Differential Equations and Boundary Value Problems
Elemenar Differenial Equaions and Boundar Value Problems Boce. & DiPrima 9 h Ediion Chaper 1: Inroducion 1006003 คณ ตศาสตร ว ศวกรรม 3 สาขาว ชาว ศวกรรมคอมพ วเตอร ป การศ กษา 1/2555 ผศ.ดร.อร ญญา ผศ.ดร.สมศ
More informationMA 366 Review - Test # 1
MA 366 Review - Tes # 1 Fall 5 () Resuls from Calculus: differeniaion formulas, implici differeniaion, Chain Rule; inegraion formulas, inegraion b pars, parial fracions, oher inegraion echniques. (1) Order
More informationMorning Time: 1 hour 30 minutes Additional materials (enclosed):
ADVANCED GCE 78/0 MATHEMATICS (MEI) Differenial Equaions THURSDAY JANUARY 008 Morning Time: hour 30 minues Addiional maerials (enclosed): None Addiional maerials (required): Answer Bookle (8 pages) Graph
More informationTheory 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 informationWeek 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)
Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =
More informationSome 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 informationwhere the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).
Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness
More informationTheory 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 informationChapter 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 informationOutline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?
of Topics wih MATLAB Shan He School for Compuaional Science Universi of Birmingham Module 6-3836: Compuaional Modelling wih MATLAB Wha will we learn from his lecure Aim of analsis: Aim of analsis. Some
More informationNumerical 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 informationSystem of Linear Differential Equations
Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y
More informationTHE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI
THE 2-BODY PROBLEM ROBERT J. VANDERBEI ABSTRACT. In his shor noe, we show ha a pair of ellipses wih a common focus is a soluion o he 2-body problem. INTRODUCTION. Solving he 2-body problem from scrach
More informationMath 2214 Solution Test 1A Spring 2016
Mah 14 Soluion Tes 1A Spring 016 sec Problem 1: Wha is he larges -inerval for which ( 4) = has a guaraneed + unique soluion for iniial value (-1) = 3 according o he Exisence Uniqueness Theorem? Soluion
More information23.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 informationOrdinary dierential equations
Chaper 5 Ordinary dierenial equaions Conens 5.1 Iniial value problem........................... 31 5. Forward Euler's mehod......................... 3 5.3 Runge-Kua mehods.......................... 36
More informationMath 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 informationOrdinary differential equations. Phys 750 Lecture 7
Ordinary differenial equaions Phys 750 Lecure 7 Ordinary Differenial Equaions Mos physical laws are expressed as differenial equaions These come in hree flavours: iniial-value problems boundary-value problems
More informationSecond-Order Differential Equations
WWW Problems and Soluions 3.1 Chaper 3 Second-Order Differenial Equaions Secion 3.1 Springs: Linear and Nonlinear Models www m Problem 3. (NonlinearSprings). A bod of mass m is aached o a wall b means
More informationu(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 information15. 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( ) 2. Review Exercise 2. cos θ 2 3 = = 2 tan. cos. 2 x = = x a. Since π π, = 2. sin = = 2+ = = cotx. 2 sin θ 2+
Review Eercise sin 5 cos sin an cos 5 5 an 5 9 co 0 a sinθ 6 + 4 6 + sin θ 4 6+ + 6 + 4 cos θ sin θ + 4 4 sin θ + an θ cos θ ( ) + + + + Since π π, < θ < anθ should be negaive. anθ ( + ) Pearson Educaion
More informationdt = 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 informationMath 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.
1 Mah 334 Tes 1 KEY Spring 21 Secion: 1 Insrucor: Sco Glasgow Daes: Ma 1 and 11. Do NOT wrie on his problem saemen bookle, excep for our indicaion of following he honor code jus below. No credi will be
More informationAnnouncements: Warm-up Exercise:
Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple
More informationLAB 6: SIMPLE HARMONIC MOTION
1 Name Dae Day/Time of Lab Parner(s) Lab TA Objecives LAB 6: SIMPLE HARMONIC MOTION To undersand oscillaion in relaion o equilibrium of conservaive forces To manipulae he independen variables of oscillaion:
More informationSecond Order Linear Differential Equations
Second Order Linear Differenial Equaions Second order linear equaions wih consan coefficiens; Fundamenal soluions; Wronskian; Exisence and Uniqueness of soluions; he characerisic equaion; soluions of homogeneous
More informationTHE 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 informationReview - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y
Review - Quiz # 1 (1) Solving Special Tpes of Firs Order Equaions I. Separable Equaions (SE). d = f() g() Mehod of Soluion : 1 g() d = f() (The soluions ma be given implicil b he above formula. Remember,
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More information23.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 information04. Kinetics of a second order reaction
4. Kineics of a second order reacion Imporan conceps Reacion rae, reacion exen, reacion rae equaion, order of a reacion, firs-order reacions, second-order reacions, differenial and inegraed rae laws, Arrhenius
More informationMA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions
MA 14 Calculus IV (Spring 016) Secion Homework Assignmen 1 Soluions 1 Boyce and DiPrima, p 40, Problem 10 (c) Soluion: In sandard form he given firs-order linear ODE is: An inegraing facor is given by
More informationExam 1 Solutions. 1 Question 1. February 10, Part (A) 1.2 Part (B) To find equilibrium solutions, set P (t) = C = dp
Exam Soluions Februar 0, 05 Quesion. Par (A) To find equilibrium soluions, se P () = C = = 0. This implies: = P ( P ) P = P P P = P P = P ( + P ) = 0 The equilibrium soluion are hus P () = 0 and P () =..
More informationPhysics 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 informationChapter 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 informationLAB # 2 - Equilibrium (static)
AB # - Equilibrium (saic) Inroducion Isaac Newon's conribuion o physics was o recognize ha despie he seeming compleiy of he Unierse, he moion of is pars is guided by surprisingly simple aws. Newon's inspiraion
More informationAP 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 informationMotion along a Straight Line
chaper 2 Moion along a Sraigh Line verage speed and average velociy (Secion 2.2) 1. Velociy versus speed Cone in he ebook: fer Eample 2. Insananeous velociy and insananeous acceleraion (Secions 2.3, 2.4)
More informationChapter 6. Systems of First Order Linear Differential Equations
Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh
More informationPhysics 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 informationFrom 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 informationMath 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 informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationHOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.
HOMEWORK # 2: MATH 2, SPRING 25 TJ HITCHMAN Noe: This is he las soluion se where I will describe he MATLAB I used o make my picures.. Exercises from he ex.. Chaper 2.. Problem 6. We are o show ha y() =
More informationarxiv:quant-ph/ v1 5 Jul 2004
Numerical Mehods for Sochasic Differenial Equaions Joshua Wilkie Deparmen of Chemisry, Simon Fraser Universiy, Burnaby, Briish Columbia V5A 1S6, Canada Sochasic differenial equaions (sdes) play an imporan
More informationIntroduction to Physical Oceanography Homework 5 - Solutions
Laure Zanna //5 Inroducion o Phsical Oceanograph Homework 5 - Soluions. Inerial oscillaions wih boom fricion non-selecive scale: The governing equaions for his problem are This ssem can be wrien as where
More informationTwo 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 informationClass 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 informationIB 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 information3.6 Derivatives as Rates of Change
3.6 Derivaives as Raes of Change Problem 1 John is walking along a sraigh pah. His posiion a he ime >0 is given by s = f(). He sars a =0from his house (f(0) = 0) and he graph of f is given below. (a) Describe
More informationWEEK-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 informationAPPM 2360 Homework Solutions, Due June 10
2.2.2: Find general soluions for he equaion APPM 2360 Homework Soluions, Due June 10 Soluion: Finding he inegraing facor, dy + 2y = 3e µ) = e 2) = e 2 Muliplying he differenial equaion by he inegraing
More informationSterilization D Values
Seriliaion D Values Seriliaion by seam consis of he simple observaion ha baceria die over ime during exposure o hea. They do no all live for a finie period of hea exposure and hen suddenly die a once,
More informationCHAPTER 12 DIRECT CURRENT CIRCUITS
CHAPTER 12 DIRECT CURRENT CIUITS DIRECT CURRENT CIUITS 257 12.1 RESISTORS IN SERIES AND IN PARALLEL When wo resisors are conneced ogeher as shown in Figure 12.1 we said ha hey are conneced in series. As
More information1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.
. Use Simpson s rule wih n 4 o esimae an () +. Soluion: Since we are using 4 seps, 4 Thus we have [ ( ) f() + 4f + f() + 4f 3 [ + 4 4 6 5 + + 4 4 3 + ] 5 [ + 6 6 5 + + 6 3 + ]. 5. Our funcion is f() +.
More informationSection 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 informationVehicle 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 informationPHYSICS 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 informationMath 2214 Solution Test 1B Fall 2017
Mah 14 Soluion Tes 1B Fall 017 Problem 1: A ank has a capaci for 500 gallons and conains 0 gallons of waer wih lbs of sal iniiall. A soluion conaining of 8 lbsgal of sal is pumped ino he ank a 10 galsmin.
More informationLecture 4 Kinetics of a particle Part 3: Impulse and Momentum
MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an
More informationAP CALCULUS AB 2003 SCORING GUIDELINES (Form B)
SCORING GUIDELINES (Form B) Quesion A blood vessel is 6 millimeers (mm) long Disance wih circular cross secions of varying diameer. x (mm) 6 8 4 6 Diameer The able above gives he measuremens of he B(x)
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationChapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis
Chaper EEE83 EEE3 Chaper # EEE83 EEE3 Linear Conroller Design and Sae Space Analysis Ordinary Differenial Equaions.... Inroducion.... Firs Order ODEs... 3. Second Order ODEs... 7 3. General Maerial...
More informationWall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.
MECHANICS APPLICATIONS OF SECOND-ORDER ODES 7 Mechanics applicaions of second-order ODEs Second-order linear ODEs wih consan coefficiens arise in many physical applicaions. One physical sysems whose behaviour
More informationStructural Dynamics and Earthquake Engineering
Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/
More informationIntegration of the equation of motion with respect to time rather than displacement leads to the equations of impulse and momentum.
Inegraion of he equaion of moion wih respec o ime raher han displacemen leads o he equaions of impulse and momenum. These equaions greal faciliae he soluion of man problems in which he applied forces ac
More information!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)
"#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5
More informationChapter 5 Kinematics
Chaper 5 Kinemaics In he firs place, wha do we mean b ime and space? I urns ou ha hese deep philosophical quesions have o be analzed ver carefull in phsics, and his is no eas o do. The heor of relaivi
More informationPROBLEMS FOR MATH 162 If a problem is starred, all subproblems are due. If only subproblems are starred, only those are due. SLOPES OF TANGENT LINES
PROBLEMS FOR MATH 6 If a problem is sarred, all subproblems are due. If onl subproblems are sarred, onl hose are due. 00. Shor answer quesions. SLOPES OF TANGENT LINES (a) A ball is hrown ino he air. Is
More informationAnswers to 1 Homework
Answers o Homework. x + and y x 5 y To eliminae he parameer, solve for x. Subsiue ino y s equaion o ge y x.. x and y, x y x To eliminae he parameer, solve for. Subsiue ino y s equaion o ge x y, x. (Noe:
More informationParametrics and Vectors (BC Only)
Paramerics and Vecors (BC Only) The following relaionships should be learned and memorized. The paricle s posiion vecor is r() x(), y(). The velociy vecor is v(),. The speed is he magniude of he velociy
More informationMath Final Exam Solutions
Mah 246 - Final Exam Soluions Friday, July h, 204 () Find explici soluions and give he inerval of definiion o he following iniial value problems (a) ( + 2 )y + 2y = e, y(0) = 0 Soluion: In normal form,
More informationMath 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 informationStochastic 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 informationLab 10: RC, RL, and RLC Circuits
Lab 10: RC, RL, and RLC Circuis In his experimen, we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors. We will sudy he way volages and currens change in
More informationt + 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 informationNote: For all questions, answer (E) NOTA means none of the above answers is correct.
Thea Logarihms & Eponens 0 ΜΑΘ Naional Convenion Noe: For all quesions, answer means none of he above answers is correc.. The elemen C 4 has a half life of 70 ears. There is grams of C 4 in a paricular
More informationCh.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC
Ch.. Group Work Unis Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Uni 2 Jusify wheher he following saemens are rue or false: a) Two sreamlines, corresponding o a same insan of ime, can never inersec
More informationMidterm Exam Review Questions Free Response Non Calculator
Name: Dae: Block: Miderm Eam Review Quesions Free Response Non Calculaor Direcions: Solve each of he following problems. Choose he BEST answer choice from hose given. A calculaor may no be used. Do no
More informationClass Notes 1: Introduction. MAE 82 Engineering Mathematics
Class Noes 1: Inroducion MAE 82 Engineering Mahemaics CHANGE Rae of Change Basic Mahemaical Models Man of he principles or laws underling he behavior of he naural world are saemens or relaions involving
More informationk 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series
Final Review A Puzzle... Consider wo massless springs wih spring consans k 1 and k and he same equilibrium lengh. 1. If hese springs ac on a mass m in parallel, hey would be equivalen o a single spring
More informationMatlab 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 information2.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 informationMEI 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 informationMath 4600: Homework 11 Solutions
Mah 46: Homework Soluions Gregory Handy [.] One of he well-known phenomenological (capuring he phenomena, bu no necessarily he mechanisms) models of cancer is represened by Gomperz equaion dn d = bn ln(n/k)
More information