MATH 3795 Lecture 18. Numerical Solution of Ordinary Differential Equations.
|
|
- Chad Kennedy
- 5 years ago
- Views:
Transcription
1 MATH 3795 Lecture 18. Numericl Solution of Ordinry Differentil Equtions. Dmitriy Leykekhmn Fll 2008 Gols Introduce ordinry differentil equtions (ODEs) nd initil vlue problems (IVPs). Exmples of IVPs. Existence nd uniqueness of solutions of IVPs. Dependence of the solution of n IVPs on prmeters. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 1
2 Systems of Ordinry Differentil Equtions. Given functions f 1,..., f n : R n+1 R nd sclrs y 0,1,..., y 0,n R, we wnt to find y 1, y 2,..., y n : R R such tht nd y 1(x) = f 1 (x, y 1 (x),..., y n (x)),. y n(x) = f n (x, y 1 (x),..., y n (x)) y 1 () = y 0,1,. y n () = y 0,n. (1) is clled system of first order ordinry differentil equtions (ODEs) nd (2) re clled initil conditions. Together it is clled n initil vlue problem (IVP). (1) (2) D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 2
3 Systems of Ordinry Differentil Equtions. If we define f : R n+1 R n, f(x, y 1,..., y n ) = we cn rewrite the system s where y 0 = (y 0,1,..., y 0,n ) T. y : R R n, y(x) = y (x) = f(x, y(x)) y() = y 0, f 1 (x, y 1..., y n ). f n (x, y 1..., y n ) y 1 (x). y n (x),, D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 3
4 n-th Order Differentil Eqution. Often, one hs to solve n-th order differentil equtions of the form z (n) (x) = g(x, z(x), z (x),..., z (n 1) (x)), x [, b] (3) with initil conditions z() = z 0, z () = z 1,..., z (n 1) () = z n 1. (3) cn be reformulted s systems of first order ODEs. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 4
5 n-th Order Differentil Eqution. If we introduce the functions y 1 (x) = z(x), y 2 (x) = z (x),. y n (x) = z (n 1) (x), then these functions stisfy the first order differentil equtions y 1(x) = y 2 (x), y 2(x) = y 3 (x),. y n 1(x) = y n (x), y n(x) = g(x, y 1 (x), y 2 (x),..., y n (x)), with initil conditions y 1 () = z 0, y 2 () = z 1,..., y n () = z n 1. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 5
6 n-th Order Differentil Eqution. If we introduce nd f(x, y 1,..., y n ) = y 0 = y 2 (x) y 3 (x). y n (x) g(x, y 1 (x),..., y n (x)) z 0 z 1. z n 1, then we rrive t the system of first order ODEs y (x) = f(x, y(x)), y() = y 0. Thus, it is sufficient to consider systems of first order ODEs. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 6
7 Exmple 1 (Predtor-Prey Model). Let s look t specil cse of n interction between two species, one of which the predtors ets the other the prey. Cndin lynx nd snowshoe hre Dt from pelt-trding records Pictures nd dt re from D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 7
8 Exmple 1 (Predtor-Prey Model). Mthemticl model. Assumptions. The predtor species is totlly dependent on the prey species s its only food supply. y(t) is the size of the predtor popultion t time t. The prey species hs n unlimited food supply; only predtor poses thret to its growth. x(t) is the size of the prey popultion t t. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 8
9 Exmple 1 (Predtor-Prey Model). Prey popultion. If there were no predtors, the second ssumption would imply tht the prey species grows exponentilly, i.e., x (t) = x(t). Since there re predtors, we must ccount for negtive component in the prey growth rte. Assumptions: The rte t which predtors encounter prey is jointly proportionl to the sizes of the two popultions. A fixed proportion of encounters leds to the deth of the prey. Predtor popultion. x (t) = x(t) bx(t)y(t). y (t) = cy(t) + px(t)y(t). D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 9
10 Exmple 1 (Predtor-Prey Model). Prey popultion. If there were no predtors, the second ssumption would imply tht the prey species grows exponentilly, i.e., x (t) = x(t). Since there re predtors, we must ccount for negtive component in the prey growth rte. Assumptions: The rte t which predtors encounter prey is jointly proportionl to the sizes of the two popultions. A fixed proportion of encounters leds to the deth of the prey. Predtor popultion. x (t) = x(t) bx(t)y(t). y (t) = cy(t) + px(t)y(t). Lotk-Volterr Predtor-Prey Model: (, b, c, p > 0 re constnts) x (t) = x(t) bx(t)y(t), y (t) = cy(t) + px(t)y(t). D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 9
11 Exmple 1 (Predtor-Prey Model). Solution of the Lotk-Volterr Predtor-Prey Model x (t) = x(t) bx(t)y(t), y (t) = cy(t) + px(t)y(t) with initil condition x(0) = 20, y(0) = 40 nd prmeters = 1, b = 0.02, c = 0.3, p = D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 10
12 Exmple 2 (Chemicl Rections). A rection involving the compounds A, B, C, D is written s σ A A + σ B B σ C C + σ D D. Here σ A, σ B, σ C, σ D re the stoichiometric coefficients. The compounds A, B re the rectnts, C, D re the products. The indictes tht the rection is irreversible. For reversible rection we use. For exmple, the reversible rection of crbon dioxide nd hydrogen to form methne plus wter is The stoichiometric coefficients re CO 2 + 4H 2 CH 4 + 2H 2 O. σ CO2 = 1, σ H2 = 4, σ CH4 = 1, σ H2O = 2. For ech rection we hve rte r (the number of rective events per second per unit volume, mesured in [mol/(sec L)]) of the rection tht together with the stoichiometric coefficients determines the chnge in concentrtions (mesured in [mol/l]) resulting from the rection. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 11
13 Exmple 2 (Chemicl Rections). For exmple, if the rte of the rection is r nd if denote the concentrtion of compound A,... by C A,..., we hve the following chnges in concentrtions: d dt C A(t) = σ A r..., d dt C B(t) = σ B r..., d dt C C(t) =... σ C r..., d dt C D(t) =... σ D r..., Usully the rection r is of the form r = kc α AC β B, where k is the rection rte constnt nd α, β re nonnegtive prmeters. The sum α + β is clled the order of the rection. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 12
14 Exmple 2 (Chemicl Rections. Autoctlytic rection.) Autoctlysis is term commonly used to describe the experimentlly observble phenomenon of homogeneous chemicl rection which shows mrked increse in rte in time, reches its pek t bout 50 percent conversion, nd the drops off. The temperture hs to remin constnt nd ll ingredients must be mixed t the strt for proper observtion. We consider the ctlytic therml decomposition of single compound A into two products B nd C, of which B is the utoctlytic gent. A cn decompose vi two routes, slow unctlyzed one (r 1 ) nd nother ctlyzed by B (r 3 ). The three essentil kinetic steps re A B + C A + B AB AB 2B + C Strt or bckground rection, Complex formtion, Autoctlytic step. Agin, we denote the concentrtion of compound A,... by C A,.... The rection rtes for the three rections re r 1 = k 1 C A, r 2 = k 2 C A C B, r 3 = k 3 C AB. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 13
15 Exmple 2 (Chemicl Rections.) The utoctlytic rection leds to system of ODEs dc A dt dc B dt dc AB dt dc C dt = k 1 C A k 2 C A C B, = k 1 C A k 2 C A C B + 2k 3 C AB, = k 2 C A C B k 3 C AB, = k 1 C A + k 3 C AB with k 1 = , k 2 = 1, k 3 = nd initil vlues C A (0) = 1, C B (0) = 0, C AB (0) = 0, C C (0) = 0. Here time is mesured in [sec] nd concentrtions re mesured in [kmol/l]. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 14
16 Exmple 3. For given A R n n consider the liner IVP y (x) = Ay(x), x [0, ) y(0) = y 0. (4) Suppose there exist V = (v 1... v n ) C n n nd Λ = dig(λ 1,..., λ n ), λ j C, such tht (eigen decomposition) or, equivlently, Insert into (4) A = V ΛV 1 Av j = λ j v j, j = 1,..., n. (V 1 y) (x) = V 1 y (x) = V 1 AV V 1 y(x) = ΛV 1 y(x), x [0, ) V 1 y(0) = V 1 y 0. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 15
17 Exmple 3. Putting z := V 1 y we cn see tht z stisfies z (x) = Λz(x), x [0, ), z(0) = V 1 y 0 := z 0. Since Λ = dig(λ 1,..., λ n ), this is equivlent to z j(x) = λ j z j (x), x [0, ), z j (0) = z 0,j, j = 1,..., n. Unique solution z j (x) = e λjx z 0,j, j = 1,..., n. Hence, the unique solution of (4) is given by n y(x) = V z(x) = v j e λjx z 0,j. j=1 D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 16
18 Exmple 3. Consider We cn show A = 1 ( A = ( ) ( ) ) 1 5 ( ) = V ΛV 1. nd y 0 = (1, 1) T. Since V 1 y 0 = 1 5 ( ) ( 1 1 ) = 1 5 ( 3 1 ) := z 0 D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 17
19 Exmple 3. This leds to z 1(x) = 100z 1 (x), z 1 (0) = 3 5 z 2(x) = 0.1z 2 (x), z 2 (0) = 1 5 D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 18
20 Existence nd Uniqueness of Solutions. Theorem (Existence Theorem of Peno) Let D R n+1 be domin, i.e. n open connected subset of R n with (, y 0 ) D. If f is continuous on D, then there exists δ > 0 such tht the IVP y (x) = f(x, y(x)), y() = y 0. hs solution on the intervl [ δ, + δ]. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 19
21 Existence nd Uniqueness of Solutions. Exmple Consider the IVP y (x) = y(x) 1/3, y(0) = 0, x 0. For rbitrry x > 0 the functions { 0, 0 x x y(x) = ± 2 3 [(x x)]3/2, x x. solve the IVP. Hence, the IVP hs infinitely mny solutions. Note tht the prtil derivtive of f(x, y) = y 1/3 with respect to y is singulr t y = 0. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 20
22 Existence nd Uniqueness of Solutions. Theorem (Existence nd Uniqueness Theorem of PicrdLindelöf) If f is continuous on D nd if there exists M > 0 such tht f(x, y) M for ll (x, y) D nd if f is Lipschitz continuous with respect to y on D = {(x, y) D : x δ, y y 0 δm}, i.e., if there exists L > 0 such tht f(x, y 1 ) f(x, y 2 ) L y 1 y 2 for ll (x, y 1 ), (x, y 2 ) D, then y (x) = f(x, y(x)), y() = y 0. hs unique solution on the intervl [ δ, + δ]. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 21
23 Existence nd Uniqueness of Solutions. The Picrd Lindelöf Theorem is bsed on the equivlence of the IVP y (x) = f(x, y(x)), y() = y 0. x [, b], nd the integrl eqution y(x) = y 0 + f(s, y(s))ds. The integrl of vector vlued function is defined component wise, i.e., f 1(s, y(s))ds f(s, y(s))ds =.. f n(s, y(s))ds D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 22
24 Dependence of the Solution of the ODE on Perturbtions of the Problem. Consider the two ODEs y (x) = f(x, y(x)), y() = y 0 (5) nd z (x) = g(x, z(x)), z() = z 0, (6) where f, g : R n+1 R n re given functions nd y 0, z 0 R n re given vectors. We view (6) s perturbtion of (5) nd we wnt to know wht the error between the exct solution y nd the solution of the perturbed problem z is. This question is lso interesting in the context of numericl solutions of IVPs D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 23
25 Dependence of the Solution of the ODE on Perturbtions of the Problem. Suppose there exist constnts ɛ 1, ɛ 2 > 0 such tht nd y 0 z 0 ɛ 1 f(x, y) g(x, y) ɛ 2 x R, y R n. Furthermore, suppose there exists L > 0 such tht f(x, y) f(x, z) L y z x R, y, z R n. Suppose tht the IVPs (5) nd (6) hve unique solutions y nd z. These solutions stisfy the integrl equtions nd y(x) = y 0 + z(x) = z 0 + f(s, y(s))ds (7) g(s, z(s))ds. (8) D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 24
26 Dependence of the Solution of the ODE on Perturbtions of the Problem. Subtrct (8) from (7) to get y(x) z(x) = y 0 z 0 + = y 0 z f(s, y(s)) g(s, z(s))ds f(s, y(s)) f(s, z(s))ds f(s, z(s)) g(s, z(s))ds. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 25
27 Dependence of the Solution of the ODE on Perturbtions of the Problem. Tking the norm we get y(x) z(x) y 0 z y 0 z 0 + }{{} ɛ 1 + f(s, y(s)) f(s, z(s))ds f(s, z(s)) g(s, z(s))ds f(s, y(s)) f(s, z(s)) ds }{{} L y(s) z(s) f(s, z(s)) g(s, z(s)) ds. }{{} ɛ 2 D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 26
28 Dependence of the Solution of the ODE on Perturbtions of the Problem. Hence if we define the error e(x) = y(x) z(x), then e(x) ɛ 1 + L e(s)ds + ɛ 2 (x ). D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 27
29 Dependence of the Solution of the ODE on Perturbtions of the Problem. Lemm (Gronwll Lemm) If h,w nd k re nonnegtive nd continuous on the intervl [, b] stisfying the inequlity h(x) w(x) + then h obeys the estimte k(s)h(s)ds x [, b], h(x) w(x) + e t k(s)ds k(t)w(t)dt. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 28
30 Dependence of the Solution of the ODE on Perturbtions of the Problem. Apply the Gronwlls lemm to the eqution e(x) ɛ 1 + ɛ 2 (x ) + L e(s)ds. Here h(x) = e(x) = y(x) z(x), w(x) = ɛ 1 + ɛ 2 (x ), nd k(x) = L. Thus, the error between the solution y of the originl IVP (5) nd the solution z of the perturbed IVP (6) obeys y(x) z(x) ɛ 1 + ɛ 2 (x ) + e L(x t) L[ɛ 1 + ɛ 2 (t )]dt = ɛ 1 e L(x ) + ɛ 2 L (el(x ) 1) x. (9) D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 29
31 Dependence of the Solution of the ODE on Perturbtions of the Problem. Exmple The solution of the differentil eqution y (x) = 3y(x), y(0) = 1 is given by y(x) = e 3x. The function f(x, y) = 3y is Lipschitz continuous with respect to y with Lipschitz constnt L = 3. The perturbed differentil eqution z (x) = 3z(x) + ɛ 2, z(0) = 1 + ɛ 1 hs the solution z(x) = (1 + ɛ 1 + ɛ2 3 )e3x ɛ2 3. Thus, y(x) z(x) = ɛ 1 e 3x + ɛ 2 3 (e3x 1) x [0, ) This shows tht the error estimte (9) is shrp. D. Leykekhmn - MATH 3795 Introduction to Computtionl Mthemtics Liner Lest Squres 30
1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation
1 1.1. Liner Constnt Coefficient Equtions Section Objective(s): Overview of Differentil Equtions. Liner Differentil Equtions. Solving Liner Differentil Equtions. The Initil Vlue Problem. 1.1.1. Overview
More informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
More informationCHEMICAL KINETICS
CHEMICAL KINETICS Long Answer Questions: 1. Explin the following terms with suitble exmples ) Averge rte of Rection b) Slow nd Fst Rections c) Order of Rection d) Moleculrity of Rection e) Activtion Energy
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationDefinite integral. Mathematics FRDIS MENDELU
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the
More informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More information(4.1) D r v(t) ω(t, v(t))
1.4. Differentil inequlities. Let D r denote the right hnd derivtive of function. If ω(t, u) is sclr function of the sclrs t, u in some open connected set Ω, we sy tht function v(t), t < b, is solution
More informationConservation Law. Chapter Goal. 5.2 Theory
Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very
More informationDefinite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová (Mendel University) Definite integrl MENDELU / Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function
More informationSturm-Liouville Eigenvalue problem: Let p(x) > 0, q(x) 0, r(x) 0 in I = (a, b). Here we assume b > a. Let X C 2 1
Ch.4. INTEGRAL EQUATIONS AND GREEN S FUNCTIONS Ronld B Guenther nd John W Lee, Prtil Differentil Equtions of Mthemticl Physics nd Integrl Equtions. Hildebrnd, Methods of Applied Mthemtics, second edition
More informationThe area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O
1 Section 5. The Definite Integrl Suppose tht function f is continuous nd positive over n intervl [, ]. y = f(x) x The re under the grph of f nd ove the x-xis etween nd is denoted y f(x) dx nd clled the
More informationOrdinary differential equations
Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview Introduction-Modelling Bsic concepts to understnd n ODE. Description nd properties
More informationAMATH 731: Applied Functional Analysis Fall Additional notes on Fréchet derivatives
AMATH 731: Applied Functionl Anlysis Fll 214 Additionl notes on Fréchet derivtives (To ccompny Section 3.1 of the AMATH 731 Course Notes) Let X,Y be normed liner spces. The Fréchet derivtive of n opertor
More informationLecture 1. Functional series. Pointwise and uniform convergence.
1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationBest Approximation in the 2-norm
Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion
More information1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.
1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More informationMath 31S. Rumbos Fall Solutions to Assignment #16
Mth 31S. Rumbos Fll 2016 1 Solutions to Assignment #16 1. Logistic Growth 1. Suppose tht the growth of certin niml popultion is governed by the differentil eqution 1000 dn N dt = 100 N, (1) where N(t)
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationThe Wave Equation I. MA 436 Kurt Bryan
1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string
More information1 1D heat and wave equations on a finite interval
1 1D het nd wve equtions on finite intervl In this section we consider generl method of seprtion of vribles nd its pplictions to solving het eqution nd wve eqution on finite intervl ( 1, 2. Since by trnsltion
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationg i fφdx dx = x i i=1 is a Hilbert space. We shall, henceforth, abuse notation and write g i f(x) = f
1. Appliction of functionl nlysis to PEs 1.1. Introduction. In this section we give little introduction to prtil differentil equtions. In prticulr we consider the problem u(x) = f(x) x, u(x) = x (1) where
More informationReview on Integration (Secs ) Review: Sec Origins of Calculus. Riemann Sums. New functions from old ones.
Mth 20B Integrl Clculus Lecture Review on Integrtion (Secs. 5. - 5.3) Remrks on the course. Slide Review: Sec. 5.-5.3 Origins of Clculus. Riemnn Sums. New functions from old ones. A mthemticl description
More informationContinuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht
More informationJim Lambers MAT 169 Fall Semester Lecture 4 Notes
Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of
More informationEntrance Exam, Real Analysis September 1, 2009 Solve exactly 6 out of the 8 problems. Compute the following and justify your computation: lim
1. Let n be positive integers. ntrnce xm, Rel Anlysis September 1, 29 Solve exctly 6 out of the 8 problems. Sketch the grph of the function f(x): f(x) = lim e x2n. Compute the following nd justify your
More informationLecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the
More informationA PREY-PREDATOR MODEL WITH COVER FOR THE PREY AND AN ALTERNATIVE FOOD FOR THE PREDATOR AND CONSTANT HARVESTING OF BOTH THE SPECIES *
Jordn Journl of Mthemtics nd Sttistics (JJMS) (), 009, pp. 43-54 A PREY-PREATOR MOEL WITH COVER FOR THE PREY A A ALTERATIVE FOO FOR THE PREATOR A COSTAT HARVESTIG OF BOTH THE SPECIES * K. LAKSHMI ARAYA.PATTABHIRAMACHARYULU
More informationAMATH 731: Applied Functional Analysis Fall Some basics of integral equations
AMATH 731: Applied Functionl Anlysis Fll 2009 1 Introduction Some bsics of integrl equtions An integrl eqution is n eqution in which the unknown function u(t) ppers under n integrl sign, e.g., K(t, s)u(s)
More informationSTEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.
STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA STEPHEN SCHECTER. The unit step function nd piecewise continuous functions The Heviside unit step function u(t) is given by if t
More informationn f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1
The Fundmentl Theorem of Clculus As we continue to study the re problem, let s think bck to wht we know bout computing res of regions enclosed by curves. If we wnt to find the re of the region below the
More informationIntroduction to the Calculus of Variations
Introduction to the Clculus of Vritions Jim Fischer Mrch 20, 1999 Abstrct This is self-contined pper which introduces fundmentl problem in the clculus of vritions, the problem of finding extreme vlues
More information1 The Lagrange interpolation formula
Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x
More informationUNIVERSITY OF MALTA DEPARTMENT OF CHEMISTRY. CH237 - Chemical Thermodynamics and Kinetics. Tutorial Sheet VIII
UNIVERSITY OF MALTA DEPARTMENT OF CHEMISTRY CH237 - Chemicl Thermodynmics nd Kinetics Tutoril Sheet VIII 1 () (i) The rte of the rection A + 2B 3C + D ws reported s 1.0 mol L -1 s -1. Stte the rtes of
More informationDefinition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim
Mth 9 Course Summry/Study Guide Fll, 2005 [1] Limits Definition of Limit: We sy tht L is the limit of f(x) s x pproches if f(x) gets closer nd closer to L s x gets closer nd closer to. We write lim f(x)
More information4.1. Probability Density Functions
STT 1 4.1-4. 4.1. Proility Density Functions Ojectives. Continuous rndom vrile - vers - discrete rndom vrile. Proility density function. Uniform distriution nd its properties. Expected vlue nd vrince of
More informationPDE Notes. Paul Carnig. January ODE s vs PDE s 1
PDE Notes Pul Crnig Jnury 2014 Contents 1 ODE s vs PDE s 1 2 Section 1.2 Het diffusion Eqution 1 2.1 Fourier s w of Het Conduction............................. 2 2.2 Energy Conservtion.....................................
More informationdifferent methods (left endpoint, right endpoint, midpoint, trapezoid, Simpson s).
Mth 1A with Professor Stnkov Worksheet, Discussion #41; Wednesdy, 12/6/217 GSI nme: Roy Zho Problems 1. Write the integrl 3 dx s limit of Riemnn sums. Write it using 2 intervls using the 1 x different
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationChapter 5 : Continuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 216 Néhémy Lim Chpter 5 : Continuous Rndom Vribles Nottions. N {, 1, 2,...}, set of nturl numbers (i.e. ll nonnegtive integers); N {1, 2,...}, set of ll
More informationReview of basic calculus
Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below
More informationOrthogonal Polynomials
Mth 4401 Gussin Qudrture Pge 1 Orthogonl Polynomils Orthogonl polynomils rise from series solutions to differentil equtions, lthough they cn be rrived t in vriety of different mnners. Orthogonl polynomils
More informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
More informationModule 2: Rate Law & Stoichiomtery (Chapter 3, Fogler)
CHE 309: Chemicl Rection Engineering Lecture-8 Module 2: Rte Lw & Stoichiomtery (Chpter 3, Fogler) Topics to be covered in tody s lecture Thermodynmics nd Kinetics Rection rtes for reversible rections
More informationA. Limits - L Hopital s Rule. x c. x c. f x. g x. x c 0 6 = 1 6. D. -1 E. nonexistent. ln ( x 1 ) 1 x 2 1. ( x 2 1) 2. 2x x 1.
A. Limits - L Hopitl s Rule Wht you re finding: L Hopitl s Rule is used to find limits of the form f ( ) lim where lim f or lim f limg. c g = c limg( ) = c = c = c How to find it: Try nd find limits by
More informationMath 1431 Section M TH 4:00 PM 6:00 PM Susan Wheeler Office Hours: Wed 6:00 7:00 PM Online ***NOTE LABS ARE MON AND WED
Mth 43 Section 4839 M TH 4: PM 6: PM Susn Wheeler swheeler@mth.uh.edu Office Hours: Wed 6: 7: PM Online ***NOTE LABS ARE MON AND WED t :3 PM to 3: pm ONLINE Approimting the re under curve given the type
More informationA. Limits - L Hopital s Rule ( ) How to find it: Try and find limits by traditional methods (plugging in). If you get 0 0 or!!, apply C.! 1 6 C.
A. Limits - L Hopitl s Rule Wht you re finding: L Hopitl s Rule is used to find limits of the form f ( x) lim where lim f x x! c g x ( ) = or lim f ( x) = limg( x) = ". ( ) x! c limg( x) = 0 x! c x! c
More informationW. We shall do so one by one, starting with I 1, and we shall do it greedily, trying
Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)
More information1.9 C 2 inner variations
46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for
More informationMath 270A: Numerical Linear Algebra
Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner
More informationLinear Inequalities. Work Sheet 1
Work Sheet 1 Liner Inequlities Rent--Hep, cr rentl compny,chrges $ 15 per week plus $ 0.0 per mile to rent one of their crs. Suppose you re limited y how much money you cn spend for the week : You cn spend
More informationSection 5.4 Fundamental Theorem of Calculus 2 Lectures. Dr. Abdulla Eid. College of Science. MATHS 101: Calculus 1
Section 5.4 Fundmentl Theorem of Clculus 2 Lectures College of Science MATHS : Clculus (University of Bhrin) Integrls / 24 Definite Integrl Recll: The integrl is used to find re under the curve over n
More informationWe divide the interval [a, b] into subintervals of equal length x = b a n
Arc Length Given curve C defined by function f(x), we wnt to find the length of this curve between nd b. We do this by using process similr to wht we did in defining the Riemnn Sum of definite integrl:
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More information( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that
Arc Length of Curves in Three Dimensionl Spce If the vector function r(t) f(t) i + g(t) j + h(t) k trces out the curve C s t vries, we cn mesure distnces long C using formul nerly identicl to one tht we
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationAdvanced Calculus: MATH 410 Uniform Convergence of Functions Professor David Levermore 11 December 2015
Advnced Clculus: MATH 410 Uniform Convergence of Functions Professor Dvid Levermore 11 December 2015 12. Sequences of Functions We now explore two notions of wht it mens for sequence of functions {f n
More informationPolynomials and Division Theory
Higher Checklist (Unit ) Higher Checklist (Unit ) Polynomils nd Division Theory Skill Achieved? Know tht polynomil (expression) is of the form: n x + n x n + n x n + + n x + x + 0 where the i R re the
More informationUnit #9 : Definite Integral Properties; Fundamental Theorem of Calculus
Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl
More informationLecture 3. Limits of Functions and Continuity
Lecture 3 Limits of Functions nd Continuity Audrey Terrs April 26, 21 1 Limits of Functions Notes I m skipping the lst section of Chpter 6 of Lng; the section bout open nd closed sets We cn probbly live
More informationPartial Derivatives. Limits. For a single variable function f (x), the limit lim
Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles
More information4 The dynamical FRW universe
4 The dynmicl FRW universe 4.1 The Einstein equtions Einstein s equtions G µν = T µν (7) relte the expnsion rte (t) to energy distribution in the universe. On the left hnd side is the Einstein tensor which
More informationMath 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of
More informationMapping the delta function and other Radon measures
Mpping the delt function nd other Rdon mesures Notes for Mth583A, Fll 2008 November 25, 2008 Rdon mesures Consider continuous function f on the rel line with sclr vlues. It is sid to hve bounded support
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationp(x) = 3x 3 + x n 3 k=0 so the right hand side of the equality we have to show is obtained for r = b 0, s = b 1 and 2n 3 b k x k, q 2n 3 (x) =
Norwegin University of Science nd Technology Deprtment of Mthemticl Sciences Pge 1 of 5 Contct during the exm: Elen Celledoni, tlf. 73593541, cell phone 48238584 PLESE NOTE: this solution is for the students
More informationFor a continuous function f : [a; b]! R we wish to define the Riemann integral
Supplementry Notes for MM509 Topology II 2. The Riemnn Integrl Andrew Swnn For continuous function f : [; b]! R we wish to define the Riemnn integrl R b f (x) dx nd estblish some of its properties. This
More informationTopic 1 Notes Jeremy Orloff
Topic 1 Notes Jerem Orloff 1 Introduction to differentil equtions 1.1 Gols 1. Know the definition of differentil eqution. 2. Know our first nd second most importnt equtions nd their solutions. 3. Be ble
More informationCS667 Lecture 6: Monte Carlo Integration 02/10/05
CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number of
More information7.2 The Definite Integral
7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where
More informationMath 120 Answers for Homework 13
Mth 12 Answers for Homework 13 1. In this problem we will use the fct tht if m f(x M on n intervl [, b] (nd if f is integrble on [, b] then (* m(b f dx M(b. ( The function f(x = 1 + x 3 is n incresing
More informationMath 42 Chapter 7 Practice Problems Set B
Mth 42 Chpter 7 Prctice Problems Set B 1. Which of the following functions is solution of the differentil eqution dy dx = 4xy? () y = e 4x (c) y = e 2x2 (e) y = e 2x (g) y = 4e2x2 (b) y = 4x (d) y = 4x
More informationA New Generalization of Lemma Gronwall-Bellman
Applied Mthemticl Sciences, Vol. 6, 212, no. 13, 621-628 A New Generliztion of Lemm Gronwll-Bellmn Younes Lourtssi LA2I, Deprtment of Electricl Engineering, Mohmmdi School Engineering Agdl, Rbt, Morocco
More informationMath Calculus with Analytic Geometry II
orem of definite Mth 5.0 with Anlytic Geometry II Jnury 4, 0 orem of definite If < b then b f (x) dx = ( under f bove x-xis) ( bove f under x-xis) Exmple 8 0 3 9 x dx = π 3 4 = 9π 4 orem of definite Problem
More informationThe Basic Functional 2 1
2 The Bsic Functionl 2 1 Chpter 2: THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 2.1 Introduction..................... 2 3 2.2 The First Vrition.................. 2 3 2.3 The Euler Eqution..................
More informationMath 115 ( ) Yum-Tong Siu 1. Lagrange Multipliers and Variational Problems with Constraints. F (x,y,y )dx
Mth 5 2006-2007) Yum-Tong Siu Lgrnge Multipliers nd Vritionl Problems with Constrints Integrl Constrints. Consider the vritionl problem of finding the extremls for the functionl J[y] = F x,y,y )dx with
More informationInterpreting Integrals and the Fundamental Theorem
Interpreting Integrls nd the Fundmentl Theorem Tody, we go further in interpreting the mening of the definite integrl. Using Units to Aid Interprettion We lredy know tht if f(t) is the rte of chnge of
More informationPresentation Problems 5
Presenttion Problems 5 21-355 A For these problems, ssume ll sets re subsets of R unless otherwise specified. 1. Let P nd Q be prtitions of [, b] such tht P Q. Then U(f, P ) U(f, Q) nd L(f, P ) L(f, Q).
More informationConsequently, the temperature must be the same at each point in the cross section at x. Let:
HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the
More informationOverview of Calculus
Overview of Clculus June 6, 2016 1 Limits Clculus begins with the notion of limit. In symbols, lim f(x) = L x c In wors, however close you emn tht the function f evlute t x, f(x), to be to the limit L
More informationConservation Law. Chapter Goal. 6.2 Theory
Chpter 6 Conservtion Lw 6.1 Gol Our long term gol is to unerstn how mthemticl moels re erive. Here, we will stuy how certin quntity chnges with time in given region (sptil omin). We then first erive the
More informationMA123, Chapter 10: Formulas for integrals: integrals, antiderivatives, and the Fundamental Theorem of Calculus (pp.
MA123, Chpter 1: Formuls for integrls: integrls, ntiderivtives, nd the Fundmentl Theorem of Clculus (pp. 27-233, Gootmn) Chpter Gols: Assignments: Understnd the sttement of the Fundmentl Theorem of Clculus.
More informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls
More informationNumerical Integration. Newton Cotes Formulas. Quadrature. Newton Cotes Formulas. To approximate the integral b
Numericl Integrtion Newton Cotes Formuls Given function f : R R nd two rel numbers, b R, < b, we clculte (pproximtely) the integrl I(f,, b) = f (x) dx K. Frischmuth (IfM UR) Numerics for CSE 08/09 8 /
More informationIntegrals - Motivation
Integrls - Motivtion When we looked t function s rte of chnge If f(x) is liner, the nswer is esy slope If f(x) is non-liner, we hd to work hrd limits derivtive A relted question is the re under f(x) (but
More informationMATHS NOTES. SUBJECT: Maths LEVEL: Higher TEACHER: Aidan Roantree. The Institute of Education Topics Covered: Powers and Logs
MATHS NOTES The Institute of Eduction 06 SUBJECT: Mths LEVEL: Higher TEACHER: Aidn Rontree Topics Covered: Powers nd Logs About Aidn: Aidn is our senior Mths techer t the Institute, where he hs been teching
More informationIntroduction to Group Theory
Introduction to Group Theory Let G be n rbitrry set of elements, typiclly denoted s, b, c,, tht is, let G = {, b, c, }. A binry opertion in G is rule tht ssocites with ech ordered pir (,b) of elements
More informationWHEN IS A FUNCTION NOT FLAT? 1. Introduction. {e 1 0, x = 0. f(x) =
WHEN IS A FUNCTION NOT FLAT? YIFEI PAN AND MEI WANG Abstrct. In this pper we prove unique continution property for vector vlued functions of one vrible stisfying certin differentil inequlity. Key words:
More informationQUADRATURE is an old-fashioned word that refers to
World Acdemy of Science Engineering nd Technology Interntionl Journl of Mthemticl nd Computtionl Sciences Vol:5 No:7 011 A New Qudrture Rule Derived from Spline Interpoltion with Error Anlysis Hdi Tghvfrd
More informationTravelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing
Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge
More informationOrdinary Differential Equations- Boundary Value Problem
Ordinry Differentil Equtions- Boundry Vlue Problem Shooting method Runge Kutt method Computer-bsed solutions o BVPFD subroutine (Fortrn IMSL subroutine tht Solves (prmeterized) system of differentil equtions
More informationMain topics for the First Midterm
Min topics for the First Midterm The Midterm will cover Section 1.8, Chpters 2-3, Sections 4.1-4.8, nd Sections 5.1-5.3 (essentilly ll of the mteril covered in clss). Be sure to know the results of the
More informationNumerical Integration. 1 Introduction. 2 Midpoint Rule, Trapezoid Rule, Simpson Rule. AMSC/CMSC 460/466 T. von Petersdorff 1
AMSC/CMSC 46/466 T. von Petersdorff 1 umericl Integrtion 1 Introduction We wnt to pproximte the integrl I := f xdx where we re given, b nd the function f s subroutine. We evlute f t points x 1,...,x n
More informationRecitation 3: More Applications of the Derivative
Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry dierentil eqution (ODE) du f(t) dt with initil condition u() : Just
More information