P-stable Exponentially fitted Obrechkoff methods
|
|
- Crystal Mason
- 5 years ago
- Views:
Transcription
1 P-stable Exponentially fitted Obrechkoff methods Marnix Van Daele, G. Vanden Berghe Vakgroep Toegepaste Wiskunde en Informatica Universiteit Gent SciCADE 7 p. /4
2 Outline Introduction on exponentially fitted (EF methods 2-step Obrechkoff methods for y = f(x, y P-stable Obrechkoff methods for y = f(x, y P-stable EF Obrechkoff methods for y = f(x, y Conclusions SciCADE 7 p. 2/4
3 Exponentially fitted methods In the past 5 years, our research group has constructed modified versions of well-known linear multistep methods Runge-Kutta methods Aim : build methods which perform very good when the solution has a known exponential of trigonometric behaviour. SciCADE 7 p. 3/4
4 Linear multistep methods A well known method to solve y = f(y y(a = y a y (a = y a is the Numerov method (order 4 y n+ 2 y n + y n = h2 2 (f(y n + f(y n + f(y n+ Construction : impose L[z(x;h] = for z(x =, x, x 2, x 3, x 4 where L[z(x;h] := z(x + h + α z(x + α z(x h h 2 ( β z (x + h + β z (x + β z (x h SciCADE 7 p. 4/4
5 Exponential fitting Consider the initial value problem y + ω 2 y = g(y y(a = y a y(a = y a. If g(y ω 2 y then y(x α cos(ω x + φ To mimic this oscillatory behaviour, one could replace polynomials by trigonometric (in the complex case : exponential functions. SciCADE 7 p. 5/4
6 EF Numerov method Construction : impose L[z(x; h] = for z(x =, x, x 2, sin(ω x, cos(ω x L[z(x;h] := z(x + h + α z(x + α z(x h h 2 ( β z (x + h + β z (x + β z (x h y n+ 2 y n + y n = h 2 (λf(y n + ( 2 λ f(y n + λf(y n+ λ = 4 sin = := ω h SciCADE 7 p. 6/4
7 EF methods Generalisation : to determine the coefficients of a method, we impose conditions on a linear functional. These conditions are related to polynomials : {x q q =,..., K} exponential or trigonometric functions, multiplied with powers of x : {x q exp(±µx q =,..., P } or, with ω = i µ, {x q cos(ωx, x q sin(ωx q =,..., P } Classical method : P = SciCADE 7 p. 7/4
8 Choice of ω local optimization based on local truncation error (lte ω is step-dependent global optimization Preservation of geometric properties (periodicity, energy,... ω is constant over the interval of integration SciCADE 7 p. 8/4
9 Obrechkoff methods Nikola Obrechkoff ( Obrechkoff methods (OM : 94 for quadrature Milne : OM for solving diff. eq. : 949 SciCADE 7 p. 9/4
10 Obrechkoff methods for y = f(x, y y n+ 2y n + y n = Two-step methods m i= ( h 2 i β i y (2i n+ + 2β i y n (2i + β i y (2i n L[z(x;h] := z(x + h 2 z(x + z(x h m ( h 2 i β i z (2i (x + h + 2β i z (2i (x + β i z (2i (x h i= symmetric method : L[z(x; h] if z(x is odd L[;h] order p L[x q ;h] =, q =,,..., p + lte = C p+2 h p+2 y (p+2 (x n + O ( h p+3 C p+2 = L[xp+2 ;h] (p + 2!h p+2 SciCADE 7 p. /4
11 Two-step OM m = : p = 4, C 6 = 24 m = 2 : p = 8, C = y n+ 2y n + y n = h 2 2 (Numerov method ( y (2 n+ + y(2 n + y (2 n y n+ 2y n + y n = h h4 52 ( y (2 n+ + 23y(2 n + y (2 n ( 3 y (4 n+ 626 y(4 n + 3y (4 n SciCADE 7 p. /4
12 Two-step OM m = 3 : p = 2, C 4 = y n+ 2y n + y n = h h h ( 229 y (2 n y(2 n + 229y (2 n ( y (4 n+ method of order p = 4m 422 y(4 n + y (4 n ( 27 y (6 n y(6 n + 27y (6 n SciCADE 7 p. 2/4
13 Stability of two-step OM y n+ 2 y n + y n = m i= ( h 2 i β i y (2i n+ + 2β i y n (2i + β i y (2i n applied to y = λ 2 y gives y n+ 2R mm (ν 2 y n + y n = ν := λh m + ( i β i ν 2i R mm (ν 2 = + i= m ( i+ β i ν 2i i= The stabilitity function R mm uniquely determines the method. A method has the interval of periodicity (, ν 2 if R mm (ν 2 < for < ν 2 < ν 2. The method is P -stable if R mm (ν 2 for all real ν. SciCADE 7 p. 3/4
14 Stability of two-step OM m = : p = 4, C 6 = 24 (Numerov ν2 = 6 y n+ 2y n + y n = h 2 2 ( y (2 n+ + y(2 n + y (2 n R 6 ν 2 SciCADE 7 p. 4/4
15 Stability of two-step OM m = 2 : p = 8, C = y n+ 2y n + y n = h h ( y (2 n+ + 5y(2 n + y (2 n ( 3 y (4 n+ 626 y(4 n + 3y (4 n+ ν 2 = 25.2 R ν 2 SciCADE 7 p. 5/4
16 Stability of two-step OM m = 3 : p = 2, C 4 = y n+ 2y n + y n = h h h ( 229 y (2 n y(2 n + 229y (2 n ( y (4 n y(4 n y (4 n ν 2 = ( 27 y (6 n y(6 n + 27y (6 n R ν 2 SciCADE 7 p. 6/4
17 P-stable 2-step OM Ananthakrishnaiah (987 idea : use some parameters to obtain P-stability. y n+ 2y n + y n = m = 3 and m = 4 e.g. m = 3 h 2 ( β y (2 n+ + 2β y (2 n + β y (2 n +h 4 ( β 2 y (4 n+ + 2β 2 y (4 n + β 2 y (4 n +h 6 ( β 3 y (6 n+ + 2β 3 y (6 n + β 3 y (6 n impose p = 6 : {x 2, x 4, x 6 }, then 3 parameters (β 2, β 3 and β 3 remain R 33 = β ν 2 + β 2 ν 4 β 3 ν 6 + β ν 2 β 2 ν 4 + β 3 ν 6 let β 3 = β 3 SciCADE 7 p. 7/4
18 Ananthakrishnaiah s idea Choose these 2 remaining parameters β 2, β 3 such that the method R 33 = N 3 D 3 becomes P-stable < (D 3 N 3 (D 3 + N 3 >.4 β 2.8. β 3 SciCADE 7 p. 8/4
19 Ananthakrishnaiah s idea Find the solution for which the phase-lag, ( 3 ν arccos R 33 = β β 2 ν 7 + O ( ν 9.4 becomes minimal. β 2.8. β 3 SciCADE 7 p. 9/4
20 Ananthakrishnaiah s idea Find the solution for which the phase-lag, ( 3 ν arccos R 33 = β β 2 ν 7 + O ( ν 9.4 becomes minimal. β 2.8. This gives (β 3, β 2 = β 3 ( 44, 6. SciCADE 7 p. 2/4
21 Ananthakrishnaiah s idea Ananthakrishnaiah m = 3 : p = 6, P-stable, C 8 = 54 y n+ 2y n + y n = h 2 2 h4 6 + h6 44 ( y (2 n+ + 8y(2 n + y (2 n ( y (4 n+ 22 y(4 n + y (4 n ( y (6 n+ + 2y(6 n + y (6 n SciCADE 7 p. 2/4
22 Ananthakrishnaiah s idea R 33 = 9 2 ν2 + 6 ν4 44 ν6 + 2 ν2 + 6 ν ν6 R R 33 = N 3 D ν < since (D 3 N 3 (D 3 + N 3 = ν2 (ν 2 2 (ν SciCADE 7 p. 22/4
23 P-stable 2-step OM We were able to generalise Ananthakrishnaiah s idea in M. Van Daele and G. Vanden Berghe, P-stable Obrechkoff methods of arbitrary order for second-order differential equations, Numerical Algorithms 44, 27, 5-3 Algoritm to construct a P-stable OM for a given m : impose order 2 m D m N m and D m + N m should both be halves of perfect squares This leads to a system of non-linear equations. Theorem : the approximant R mm (ν 2, obtained by generalising Ananthakrishnaiah s approach, is given by the real part of the (m, m-padé approximant of exp(i ν. This leads to a system of linear equations. SciCADE 7 p. 23/4
24 Example R 33 (ν 2 = 9 2 ν2 + 6 ν4 44 ν6 + 2 ν2 + 6 ν ( ν6 + 2 = R i ν ν2 2 i ν3 2 i ν = ( ν2 + 2 i ν3 ν2 2 ( 2 ν 2 ν3 2 ( ν2 2 ( + 2 ν 2 ν3 2 where + 2 ν + ν2 + 2 ν3 2 ν + ν2 = exp(ν + O ( ν 7 2 ν3 ( R 33 (ν 2 = ν2 2 ( 2 ν 2 ν3 2 ( ν2 2 + ( 2 ν 2 ν3 2 SciCADE 7 p. 24/4
25 Conclusion How does the P-stable Obrechkoff method m ( y n+ 2y n + y n = h 2 i β i y (2i n+ + 2β i y n (2i + β i y (2i n i= of order 2 m look like? i β i = ( i+ a 2 i + 2 ( j+ a j a 2i j β i = a 2 i + 2 i j= where a j = j= a j a 2i j ( m j ( 2 m j for j m for j > m i =..., m SciCADE 7 p. 25/4
26 P-stable Expon. fitted OM How to obtain P-stable exponentially-fitted Obrechkoff methods? m ( y n+ 2 y n + y n = h 2 i β i y (2i n+ + 2β i y n (2i + β i y (2i n i= applied to y = λ 2 y gives y n+ 2R mm (, ν 2 y n + y n = with := ω h and ν := λh Padé-approximants = exponentially-fitted Padé approximants Construction of Polynomially-fitted (m, m Padé approximants : P m (x P m ( x = exp(x + O ( x 2 m+ exp(x P m ( x P m (x = O ( x 2 m+ d2 q dx 2 q (exp(x P m( x P m (x x= = q =,..., m SciCADE 7 p. 26/4
27 EF Padé approximants F(x, t = exp(t x V m ( t x V m (t x V m (x = + 2q x 2q F(x, t (x,t=(, = q =,..., K ( R q t q F(x, t = q =,..., P (x,t=(i, where K m and P + K + = m. m a i x i j= This leads to a system of m linear equations in the unknowns a i, i =,..., m. The EF (K, P Padé approximant to exp(ν is then given by (K,PˆP m m (ν = V m (ν/v m ( ν. SciCADE 7 p. 27/4
28 EF Padé approximants : m = F(x, t = exp(t x V ( t x V (t x V m (x = + a x (K =, P = 2 x 2 F(x, t (x,t=(, = 2 ( 2 a = a = 2 (, ˆP (ν = + 2 x (K =, P = 2 x ( R (F(i, = R e i ( i a ( + i a = a = sin (cos + = 2 tan(/2 /2 (,ˆP (ν = tan(/2 /2 x tan(/2 /2 x SciCADE 7 p. 28/4
29 EF Padé approximants : m = m = (K, P a (, (, 2 tan 2 SciCADE 7 p. 29/4
30 EF Padé approximants : m = m = (K, P a (, (, O(4 SciCADE 7 p. 3/4
31 EF Padé approximants : m = m = (K, P a (, a.5.5 (, a 2π 4π 6π 8π π.5.5 2π 4π 6π 8π π SciCADE 7 p. 3/4
32 EF Padé approximants : m = 2 m = 2 (K, P a a 2 (2, (, (, ( cos ( + sin 2 2 tan tan 2 sin ( + sin 2 SciCADE 7 p. 32/4
33 EF Padé approximants : m = 2 m = 2 (K, P a a 2 (2, (, (, O( O ( O(4 SciCADE 7 p. 33/4
34 EF Padé approximants : m = 2 m = 2 (K, P a a 2 (2, a.5 a π 4π 6π 8π π. 2π 4π 6π 8π π (, a.5 a (, a 2π 4π 6π 8π π π 4π 6π 8π π a 2 2π 4π 6π 8π π π 4π 6π 8π π SciCADE 7 p. 34/4
35 m = 3 (K, P a a 2 a 3 (3, a.5 a 2..5 a π 4π 6π 8π π.5. 2π 4π 6π 8π π π 4π 6π 8π π.5 (2, a.5 a 2..5 a (, a 2π 4π 6π 8π π.5 a π 4π 6π 8π π a 3.5 2π 4π 6π 8π π π 4π 6π 8π π 2π 4π 6π 8π π.5 2π 4π 6π 8π π.5 (,2 a.5 a 2..5 a π 4π 6π 8π π 2π 4π 6π 8π π.5 2π 4π 6π 8π π SciCADE 7 p. 35/4
36 The Stiefel-Bettis problem z + z =. exp(ix z( = z ( =.995i x 4π equivalent real form u + u =. cos(x, u( =, u ( =, v + v =. sin(x, v( =, v ( = z(x = u(x + i v(x u(x = cos(x +.5x sin(x v(x = sin(x.5x cos(x γ(x = u 2 (x + v 2 (x = + (.5x 2 γ(4 π h = 2 i π ω = SciCADE 7 p. 36/4
37 The Stiefel-Bettis problem m = log(abs. err. in γ(4π log(h : P = : P = SciCADE 7 p. 37/4
38 The Stiefel-Bettis problem m = 2 log(abs. err. in γ(4π log(h : P = : P = : P = SciCADE 7 p. 38/4
39 The Stiefel-Bettis problem m = log(abs. err. in γ(4π log(h : P = : P = : P = : P = 2 SciCADE 7 p. 39/4
40 Conclusion two-step P-stable exponentially fitted Obrechkoff methods : for any given m, P-stable EF (K, P-methods of order 2 m exist the construction is based on EF Padé approximants to the exponential function the coefficients depend on a parameter the coefficients are continous functions of iff P is odd numerical example was given to illustrate the theory SciCADE 7 p. 4/4
Deferred correction based on exponentially-fitted mono-implicit Runge-Kutta methods
Deferred correction based on exponentially-fitted mono-implicit Runge-Kutta methods M. Van Daele Department of Applied Mathematics, Computer Science and Statistics Ghent University 3th International Conference
More informationFourth-order symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: a review
Fourth-order symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: a review G. Vanden Berghe and M. Van Daele Vakgroep Toegepaste Wiskunde en Informatica, Universiteit Gent, Krijgslaan
More informationTrigonometrically tted predictor corrector methods for IVPs with oscillating solutions
Journal of Computational and Applied Mathematics 58 (23) 35 44 www.elsevier.com/locate/cam Trigonometrically tted predictor corrector methods for IVPs with oscillating solutions G.Psihoyios a, T.E. Simos
More information1 Arithmetic calculations (calculator is not allowed)
1 ARITHMETIC CALCULATIONS (CALCULATOR IS NOT ALLOWED) 1 Arithmetic calculations (calculator is not allowed) 1.1 Check the result Problem 1.1. Problem 1.2. Problem 1.3. Problem 1.4. 78 5 6 + 24 3 4 99 1
More informationSymplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: revisited
Symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: revisited G. Vanden Berghe and M. Van Daele Vakgroep Toegepaste Wiskunde en Informatica, Universiteit Gent, Krijgslaan 8
More informationComplex numbers. Reference note to lecture 9 ECON 5101 Time Series Econometrics
Complex numbers. Reference note to lecture 9 ECON 5101 Time Series Econometrics Ragnar Nymoen March 31 2014 1 Introduction This note gives the definitions and theorems related to complex numbers that are
More informationA New Embedded Phase-Fitted Modified Runge-Kutta Method for the Numerical Solution of Oscillatory Problems
Applied Mathematical Sciences, Vol. 1, 16, no. 44, 157-178 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/ams.16.64146 A New Embedded Phase-Fitted Modified Runge-Kutta Method for the Numerical Solution
More informationTest #2 Math 2250 Summer 2003
Test #2 Math 225 Summer 23 Name: Score: There are six problems on the front and back of the pages. Each subpart is worth 5 points. Show all of your work where appropriate for full credit. ) Show the following
More informationAn Optimized Symmetric 8-Step Semi-Embedded Predictor-Corrector Method for IVPs with Oscillating Solutions
Appl. Math. Inf. Sci., No., - () Applied Mathematics & Information Sciences An International Journal An Optimized Symmetric -Step Semi-Embedded Predictor-Corrector Method for IVPs with Oscillating Solutions
More informationfunction independent dependent domain range graph of the function The Vertical Line Test
Functions A quantity y is a function of another quantity x if there is some rule (an algebraic equation, a graph, a table, or as an English description) by which a unique value is assigned to y by a corresponding
More informationSymplectic Runge-Kutta methods satisfying effective order conditions
/30 Symplectic Runge-Kutta methods satisfying effective order conditions Gulshad Imran John Butcher (supervisor) University of Auckland 5 July, 20 International Conference on Scientific Computation and
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 informationFirst and Second Order ODEs
Civil Engineering 2 Mathematics Autumn 211 M. Ottobre First and Second Order ODEs Warning: all the handouts that I will provide during the course are in no way exhaustive, they are just short recaps. Notation
More informationNotes for the Physics-Based Calculus workshop
Notes for the hysics-based Calculus workshop Adam Coffman June 9, 24 Trigonometric Functions We recall the following identities for trigonometric functions. Theorem.. For all x R, cos( x = cos(x and sin(
More information1. Use the properties of exponents to simplify the following expression, writing your answer with only positive exponents.
Math120 - Precalculus. Final Review. Fall, 2011 Prepared by Dr. P. Babaali 1 Algebra 1. Use the properties of exponents to simplify the following expression, writing your answer with only positive exponents.
More informationThe accurate numerical solution of the Schrödinger equation with an explicitly time-dependent Hamiltonian
The accurate numerical solution of the Schrödinger equation with an explicitly time-dependent Hamiltonian Marnix Van Daele Department of Applied Mathematics, Computer Science and Statistics Ghent University
More informationSolutions 2: Simple Harmonic Oscillator and General Oscillations
Massachusetts Institute of Technology MITES 2017 Physics III Solutions 2: Simple Harmonic Oscillator and General Oscillations Due Wednesday June 21, at 9AM under Rene García s door Preface: This problem
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 informationDenition and some Properties of Generalized Elementary Functions of a Real Variable
Denition and some Properties of Generalized Elementary Functions of a Real Variable I. Introduction The term elementary function is very often mentioned in many math classes and in books, e.g. Calculus
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 information1 + lim. n n+1. f(x) = x + 1, x 1. and we check that f is increasing, instead. Using the quotient rule, we easily find that. 1 (x + 1) 1 x (x + 1) 2 =
Chapter 5 Sequences and series 5. Sequences Definition 5. (Sequence). A sequence is a function which is defined on the set N of natural numbers. Since such a function is uniquely determined by its values
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 informationTwo Optimized Runge-Kutta Methods for the Solution of the Schrödinger Equation
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. (8) 7-77 ISSN - Two Optimized Runge-Kutta Methods for the Solution of the Schrödinger Equation T.V. Triantafyllidis,
More informationSecond-Order Homogeneous Linear Equations with Constant Coefficients
15 Second-Order Homogeneous Linear Equations with Constant Coefficients A very important class of second-order homogeneous linear equations consists of those with constant coefficients; that is, those
More informationAn Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations
An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical
More informationHarmonic oscillator. U(x) = 1 2 bx2
Harmonic oscillator The harmonic oscillator is a familiar problem from classical mechanics. The situation is described by a force which depends linearly on distance as happens with the restoring force
More informationFind all of the real numbers x that satisfy the algebraic equation:
Appendix C: Factoring Algebraic Expressions Factoring algebraic equations is the reverse of expanding algebraic expressions discussed in Appendix B. Factoring algebraic equations can be a great help when
More informationFirst order differential equations
First order differential equations Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T. First
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 information1MA6 Partial Differentiation and Multiple Integrals: I
1MA6/1 1MA6 Partial Differentiation and Multiple Integrals: I Dr D W Murray Michaelmas Term 1994 1. Total differential. (a) State the conditions for the expression P (x, y)dx+q(x, y)dy to be the perfect
More informationMATH 409 Advanced Calculus I Lecture 11: More on continuous functions.
MATH 409 Advanced Calculus I Lecture 11: More on continuous functions. Continuity Definition. Given a set E R, a function f : E R, and a point c E, the function f is continuous at c if for any ε > 0 there
More information(c) Find the equation of the degree 3 polynomial that has the same y-value, slope, curvature, and third derivative as ln(x + 1) at x = 0.
Chapter 7 Challenge problems Example. (a) Find the equation of the tangent line for ln(x + ) at x = 0. (b) Find the equation of the parabola that is tangent to ln(x + ) at x = 0 (i.e. the parabola has
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 informationStep 1: Greatest Common Factor Step 2: Count the number of terms If there are: 2 Terms: Difference of 2 Perfect Squares ( + )( - )
Review for Algebra 2 CC Radicals: r x p 1 r x p p r = x p r = x Imaginary Numbers: i = 1 Polynomials (to Solve) Try Factoring: i 2 = 1 Step 1: Greatest Common Factor Step 2: Count the number of terms If
More informationOn second order sufficient optimality conditions for quasilinear elliptic boundary control problems
On second order sufficient optimality conditions for quasilinear elliptic boundary control problems Vili Dhamo Technische Universität Berlin Joint work with Eduardo Casas Workshop on PDE Constrained Optimization
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 information1. Use the properties of exponents to simplify the following expression, writing your answer with only positive exponents.
Math120 - Precalculus. Final Review Prepared by Dr. P. Babaali 1 Algebra 1. Use the properties of exponents to simplify the following expression, writing your answer with only positive exponents. (a) 5
More informationReview session Midterm 1
AS.110.109: Calculus II (Eng) Review session Midterm 1 Yi Wang, Johns Hopkins University Fall 2018 7.1: Integration by parts Basic integration method: u-sub, integration table Integration By Parts formula
More informationHomogeneous Linear ODEs of Second Order Notes for Math 331
Homogeneous Linear ODEs of Second Order Notes for Math 331 Richard S. Ellis Department of Mathematics and Statistics University of Massachusetts Amherst, MA 01003 In this document we consider homogeneous
More informationChapter 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 informationFLAP M6.1 Introducing differential equations COPYRIGHT 1998 THE OPEN UNIVERSITY S570 V1.1
F1 (a) This equation is of first order and of first degree. It is a linear equation, since the dependent variable y and its derivative only appear raised to the first power, and there are no products of
More informationMath Practice Exam 3 - solutions
Math 181 - Practice Exam 3 - solutions Problem 1 Consider the function h(x) = (9x 2 33x 25)e 3x+1. a) Find h (x). b) Find all values of x where h (x) is zero ( critical values ). c) Using the sign pattern
More informationExponentially Fitted Symplectic Runge-Kutta-Nyström methods
Appl. Math. Inf. Sci. 7, No. 1, 81-85 (2013) 81 Applied Mathematics & Information Sciences An International Journal Exponentially Fitted Symplectic Runge-Kutta-Nyström methods Th. Monovasilis 1, Z. Kalogiratou
More informationPolytechnic Institute of NYU MA 2132 Final Practice Answers Fall 2012
Polytechnic Institute of NYU MA Final Practice Answers Fall Studying from past or sample exams is NOT recommended. If you do, it should be only AFTER you know how to do all of the homework and worksheet
More informationResearch Article A Family of Trigonometrically Fitted Enright Second Derivative Methods for Stiff and Oscillatory Initial Value Problems
Journal of Applied Mathematics Volume 2015, Article ID 343295, 17 pages http://dx.doi.org/10.1155/2015/343295 Research Article A Family of Trigonometrically Fitted Enright Second Derivative Methods for
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 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 informationMath 106: Review for Exam II - SOLUTIONS
Math 6: Review for Exam II - SOLUTIONS INTEGRATION TIPS Substitution: usually let u a function that s inside another function, especially if du (possibly off by a multiplying constant) is also present
More information6.1 The Inverse Sine, Cosine, and Tangent Functions Objectives
Objectives 1. Find the Exact Value of an Inverse Sine, Cosine, or Tangent Function. 2. Find an Approximate Value of an Inverse Sine Function. 3. Use Properties of Inverse Functions to Find Exact Values
More informationSolution of Fourth Order Obstacle Problems Using Quintic B-Splines
Applied Mathematical Sciences, Vol. 6, 202, no. 94, 465-4662 Solution of Fourth Order Obstacle Problems Using Quintic B-Splines Shahid S. Siddiqi, Ghazala Akram and Kalsoom Arshad Department of Mathematics
More informationECM Calculus and Geometry. Revision Notes
ECM1702 - Calculus and Geometry Revision Notes Joshua Byrne Autumn 2011 Contents 1 The Real Numbers 1 1.1 Notation.................................................. 1 1.2 Set Notation...............................................
More informationMath 106: Review for Exam II - SOLUTIONS
Math 6: Review for Exam II - SOLUTIONS INTEGRATION TIPS Substitution: usually let u a function that s inside another function, especially if du (possibly off by a multiplying constant) is also present
More informationFinal Exam Solutions
Final Exam Solutions Laurence Field Math, Section March, Name: Solutions Instructions: This exam has 8 questions for a total of points. The value of each part of each question is stated. The time allowed
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More informationCalculus I Announcements
Slide 1 Calculus I Announcements Read sections 3.9-3.10 Do all the homework for section 3.9 and problems 1,3,5,7 from section 3.10. The exam is in Thursday, October 22nd. The exam will cover sections 3.2-3.10,
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More informationMATH 150 TOPIC 16 TRIGONOMETRIC EQUATIONS
Math 150 T16-Trigonometric Equations Page 1 MATH 150 TOPIC 16 TRIGONOMETRIC EQUATIONS In calculus, you will often have to find the zeros or x-intercepts of a function. That is, you will have to determine
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 information17 Source Problems for Heat and Wave IB- VPs
17 Source Problems for Heat and Wave IB- VPs We have mostly dealt with homogeneous equations, homogeneous b.c.s in this course so far. Recall that if we have non-homogeneous b.c.s, then we want to first
More informationEXAMPLE OF ONE-STEP METHOD
EXAMPLE OF ONE-STEP METHOD Consider solving y = y cos x, y(0) = 1 Imagine writing a Taylor series for the solution Y (x), say initially about x = 0. Then Y (h) = Y (0) + hy (0) + h2 2 Y (0) + h3 6 Y (0)
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Matematics National Institute of Tecnology Durgapur Durgapur-7109 email: anita.buie@gmail.com 1 . Capter 8 Numerical Solution of Ordinary
More informationSome functions and their derivatives
Chapter Some functions an their erivatives. Derivative of x n for integer n Recall, from eqn (.6), for y = f (x), Also recall that, for integer n, Hence, if y = x n then y x = lim δx 0 (a + b) n = a n
More informationResearch Article P-Stable Higher Derivative Methods with Minimal Phase-Lag for Solving Second Order Differential Equations
Hindawi Publishing Corporation Journal of Applied Mathematics Volume 2011, Article ID 407151, 15 pages doi:10.1155/2011/407151 Research Article P-Stable Higher Derivative Methods with Minimal Phase-Lag
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationProblem Set 1 October 30, 2017
1. e π can be calculated from e x = x n. But that s not a good approximation method. The n! reason is that π is not small compared to 1. If you want O(0.1) accuracy, then the last term you need to include
More informationDifferential Equations
Differential Equations Definitions Finite Differences Taylor Series based Methods: Euler Method Runge-Kutta Methods Improved Euler, Midpoint methods Runge Kutta (2nd, 4th order) methods Predictor-Corrector
More informationHighly Oscillatory Problems: Computation, Theory and Applications
Highly Oscillatory Problems: Computation, Theory and Applications Isaac Newton Institute for Mathematical Sciences Edited by 1 Rapid function approximation by modified Fourier series Daan Huybrechs and
More informationAP Calculus Summer Prep
AP Calculus Summer Prep Topics from Algebra and Pre-Calculus (Solutions are on the Answer Key on the Last Pages) The purpose of this packet is to give you a review of basic skills. You are asked to have
More informationPHYS 404 Lecture 1: Legendre Functions
PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions
More informationAnswers to Problem Set Number MIT (Fall 2005).
Answers to Problem Set Number 5. 18.305 MIT Fall 2005). D. Margetis and R. Rosales MIT, Math. Dept., Cambridge, MA 02139). November 23, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,
More informationComputational Physics (6810): Session 8
Computational Physics (6810): Session 8 Dick Furnstahl Nuclear Theory Group OSU Physics Department February 24, 2014 Differential equation solving Session 7 Preview Session 8 Stuff Solving differential
More informationPartial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators
PDE-SEP-HEAT-1 Partial Differential Equations Separation of Variables 1 Partial Differential Equations and Operators et C = C(R 2 ) be the collection of infinitely differentiable functions from the plane
More information1) [16 points] The light intensity of a particular light source (not an ideal blackbody) is given by the equation. p = - iħ d (2.
CHM 3411 First hour exam February 8, 2019 There are 5 problems on the exam. Do all of the problems. Show your work. 1) [16 points] The light intensity of a particular light source (not an ideal blackbody)
More informationWeek 5: Functions and graphs
Calculus and Linear Algebra for Biomedical Engineering Week 5: Functions and graphs H. Führ, Lehrstuhl A für Mathematik, RWTH Aachen, WS 07 Motivation: Measurements at fixed intervals 1 Consider a sequence
More informationSolved Examples. (Highest power of x in numerator and denominator is ½. Dividing numerator and denominator by x)
Solved Examples Example 1: (i) (ii) lim x (x 4 + 2x 3 +3) / (2x 4 -x+2) lim x x ( (x+c)- x) (iii) lim n (1-2+3-4+...(2n-1)-2n)/ (n 2 +1) (iv) lim x 0 ((1+x) 5-1)/3x+5x 2 (v) lim x 2 ( (x+7)-3 (2x-3))/((x+6)
More informationTopics from Algebra and Pre-Calculus. (Key contains solved problems)
Topics from Algebra and Pre-Calculus (Key contains solved problems) Note: The purpose of this packet is to give you a review of basic skills. You are asked not to use the calculator, except on p. (8) and
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 informationFunctions. Remark 1.2 The objective of our course Calculus is to study functions.
Functions 1.1 Functions and their Graphs Definition 1.1 A function f is a rule assigning a number to each of the numbers. The number assigned to the number x via the rule f is usually denoted by f(x).
More informationAdebayo O. Adeniran1, Saheed O. Akindeinde2 and Babatunde S. Ogundare2 * Contents. 1. Introduction
Malaya Journal of Matematik Vol. No. 73-73 8 https://doi.org/.37/mjm/ An accurate five-step trigonometrically-fitted numerical scheme for approximating solutions of second order ordinary differential equations
More informationEigenvalue Problem. 1 The First (Dirichlet) Eigenvalue Problem
Eigenvalue Problem A. Salih Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram July 06 The method of separation variables for solving the heat 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 informationPhasor Young Won Lim 05/19/2015
Phasor Copyright (c) 2009-2015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version
More informationMathematical Economics: Lecture 2
Mathematical Economics: Lecture 2 Yu Ren WISE, Xiamen University September 25, 2012 Outline 1 Number Line The number line, origin (Figure 2.1 Page 11) Number Line Interval (a, b) = {x R 1 : a < x < b}
More informationExponential Integrators
Exponential Integrators John C. Bowman (University of Alberta) May 22, 2007 www.math.ualberta.ca/ bowman/talks 1 Exponential Integrators Outline Exponential Euler History Generalizations Stationary Green
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 informationDr. Sophie Marques. MAM1020S Tutorial 8 August Divide. 1. 6x 2 + x 15 by 3x + 5. Solution: Do a long division show your work.
Dr. Sophie Marques MAM100S Tutorial 8 August 017 1. Divide 1. 6x + x 15 by 3x + 5. 6x + x 15 = (x 3)(3x + 5) + 0. 1a 4 17a 3 + 9a + 7a 6 by 3a 1a 4 17a 3 + 9a + 7a 6 = (4a 3 3a + a + 3)(3a ) + 0 3. 1a
More informationAn Analysis of Five Numerical Methods for Approximating Certain Hypergeometric Functions in Domains within their Radii of Convergence
An Analysis of Five Numerical Methods for Approximating Certain Hypergeometric Functions in Domains within their Radii of Convergence John Pearson MSc Special Topic Abstract Numerical approximations of
More informationUniversity Calculus I. Worksheet # 8 Mar b. sin tan e. sin 2 sin 1 5. b. tan. c. sec sin 1 ( x )) cos 1 ( x )) f. csc. c.
MATH 6 WINTER 06 University Calculus I Worksheet # 8 Mar. 06-0 The topic covered by this worksheet is: Derivative of Inverse Functions and the Inverse Trigonometric functions. SamplesolutionstoallproblemswillbeavailableonDL,
More informationDepartment of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
Department of Applied Mathematics and Theoretical Physics AMA 204 Numerical analysis Exam Winter 2004 The best six answers will be credited All questions carry equal marks Answer all parts of each question
More informationTaylor series. Chapter Introduction From geometric series to Taylor polynomials
Chapter 2 Taylor series 2. Introduction The topic of this chapter is find approximations of functions in terms of power series, also called Taylor series. Such series can be described informally as infinite
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 informationThe Area bounded by Two Functions
The Area bounded by Two Functions The graph below shows 2 functions f(x) and g(x) that are continuous between x = a and x = b and f(x) g(x). The area shaded in green is the area between the 2 curves. We
More informationLinear Multistep Methods
Linear Multistep Methods Linear Multistep Methods (LMM) A LMM has the form α j x i+j = h β j f i+j, α k = 1 i 0 for the approximate solution of the IVP x = f (t, x), x(a) = x a. We approximate x(t) on
More informationSolutions to Exam 1, Math Solution. Because f(x) is one-to-one, we know the inverse function exists. Recall that (f 1 ) (a) =
Solutions to Exam, Math 56 The function f(x) e x + x 3 + x is one-to-one (there is no need to check this) What is (f ) ( + e )? Solution Because f(x) is one-to-one, we know the inverse function exists
More informationCalculus II Practice Test 1 Problems: , 6.5, Page 1 of 10
Calculus II Practice Test Problems: 6.-6.3, 6.5, 7.-7.3 Page of This is in no way an inclusive set of problems there can be other types of problems on the actual test. To prepare for the test: review homework,
More information1 Separation of Variables
Jim ambers ENERGY 281 Spring Quarter 27-8 ecture 2 Notes 1 Separation of Variables In the previous lecture, we learned how to derive a PDE that describes fluid flow. Now, we will learn a number of analytical
More information= x iy. Also the. DEFINITION #1. If z=x+iy, then the complex conjugate of z is given by magnitude or absolute value of z is z =.
Handout # 4 COMPLEX FUNCTIONS OF A COMPLEX VARIABLE Prof. Moseley Chap. 3 To solve z 2 + 1 = 0 we "invent" the number i with the defining property i 2 = ) 1. We then define the set of complex numbers as
More information(e) 2 (f) 2. (c) + (d). Limits at Infinity. 2.5) 9-14,25-34,41-43,46-47,56-57, (c) (d) 2
Math 150A. Final Review Answers, Spring 2018. Limits. 2.2) 7-10, 21-24, 28-1, 6-8, 4-44. 1. Find the values, or state they do not exist. (a) (b) 1 (c) DNE (d) 1 (e) 2 (f) 2 (g) 2 (h) 4 2. lim f(x) = 2,
More informationA( x) B( x) C( x) y( x) 0, A( x) 0
3.1 Lexicon Revisited The nonhomogeneous nd Order ODE has the form: d y dy A( x) B( x) C( x) y( x) F( x), A( x) dx dx The homogeneous nd Order ODE has the form: d y dy A( x) B( x) C( x) y( x), A( x) dx
More informationSOLUTIONS TO PROBLEMS FROM ASSIGNMENT 5. Problems 3.1:6bd
SOLUTIONS TO PROBLEMS FROM ASSIGNMENT 5 Statement. Solve the problem Problems 3.1:6bd u t = u xx, (t ), with the initial condition u(x, ) = f(x), where the functions u(x, t) and f(x) are assumed to be
More information