Math Week 5 concepts and homework, due Friday February 10

Similar documents
Chapter 2. Numerical Integration also called quadrature. 2.2 Trapezoidal Rule. 2.1 A basic principle Extending the Trapezoidal Rule DRAWINGS

Lecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University

NUMERICAL INTEGRATION

Math& 152 Section Integration by Parts

University of Houston, Department of Mathematics Numerical Analysis II

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

Z b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...

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

Key words. Numerical quadrature, piecewise polynomial, convergence rate, trapezoidal rule, midpoint rule, Simpson s rule, spectral accuracy.

Topic 6b Finite Difference Approximations

Math 20C Multivariable Calculus Lecture 5 1. Lines and planes. Equations of lines (Vector, parametric, and symmetric eqs.). Equations of lines

1 The Riemann Integral

Lecture 14: Quadrature

Best Approximation. Chapter The General Case

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

1 Part II: Numerical Integration

Chapter 7 Notes, Stewart 8e. 7.1 Integration by Parts Trigonometric Integrals Evaluating sin m x cos n (x) dx...

Numerical Integration

ODE: Existence and Uniqueness of a Solution

Numerical Analysis: Trapezoidal and Simpson s Rule

12 Basic Integration in R

Math 100 Review Sheet

B.Sc. in Mathematics (Ordinary)

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

Midpoint Approximation

Math 1B, lecture 4: Error bounds for numerical methods

Best Approximation in the 2-norm

Anti-derivatives/Indefinite Integrals of Basic Functions

Math 113 Exam 1-Review

Definite integral. Mathematics FRDIS MENDELU

Section 6.1 INTRO to LAPLACE TRANSFORMS

Math 360: A primitive integral and elementary functions

Chapters 4 & 5 Integrals & Applications

The Regulated and Riemann Integrals

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Review of basic calculus

Review of Calculus, cont d

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Numerical integration

CAAM 453 NUMERICAL ANALYSIS I Examination There are four questions, plus a bonus. Do not look at them until you begin the exam.

Orthogonal Polynomials

If u = g(x) is a differentiable function whose range is an interval I and f is continuous on I, then f(g(x))g (x) dx = f(u) du

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

1 Review: Volumes of Solids (Stewart )

f(a+h) f(a) x a h 0. This is the rate at which

6.5 Numerical Approximations of Definite Integrals

Integration Techniques

If deg(num) deg(denom), then we should use long-division of polynomials to rewrite: p(x) = s(x) + r(x) q(x), q(x)

3.4 Numerical integration

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Chapter 2 Differentiation

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Improper Integrals, and Differential Equations

Ordinary Differential Equations- Boundary Value Problem

Abstract inner product spaces

Math 113 Exam 2 Practice

1 The Lagrange interpolation formula

Math 113 Fall Final Exam Review. 2. Applications of Integration Chapter 6 including sections and section 6.8

LECTURE 19. Numerical Integration. Z b. is generally thought of as representing the area under the graph of fèxè between the points x = a and

Math 113 Exam 2 Practice

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

1 Techniques of Integration

Chapter 6 Notes, Larson/Hostetler 3e

Lecture 19: Continuous Least Squares Approximation

Multiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution

4181H Problem Set 11 Selected Solutions. Chapter 19. n(log x) n 1 1 x x dx,

Lecture 20: Numerical Integration III

Theoretical foundations of Gaussian quadrature

First midterm topics Second midterm topics End of quarter topics. Math 3B Review. Steve. 18 March 2009

New Expansion and Infinite Series

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Part IB Numerical Analysis

Numerical Integration. 1 Introduction. 2 Midpoint Rule, Trapezoid Rule, Simpson Rule. AMSC/CMSC 460/466 T. von Petersdorff 1

Chapter 3 Solving Nonlinear Equations

How can we approximate the area of a region in the plane? What is an interpretation of the area under the graph of a velocity function?

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

Chapter 0. What is the Lebesgue integral about?

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

Math 31S. Rumbos Fall Solutions to Assignment #16

Calculus and linear algebra for biomedical engineering Week 11: The Riemann integral and its properties

Chapter 5 : Continuous Random Variables

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

A P P E N D I X POWERS OF TEN AND SCIENTIFIC NOTATION A P P E N D I X SIGNIFICANT FIGURES

We know that if f is a continuous nonnegative function on the interval [a, b], then b

Math 116 Calculus II

Quantum Physics II (8.05) Fall 2013 Assignment 2

Linearity, linear operators, and self adjoint eigenvalue problems

1 Functions Defined in Terms of Integrals

COT4501 Spring Homework VII

Summer MTH142 College Calculus 2. Section J. Lecture Notes. Yin Su University at Buffalo

APPROXIMATE INTEGRATION

20 MATHEMATICS POLYNOMIALS

THE HANKEL MATRIX METHOD FOR GAUSSIAN QUADRATURE IN 1 AND 2 DIMENSIONS

DOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES

Notes on Calculus II Integral Calculus. Miguel A. Lerma

Chapter 8: Methods of Integration

Transcription:

Mt 2280-00 Week 5 concepts nd omework, due Fridy Februry 0 Recll tt ll problems re good for seeing if you cn work wit te underlying concepts; tt te underlined problems re to be nded in; nd tt te Fridy quiz will be drwn from ll of tese concepts nd from tese or relted problems 3 Solving initil vlue problems for liner omogeneous second order differentil equtions, given bsis for te solution spce Finding generl solutions for constnt coefficient omogeneous DE's by sercing for exponentil or oter functions Superposition for liner differentil equtions, nd its filure for nonliner DE's,, (in use initil vlues y 0 = 0, y 0 = 5 rter tn te ones in te text), 0,, 2, 4 (In 4 use te initil vlues y =, y = 7 rter tn te ones in te text), 7, 8, 27, 33, 39 w53) In 3 bove, te text tells you tt y x = e 2 x x = e 3 x re two independent solutions to te second order omogeneous differentil eqution y y y = 0 Verify tt you could ve found tese two exponentil solutions vi te following guessing lgoritm: Try y x = e r x were te constnt r is to be determined Substitute tis possible solution into te omogeneous differentil eqution nd find te only two vlues of r for wic y x will stisfy te DE (See Teorem 5 in te text) w53b) In 30 bove, te text tells you tt y x = e 5 x x = x e 5 x re two independent solutions to te second order omogeneous differentil eqution y 0 y 25 y = 0 Follow te procedure in prt of trying for solutions of te form y x = e r x, nd ten use te repeted roots Teorem in te text, to recover tese two solutions 32 Testing collections of functions for dependence nd independence Solving IVP's for omogeneous nd non-omogeneous differentil equtions Superposition, 2, 5, 8,,3,, 2, 25, 2 Here re two problems tt explicitly connect ides from sections 3-32 wit liner lgebr concepts from Mt 2270 w54) Consider te 3 rd order omogeneous liner differentil eqution for y x y x = 0 nd let W be te solution spce w54) Use successive ntidifferentition to solve tis differentil eqution Interpret your results using vector spce concepts to sow tt te functions y 0 x =, y x = x x = x 2 re bsis for W Tus te dimension of W is 3 w54b) Sow tt te functions z 0 x =, z x = x 2, z 2 x = 2 x 2 2 re lso bsis for W Hint: If you verify tt tey solve te differentil eqution nd tt tey're linerly independent, tey will utomticlly spn te 3-dimensionl solution spce nd terefore be bsis w54c) Use liner combintion of te solution bsis from prt b, in order to solve te initil vlue

problem below Notice ow tis bsis is dpted to initil vlue problems t x 0 = 2, weres for n IVP t x 0 = 0 te bsis in would ve been esier to use y x = 0 y 2 = 7 y 2 = 3 y 2 = 5 w55) Consider te tree functions y x = cos 2 x x = sin 2 x, y 3 x = sin 2 x w55) Sow tt ll tree functions solve te differentil eqution y 4 y = 0 w5b) Te differentil eqution bove is second order liner omogeneous DE, so te solution spce is 2 dimensionl Tus te tree functions y, y 3 bove must be linerly dependent Find liner dependency (Hint: use trigonometry ddition ngle formul) w55c) Explicitly verify tt every initil vlue problem y 4 y = 0 y 0 = b y 0 = b 2 s solution of te form y x = c cos 2 x c 2 sin 2 x, nd tt c, c 2 re uniquely determined by b, b 2 (Tus cos 2 x, sin 2 x re bsis for te solution spce of y 4 y = 0: every solution y x s initil vlues tt cn be mtced wit liner combintion of y, but once te initil vlues mtc te solutions must gree by te uniqueness teorem, so y spn te solution spce; y re linerly independent becuse if c cos 2 x c 2 sin 2 x = y x 0 ten y 0 = y 0 = 0 so lso c = c 2 = 0) w55d) Find by inspection, prticulr solutions y x to te two non-omogeneous differentil equtions y 4 y = 28, y 4 y = x Hint: one of tem could be constnt, te oter could be multiple of x w55e) Use superposition (linerity) nd your work from c,d to find te generl solution to te nonomogeneous differentil eqution y 4 y = 28 x w55f) Solve te initil vlue problem, using your work bove: y 4 y = 28 x y 0 = 0 y 0 = 0

