n 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
|
|
- Eugenia Karen Barker
- 5 years ago
- Views:
Transcription
1 Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and β gives x n+1 = x n 1f n x n f n 1 x = x (1) f n f n f n x n 1 n n 1 f n f n 1 b Denote the root ξ, ie x k = + k, k =, 1,, Taylor expanding around the root gives f k = f( + k ) = k f ( ) + 1 k f ( ) + Substitution into (1) gives, after simple algebra n 1 n+1 f = 1 n 1 n f + n n+1 f + 1 n n+1 f 1 n 1 n f 1 n 1 n+1 f + Given that the successive errors are getting smaller, only the first term on the right can balance the one on the left Hence n+1 = 1 () f f n n 1 + {smaller terms} We drop the {smaller terms} and set c 1 = 1 Equation () can then be written c 1 n+1 = c 1 n $ c 1 n 1 After taking logs, this becomes f f d n+1 = d n + d n 1 (3) where d k = ln(c k ) The recursion relation (3) is easiest solved by its characteristic equation r r 1 =, with roots r 1, = 1! n The growing solution to (3) is therefore d n = c, and we get ln (c 1 n ) = c 5 +1 n = 5 +1 ln(c 1 n 1 ) = ln ((c 1 n 1 ) ( 5 +1)/ ), and the result follows c Given that the points {x n+1, y n+1, } and {x k, y k, f k }, k = n, n 1, n all lie on the first plane, we have x n+1 + y n+1 + = x n + y n +f n = x n 1 + y n 1 +f n 1 = x n + y n +f n = Subtracting the top equation from the remaining ones gives x n x n+1 x n 1 x n+1 x n x n+1 + y n y n+1 y n 1 y n+1 y n y n+1 + f n f n 1 f n = (4) Similarly for the second plane x n x n+1 x n 1 x n+1 x n x n+1 + y n y n+1 y n 1 y n+1 y n y n+1 + g n g n 1 g n = (5) The relations (4) and (5) show that only two of the four different vectors that appear in their left hand sides can be independent Hence, the two determinants in the hint for the problem must both be zero
2 Numerical Quadrature Consider the singular integral I := 1 r α dr, for < α < 1 a Since the integrand is concave, the right-hand rule satisfies the following bound, I R := N h(jh) α I j=1 where h = 1/N This same rule can be viewed as a left-hand rule for the integral on the interval (h, 1 + h) which yields the bound 1+h h r α dr I R Use this information to develop a bound on the error E R = I I R Find a simple lower bound on the error The two bounds above yield I I R I 1+h h r α dr h r α dr 1+h 1 1 α (h1 α ((1 + h) 1 α 1)) 1 1 α h1 α 1 r α dr For a lower bound on the error, we have the error in the first interval I I R h r α h 1 α = α 1 α h1 α b Use the same technique to develop an upper bound on the error for the midpoint rule I M := N h((j 1 )h) α j=1 Shift to the left by h origin to get and subtract the part of the rule to the left of the I M ( h )1 α I Shift to the right by h to get the bound 1 1+ h h I M
3 Putting this together we have ( h )1 α I I M h 1 1 α (h )1 α Choosing the larger of the two extremes yields 1+ h α ((h )1 α ((1 + h )1 α 1)) I I M 1 1 α (h )1 α c Let f(r) C 4 Find a transformation on the problem that allows the use of standard quadrature rules for the integral 1 r α f(r)dr Write the error bounds for the composite trapezoidal rule on the transformed integral What transformation assures O(h ) convergence Consider the transformation r = s k This yields 1 r α f(r)dr = 1 s ((1 α)k 1) f(s k )ds For k > 1 1 α the integrand, F (s) = s((1 α)k 1) f(s k ) is bounded The composite trapezoidal rule has error bound I I T h 1 F If we choose k = 3 1 α, then F (s) = s f(s k ) This yields F (s) = sf(s k ) + ks 1+k f (s k ) F (s) = f(s k ) + ks k f (s k ) + k(k + 1)s k f (s k ) + k s k f (s k ) Since f C 4, we know that F is bounded
4 Interpolation/approximation a) Consider a function f(x) sampled at x =, 1,, 3 with the values f() =, f(1) = 1, f() = and, f(3) = 1 Find the 3rd degree polynomial that fits the data Solution Using the Lagrange interpolants l i (x) = N j=,j i we can write the interpolating polynomial as (x x j ) (x i x j ) p(x) = = 3 f(x i )l i (x) = l 1 (x) l 3 (x) = i= x(x )(x 3) x(x 1)(x ) 3 = x3 6x + 8x 3 Using B-splines we can piecewise interpolate the data s,, s N given at x =,, N by s(x) = N B (m) (x k)s k (1) k= with the B-splines of order m + 1 given by m+1 ( m + 1 B (m) (x) = ( 1) l l l= ) (x + (m + 1)/ l) m + and x m + is defined as xm for x > and zero for x m! () b) Using the properties of B-splines, relate the smoothness of s(x) in (1) 1
5 to the integer m in the () In other words, for a given m, what s the differentiability of s(x)? Solution Using the smoothness at the knots imposed by spline interpolation, we get the following table: m Spline interpolation Smoothness of s(x) m = constant Discontinuous m = 1 linear Continuous m = quadratic Differentiable m = 3 cubic Twice differentiable m = k s(x) C k 1 ( ) c) Write down the expression and draw a picture of the B-spline used for piecewise linear interpolation of a data set Solution We have that B (1) (x) = ( ( 1) l l l= ) (x + 1 l) + where the first term is given by the second term is given by (x + 1) + = x + = and the third term is given by (x 1) + = {, x x x > 1, {, x x x >, {, x 1 x 1 x > 1
6 Hence,, x 1 B (1) 1 + x, 1 < x (x) = 1 + x x = 1 x, < x x x + x 1 =, x > 1 that is, the tent function = { 1 x, x < 1 x 1, The tent function d) Compute the Fourier transform of a function obtained by piecewise linear spline interpolation of a data set, and show that it decays as O(ω ) Hint: You can use that if g(x) is real and even, then its Fourier transform ( ) ĝ(ω) = g(x)e πiωx dx = Re g(x)e πiωx dx Solution We have that N ŝ(ω) = s k B (1) (x k)e πiωx dx = k= N s k e πiωk B (1) (t)e πiωt dt k= 3
7 Using the hint and the expression for B (1) (x) we have that where B (1) (ω) ( 1 ) B (1) (t)e πiωt dt = Re (1 t)e πiωt dt 1 [ (1 t)e (1 t)e πiωt πiωt dt = πiω The real part of (3) is given by ] 1 = 1 πiω + 1 πiω + 1 πiω [ e πiωt ] 1 1 πiω = 1 πiω 1 4π ω (eπiω 1) e πiωt dt (3) 1 4π ω (cos πω 1) = sin πω π ω and hence B (1) (ω) = sin πω π ω which decays as 1 ω 4
8 Linear algebra: 4 a We modify the system Ax = b to (U & AU)U & x = U & b where U is the DFT (discrete Fourier transform) matrix The matrix U & AU is then diagonal, and it will cost one FFT to obtain its elements Likewise, it will take one FFT to compute U & b and finally one more FFT to extract x from U & x With each FFT costing O(N log N) operations, this will also become the estimate for the total cost b Unpivoted Gaussian elimination for a banded system will produce banded L and U factors at a cost of O(1) operations for each row/column of U and L For each of the two beck substitutions, the cost will again be O(1) operations for each computed entry The total cost will therefore be O(N) operations c We can interchange rows and columns in A so that its structure instead becomes A = * * * * By elementary row operations, we can modify each of the diagonal % blocks to purely diagonal form Using the top N- equations, we can next eliminate all entries in the last two rows - up to their last two columns We then diagonalize the bottom right % block, and then use the last two rows to eliminate all remaining entries of the last two columns We have after that a purely diagonal system to solve Each of the major steps described (diagonalizing the top left (N )%(N ) block, eliminating the last two rows and eliminating the last two columns) takes only O(N) operations Hence, this will also be the overall cost Alternatively, one can note that the system takes the form of a tridiagonal matrix with two extra full rows and columns appended The system can be solved as a tridiagonal system plus use of Woodbury's formula, each requiring O(N) operations
9 ODEs a) Derive the coefficients for the backward differentiation formula c u n+1 + c 1 u n + c u n 1 = f(t n+1, u n+1 ) that result in a nd order method for solving the ODE u = f(t, u) Solution The formula should hold for {u = 1, f = }, {u = t, f = 1}, and {u = t, f = t} which, for t n 1 = k, t n = and t n+1 = k, gives the system c + c 1 + c = k(c c ) = 1 k (c + c ) = k with the solution c = 3 k, c 1 = /k and c = 1 k b) Show that the scheme derived in a) is convergent Solution To establish convergence we need to confirm consistency (OK, since our scheme is nd order accurate), and stability To show stability, we need to check the root condition The characteristic equation is given by 3r r + 1 = with the roots r = 1 and r = 1/3 These roots are within the unit circe, and since there is no double root on the unit circle, the root condition is satisfied and the scheme is stable 1
10 Consider a stiff linear system of ODEs of the form where A is a negative definite matrix c) Explain what is meant by a stiff system du dt = Au, u() = u (1) Solution There is no strict definition, but we can think of stiff ODEs as something like the numerical solution requires a significant reduction in step size to achieve stability, or the ratio of the largest and smallest eigenvalues of the linear system is large, or the system contains both rapidly and slowly decaying modes d) By studying the stability domain of the scheme derived in a), determine if this scheme is suitable for solving (1) Hint: You don t have to derive the entire stability domain It suffices to establish stability for the part of the complex plane that is relevant for the system in (1) Solution Since A is negative definite and therefore has only negative eigenvalues, it suffices to investigate stability for the test equation u = λu where λ > Applying our scheme to this test equation with time step k gives us 3 k u n+1 k u n + u n 1 k = λu n+1
11 We set µ kλ which gives the characteristic equation with the solution r 4r 3 + µ µ = r = ± 1 µ 3 + µ In order to see if r 1, we distinguish two cases: µ 1/: Then the numerator is real and bounded by 3, and hence r 3 3+µ 1 µ > 1/: The numerator is complex valued and r can be written as r = ± i µ µ Then r = 4 + µ 1 (3 + µ) = µ < 1 Hence, since r 1 for the test equation, the scheme is stable for all negative eigenvalues, and therefore the scheme is suitable for the problem 3
12 PDEs Consider the wave equation u tt = cu xx, x [, 1], t >, with initial conditions u(x, ) = f(x), u t (x, ) = g(x), and boundary conditions u(, t) = u(1, t) =, for t > Let x = 1/N and define the mesh (x k, t l ) = (k x, l t) Let u l k approximate u(x k, t l ) Use von Neumann analysis to determine the relationship between t and x for which the following scheme is stable: Justify your answer u l+1 k u l k + u l 1 k t = c ul k+1 u l k + u l k 1 x Solution Using von Neumann analysis we set u l k = a l e ikθ and plug into the difference formula to determine the amplification factor, a This yields ( ) t a + 1/a = c (e iθ + e iθ ) x ( ) t = 4c sin(θ/) x Let s = sin(θ/) and λ = ( ) t x This yields the following equation for a: Solving for a yields: a ( 4cλ s )a + 1 = a = (1 cλ s ) ± (1 cλ s ) 1 i For < cλ s < the roots are complex with modulous 1 ii For = cλ s, a = 1 is the only root while for cλ s = a = 1 is the only root iii For cλ s > there is a root with a > 1 Since we want a 1 for all θ, and if a = 1 we want a simple root, we have the condition: λ < 1 c or t < x c 1
Preliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationCubic Splines MATH 375. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Cubic Splines
Cubic Splines MATH 375 J. Robert Buchanan Department of Mathematics Fall 2006 Introduction Given data {(x 0, f(x 0 )), (x 1, f(x 1 )),...,(x n, f(x n ))} which we wish to interpolate using a polynomial...
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationNumerical Analysis Preliminary Exam 10.00am 1.00pm, January 19, 2018
Numerical Analysis Preliminary Exam 0.00am.00pm, January 9, 208 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More informationx n+1 = x n f(x n) f (x n ), n 0.
1. Nonlinear Equations Given scalar equation, f(x) = 0, (a) Describe I) Newtons Method, II) Secant Method for approximating the solution. (b) State sufficient conditions for Newton and Secant to converge.
More informationCS 257: Numerical Methods
CS 57: Numerical Methods Final Exam Study Guide Version 1.00 Created by Charles Feng http://www.fenguin.net CS 57: Numerical Methods Final Exam Study Guide 1 Contents 1 Introductory Matter 3 1.1 Calculus
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More informationExam 2. Average: 85.6 Median: 87.0 Maximum: Minimum: 55.0 Standard Deviation: Numerical Methods Fall 2011 Lecture 20
Exam 2 Average: 85.6 Median: 87.0 Maximum: 100.0 Minimum: 55.0 Standard Deviation: 10.42 Fall 2011 1 Today s class Multiple Variable Linear Regression Polynomial Interpolation Lagrange Interpolation Newton
More informationSolutions Preliminary Examination in Numerical Analysis January, 2017
Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationDepartment of Applied Mathematics Preliminary Examination in Numerical Analysis August, 2013
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August, 013 August 8, 013 Solutions: 1 Root Finding (a) Let the root be x = α We subtract α from both sides of x n+1 = x
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationEngg. Math. II (Unit-IV) Numerical Analysis
Dr. Satish Shukla of 33 Engg. Math. II (Unit-IV) Numerical Analysis Syllabus. Interpolation and Curve Fitting: Introduction to Interpolation; Calculus of Finite Differences; Finite Difference and Divided
More informationNumerical integration and differentiation. Unit IV. Numerical Integration and Differentiation. Plan of attack. Numerical integration.
Unit IV Numerical Integration and Differentiation Numerical integration and differentiation quadrature classical formulas for equally spaced nodes improper integrals Gaussian quadrature and orthogonal
More informationScientific Computing
2301678 Scientific Computing Chapter 2 Interpolation and Approximation Paisan Nakmahachalasint Paisan.N@chula.ac.th Chapter 2 Interpolation and Approximation p. 1/66 Contents 1. Polynomial interpolation
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More information5 Numerical Integration & Dierentiation
5 Numerical Integration & Dierentiation Department of Mathematics & Statistics ASU Outline of Chapter 5 1 The Trapezoidal and Simpson Rules 2 Error Formulas 3 Gaussian Numerical Integration 4 Numerical
More information1 Piecewise Cubic Interpolation
Piecewise Cubic Interpolation Typically the problem with piecewise linear interpolation is the interpolant is not differentiable as the interpolation points (it has a kinks at every interpolation point)
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 informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationBindel, Spring 2012 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 16. f(x) dx,
Panel integration Week 12: Monday, Apr 16 Suppose we want to compute the integral b a f(x) dx In estimating a derivative, it makes sense to use a locally accurate approximation to the function around the
More informationLösning: Tenta Numerical Analysis för D, L. FMN011,
Lösning: Tenta Numerical Analysis för D, L. FMN011, 090527 This exam starts at 8:00 and ends at 12:00. To get a passing grade for the course you need 35 points in this exam and an accumulated total (this
More informationChapter 2 Interpolation
Chapter 2 Interpolation Experiments usually produce a discrete set of data points (x i, f i ) which represent the value of a function f (x) for a finite set of arguments {x 0...x n }. If additional data
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationGENG2140, S2, 2012 Week 7: Curve fitting
GENG2140, S2, 2012 Week 7: Curve fitting Curve fitting is the process of constructing a curve, or mathematical function, f(x) that has the best fit to a series of data points Involves fitting lines and
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationNumerical interpolation, extrapolation and fi tting of data
Chapter 6 Numerical interpolation, extrapolation and fi tting of data 6.1 Introduction Numerical interpolation and extrapolation is perhaps one of the most used tools in numerical applications to physics.
More information5. Hand in the entire exam booklet and your computer score sheet.
WINTER 2016 MATH*2130 Final Exam Last name: (PRINT) First name: Student #: Instructor: M. R. Garvie 19 April, 2016 INSTRUCTIONS: 1. This is a closed book examination, but a calculator is allowed. The test
More informationNUMERICAL METHODS. x n+1 = 2x n x 2 n. In particular: which of them gives faster convergence, and why? [Work to four decimal places.
NUMERICAL METHODS 1. Rearranging the equation x 3 =.5 gives the iterative formula x n+1 = g(x n ), where g(x) = (2x 2 ) 1. (a) Starting with x = 1, compute the x n up to n = 6, and describe what is happening.
More informationJim Lambers ENERGY 281 Spring Quarter Lecture 5 Notes
Jim ambers ENERGY 28 Spring Quarter 27-8 ecture 5 Notes These notes are based on Rosalind Archer s PE28 lecture notes, with some revisions by Jim ambers. Fourier Series Recall that in ecture 2, when we
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationWe consider the problem of finding a polynomial that interpolates a given set of values:
Chapter 5 Interpolation 5. Polynomial Interpolation We consider the problem of finding a polynomial that interpolates a given set of values: x x 0 x... x n y y 0 y... y n where the x i are all distinct.
More informationLecture 1 INF-MAT3350/ : Some Tridiagonal Matrix Problems
Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems Tom Lyche University of Oslo Norway Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems p.1/33 Plan for the day 1. Notation
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationExtra Problems for Math 2050 Linear Algebra I
Extra Problems for Math 5 Linear Algebra I Find the vector AB and illustrate with a picture if A = (,) and B = (,4) Find B, given A = (,4) and [ AB = A = (,4) and [ AB = 8 If possible, express x = 7 as
More information8.3 Partial Fraction Decomposition
8.3 partial fraction decomposition 575 8.3 Partial Fraction Decomposition Rational functions (polynomials divided by polynomials) and their integrals play important roles in mathematics and applications,
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More information1 Solutions to selected problems
Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p
More informationMODULE 7. where A is an m n real (or complex) matrix. 2) Let K(t, s) be a function of two variables which is continuous on the square [0, 1] [0, 1].
Topics: Linear operators MODULE 7 We are going to discuss functions = mappings = transformations = operators from one vector space V 1 into another vector space V 2. However, we shall restrict our sights
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 informationSPLINE INTERPOLATION
Spline Background SPLINE INTERPOLATION Problem: high degree interpolating polynomials often have extra oscillations. Example: Runge function f(x = 1 1+4x 2, x [ 1, 1]. 1 1/(1+4x 2 and P 8 (x and P 16 (x
More informationLecture Note 3: Interpolation and Polynomial Approximation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Interpolation and Polynomial Approximation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 10, 2015 2 Contents 1.1 Introduction................................ 3 1.1.1
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 informationTABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1. Chapter Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9
TABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1 Chapter 01.01 Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9 Chapter 01.02 Measuring errors 11 True error 11 Relative
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationIntroductory Numerical Analysis
Introductory Numerical Analysis Lecture Notes December 16, 017 Contents 1 Introduction to 1 11 Floating Point Numbers 1 1 Computational Errors 13 Algorithm 3 14 Calculus Review 3 Root Finding 5 1 Bisection
More informationUNIVERSITY OF SOUTHAMPTON. A foreign language dictionary (paper version) is permitted provided it contains no notes, additions or annotations.
UNIVERSITY OF SOUTHAMPTON MATH055W SEMESTER EXAMINATION 03/4 MATHEMATICS FOR ELECTRONIC & ELECTRICAL ENGINEERING Duration: 0 min Solutions Only University approved calculators may be used. A foreign language
More informationMath 10b Ch. 8 Reading 1: Introduction to Taylor Polynomials
Math 10b Ch. 8 Reading 1: Introduction to Taylor Polynomials Introduction: In applications, it often turns out that one cannot solve the differential equations or antiderivatives that show up in the real
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 7 Interpolation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationDifferential Equations
Differential Equations Problem Sheet 1 3 rd November 2011 First-Order Ordinary Differential Equations 1. Find the general solutions of the following separable differential equations. Which equations are
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationON APPROXIMATION OF LAPLACIAN EIGENPROBLEM OVER A REGULAR HEXAGON WITH ZERO BOUNDARY CONDITIONS 1) 1. Introduction
Journal of Computational Mathematics, Vol., No., 4, 75 86. ON APPROXIMATION OF LAPLACIAN EIGENPROBLEM OVER A REGULAR HEXAGON WITH ZERO BOUNDARY CONDITIONS ) Jia-chang Sun (Parallel Computing Division,
More informationCalculus I Review Solutions
Calculus I Review Solutions. Compare and contrast the three Value Theorems of the course. When you would typically use each. The three value theorems are the Intermediate, Mean and Extreme value theorems.
More informationLectures 9-10: Polynomial and piecewise polynomial interpolation
Lectures 9-1: Polynomial and piecewise polynomial interpolation Let f be a function, which is only known at the nodes x 1, x,, x n, ie, all we know about the function f are its values y j = f(x j ), j
More informationQuaternion Cubic Spline
Quaternion Cubic Spline James McEnnan jmcennan@mailaps.org May 28, 23 1. INTRODUCTION A quaternion spline is an interpolation which matches quaternion values at specified times such that the quaternion
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationMath 353 Lecture Notes Week 6 Laplace Transform: Fundamentals
Math 353 Lecture Notes Week 6 Laplace Transform: Fundamentals J. Wong (Fall 217) October 7, 217 What did we cover this week? Introduction to the Laplace transform Basic theory Domain and range of L Key
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationCubic Splines; Bézier Curves
Cubic Splines; Bézier Curves 1 Cubic Splines piecewise approximation with cubic polynomials conditions on the coefficients of the splines 2 Bézier Curves computer-aided design and manufacturing MCS 471
More informationLecture Note 3: Polynomial Interpolation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Polynomial Interpolation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 24, 2013 1.1 Introduction We first look at some examples. Lookup table for f(x) = 2 π x 0 e x2
More informationFinite Element Method for Ordinary Differential Equations
52 Chapter 4 Finite Element Method for Ordinary Differential Equations In this chapter we consider some simple examples of the finite element method for the approximate solution of ordinary differential
More information12.0 Properties of orthogonal polynomials
12.0 Properties of orthogonal polynomials In this section we study orthogonal polynomials to use them for the construction of quadrature formulas investigate projections on polynomial spaces and their
More informationThe use of exact values at quadrature points in the boundary element method
The use of exact values at quadrature points in the boundary element method N.D. Stringfellow, R.N.L. Smith Applied Mathematics & OR Group, Cranfield University, Shrivenham, Swindon, U.K., SN6 SLA E-mail:
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationName: INSERT YOUR NAME HERE. Due to dropbox by 6pm PDT, Wednesday, December 14, 2011
AMath 584 Name: INSERT YOUR NAME HERE Take-home Final UWNetID: INSERT YOUR NETID Due to dropbox by 6pm PDT, Wednesday, December 14, 2011 The main part of the assignment (Problems 1 3) is worth 80 points.
More informationMA3232 Summary 5. d y1 dy1. MATLAB has a number of built-in functions for solving stiff systems of ODEs. There are ode15s, ode23s, ode23t, ode23tb.
umerical solutions of higher order ODE We can convert a high order ODE into a system of first order ODEs and then apply RK method to solve it. Stiff ODEs Stiffness is a special problem that can arise in
More informationLecture 10: Eigenvalue Estimation
CS 880: Quantum Information Processing 9/7/010 Lecture 10: Eigenvalue Estimation Instructor: Dieter van Melkebeek Scribe: Dalibor Zelený Last time we discussed the quantum Fourier transform, and introduced
More informationThe Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods
The Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods by Hae-Soo Oh Department of Mathematics, University of North Carolina at Charlotte, Charlotte, NC 28223 June
More informationMA 3021: Numerical Analysis I Numerical Differentiation and Integration
MA 3021: Numerical Analysis I Numerical Differentiation and Integration Suh-Yuh Yang ( 楊肅煜 ) Department of Mathematics, National Central University Jhongli District, Taoyuan City 32001, Taiwan syyang@math.ncu.edu.tw
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationNumerical Integration (Quadrature) Another application for our interpolation tools!
Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
More informationPoisson Solvers. William McLean. April 21, Return to Math3301/Math5315 Common Material.
Poisson Solvers William McLean April 21, 2004 Return to Math3301/Math5315 Common Material 1 Introduction Many problems in applied mathematics lead to a partial differential equation of the form a 2 u +
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationSpan and Linear Independence
Span and Linear Independence It is common to confuse span and linear independence, because although they are different concepts, they are related. To see their relationship, let s revisit the previous
More informationQ1 Q2 Q3 Q4 Tot Letr Xtra
Mathematics 54.1 Final Exam, 12 May 2011 180 minutes, 90 points NAME: ID: GSI: INSTRUCTIONS: You must justify your answers, except when told otherwise. All the work for a question should be on the respective
More informationPolynomial Interpolation
Polynomial Interpolation (Com S 477/577 Notes) Yan-Bin Jia Sep 1, 017 1 Interpolation Problem In practice, often we can measure a physical process or quantity (e.g., temperature) at a number of points
More informationMath 251 December 14, 2005 Answer Key to Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Name Section Math 51 December 14, 5 Answer Key to Final Exam There are 1 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning
More informationMATH 4330/5330, Fourier Analysis Section 8, The Fourier Transform on the Line
MATH 4330/5330, Fourier Analysis Section 8, The Fourier Transform on the Line What makes the Fourier transform on the circle work? What is it about the functions φ n (x) e 2πinx that underlies their importance
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More information2. Interpolation. Closing the Gaps of Discretization...
2. Interpolation Closing the Gaps of Discretization... Numerisches Programmieren, Hans-Joachim Bungartz page 1 of 48 2.1. Preliminary Remark The Approximation Problem Problem: Approximate a (fully or partially
More information