Taylor Series Method for Second-Order Polynomial Differential Equations
|
|
- Alexia Young
- 5 years ago
- Views:
Transcription
1 Taylor Series Method for Second-Order Polynomial Differential Equations Viktor N. Latypov Sergey V. Sokolov Saint-Petersburg State University SCP-2015
2 Introduction New numerical algorithm to solve second-order polynomial ODEs Based on series expansion
3 Introduction New numerical algorithm to solve second-order polynomial ODEs Based on series expansion Uses only basic linear algebra
4 Introduction New numerical algorithm to solve second-order polynomial ODEs Based on series expansion Uses only basic linear algebra Easy to implement
5 Lagrange Equations Mechanics/Lagrange equations D(q) q + C(q, q) q + K(q) = 0
6 Lagrange Equations Mechanics/Lagrange equations D(q) q + C(q, q) q + K(q) = 0 General second-order system A(q) q + b(q, q) = 0
7 Lagrange Equations Mechanics/Lagrange equations General second-order system D(q) q + C(q, q) q + K(q) = 0 A(q) q + b(q, q) = 0 Goal: Derive series expansion for the solution q k (t) = + p=0 q p k (t t 0) p
8 First-Order systems For the first-order Cauchy problem i N n j 1 ẋ k = a k +... a k;i1,i 2,...,i j x i1... x ij. j=1 i 1=1 i j =1 x k (t 0 ) = x k,0, k = 1,..., n.
9 First-Order systems For the first-order Cauchy problem i N n j 1 ẋ k = a k +... a k;i1,i 2,...,i j x i1... x ij. j=1 i 1=1 i j =1 x k (t 0 ) = x k,0, k = 1,..., n. we can directly substitute x k (t) = + q=0 x q k (t t 0) q
10 First-Order systems For the first-order Cauchy problem i N n j 1 ẋ k = a k +... a k;i1,i 2,...,i j x i1... x ij. j=1 i 1=1 i j =1 x k (t 0 ) = x k,0, k = 1,..., n. we can directly substitute to get (for N = 2) x q+1 k = 1 q + 1 x k (t) = + q=0 x q k (t t 0) q δ q a k + n i=1 a k;ix q i + + n i i=1 j=1 a ] q k;ij l=0 x i lx q l j. (1)
11 First-Order system: an example For the non-polynomial equation ẋ 1 = sin(x 1 ) + e x1
12 First-Order system: an example For the non-polynomial equation ẋ 1 = sin(x 1 ) + e x1 one may find a set of auxiliary variables x 2 = sin(x 1 ) x 3 = cos(x 1 ) x 4 = e x1
13 First-Order system: an example For the non-polynomial equation ẋ 1 = sin(x 1 ) + e x1 one may find a set of auxiliary variables x 2 = sin(x 1 ) x 3 = cos(x 1 ) x 4 = e x1 such that the new system becomes polynomial: ẋ 1 = x 2 + x 4 ẋ 2 = x 3 (x 2 + x 4 ) ẋ 3 = x 2 (x 2 + x 4 ) ẋ 4 = x 4 (x 2 + x 4 )
14 Second-Order system: an example For the double pendulum we have a non-polynomial equation with ( ) (m1 + m D = 2 ) L 2 1 m 2 L 1 L 2 c 1 2 m 2 L 1 L 2 c 1 2 m 2 L 2, 2
15 Second-Order system: an example For the double pendulum we have a non-polynomial equation with ( ) (m1 + m D = 2 ) L 2 1 m 2 L 1 L 2 c 1 2 m 2 L 1 L 2 c 1 2 m 2 L 2, 2 ( ) m2 L C = 1 L 2 s 1 2 θ2 0, m 2 L 1 L 2 s 1 2 θ1 0
16 Second-Order system: an example For the double pendulum we have a non-polynomial equation with ( ) (m1 + m D = 2 ) L 2 1 m 2 L 1 L 2 c 1 2 m 2 L 1 L 2 c 1 2 m 2 L 2, 2 ( ) m2 L C = 1 L 2 s 1 2 θ2 0, m 2 L 1 L 2 s 1 2 θ1 0 ( ) L1 (m K = g 1 + m 2 ) s 1. m 2 L 2 s 2
17 Second-Order system: an example For the double pendulum we have a non-polynomial equation with ( ) (m1 + m D = 2 ) L 2 1 m 2 L 1 L 2 c 1 2 m 2 L 1 L 2 c 1 2 m 2 L 2, 2 ( ) m2 L C = 1 L 2 s 1 2 θ2 0, m 2 L 1 L 2 s 1 2 θ1 0 ( ) L1 (m K = g 1 + m 2 ) s 1. m 2 L 2 s 2 Here q 1 = θ 1 and q 2 = θ 2, s i = sin θ i, c i = cos θ i, s 1 2 = sin (θ 1 θ 2 ).
18 Second-Order system: polynomial form Introducing q 3 = s 1, q 4 = c 1, q 5 = s 2 and q 6 = c 2 we reduce original equations to the polynomial one.
19 Second-Order system: polynomial form Introducing q 3 = s 1, q 4 = c 1, q 5 = s 2 and q 6 = c 2 we reduce original equations to the polynomial one. First derivative with respect to t: q 3 = q 4 q 1 q 4 = q 3 q 1 q 5 = q 6 q 2 q 6 = q 5 q 2 (2)
20 Second-Order system: polynomial form Introducing q 3 = s 1, q 4 = c 1, q 5 = s 2 and q 6 = c 2 we reduce original equations to the polynomial one. First derivative with respect to t: Second derivative: q 3 = q 4 q 1 + q 4 q 1 q 4 = q 3 q 1 q 3 q 1 q 5 = q 6 q 2 + q 6 q 2 q 6 = q 5 q 2 q 3 q 2 q 3 = q 4 q 1 q 4 = q 3 q 1 q 5 = q 6 q 2 q 6 = q 5 q 2 (2)
21 Second-Order system: polynomial form Introducing q 3 = s 1, q 4 = c 1, q 5 = s 2 and q 6 = c 2 we reduce original equations to the polynomial one. First derivative with respect to t: Second derivative: q 3 = q 4 q 1 + q 4 q 1 q 4 = q 3 q 1 q 3 q 1 q 5 = q 6 q 2 + q 6 q 2 q 6 = q 5 q 2 q 3 q 2 q 3 = q 4 q 1 q 4 = q 3 q 1 q 5 = q 6 q 2 q 6 = q 5 q 2 (2)
22 Second-Order system: polynomial form Introducing q 3 = s 1, q 4 = c 1, q 5 = s 2 and q 6 = c 2 we reduce original equations to the polynomial one. First derivative with respect to t: q 3 = q 4 q 1 q 4 = q 3 q 1 q 5 = q 6 q 2 q 6 = q 5 q 2 (2) Second derivative: q 3 = q 4 q 1 + q 4 q 1 q 4 = q 3 q 1 q 3 q 1 q 5 = q 6 q 2 + q 6 q 2 q 6 = q 5 q 2 q 3 q 2 q 4 q 1 + q 3 = q 4 q 1 q 3 q 1 + q 4 = q 3 q 1 q 6 q 2 + q 5 = q 6 q 2 q 3 q 2 + q 6 = q 5 q 2
23 Second-Order system: polynomial form D = C = d 11 d d 21 d q q q q c 11 c c 21 c q q q q ,.
24 Polynomial system: notations Equation A(q) q + b(q, q) = 0 where A = (a ij ) n i,j=1 and b = (b 1,..., b n ) T.
25 Polynomial system: notations Equation A(q) q + b(q, q) = 0 where A = (a ij ) n i,j=1 and b = (b 1,..., b n ) T. Multi-indices q p = q p qpn n, p = p p n, r = r r n, a p1p 2...,p n,r 1,r 2...,r n;ij = a pr;ij, b p1p 2...,p n,r 1,r 2...,r n;i = b pr;i. (3)
26 Polynomial system: notations Equation A(q) q + b(q, q) = 0 where A = (a ij ) n i,j=1 and b = (b 1,..., b n ) T. Multi-indices q p = q p qpn n, p = p p n, r = r r n, a p1p 2...,p n,r 1,r 2...,r n;ij = a pr;ij, b p1p 2...,p n,r 1,r 2...,r n;i = b pr;i. (3) Components a ij, b i polynomials of degree α ij and β i respectively. a ij (q, q) = b i (q, q) = α ij s=0 p + r =s β i s=0 p + r =s a pr;ij q p q r, b pr;i q p q r
27 Series expansion: shortcuts Substituting q i (t) = + m=0 q i m(t t 0 ) m into A(q) q + b(q, q) = 0 requires
28 Series expansion: shortcuts Substituting q i (t) = + m=0 q i m(t t 0 ) m into A(q) q + b(q, q) = 0 requires 1. Series products h = f g = + m=0 k=0 m (f k g m k ) (t t 0 ) m = h m = (f m g m ) = + m=0 m f k g m k. k=0 h m (t t 0 ) m.
29 Series expansion: shortcuts 2. Series powers Arbitrary power of q i (t) is q p i (t) = + q p m=0 i m = m k=0 q p i m (t t 0) m, ( q p 1 i k q i m k ).
30 Series expansion: shortcuts 3. Derivatives Derivatives of q i (t): + q i r (t) = (m + 1)qi m+1 r (t t 0) m, q i (t) = m=0 + m=0 (m + 1)(m + 2)q i m+2 (t t 0 ) m.
31 Series expansion: equations for coefficients Components of A(q) q + b(q, q) = 0 are n (a ij m ((m + 1)(m + 2)q j m+2 )) + b i m = 0. j=1
32 Series expansion: equations for coefficients Components of A(q) q + b(q, q) = 0 are n (a ij m ((m + 1)(m + 2)q j m+2 )) + b i m = 0. j=1 a ij m = b i m = α ij s=0 p + r =s β i s=0 p + r =s a pr;ij q p m [ (m + 1) q m+1 ] r. (4) b pr;i q p m [ (m + 1) q m+1 ] r. (5)
33 Series expansion: equations for coefficients Components of A(q) q + b(q, q) = 0 are n (a ij m ((m + 1)(m + 2)q j m+2 )) + b i m = 0. j=1 a ij m = b i m = α ij s=0 p + r =s β i s=0 p + r =s a pr;ij q p m [ (m + 1) q m+1 ] r. (4) b pr;i q p m [ (m + 1) q m+1 ] r. (5) q p m = q p1 1 m q p2 2 m... q pn n m (6)
34 Series expansion: equations for coefficients Components of A(q) q + b(q, q) = 0 are n (a ij m ((m + 1)(m + 2)q j m+2 )) + b i m = 0. j=1 a ij m = b i m = q p α ij s=0 p + r =s β i s=0 p + r =s i m = m a pr;ij q p m [ (m + 1) q m+1 ] r. (4) b pr;i q p m [ (m + 1) q m+1 ] r. (5) q p m = q p1 1 m q p2 2 m... q pn n m (6) k=0 q p 1 i k q i m k = q i m... q i m. }{{} p
35 Series expansion: equations for coefficients For q m+2 we get (m + 1)(m + 2)A 0 q m+2 + c m = 0,
36 Series expansion: equations for coefficients For q m+2 we get (m + 1)(m + 2)A 0 q m+2 + c m = 0, where A 0 = { a ij 0 } n i,j=1, c i;0 = b i 0, and for m > 0
37 Series expansion: equations for coefficients For q m+2 we get (m + 1)(m + 2)A 0 q m+2 + c m = 0, where A 0 = { a ij 0 } n i,j=1, c i;0 = b i 0, and for m > 0 c i;m = m 1 k=0 j=1 n (a ij m k ((k + 1)(k + 2)q j k+2 )) + b i m. (7)
38 Summary Given initial conditions q 0 and q 1 for each m q m+2 = A 1 0 c m (m + 1) (m + 2), A 0 = { a ij 0 } n i,j.
Chapter 8: Taylor s theorem and L Hospital s rule
Chapter 8: Taylor s theorem and L Hospital s rule Theorem: [Inverse Mapping Theorem] Suppose that a < b and f : [a, b] R. Given that f (x) > 0 for all x (a, b) then f 1 is differentiable on (f(a), f(b))
More informationCHALLENGE! (0) = 5. Construct a polynomial with the following behavior at x = 0:
TAYLOR SERIES Construct a polynomial with the following behavior at x = 0: CHALLENGE! P( x) = a + ax+ ax + ax + ax 2 3 4 0 1 2 3 4 P(0) = 1 P (0) = 2 P (0) = 3 P (0) = 4 P (4) (0) = 5 Sounds hard right?
More informationTaylor Series. Math114. March 1, Department of Mathematics, University of Kentucky. Math114 Lecture 18 1/ 13
Taylor Series Math114 Department of Mathematics, University of Kentucky March 1, 2017 Math114 Lecture 18 1/ 13 Given a function, can we find a power series representation? Math114 Lecture 18 2/ 13 Given
More informationSection 8.7. Taylor and MacLaurin Series. (1) Definitions, (2) Common Maclaurin Series, (3) Taylor Polynomials, (4) Applications.
Section 8.7 Taylor and MacLaurin Series (1) Definitions, (2) Common Maclaurin Series, (3) Taylor Polynomials, (4) Applications. MATH 126 (Section 8.7) Taylor and MacLaurin Series The University of Kansas
More informationTransverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization)
Transverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization) Anton Shiriaev 1,2, Leonid Freidovich 1, Sergey Gusev 3 1 Department
More information8.5 Taylor Polynomials and Taylor Series
8.5. TAYLOR POLYNOMIALS AND TAYLOR SERIES 50 8.5 Taylor Polynomials and Taylor Series Motivating Questions In this section, we strive to understand the ideas generated by the following important questions:
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationReview: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function.
Taylor Series (Sect. 10.8) Review: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function. Review: Power series define functions Remarks:
More informationPower Series Solutions We use power series to solve second order differential equations
Objectives Power Series Solutions We use power series to solve second order differential equations We use power series expansions to find solutions to second order, linear, variable coefficient equations
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationProblem 1 Kaplan, p. 436: 2c,i
Mathematical Methods I - Fall 03 Homework Solutions Page Problem Kaplan, p 36: c,i Find the first three nonzero terms of the following Taylor series: c) ln + x) about x = 0 i) arctanh x about x = 0 The
More informationReview for Exam 2. Review for Exam 2.
Review for Exam 2. 5 or 6 problems. No multiple choice questions. No notes, no books, no calculators. Problems similar to homeworks. Exam covers: Regular-singular points (5.5). Euler differential equation
More information1 Lecture 24: Linearization
1 Lecture 24: Linearization 1.1 Outline The linearization of a function at a point a. Linear approximation of the change in f. Error, absolute error. Examples 1.2 Linearization Functions can be complicated.
More informationf (x) = k=0 f (0) = k=0 k=0 a k k(0) k 1 = a 1 a 1 = f (0). a k k(k 1)x k 2, k=2 a k k(k 1)(0) k 2 = 2a 2 a 2 = f (0) 2 a k k(k 1)(k 2)x k 3, k=3
1 M 13-Lecture Contents: 1) Taylor Polynomials 2) Taylor Series Centered at x a 3) Applications of Taylor Polynomials Taylor Series The previous section served as motivation and gave some useful expansion.
More informationRepresent this system in terms of a block diagram consisting only of. g From Newton s law: 2 : θ sin θ 9 θ ` T
Exercise (Block diagram decomposition). Consider a system P that maps each input to the solutions of 9 4 ` 3 9 Represent this system in terms of a block diagram consisting only of integrator systems, represented
More informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More informationMath 113 (Calculus 2) Exam 4
Math 3 (Calculus ) Exam 4 November 0 November, 009 Sections 0, 3 7 Name Student ID Section Instructor In some cases a series may be seen to converge or diverge for more than one reason. For such problems
More informationSection 5.8. Taylor Series
Difference Equations to Differential Equations Section 5.8 Taylor Series In this section we will put together much of the work of Sections 5.-5.7 in the context of a discussion of Taylor series. We begin
More informationMathematics 104 Fall Term 2006 Solutions to Final Exam. sin(ln t) dt = e x sin(x) dx.
Mathematics 14 Fall Term 26 Solutions to Final Exam 1. Evaluate sin(ln t) dt. Solution. We first make the substitution t = e x, for which dt = e x. This gives sin(ln t) dt = e x sin(x). To evaluate the
More informationMain topics and some repetition exercises for the course MMG511/MVE161 ODE and mathematical modeling in year Main topics in the course:
Main topics and some repetition exercises for the course MMG5/MVE6 ODE and mathematical modeling in year 04. Main topics in the course:. Banach fixed point principle. Picard- Lindelöf theorem. Lipschitz
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More informationDifferentiation. Table of contents Definition Arithmetics Composite and inverse functions... 5
Differentiation Table of contents. Derivatives................................................. 2.. Definition................................................ 2.2. Arithmetics...............................................
More informationChapter 1 Numerical approximation of data : interpolation, least squares method
Chapter 1 Numerical approximation of data : interpolation, least squares method I. Motivation 1 Approximation of functions Evaluation of a function Which functions (f : R R) can be effectively evaluated
More informationENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations
ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential
More informationChapter 4: Partial differentiation
Chapter 4: Partial differentiation It is generally the case that derivatives are introduced in terms of functions of a single variable. For example, y = f (x), then dy dx = df dx = f. However, most of
More informationSecond quantization: where quantization and particles come from?
110 Phys460.nb 7 Second quantization: where quantization and particles come from? 7.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system? 7.1.1.Lagrangian Lagrangian
More informationCalculus I Exam 1 Review Fall 2016
Problem 1: Decide whether the following statements are true or false: (a) If f, g are differentiable, then d d x (f g) = f g. (b) If a function is continuous, then it is differentiable. (c) If a function
More information8.7 Taylor s Inequality Math 2300 Section 005 Calculus II. f(x) = ln(1 + x) f(0) = 0
8.7 Taylor s Inequality Math 00 Section 005 Calculus II Name: ANSWER KEY Taylor s Inequality: If f (n+) is continuous and f (n+) < M between the center a and some point x, then f(x) T n (x) M x a n+ (n
More informationReview of Linear Algebra Definitions, Change of Basis, Trace, Spectral Theorem
Review of Linear Algebra Definitions, Change of Basis, Trace, Spectral Theorem Steven J. Miller June 19, 2004 Abstract Matrices can be thought of as rectangular (often square) arrays of numbers, or as
More informationMath 106 Fall 2014 Exam 2.1 October 31, ln(x) x 3 dx = 1. 2 x 2 ln(x) + = 1 2 x 2 ln(x) + 1. = 1 2 x 2 ln(x) 1 4 x 2 + C
Math 6 Fall 4 Exam. October 3, 4. The following questions have to do with the integral (a) Evaluate dx. Use integration by parts (x 3 dx = ) ( dx = ) x3 x dx = x x () dx = x + x x dx = x + x 3 dx dx =
More informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
More informationA Sumudu based algorithm for solving differential equations
Computer Science Journal of Moldova, vol.15, no.3(45), 2007 A Sumudu based algorithm for solving differential equations Jun Zhang Abstract An algorithm based on Sumudu transform is developed. The algorithm
More informationMSMS Matlab Problem 02
MSMS 2014-2015 Matlab Problem 02 Basilio Bona DAUIN PoliTo Problem formulation The planar system illustrated in Figure 1 consists of a cart C sliding with friction along the horizontal rail; the cart supports
More informationSection 5.2 Series Solution Near Ordinary Point
DE Section 5.2 Series Solution Near Ordinary Point Page 1 of 5 Section 5.2 Series Solution Near Ordinary Point We are interested in second order homogeneous linear differential equations with variable
More informationSection 10.7 Taylor series
Section 10.7 Taylor series 1. Common Maclaurin series 2. s and approximations with Taylor polynomials 3. Multiplication and division of power series Math 126 Enhanced 10.7 Taylor Series The University
More informationThe geometry of polynomial diagrams
Theoretical Mathematics & Applications, vol.4, no.1, 014, 71-76 ISSN: 179-9687 (print), 179-9709 (online) Scienpress Ltd, 014 The geometry of polynomial diagrams Alennikov A. Maksim 1 Abstract In this
More informationAnnouncements. Related Rates (last week), Linear approximations (today) l Hôpital s Rule (today) Newton s Method Curve sketching Optimization problems
Announcements Assignment 4 is now posted. Midterm results should be available by the end of the week (assuming the scantron results are back in time). Today: Continuation of applications of derivatives:
More informationMath 181, Exam 2, Study Guide 2 Problem 1 Solution. 1 + dx. 1 + (cos x)2 dx. 1 + cos2 xdx. = π ( 1 + cos π 2
Math 8, Exam, Study Guide Problem Solution. Use the trapezoid rule with n to estimate the arc-length of the curve y sin x between x and x π. Solution: The arclength is: L b a π π + ( ) dy + (cos x) + cos
More informationMATH 31B: MIDTERM 2 REVIEW. sin 2 x = 1 cos(2x) dx = x 2 sin(2x) 4. + C = x 2. dx = x sin(2x) + C = x sin x cos x
MATH 3B: MIDTERM REVIEW JOE HUGHES. Evaluate sin x and cos x. Solution: Recall the identities cos x = + cos(x) Using these formulas gives cos(x) sin x =. Trigonometric Integrals = x sin(x) sin x = cos(x)
More informationTopics Covered in Calculus BC
Topics Covered in Calculus BC Calculus BC Correlation 5 A Functions, Graphs, and Limits 1. Analysis of graphs 2. Limits or functions (including one sides limits) a. An intuitive understanding of the limiting
More informationMath Review for Exam Answer each of the following questions as either True or False. Circle the correct answer.
Math 22 - Review for Exam 3. Answer each of the following questions as either True or False. Circle the correct answer. (a) True/False: If a n > 0 and a n 0, the series a n converges. Soln: False: Let
More informationAdvanced Eng. Mathematics
Koya University Faculty of Engineering Petroleum Engineering Department Advanced Eng. Mathematics Lecture 6 Prepared by: Haval Hawez E-mail: haval.hawez@koyauniversity.org 1 Second Order Linear Ordinary
More informationExamples of the Fourier Theorem (Sect. 10.3). The Fourier Theorem: Continuous case.
s of the Fourier Theorem (Sect. 1.3. The Fourier Theorem: Continuous case. : Using the Fourier Theorem. The Fourier Theorem: Piecewise continuous case. : Using the Fourier Theorem. The Fourier Theorem:
More informationSample Quiz 8, Problem 1. Solving Higher Order Constant-Coefficient Equations
Sample Quiz 8 Sample Quiz 8, Problem. Solving Higher Order Constant-Coefficient Equations The Algorithm applies to constant-coefficient homogeneous linear differential equations of order N, for example
More informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationTaylor Series and Maclaurin Series
Taylor Series and Maclaurin Series Definition (Taylor Series) Suppose the function f is infinitely di erentiable at a. The Taylor series of f about a (or at a or centered at a) isthepowerseries f (n) (a)
More informationb n x n + b n 1 x n b 1 x + b 0
Math Partial Fractions Stewart 7.4 Integrating basic rational functions. For a function f(x), we have examined several algebraic methods for finding its indefinite integral (antiderivative) F (x) = f(x)
More informationChapter 4 Sequences and Series
Chapter 4 Sequences and Series 4.1 Sequence Review Sequence: a set of elements (numbers or letters or a combination of both). The elements of the set all follow the same rule (logical progression). The
More informationPhysics GRE Practice
Physics GRE Practice Chapter 3: Harmonic Motion in 1-D Harmonic Motion occurs when the acceleration of a system is proportional to the negative of its displacement. or a x (1) ẍ x (2) Examples of Harmonic
More informationCSL105: Discrete Mathematical Structures. Ragesh Jaiswal, CSE, IIT Delhi
Definition (Linear homogeneous recurrence) A linear homogeneous recurrence relation of degree k with constant coefficients is a recurrence relation of the form a n = c 1 a n 1 + c 2 a n 2 +... + c k a
More informationSESSION 6 Trig. Equations and Identities. Math 30-1 R 3. (Revisit, Review and Revive)
SESSION 6 Trig. Equations and Identities Math 30-1 R 3 (Revisit, Review and Revive) 1 P a g e 2 P a g e Mathematics 30-1 Learning Outcomes Specific Outcome 5: Solve, algebraically and graphically, first
More informationMATH 167: APPLIED LINEAR ALGEBRA Chapter 3
MATH 167: APPLIED LINEAR ALGEBRA Chapter 3 Jesús De Loera, UC Davis February 18, 2012 Orthogonal Vectors and Subspaces (3.1). In real life vector spaces come with additional METRIC properties!! We have
More informationTransformation Matrices; Geometric and Otherwise As examples, consider the transformation matrices of the C 3v
Transformation Matrices; Geometric and Otherwise As examples, consider the transformation matrices of the v group. The form of these matrices depends on the basis we choose. Examples: Cartesian vectors:
More informationLecture 5: Function Approximation: Taylor Series
1 / 10 Lecture 5: Function Approximation: Taylor Series MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Better
More information2003 Mathematics. Advanced Higher. Finalised Marking Instructions
2003 Mathematics Advanced Higher Finalised Marking Instructions 2003 Mathematics Advanced Higher Section A Finalised Marking Instructions Advanced Higher 2003: Section A Solutions and marks A. (a) Given
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationPHYS 408, Optics. Problem Set 4 - Spring Posted: Fri, March 4, 2016 Due: 5pm Thu, March 17, 2016
PHYS 408, Optics Problem Set 4 - Spring 06 Posted: Fri, March 4, 06 Due: 5pm Thu, March 7, 06. Refraction at a Spherical Boundary. Derive the M matrix of.4-6 in the textbook. You may use Snell s Law directly..
More informationSeries Solutions of Linear ODEs
Chapter 2 Series Solutions of Linear ODEs This Chapter is concerned with solutions of linear Ordinary Differential Equations (ODE). We will start by reviewing some basic concepts and solution methods for
More information7.2 Parts (Undoing the Product Rule)
7.2 Parts (Undoing the Product Rule) 7.1 Substitution (Undoing the Chain Rule) Often w inside of some other function, dw up to constant What I want you to show me... w, dw, with respect to w 7.2 Parts
More informationMTH30 Review Sheet. y = g(x) BRONX COMMUNITY COLLEGE of the City University of New York DEPARTMENT OF MATHEMATICS & COMPUTER SCIENCE
BRONX COMMUNITY COLLEGE of the City University of New York DEPARTMENT OF MATHEMATICS & COMPUTER SCIENCE MTH0 Review Sheet. Given the functions f and g described by the graphs below: y = f(x) y = g(x) (a)
More informationy= sin3 x+sin6x x 1 1 cos(2x + 4 ) = cos x + 2 = C(x) (M2) Therefore, C(x) is periodic with period 2.
. (a).5 0.5 y sin x+sin6x 0.5.5 (A) (C) (b) Period (C) []. (a) y x 0 x O x Notes: Award for end points Award for a maximum of.5 Award for a local maximum of 0.5 Award for a minimum of 0.75 Award for the
More informationComplex Analysis Prelim Written Exam Spring 2015
Prelim Written Exam Spring 2015 Questions are equally weighted. Give essential explanations and justifications: a large part of each question is demonstration that you understand the context and understand
More informationCHANGE OF BASIS FOR LINEAR SPACES
CHANGE OF ASIS FOR LINEAR SPACES. Motivating Example Consider P 2 the space of polynomials of degree at most 2. Let : P 2 P 2 be defined by (p)(x) = p (x) p(x). his is a linear transformation. We would
More informationRESEARCH ARTICLE. Spectral Method for Solving The General Form Linear Fredholm Volterra Integro Differential Equations Based on Chebyshev Polynomials
Journal of Modern Methods in Numerical Mathematics ISSN: 29-477 online Vol, No, 2, c 2 Modern Science Publishers wwwm-sciencescom RESEARCH ARTICLE Spectral Method for Solving The General Form Linear Fredholm
More informationDistribution results for additive functions
Distribution results for additive functions Manfred G. Madritsch Institute of Statistics Graz University of Technology madritsch@finanz.math.tugraz.at Universiteit Stellenbosch Supported by the Austrian
More informationCombinatorial Analysis of the Geometric Series
Combinatorial Analysis of the Geometric Series David P. Little April 7, 205 www.math.psu.edu/dlittle Analytic Convergence of a Series The series converges analytically if and only if the sequence of partial
More informationQUANTUM MECHANICS I PHYS 516. Solutions to Problem Set # 5
QUANTUM MECHANICS I PHYS 56 Solutions to Problem Set # 5. Crossed E and B fields: A hydrogen atom in the N 2 level is subject to crossed electric and magnetic fields. Choose your coordinate axes to make
More information1.2 Finite Precision Arithmetic
MACM Assignment Solutions.2 Finite Precision Arithmetic.2:e Rounding Arithmetic Use four-digit rounding arithmetic to perform the following calculation. Compute the absolute error and relative error with
More informationMath 12 Final Exam Review 1
Math 12 Final Exam Review 1 Part One Calculators are NOT PERMITTED for this part of the exam. 1. a) The sine of angle θ is 1 What are the 2 possible values of θ in the domain 0 θ 2π? 2 b) Draw these angles
More informationThis practice exam is intended to help you prepare for the final exam for MTH 142 Calculus II.
MTH 142 Practice Exam Chapters 9-11 Calculus II With Analytic Geometry Fall 2011 - University of Rhode Island This practice exam is intended to help you prepare for the final exam for MTH 142 Calculus
More informationPRINCIPLE OF MATHEMATICAL INDUCTION
Chapter 4 PRINCIPLE OF MATHEMATICAL INDUCTION 4.1 Overview Mathematical induction is one of the techniques which can be used to prove variety of mathematical statements which are formulated in terms of
More informationMath 2233 Homework Set 7
Math 33 Homework Set 7 1. Find the general solution to the following differential equations. If initial conditions are specified, also determine the solution satisfying those initial conditions. a y 4
More informationf (r) (a) r! (x a) r, r=0
Part 3.3 Differentiation v1 2018 Taylor Polynomials Definition 3.3.1 Taylor 1715 and Maclaurin 1742) If a is a fixed number, and f is a function whose first n derivatives exist at a then the Taylor polynomial
More informationFunctional Analysis Exercise Class
Functional Analysis Exercise Class Week: December 4 8 Deadline to hand in the homework: your exercise class on week January 5. Exercises with solutions ) Let H, K be Hilbert spaces, and A : H K be a linear
More information2. Second-order Linear Ordinary Differential Equations
Advanced Engineering Mathematics 2. Second-order Linear ODEs 1 2. Second-order Linear Ordinary Differential Equations 2.1 Homogeneous linear ODEs 2.2 Homogeneous linear ODEs with constant coefficients
More informationMa 530 Power Series II
Ma 530 Power Series II Please note that there is material on power series at Visual Calculus. Some of this material was used as part of the presentation of the topics that follow. Operations on Power Series
More informationSolutions to homework assignment #3 Math 119B UC Davis, Spring = 12x. The Euler-Lagrange equation is. 2q (x) = 12x. q(x) = x 3 + x.
1. Find the stationary points of the following functionals: (a) 1 (q (x) + 1xq(x)) dx, q() =, q(1) = Solution. The functional is 1 L(q, q, x) where L(q, q, x) = q + 1xq. We have q = q, q = 1x. The Euler-Lagrange
More informationStudy # 1 11, 15, 19
Goals: 1. Recognize Taylor Series. 2. Recognize the Maclaurin Series. 3. Derive Taylor series and Maclaurin series representations for known functions. Study 11.10 # 1 11, 15, 19 f (n) (c)(x c) n f(c)+
More informationx 2 y = 1 2. Problem 2. Compute the Taylor series (at the base point 0) for the function 1 (1 x) 3.
MATH 8.0 - FINAL EXAM - SOME REVIEW PROBLEMS WITH SOLUTIONS 8.0 Calculus, Fall 207 Professor: Jared Speck Problem. Consider the following curve in the plane: x 2 y = 2. Let a be a number. The portion of
More informationHow to Solve Linear Differential Equations
How to Solve Linear Differential Equations Definition: Euler Base Atom, Euler Solution Atom Independence of Atoms Construction of the General Solution from a List of Distinct Atoms Euler s Theorems Euler
More informationx 2 = 1 Clearly, this equation is not true for all real values of x. Nevertheless, we can solve it by taking careful steps:
Sec. 01 notes Solving Trig Equations: The Easy Ones Main Idea We are now ready to discuss the solving of trigonometric equations. Recall that, generally speaking, identities are equations which hold true
More informationWinter 2014 Practice Final 3/21/14 Student ID
Math 4C Winter 2014 Practice Final 3/21/14 Name (Print): Student ID This exam contains 5 pages (including this cover page) and 20 problems. Check to see if any pages are missing. Enter all requested information
More informationIntroduction Derivation General formula List of series Convergence Applications Test SERIES 4 INU0114/514 (MATHS 1)
MACLAURIN SERIES SERIES 4 INU0114/514 (MATHS 1) Dr Adrian Jannetta MIMA CMath FRAS Maclaurin Series 1/ 21 Adrian Jannetta Recap: Binomial Series Recall that some functions can be rewritten as a power series
More informationInterpolation and Approximation
Interpolation and Approximation The Basic Problem: Approximate a continuous function f(x), by a polynomial p(x), over [a, b]. f(x) may only be known in tabular form. f(x) may be expensive to compute. Definition:
More information2005 Mathematics. Advanced Higher. Finalised Marking Instructions
2005 Mathematics Advanced Higher Finalised Marking Instructions These Marking Instructions have been prepared by Examination Teams for use by SQA Appointed Markers when marking External Course Assessments.
More informationHermite Interpolation
Jim Lambers MAT 77 Fall Semester 010-11 Lecture Notes These notes correspond to Sections 4 and 5 in the text Hermite Interpolation Suppose that the interpolation points are perturbed so that two neighboring
More informationChapter 4: More Applications of Differentiation
Chapter 4: More Applications of Differentiation Autumn 2017 Department of Mathematics Hong Kong Baptist University 1 / 68 In the fall of 1972, President Nixon announced that, the rate of increase of inflation
More informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
More informationAdvanced Calculus Math 127B, Winter 2005 Solutions: Final. nx2 1 + n 2 x, g n(x) = n2 x
. Define f n, g n : [, ] R by f n (x) = Advanced Calculus Math 27B, Winter 25 Solutions: Final nx2 + n 2 x, g n(x) = n2 x 2 + n 2 x. 2 Show that the sequences (f n ), (g n ) converge pointwise on [, ],
More informationThe Derivative. Appendix B. B.1 The Derivative of f. Mappings from IR to IR
Appendix B The Derivative B.1 The Derivative of f In this chapter, we give a short summary of the derivative. Specifically, we want to compare/contrast how the derivative appears for functions whose domain
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationBuckingham s magical pi theorem
Buckingham s magical pi theorem and the Lie symmetries of nature Harald Hanche-Olsen Theoretical physics colloquium 2002 11 05 p.1/31 A simple example θ λ The period t depends on λ, g, θ max. But how?
More information13 Maximum Modulus Principle
3 Maximum Modulus Principle Theorem 3. (maximum modulus principle). If f is non-constant and analytic on an open connected set Ω, then there is no point z 0 Ω such that f(z) f(z 0 ) for all z Ω. Remark
More informationdx n a 1(x) dy
HIGHER ORDER DIFFERENTIAL EQUATIONS Theory of linear equations Initial-value and boundary-value problem nth-order initial value problem is Solve: a n (x) dn y dx n + a n 1(x) dn 1 y dx n 1 +... + a 1(x)
More informationTrigonometry LESSON SIX - Trigonometric Identities I Lesson Notes
LESSON SIX - Trigonometric Identities I Example Understanding Trigonometric Identities. a) Why are trigonometric identities considered to be a special type of trigonometric equation? Trigonometric Identities
More informationLecture Notes for Math 251: ODE and PDE. Lecture 12: 3.3 Complex Roots of the Characteristic Equation
Lecture Notes for Math 21: ODE and PDE. Lecture 12: 3.3 Complex Roots of the Characteristic Equation Shawn D. Ryan Spring 2012 1 Complex Roots of the Characteristic Equation Last Time: We considered the
More informationSection x7 +
Difference Equations to Differential Equations Section 5. Polynomial Approximations In Chapter 3 we discussed the problem of finding the affine function which best approximates a given function about some
More informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
More informationMultiple Choice Answers. MA 114 Calculus II Spring 2013 Final Exam 1 May Question
MA 114 Calculus II Spring 2013 Final Exam 1 May 2013 Name: Section: Last 4 digits of student ID #: This exam has six multiple choice questions (six points each) and five free response questions with points
More informationWeek #6 - Taylor Series, Derivatives and Graphs Section 10.1
Week #6 - Taylor Series, Derivatives and Graphs Section 10.1 From Calculus, Single Variable by Hughes-Hallett, Gleason, McCallum et. al. Copyright 005 by John Wiley & Sons, Inc. This material is used by
More information