33: using te lgoritm for finding te generl solution to constnt coefficient liner omogeneous differentil equtions: rel roots, Euler's formul nd complex roots, repeted roots; solving ssocited initil vlue problems 33: 3, 9,, 23, 27 NOTE: w5 nd w57 re postponed until next week (Fridy) We will tlk bout section 33 tis Fridy Feb 0 w5) Do te following problems for omogeneous liner differentil equtions by nd Tey re testing your bility to use te lgoritm for finding bses for te solution spces, bsed on te crcteristic polynomil Ceck your work wit Mple (or oter softwre) In Mple you will wnt to use te "dsolve" commnd to ceck differentil eqution solutions, nd my wnt to use te "fctor" commnd to ceck your fctoriztions of te crcteristic polynomils Hnd in printout of your computer verifictions for te differentil eqution solutions, long wit your written work ) Find te generl solution to te differentil eqution for y x y 3 5 y 3 y 9 y = 0 Hint: Find root r of te cubic crcteristic polynomil, ten divide it by r r to get quotient qudrtic polynomil b) Find te generl solution to te differentil eqution for x t x 4 x 3 x = 0 Hint: completing te squre works well ere - probbly better tn te qudrtic formul c) Solve te initil vlue problem for te differentil eqution in b, wit x 0 = 0, x 0 = 9 d) Find te generl solution to te differentil eqution for y x y 4 8 y = 0 e) Find te generl solution to te differentil eqution for y x y 5 y 3 9 y = 0 w57) Euler's formul e i = cos i sin is extremely useful in iger mtemtics/science/engineering In clss we discuss ow tis definition is motivted by Tylor series Amzingly, te rule of exponents is true for suc expressions In oter words e i i = e i e i Ceck tis identity by rewriting te left side s e i nd ten using Euler's formul to expnd bot sides You'll notice tt te identity is true becuse of te ddition ngle formuls for cos nd sin (And so tis gives good wy to recover te trig identities if you ppen to forget tem) 24-2 Numericl metods problems re t te end of tis ssignment w5) Runge-Kutt is bsed on Simpson's rule for numericl integrtion Simpson's rule is bsed on te fct tt for subintervl d, d of lengt, te prbol y = p x wic psses troug te points

d, y 0, d 2, y, d s integrl d d p x dx = y 0 4 y y 2 w5) Te integrl pproximtion bove follows from one on te intervl, by n ffine cnge of vribles So first consider te intervl, We wis to find te prbolic function q x = x 2 b x c wit unknown prmeters, b, c We wnt q = y 0, q 0 = y, q = y 2 Tis gives 3 equtions in 3 unknowns, to find, b, c in terms of y 0, y Write down tese liner equtions nd find, b, c w5b) Compute q x dx for tes vlues of, b, c you find in prt, nd verify te identity q x dx = 3 y 0 4y y 2 Note, te formul for generl intervl follows from cnge of vribles, s indicted below: Remrk: If you've forgotten, or if you never tlked bout Simpson's rule in your Clculus clss, ere's ow it goes: In order to pproximte te definite integrl of f x on te intervl, b, you subdivide, b into n subintervls of widt x = b n = Ten dd te midpoints of ec subintervl Lbel tese x vlues (incluidng midpoints) s x 0 =, x = 2, x 2 =, x 3 = x 3 0 2, x 4 = x 0 2, x 2 n = x 0 2n = b, wit corresponding y vlues y i = f x i, i = 0,2 n On ec successive pir of intervls x 2 k, x 2 k use te prbolic estimte x 2 k x 2 k f u du f x 2 k 4 f x 2 k f x 2 k 2 = y 2 k 4 y 2 k y 2 k 2 bove, estimting te integrl of f by te integrl of te interpolting prbol on te subintervl Tis

yields te very ccurte (for lrge enoug n) Simpson's rule formul b f x dx y 4 y y y 4 y y y 4 y y, 0 2 2 3 4 2 n 2 2 n 2 n ie b f x dx b n y 0 4 y 2 y 2 4 y 3 2 y 4 2 y 2 n 2 4 y 2 n y 2 n ttp://enwikipediorg/wiki/simpson's_rule w52) (Fmous numbers revisited, section 2, pge 35, of text) Te mty numbers e, ln 2, cn be well-pproximted using pproximte solutions to differentil equtions We illustrte tis on Wednesdy Feb 4 for e, wic is y for te solution to te IVP y x = y y 0 = Apply Runge-Kutt wit n = 0, 20, 40 subintervls, successively doubling te number of subintervls until you obtin te trget number below - rounded to 9 deciml digits - twice in succession We will do tis in clss for e, nd you cn modify tt code if you wis w52) ln 2 is y 2, were y x solves te IVP (since y x = ln x ) y x = x y = 0 w52b) is y, were y x solves te IVP 4 y x = x 2 y 0 = 0 (since y x = 4 rctn x ) Note tt in,b you re ctully "just" using Simpson's rule from Clculus, since te rigt sides of tese DE's only depend on te vrible x nd not on te vlue of te function y x For reference: > Digits : #ow mny digits to use in floting point numbers nd clculultions evlf e ; #evlute te floting point of e evlf ; 27828828459045 3459253589793 ()