Maple Tools for Differential Equations A. J. Meir
|
|
- Gillian Woods
- 5 years ago
- Views:
Transcription
1 Maple Tools for Differenial Equaions A. J. Meir Copyrigh (C) A. J. Meir. All righs reserved. This workshee is for educaional use only. No par of his publicaion may be reproduced or ransmied for profi in any form or by any means, elecronic or mechanical, including phoocopy, recording, or any informaion sorage and rerieval sysem wihou prior wrien permission from he auhor. No for profi disribuion of he sofware is allowed wihou prior wrien permission, providing ha he workshee is no modified in any way and full credi o he auhor is acknowledged. (This documen is largely based on a similar documen originally wrien by P. G. Schmid.) Maple will find symbolic soluions (ha is, closed-form soluions or soluion formulas) for many ypes of d.e.s and i.v.p.s, including separable and linear firs-order equaions. Moreover, Maple can compue and graph numerical soluions (ha is, approximae soluions) for virually all i.v.p.s of pracical ineres. 1. Finding and Ploing Symbolic Soluions. To find a general soluion for a given d.e., use Maple's dsolve command, whose general forma is dsolve(de,unknown_funcion). The following command, for example, finds a general soluion for he d.e. x' = Kk x (he radioacivedecay equaion), where he unknown funcion is x = x and k is a parameer: dsolve(d(x)()=-k*,); x = _C1 e Kk (1) Noe ha inside he dsolve command he derivaive, x', should be wrien as D x and ha he independen variable,, mus be carried along (so ha he d.e. becomes D x = Kk x insead of x' = Kk x ). In Maple's response, _C1 denoes an arbirary consan. (In general, Maple renders arbirary consans as _C1, _C, _C3, and so forh.) In sandard calculus noaion, we would wrie he soluion formula as x = c e Kk. The nex command produces a general soluion for he d.e. x'' + a x = 0 (he elasic-spring equaion wih a = k ), where he unknown funcion is again x = x and a is a parameer: m dsolve(d(d(x))()+a^*,); x = _C1 sin a C_C cos a () Noe ha he second derivaive, x'', should be wrien as D D x and ha he independen variable mus be carried along, jus as in he previous example. Wih a minor modificaion, we can use he dsolve command o find symbolic soluions of i.v.p.s; he modified forma is dsolve({de,iniial_condiion(s)},unknown_funcion). Noe he "curly braces" around he d.e. and iniial condiion(s). For example, le's solve he firs-order i. v.p. x' = Kk x, x 0 = x0 : dsolve({d(x)()=-k*,x(0)=x0},); x = x0 e Kk (3) Nex, we solve he second-order i.v.p. x'' + a x = 0, x 0 = x0, x'(0) = v0 : 1
2 dsolve({d(d(x))()+a^*=0,x(0)=x0,d(x)(0)=v0},); v0 sin a x = Cx0 cos a a () In order o plo a soluion generaed by he dsolve command, we firs have o give i a name, say, s. The following command ses s equal o he righ-hand side of he immediaely preceding soluion formula: s:=rhs(%); s := v0 sin a a Cx0 cos a (5) Noe: The percen (%) is jus a placeholder for he las oupu generaed before execuion of he presen command (which, in his case, happens o be he soluion formula produced by he preceding dsolve command). Insead of using he percen, we could have given a name (say, sol ) o he soluion formula and used ha name in he rhs command: sol:=dsolve({d(d(x))()+a^*=0,x(0)=x0,d(x)(0)=v0},); v0 sin a sol := x = Cx0 cos a a (6) s:=rhs(sol); s := v0 sin a a Cx0 cos a (7) If we now specify he value of he parameer a and he iniial values x0 and v0, we can plo s jus like any oher funcion. Se, for example, a = 1, x0 =, and v0 = K3 : a:=1; x0:=; v0:=-3; a := 1 x0 := v0 := K3 (8) In his case, he soluion of he i.v.p. is s; K3 sin C cos (9) and a graph can be generaed wih he usual plo command: plo(s,=-10..0);
3 3 1 K K1 K K3 To obain he soluion and is graph for oher values of he parameer a or oher iniial values x0 and v0, jus re-execue he las hree commands, changing only he firs one as needed. For furher deails on ploing, especially on simulaneous ploing of several graphs, see he handou "More on Basic Maple Commands".. Ploing Direcion Fields and Numerical Soluions. The basic Maple command for he graphical represenaion of direcion fields and approximae soluion curves is DEplo ; i is par of a package of rouines called DEools, which mus be explicily loaded wih he command wih(deools): To plo a direcion field for a firs-order d.e. of he form x' = f, x, in he recangle = a..b, x = c..d, use he DEplo command in he form DEplo(D(x)()=f(,),,=a..b,x=c..d). The following command, for example, generaes a direcion field for x' = x in he recangle = K3..3, x =..6 : DEplo(D(x)()=^,,=-3..3,x=-6..6); 3
4 6 K3 K K K K The defaul number of grid poins used is 0 in eiher coordinae direcion. You can change ha, if needed, by adding he opion dirgrid=[m,n] inside he DEplo command, where m and n denoe he desired number of grid poins in he horizonal and verical direcion, respecively. Example: DEplo(D(x)()=^,,=-3..3,x=-6..6,dirgrid=[15,30]);
5 6 K3 K K K K In addiion o a direcion field, Maple will compue and display one or more approximae soluion curves of a given d.e. if suiable iniial condiions are specified in he form {[x(0)=x0],[x(1)=x1],[x()=x],...}, immediaely following he range specificaions in he DEplo command. Noe ha he se of iniial condiions is enclosed in curly braces, while each individual condiion appears in brackes. Boh, brackes and braces, mus be presen even if only a single iniial condiion is given. For example, he following command produces a direcion field for x' = x, as before, along wih he (approximae) graph of a paricular soluion, saisfying he iniial condiion x 0 =1 : DEplo(D(x)()=^,,=-3..3,x=-6..6,{[x(0)=1]}); 5
6 6 K3 K K K K Now we will add a second soluion curve, saisfying he iniial condiion x 0 = K1 : DEplo(D(x)()=^,,=-3..3,x=-6..6,{[x(0)=1],[x(0)=-1]}) ; 6 K3 K K K K 6
7 In general, he given iniial poins do no have o lie inside he plo window. If iniial condiions are given, specificaion of he verical range (x-range) is opional; if omied, DEplo will choose he x-range auomaically, depending on where he soluion curves "go". This is someimes useful, bu no always. Le's ry i wih he above example (same command as before, bu wihou x-range specificaion): DEplo(D(x)()=^,,=-3..3,x=-6..6,{[x(0)=1],[x(0)=-1]}) ; 6 K3 K K K K Can you explain he poor oucome? (Hin: The exac soluion wih iniial condiion x 0 = 1 has a verical asympoe a = 1.) If only he soluion curves are of ineres, display of he direcion field arrows can be suppressed by adding he opion arrows=none inside he DEplo command, following he specificaion of iniial condiions. Example: DEplo(D(x)()=^,,=-3..3,x=-6..6,{[x(0)=1],[x(0)=-1]}, arrows=none); 7
8 6 K3 K K K K Many furher opions are available. In fac, mos of he opions available wih he basic plo command will also work wih DEplo. Very useful for prining purposes is he opion linecolor=black, which causes he soluion curves o be drawn in black raher han yellow (he defaul, which looks piiful on a black-and-whie prinou). The color and appearance of he direcion field arrows can be cusomized, oo, bu his will hardly ever be necessary. For deails, consul he relevan help pages (by yping?deplo or?plo,opions). Noe: pions may be added o he DEplo command in any order, bu mus be preceded by he range specificaion(s) and iniial condiion(s) (if any). Example: DEplo(D(x)()=^,,=-3..3,x=-6..6,{[x(0)=1],[x(0)=-1]}, dirgrid=[10,10],linecolor=black,color=green,ile=`dx/d=x^`); 8
9 dx/d=x^ 6 K3 K K K K The defaul mehod ha DEplo uses o compue approximae soluions of he given d.e. is "rk", a socalled h order Runge-Kua mehod, much more sophisicaed han he classical Euler mehod. You can force Maple o use he Euler mehod by specifying he opion mehod=classical. Boh "rk" and "Euler" employ a fixed sepsize; he defaul sepsize is (b Ka)/0, where a..b is he horizonal range (-range) of he plo. Frequenly, i will be necessary o reduce he sepsize, for improved numerical accuracy. This can be done by specifying he opion sepsize=h, where h denoes he desired sepsize, smaller han (b Ka)/0. Le's ry he following example, where we approximae he same soluions of x' = x as before, bu wih a larger verical range (K10..10) and using he classical Euler mehod (insead of "rk"): DEplo(D(x)()=^,,=-3..3,x=-9..9,{[x(0)=1],[x(0)=-1]}, mehod=classical,linecolor=black); 9
10 8 6 K3 K K K K K8 The horizonal par of he lower soluion curve is clearly erroneous -- i cerainly doesn' mach he direcion field! The (defaul) sepsize used in he plo is.3 (namely, (b Ka)/0 wih b = 3 and a = K3 ), and apparenly his is no accurae enough. Reducing he sepsize should remedy he problem: DEplo(D(x)()=^,,=-3..3,x=-9..9,{[x(0)=1],[x(0)=-1]}, mehod=classical,sepsize=.1,linecolor=black); 10
11 8 6 K3 K K K K K8 To accuraely rack a rapidly growing or highly oscillaory soluion, i may be necessary o employ a very small sepsize and compue millions of daa poins. No only will his consume a lo of ime -- soring millions of daa poins may also ax Maple's memory. Forunaely, i is no nessary o sore all he compued daa poins (he number of poins used in creaing he plo is anyway limied by he resoluion of compuer screen and laser priner). Wih he opion ieraions=n, you can insruc Maple o sore and plo only every nh of he compued daa poins. This opion, if applied judiciously, can save a lo of ime and memory, wih no loss of numerical accuracy. As a drasic example, we will compue and display an approximae soluion of he i.v.p. dx =1Ke x, x 0 =0, d on he inerval = 0..15, using "rk" (he defaul mehod) wih a sepsize of.001 and 100 ieraions beween ploed poins. (Thus, Maple will compue 15,000 daa poins, bu sore and plo only 150.) DEplo(D(x)()=1-exp()*^,,=0..15,{[x(0)=0]}, sepsize=.001,ieraions=100,arrows=none,linecolor=black); 11
12 Generaing his plo sill akes quie some ime, bu wihou he opion ieraions=100, i would ake abou en imes as long, wihou any visible improvemen! n he oher hand, if you ried o speed hings up by using a sepsize significanly greaer han.001, he plo would no be reliable. Check i ou! Beyond "rk" (he defaul mehod) and "Euler" ( mehod=classical ), you have a choice of several, much more sophisicaed d.e. solvers. ne of hem is "rkf5", a so-called h/5h order Runge-Kua-Fehlberg mehod, which can be invoked by specifying mehod=rkf5 in he DEplo command. While "rk" and "Euler" use a fixed sepsize, "rkf5" employs a variable sepsize ha's auomaically adjused in accordance wih he "siffness" of he problem. As an adapive (self-correcing) mehod, "rkf5" is more reliable han "rk" or "Euler", bu usually slower. Bu in "siff" problems (like he previous example), "rkf5" offers a double advanage: You don' have o worry abou he sepsize, and i is faser han "rk" or "Euler"! (In general, conrolling and adjusing he sepsize as needed is more efficien han using a very small sepsize hroughou.) Since "rkf5" conrols he (numerical) sepsize auomaically, specifying a sepsize when using his mehod has no effec on he numerical accuracy of he procedure; i does, however, affec he number of poins ploed and, hereby, he apparen smoohness of he graph. More precisely, in conjuncion wih "rkf5", he specified sepsize becomes simply he horizonal disance beween consecuive poins ploed. For example, if we solve he previous problem wih "rkf5", we should specify a sepsize of.1 o obain he same number of poins ploed as before, when we used "rk" wih a (numerical) sepsize of.001 and 100 ieraions beween poins ploed. Le's ry i: DEplo(D(x)()=1-exp()*^,,=0..15,{[x(0)=0]}, mehod=rkf5,sepsize=.1,arrows=none,linecolor=black); 1
13 DEplo will no only solve firs-order, bu also higher-order i.v.p.s (alhough no direcion fields can be generaed for he laer). The following command, for example, produces an approximae soluion of he i.v.p. x'' + x = 0, x 0 =, x'(0) = -3, whose symbolic soluion is cos K3 sin : DEplo(D(D(x))()+=0,,=-10..0,{[x(0)=,D(x)(0)=-3]}, sepsize=.,linecolor=black); 13
14 3 1 K K1 K K3 Noe ha he wo iniial condiions appear inside a single pair of brackes (because hey deermine only one soluion curve). The sepsize specificaion was necessary, due o he oscillaory naure of he graph and he lengh of he plo inerval; he defaul sepsize would have been 30/0 = 1.5! (Check ou wha you would ge wihou sepsize specificaion.) Basically, you have he same opions as in he firs-order case (excep for hose relaing o direcion fields, of course), and jus as in he firs-order case, several approximae soluions can be generaed simulaneously, by specifying several ses of iniial daa in he form {[x(0)=x0,d(x)(0)=v0],[x(1)=x1,d(x)(1)=v1],[x(3)=x3,d(x)(3)=v3],...}. Noe ha, in analogy o he firs-order case, he enire se of iniial condiions is enclosed in curly braces, while he iniial condiions for each individual soluion curve appear in brackes. Example: DEplo(D(D(x))()+=0,,=-10..0,{[x(0)=,D(x)(0)=-3],[x (0)=-3,D(x)(0)=]},linecolor=black); 1
15 3 1 K K1 K K3 Prey lousy! bviously, we should reduce he sepsize or use a more flexible solver. Le's do boh and also add some perks: DEplo(D(D(x))()+=0,,=-10..0,{[x(0)=,D(x)(0)=-3],[x (0)=-3,D(x)(0)=]}, mehod=rkf5,sepsize=.,linecolor=[blue,green], ile=`x(0)=, x'(0)=-3 (blue); x(0)=-3, x'(0)= (green)`, ilefon=[times,rman,1]); 15
16 x(0)=, x'(0)=-3 (blue); x(0)=-3, x'(0)= (green) 3 1 K K1 K K3 16
Math Linear Differential Equations
Mah 65 - Linear Differenial Equaions A. J. Meir Copyrigh (C) A. J. Meir. All righs reserved. This workshee is for educaional use only. No par of his publicaion may be reproduced or ransmied for profi in
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 informationdy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page
Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,
More 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 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 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 informationChapter 7: Solving Trig Equations
Haberman MTH Secion I: The Trigonomeric Funcions Chaper 7: Solving Trig Equaions Le s sar by solving a couple of equaions ha involve he sine funcion EXAMPLE a: Solve he equaion sin( ) The inverse funcions
More information!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)
"#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5
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 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 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 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 informationSingle and Double Pendulum Models
Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double
More informationKEY. 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 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 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 informationThe equation to any straight line can be expressed in the form:
Sring Graphs Par 1 Answers 1 TI-Nspire Invesigaion Suden min Aims Deermine a series of equaions of sraigh lines o form a paern similar o ha formed by he cables on he Jerusalem Chords Bridge. Deermine he
More informationMATH 128A, SUMMER 2009, FINAL EXAM SOLUTION
MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange
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 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 informationKEY. 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 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 informationEE 301 Lab 2 Convolution
EE 301 Lab 2 Convoluion 1 Inroducion In his lab we will gain some more experience wih he convoluion inegral and creae a scrip ha shows he graphical mehod of convoluion. 2 Wha you will learn This lab will
More information2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes
Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion
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 informationSOLUTIONS 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 informationIn this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should
Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion
More information2. 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 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 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 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 informationODEs 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 informationSMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.
SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a
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 informationPredator - 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 informationDifferential 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 informationHamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t
M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n
More informationChapter 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 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 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 informationMath 23 Spring Differential Equations. Final Exam Due Date: Tuesday, June 6, 5pm
Mah Spring 6 Differenial Equaions Final Exam Due Dae: Tuesday, June 6, 5pm Your name (please prin): Insrucions: This is an open book, open noes exam. You are free o use a calculaor or compuer o check your
More informationMathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation
Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion
More informationChallenge 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 informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More informationProblem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims
Problem Se 5 Graduae Macro II, Spring 2017 The Universiy of Nore Dame Professor Sims Insrucions: You may consul wih oher members of he class, bu please make sure o urn in your own work. Where applicable,
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More informationSolutions 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 informationBias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé
Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070
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 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 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 informationIntroduction to Probability and Statistics Slides 4 Chapter 4
Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random
More informationThis document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC
his documen was generaed a 1:4 PM, 9/1/13 Copyrigh 213 Richard. Woodward 4. End poins and ransversaliy condiions AGEC 637-213 F z d Recall from Lecure 3 ha a ypical opimal conrol problem is o maimize (,,
More information5.2. The Natural Logarithm. Solution
5.2 The Naural Logarihm The number e is an irraional number, similar in naure o π. Is non-erminaing, non-repeaing value is e 2.718 281 828 59. Like π, e also occurs frequenly in naural phenomena. In fac,
More informationd 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 informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More 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 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 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 informationMATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence
MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,
More 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 informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion
More information4.5 Constant Acceleration
4.5 Consan Acceleraion v() v() = v 0 + a a() a a() = a v 0 Area = a (a) (b) Figure 4.8 Consan acceleraion: (a) velociy, (b) acceleraion When he x -componen of he velociy is a linear funcion (Figure 4.8(a)),
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More 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 informationFishing 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 informationA First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18
A Firs ourse on Kineics and Reacion Engineering lass 19 on Uni 18 Par I - hemical Reacions Par II - hemical Reacion Kineics Where We re Going Par III - hemical Reacion Engineering A. Ideal Reacors B. Perfecly
More informationThis document was generated at 7:34 PM, 07/27/09 Copyright 2009 Richard T. Woodward
his documen was generaed a 7:34 PM, 07/27/09 Copyrigh 2009 Richard. Woodward 15. Bang-bang and mos rapid approach problems AGEC 637 - Summer 2009 here are some problems for which he opimal pah does no
More informationLab #2: Kinematics in 1-Dimension
Reading Assignmen: Chaper 2, Secions 2-1 hrough 2-8 Lab #2: Kinemaics in 1-Dimension Inroducion: The sudy of moion is broken ino wo main areas of sudy kinemaics and dynamics. Kinemaics is he descripion
More 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 informationOnline 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( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
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 informationR t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t
Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,
More informationThe 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 informationChapter Floating Point Representation
Chaper 01.05 Floaing Poin Represenaion Afer reading his chaper, you should be able o: 1. conver a base- number o a binary floaing poin represenaion,. conver a binary floaing poin number o is equivalen
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 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 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 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 information= ( ) ) or a system of differential equations with continuous parametrization (T = R
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
More informationDecimal moved after first digit = 4.6 x Decimal moves five places left SCIENTIFIC > POSITIONAL. a) g) 5.31 x b) 0.
PHYSICS 20 UNIT 1 SCIENCE MATH WORKSHEET NAME: A. Sandard Noaion Very large and very small numbers are easily wrien using scienific (or sandard) noaion, raher han decimal (or posiional) noaion. Sandard
More information5.1 - Logarithms and Their Properties
Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We
More informationY 0.4Y 0.45Y Y to a proper ARMA specification.
HG Jan 04 ECON 50 Exercises II - 0 Feb 04 (wih answers Exercise. Read secion 8 in lecure noes 3 (LN3 on he common facor problem in ARMA-processes. Consider he following process Y 0.4Y 0.45Y 0.5 ( where
More informationEchocardiography Project and Finite Fourier Series
Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every
More informationConstant Acceleration
Objecive Consan Acceleraion To deermine he acceleraion of objecs moving along a sraigh line wih consan acceleraion. Inroducion The posiion y of a paricle moving along a sraigh line wih a consan acceleraion
More informationGuest Lectures for Dr. MacFarlane s EE3350 Part Deux
Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A
More informationt 2 B F x,t n dsdt t u x,t dxdt
Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.
More informationModule 2 F c i k c s la l w a s o s f dif di fusi s o i n
Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms
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 informationCircuit Variables. AP 1.1 Use a product of ratios to convert two-thirds the speed of light from meters per second to miles per second: 1 ft 12 in
Circui Variables 1 Assessmen Problems AP 1.1 Use a produc of raios o conver wo-hirds he speed of ligh from meers per second o miles per second: ( ) 2 3 1 8 m 3 1 s 1 cm 1 m 1 in 2.54 cm 1 f 12 in 1 mile
More informationModule 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 informationSummary of shear rate kinematics (part 1)
InroToMaFuncions.pdf 4 CM465 To proceed o beer-designed consiuive equaions, we need o know more abou maerial behavior, i.e. we need more maerial funcions o predic, and we need measuremens of hese maerial
More information1 Differential Equation Investigations using Customizable
Differenial Equaion Invesigaions using Cusomizable Mahles Rober Decker The Universiy of Harford Absrac. The auhor has developed some plaform independen, freely available, ineracive programs (mahles) for
More information1. An introduction to dynamic optimization -- Optimal Control and Dynamic Programming AGEC
This documen was generaed a :45 PM 8/8/04 Copyrigh 04 Richard T. Woodward. An inroducion o dynamic opimizaion -- Opimal Conrol and Dynamic Programming AGEC 637-04 I. Overview of opimizaion Opimizaion is
More informationDIFFERENTIAL GEOMETRY HW 5
DIFFERENTIAL GEOMETRY HW 5 CLAY SHONKWILER 3. Le M be a complee Riemannian manifold wih non-posiive secional curvaure. Prove ha d exp p v w w, for all p M, all v T p M and all w T v T p M. Proof. Le γ
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 informationarxiv: v1 [math.ca] 15 Nov 2016
arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy
More information( ) ( ) 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