(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
|
|
- Leslie Garrett
- 5 years ago
- Views:
Transcription
1 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 would you obtain an upper bound value to the magnitude of the error? (b) Given f(x) = (x + 1) log e (x), obtain the cubic Taylor polynomial P 3 (x) about x 0 = 1, and hence obtain an upper bound to the error, as a function of x, if x x 0 < 1. Use this information to estimate an upper bound to SOLUTION 1.5 (f(x) P 3 (x)) dx 0.5 (a) The Lagrange formula for the error is given by E(x) = (x x 0) n+1 f (n+1) (α), (n + 1)! where α is an unknown point between x and x 0 : either x 0 < α < x or x < α < x 0. If you can find upper and lower bounds for f (n+1) (α), then you can find an upper bound value for the error. (b) The cubic Taylor polynomial is given by P 3 (x) = f(x 0 ) + (x x 0 )f (x 0 ) + (x x 0) f (x 0 ) + (x x 0) 3 f (x 0 ). 6 Since our function f(x) = (x + 1) log e (x), it follows that f (x) = x log e (x) + x + 1/x, f (x) = log e (x) + 3 1/x, f (x) = /x + /x 3. Substituting our choice for x 0 = 1, we thus find that f(1) = 0 f (1) = f (1) = f (1) = 4. The Taylor polynomial P 3 (x) is thus given by P 3 (x) = (x 1) + (x 1) + (x 1)3. 3 In order to find a bound for the error, we need to calculate the fourth derivative of f(x) f (4) (x) = /x 6/x 4.
2 We need to find a bound for the case that x x 0 < 1/. The largest value for f (4) (x) on the interval 1/ < x < 3/ is found at x = 1/ : The upper bound for the error in our Taylor polynomial is thus given by E(x) = (x x 0) 4 ( 104). 4 For x = 1/, this error gives E = (f(x) P 3 (x)) dx 1.5 (x 1) 4 [ ] ( 104)dx = (x 1) ( 104) (f(x) P 3 (x)) dx
3 QUESTION. (a) Derive the Newton-Raphson method for solving equations in one variable. Compared to the false-position method, name an advantage and two disadvantages of the Newton-Raphson method. (b) Use the Newton-Raphson method to find an estimate for 3 with 6 figures accuracy by solving the equation f(x) = (x 3)(x 15) = 0 using an initial guess p 0 =. Why does the Newton-Raphson method converge faster to 3 for solving f(x) than for solving g(x) = (x 3) = 0? SOLUTION (a) In the Newton-Raphson method for finding the root of a function f(x) = 0, we start with an initial guess for the root p 0. We approximate the function f(x) by its Taylor polynomial of degree 1 P 1 (x) around p 0 and assume that the root of this Taylor polynomial is an improved guess for the root of the function f(x). We then iterate until we achieve convergence. Around the point p 0, the Taylor polynomial of degree 1 is given by P 1 (x) = f(p 0 ) + (x p 0 )f (p 0 ). The root of this Taylor polynomial is then given by or x p 0 = f(p 0) f (p 0 ) x = p 0 f(p 0) f (p 0 ). Compared to the false-position method, the Newton-Raphson method converges significantly faster. The drawbacks of the method are that convergence is no longer guaranteed, and that we need to be able to evaluate the derivative of the function we re solving for. (b) To apply the Newton-Raphson method, we need the derivative of f(x): f (x) = x(x 15) + x(x 3) = 4x(x 9) Our initial guess is p 0 = The next guess thus becomes The next guess becomes p = p 1 f(p 1) f (p 1 ) p 1 = p 0 f(p 0) f (p 0 ) = = 1.75 = =
4 The next guess becomes p 3 = p f(p ) f (p ) = = The last two guesses differ by , and 6 figures accuracy has thus been achieved. The Newton-Raphson method is based on Taylor polynomials. This means that the closer the function resembles a Taylor polynomial of degree 1 the faster convergence will be. This closeness can be estimated by looking at the Taylor polynomial of degree. We thus need to look at the second derivative of the functions f(x) and g(x) at 3. g (x) = f (x) = 1x 36 For f(x), the second derivative is 0 at 3. Close to 3, f(x) resembles a Taylor polynomial of degree 1 much better than g(x).
5 QUESTION 3. (a) What is meant by LU decomposition of a matrix? When solving a system of linear equations using LU decomposition, what advantages does LU decomposition have over Gauss Elimination? (b) Use LU decomposition to solve SOLUTION x 1 x x 3 = (a) LU decomposition is the decomposition of a matrix A into two matrices L and U, which are lower-triangular and upper triangular, respectively. Since the matrices L and U have n + n unknowns and there are n equations, we have the freedom to choose n coefficients. The most common choices are : the diagonal elements of L are all equal to 1 the diagonal elements of U are all equal to 1 the diagonal elements of L and U are equal to each other The advantage of LU decomposition over Gauss Elimination for solving a system of linear equations is that LU decomposition needs to be performed only once. Thus, if we need to solve Ax = b for many different b at different times, then we need to do the decomposition once, and all solutions for any b vector can then be obtained directly using forward- and back-substitutions. One calculation is of O(n 3 ), while all othe calculations are O(n ). For Gauss Elimination, we need to repeat the entire procedure whenever a new vector b is presented. The previous calculations cannot be used to simplify the calculations for a new b. Every calculation is thus of order O(n 3 ). (b) First we need to decompose A into the separate matrices L and U. We have to choose n = 3 coefficients: the diagonal elements of U are equal to 1. We must then solve the following equation: = α α 1 α 0 α 31 α 3 α 33 1 The first column of A gives the following equations α 11 = 0.5, α 1 = 1, and 1 β 1 β β
6 α 31 = 1. The remaining elements on the first row of A then give 1 = 0.5 β 1 or β 1 =, and 1 = 0.5 β 13 or β 13 =. Focussing our attention on the remaining elements of the second column of A: 4 = α 1 β 1 + α = + α. Thus α = 1 = α 31 β 1 + α 3 = + α 3. Thus α 3 = 1 The last element of the second row : = α 1 β 13 + α β 3 = + β 3. Thus β 3 = The last element of the third row : 3.5 = α 31 β 13 + α 3 β 3 + α 33 = + α 33. Thus α 33 = 0.5. So we have = Now we have to solve the equation Ly = b. Or y 1 y y 3 = Solving the equations using forward-substitution, we find 0.5y 1 = or y 1 = 4 y 1 + y = 1 or y =.5 y 1 y + 0.5y 3 = or y 3 = 9. Now, we have to use backsubstition to solve Ux = y and find x x 1 x x 3 = This gives the following equations: x 3 = 9 x + x 3 =.5 or x = 15.5, x 1 x + x 3 = 4 or x 1 = = 45. If we substitute this answer back into the original equation, we obtain and our answer is thus correct =. 1.
7 4. Question (a) Why are polynomials often used in interpolation? Why are Taylor polynomials not useful for interpolation? What is meant by Lagrangian interpolation? (b) Employ Neville s algorithm to obtain f(1) for the function tabulated below. Use the four most appropriate data points. Solution x f(x) (a) Polynomials are often used in interpolation because of Weierstrass theorem, which states that for every defined and continuous function f(x) on [a, b] and given ε > 0, there exists a polynomial P such that f(x) P (x) < ε for all x in [a, b]. Moreover, polynomials are easy to evaluate, differentiate and integrate. Taylor polynomials are not useful for interpolation, since all the information is concentrated at a single point x 0. Therefore it is an excellent approximation near x 0, but the error increases rapidly when you move away from x 0. In Lagrangian interpolation, we approximate a function f(x) by a polynomial P (x) whose values agree with f(x) at a set of distinct points {x k }: f(x k ) = P (x k ). (b) Neville s algorithm involves the repeated linear interpolation of two polynomials: P mab = (x x a)p mb (x x b )P ma x b x a where subscripts indicate the k-values where P (x k ) = f(x k ). m indicates a set of x k s common to both polynomials involved in the interpolation. For the present purpose, we need the 4 values closest to 1: 0.7, 0.9, 1.1, 1.3. In Neville s algorithm, we combine nearest neighbours. Thus we need to find P 134 (1). P 1 (1) = (1 0.7) (1 0.9) = P 3 (1) = P 34 (1) = (1 0.9) (1 1.1) (1 1.1)0.64 (1 1.3) = =
8 P 13 (1) = P 34 (1) = P 134 (1) = (1 0.7) (1 1.1)( 0.003) (1 0.9) (1 1.3)(0.1004) (1 0.7) (1 1.3)(0.070) = = =
9 5. Question (a) Show that the error e (k) in the k th iteration of the Jacobi method for solving the set of simultaneous linear equations Ax = b satisfies: e (k) = P k e (0) where the iteration matrix P is independent of k, and e (0) is the initial error vector. Show that the method converges for any e (0) if P ij < 1, all i. j (b) Apply three iterations of the Jacobi method to the following system of equations: 3x 1 + x x 3 = 1 x 3x 3 = 5 4x 1 + 6x x 3 = 6 starting from (0,0,0) using 4 significant figures. Identify the matrix P and show that convergence is guaranteed. Solution (a) In the Jacobi method, the matrix equation Ax = b is transformed to an equation x = P x + d, where P = D 1 (L + U) and d = D 1 b. D contains the diagonal elements of A, L and U are the strictly lower-triangular and upper-triangular parts of A respectively. In the Jacobi method, we iteratively calculate new guesses by while the exact solution is given by x (k) = P x (k 1) + d x = P x + d Subtracting the two equations gives for the error term x (k) x = P x (k 1) P x = P (x (k 1) x) Thus, e (k) = P x (k 1) P x = P e (k 1) = P k e (0)
10 The Jacobi method converges, if the norm of the error vector goes to zero when k goes to infinity. This will happen if the spectral radius of the matrix P is smaller than 1. The spectral radius is a lower bound to the norms. Thus if we find a norm of P that is smaller than 1, then we know that the spectral radius is smaller than 1. For example, the Chebyshev norm is given by A = A ij. j If this sum is less than 1, the spectral radius is less than 1 and the method converges. (b) The system of equations is first transformed to x (k+1) 1 = 1 xk + x k 3 3 x (k+1) = 6 + 4xk 1 + x k 3 6 x (k+1) 3 = 5 xk 3 Our initial guess is x 0 = (0, 0, 0). The first guess is x 1 = ( , 1, 1.667). The second guess is x = ( , 1.333, ). The third guess is x 3 = ( , 1.593,.111). P = 0 1/3 1/3 /3 0 1/3 0 1/3 0 The norms P 1 and P e are < 1. The spectral radius is thus < 1 and convergence is thus assured.
11 6. QUESTION (a) The spectrum of an isolated narrow resonance can be regarded as a smooth background signal plus a highly-localized disturbance on the signal. Suppose we make a finite-difference table of such a spectrum, and that only one point is affected by the localized disturbance. Explain how we can use the finite-difference table to locate the element in the table that is influenced by the disturbance, and how we can estimate the size of the disturbance. (b) Use finite-difference methods up to third differences to estimate f(1.1) and f(1.95) given the following values of f(x) : x f(x) SOLUTION (a) A finite-difference table starts of an ordered list of data points separated by a fixed step-size h. We then add another list to the table by calculating the difference between consecutive points in our first list. We can then continue by iteratively calculating differences between consecutive points for each new list created. When a function is smooth, the numbers will decrease with increasing order of differences. A localized distrubance can be considered an error on a data point. In the case of an error, the numbers will increase with increasing order of differences, and in fact they increase according to the binomial coefficients. In our table, we thus have to look for a set of numbers that increases with increasing of order of differences according to the binomial coefficients. The center of the cone gives the position of the error, while we can get an estimate of the error by dividing the n th difference by the appropriate binomial coefficient. (b) First of all, we need to make a finite-difference table. Since we want to go up to third differences, we need to determine 3 difference columns.
12 For f(1.), we need to use Newton s forward formula: f(1 + ph) = (1 + ) p f 1 = Thus, ( 1 + p + p(p 1) + f(1.) = f ) p(p 1)(p ) 3 f 1 6 f(1.) = ( 0.034) = For f(1.95), we need to use Newton s backward formula: f(+ph) = (1 ) p f 6 = Thus ( 1 + p + f(1.95) = f 6 + ( 1/4) 6 + ( 1/4)(3/4) p(p + 1) + ) p(p + 1)(p + ) 3 f ( 1/4)(3/4)(7/4) f(1.95) = ( ) =
13 7. QUESTION (a) What is meant by the L p and Chebyshev norms of a function? Use these to define the least squares and minimax approximations to a function, and discuss the advantages and disadvantages of the two approximations. (b) The 6 th order Taylor polynomial for the function f(x) = e (cos(x) 1) around 0 is given by P 6 (x) = 1 x + x4 6 31x6 70. Find a minimax quadratic approximation Q(x) on the interval [-0.5,0.5] to the fourth-order Taylor polynomial P 4 (x) of f(x) = e (cos(x) 1) around 0 using the Chebyshev polynomial T 4 (x) = 8x 4 8x + 1. Can Q(x) be considered a near-minimax quadratic approximation to f(x), and why? SOLUTION (a) The L p norm of a function f(x) on an interval [a, b] is given by [ b 1/p f p = f(x) dx] p. a The Chebyshev norm of a function f(x) on an interval [a, b] is given by f = max a x b f(x). The least squares approximation to a function f(x) is a function F (x) such that the L norm of f(x) F (x) is a minimum. The Chebyshev approximationto a function f(x) is a function F (x) such that the Chebyshev norm of f(x) F (x) is a minimum. The advantage of the least-squares approximation is that it is easy to calculate. On the other hand, a disadvantage is that f(x) F (x) may be large for some x. The advantage of the Chebyshev approximation is that you know the maximum error in the approximation F (x). The drawback is that the approximation is difficult to calculate. (b) We need to economise P 4 (x) using T 4 (x). P 4 (x) = 1 x + x4 6
14 In economizing, we first need to transform our interval [-0.5,0.5] to the interval [-1,1]. We can achieve this by the transformation t = x. We find a new polynomial P 4 (t), P 4 (t) = P 4 (t/) = 1 t 8 + t4 96 We can now economise P 4 (t) to give us a polynomial Q(t) by subtracting (1/96)(1/8)T 4 (x). Thus Q(t) = P 4 (t) T 4(x) = 1 t 8 + t4 96 t t = t 96 To get our minimax approximation, we need to transform Q(t) back to Q(x) using x = t. Thus Q(x) = Q(x) = x 4. In order to estimate whether Q(x) is a near-minimax approximation to f(x), we need to compare the errors made in the Taylor approximation and in the economization. The error in the economization is The error in the Taylor approximation can be estimated from the neglected term in P 6 (x). This last term is largest at the boundaries, x = 1/. The size of the error is approximately , while the actual error at 0.5 is The error in the economization is thus larger than the error of Taylor series, and Q(x) can be regarded as a near-minimax approximation to f(x) on [-0.5,0.5].
15 Question 8. (a) What is composite quadrature? Write down the composite Trapezoidal rule, indicating how the error varies with the step size h. (b) Evaluate 1 1 cos(e x )dx using the composite Trapezoidal rule and four significant figures. Use four different stepsizes :, 1, 0.5 and end with a stepsize of 0.5. Use these results to eliminate the leading term in the error expansion. Solution (a) In composite quadrature the interval of integration is split into subintervals. We then use the same low-order quadrature rule, eg. the trapezoidal rule, to each subinterval. The composite trapezoidal rule is b f(x)dx h [ ] n f 0 + f k + f n a k=1 (b a) h f (µ), 1 where h = (b a)/n and µ an unknown point in the interval (a, b). The error term has the form E(h) = a n h n, n=1 ie. it contains terms containing even powers of h. The knowledge of the error dependence can be exploited using Richardson s extrapolation to give the Romberg method. (b) We have to calculate I = 1 1 cos(e x )dx. We will approximate this integral by I(h) where h is the stepsize. I() = [f(1) + f( 1)] = ( ) = I(1) = I() + 1 [f(0)] = = I(0.5) = I(1) [f( 0.5) + f(0.5)] = = I(0.5) = I(0.5) [f( 0.75) + f( 0.5) + f(0.5) + f(0.75)]
16 I(0.5) = = The leading term in the error is a 1 h. Richardson s extrapolation tells us that the leading term ah n can be eleminated using stepsizes h and αh following I 1 (αh) = I 0(αh) α n I 0 (h) 1 α n In the present example α = 1/ and n =. Thus I 1 (1) = (0.014) = I 1 (0.5) = I 1 (0.5) = (0.5510) (0.6474) = =
17 Question 9. (a) Given the first order ordinary differential equation dy dt = f(t, y(t)), derive the general nd order Runge-Kutta formula, and indicate how it can be used to derive the mid-point method. (b) Given dy dx = y + x + 1 and y(1) = 0, use the mid-point method to obtain a value for y() using a stepsize of 0.5. Solution (a) In the second-order Runge Kutta method, we approximate the function value y(t + h) using a second-order Taylor polynomial expanded around y(t): y(t + h) = y(t) + hy (t) + h y (t) The first and second term on the right-hand side are trivial, since y(t) is known and y (t) = f(t, y(t)) is our differential equation. For the last term, y = df dt = df dt + df dy dy dt = f t + f y f If we now approximate f(t + αh, y + αhf(t, y)) by a first-order Taylor polynomial around f(t, y), we find Thus, f(t + αh, y + αhf(t, y)) = f(t, y) + αhf t + αhff y = f(t, y) + αhy y = 1 (f(t + αh, y + αhf(t, y)) f(t, y)) αh If we substitute this expression in the second-order Taylor polynomial, we find y(t + h) = y(t) + hf(t, y) + h 1 (f(t + αh, y + αhf(t, y)) f(t, y)) αh
18 or ( y(t + h) = y(t) + h 1 1 ) f(t, y) + h f(t + αh, y + αhf(t, y)) α α (b) In this last expression, we can choose any value for α. method corresponds to choosing α = 1/. dy dx = y + x + 1, The midpoint y(1) = 0, a stepsize of 0.5 and we need to use the midpoint method : y(t + h) = y(t) + hf(t + h/, y + hf(t, y)/). So, we start with (t, y) = (1, 0) y(1.5) = y(1) + 0.5f(1.15, 0.15f(1, 0)) f(1, 0) = 1.5 y(1.5) = f(1.15, ) = = y(1.5) = y(1.5) + 0.5f(1.375, f(1.5, )) f(1.5, ) = y(1.5) = f(1.375, ) = = y(1.75) = y(1.5) + 0.5f(1.65, f(1.5, 0.734)) f(1.5, 0.734) = 1.7 y(1.75) = f(1.65, ) = = y() = y(1.75) + 0.5f(1.875, f(1.75, )) f(1.75, ) = y() = f(1.875, 1.098) = = 1.167
Department 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 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 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 informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
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 informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation
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 informationNumerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point
Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error In this method we assume initial value of x, and substitute in the equation. Then modify x and continue till we
More informationIterative Methods. Splitting Methods
Iterative Methods Splitting Methods 1 Direct Methods Solving Ax = b using direct methods. Gaussian elimination (using LU decomposition) Variants of LU, including Crout and Doolittle Other decomposition
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 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 informationVirtual University of Pakistan
Virtual University of Pakistan File Version v.0.0 Prepared For: Final Term Note: Use Table Of Content to view the Topics, In PDF(Portable Document Format) format, you can check Bookmarks menu Disclaimer:
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 informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
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 informationHence a root lies between 1 and 2. Since f a is negative and f(x 0 ) is positive The root lies between a and x 0 i.e. 1 and 1.
The Bisection method or BOLZANO s method or Interval halving method: Find the positive root of x 3 x = 1 correct to four decimal places by bisection method Let f x = x 3 x 1 Here f 0 = 1 = ve, f 1 = ve,
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 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 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 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 informationPart IB Numerical Analysis
Part IB Numerical Analysis Definitions Based on lectures by G. Moore Notes taken by Dexter Chua Lent 206 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after
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 informationM.SC. PHYSICS - II YEAR
MANONMANIAM SUNDARANAR UNIVERSITY DIRECTORATE OF DISTANCE & CONTINUING EDUCATION TIRUNELVELI 627012, TAMIL NADU M.SC. PHYSICS - II YEAR DKP26 - NUMERICAL METHODS (From the academic year 2016-17) Most Student
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 informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More information3.1 Interpolation and the Lagrange Polynomial
MATH 4073 Chapter 3 Interpolation and Polynomial Approximation Fall 2003 1 Consider a sample x x 0 x 1 x n y y 0 y 1 y n. Can we get a function out of discrete data above that gives a reasonable estimate
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 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 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 informationMTH603 FAQ + Short Questions Answers.
Absolute Error : Accuracy : The absolute error is used to denote the actual value of a quantity less it s rounded value if x and x* are respectively the rounded and actual values of a quantity, then absolute
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 informationPage No.1. MTH603-Numerical Analysis_ Muhammad Ishfaq
Page No.1 File Version v1.5.3 Update: (Dated: 3-May-011) This version of file contains: Content of the Course (Done) FAQ updated version.(these must be read once because some very basic definition and
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 informationNumerical methods. Examples with solution
Numerical methods Examples with solution CONTENTS Contents. Nonlinear Equations 3 The bisection method............................ 4 Newton s method.............................. 8. Linear Systems LU-factorization..............................
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 informationOrdinary Differential Equations
CHAPTER 8 Ordinary Differential Equations 8.1. Introduction My section 8.1 will cover the material in sections 8.1 and 8.2 in the book. Read the book sections on your own. I don t like the order of things
More informationExamination paper for TMA4130 Matematikk 4N: SOLUTION
Department of Mathematical Sciences Examination paper for TMA4 Matematikk 4N: SOLUTION Academic contact during examination: Morten Nome Phone: 9 84 97 8 Examination date: December 7 Examination time (from
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 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 informationSOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS BISECTION METHOD
BISECTION METHOD If a function f(x) is continuous between a and b, and f(a) and f(b) are of opposite signs, then there exists at least one root between a and b. It is shown graphically as, Let f a be negative
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE712B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationSRI RAMAKRISHNA INSTITUTE OF TECHNOLOGY DEPARTMENT OF SCIENCE & HUMANITIES STATISTICS & NUMERICAL METHODS TWO MARKS
SRI RAMAKRISHNA INSTITUTE OF TECHNOLOGY DEPARTMENT OF SCIENCE & HUMANITIES STATISTICS & NUMERICAL METHODS TWO MARKS UNIT-I HYPOTHESIS TESTING 1. What are the applications of distributions? * Test the hypothesis
More informationSBAME CALCULUS OF FINITE DIFFERENCES AND NUMERICAL ANLAYSIS-I Units : I-V
SBAME CALCULUS OF FINITE DIFFERENCES AND NUMERICAL ANLAYSIS-I Units : I-V Unit I-Syllabus Solutions of Algebraic and Transcendental equations, Bisection method, Iteration Method, Regula Falsi method, Newton
More informationMath Introduction to Numerical Analysis - Class Notes. Fernando Guevara Vasquez. Version Date: January 17, 2012.
Math 5620 - Introduction to Numerical Analysis - Class Notes Fernando Guevara Vasquez Version 1990. Date: January 17, 2012. 3 Contents 1. Disclaimer 4 Chapter 1. Iterative methods for solving linear systems
More informationNUMERICAL MATHEMATICS AND COMPUTING
NUMERICAL MATHEMATICS AND COMPUTING Fourth Edition Ward Cheney David Kincaid The University of Texas at Austin 9 Brooks/Cole Publishing Company I(T)P An International Thomson Publishing Company Pacific
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 informationMath 5630: Iterative Methods for Systems of Equations Hung Phan, UMass Lowell March 22, 2018
1 Linear Systems Math 5630: Iterative Methods for Systems of Equations Hung Phan, UMass Lowell March, 018 Consider the system 4x y + z = 7 4x 8y + z = 1 x + y + 5z = 15. We then obtain x = 1 4 (7 + y z)
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 informationUnit I (Testing of Hypothesis)
SUBJECT NAME : Statistics and Numerical Methods SUBJECT CODE : MA645 MATERIAL NAME : Part A questions REGULATION : R03 UPDATED ON : November 07 (Upto N/D 07 Q.P) Unit I (Testing of Hypothesis). State level
More informationChapter 5: Numerical Integration and Differentiation
Chapter 5: Numerical Integration and Differentiation PART I: Numerical Integration Newton-Cotes Integration Formulas The idea of Newton-Cotes formulas is to replace a complicated function or tabulated
More informationINSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad MECHANICAL ENGINEERING TUTORIAL QUESTION BANK
Course Name Course Code Class Branch INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad - 500 043 Mathematics-II A30006 II-I B. Tech Freshman Engineering Year 016 017 Course Faculty MECHANICAL ENGINEERING
More informationINSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad
1 P a g e INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad - 500 04 Name : Mathematics-II Code : A0006 Class : II B. Tech I Semester Branch : CIVIL Year : 016 017 FRESHMAN ENGINEERING
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
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 informationChapter 11 ORDINARY DIFFERENTIAL EQUATIONS
Chapter 11 ORDINARY DIFFERENTIAL EQUATIONS The general form of a first order differential equations is = f(x, y) with initial condition y(a) = y a We seek the solution y = y(x) for x > a This is shown
More informationNumerical Analysis Exam with Solutions
Numerical Analysis Exam with Solutions Richard T. Bumby Fall 000 June 13, 001 You are expected to have books, notes and calculators available, but computers of telephones are not to be used during the
More informationPreliminary 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 informationApplied Numerical Analysis (AE2220-I) R. Klees and R.P. Dwight
Applied Numerical Analysis (AE0-I) R. Klees and R.P. Dwight February 018 Contents 1 Preliminaries: Motivation, Computer arithmetic, Taylor series 1 1.1 Numerical Analysis Motivation..........................
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 21 Numerical Differentiation PowerPoints organized by Dr. Michael R. Gustafson II, Duke University 1 All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
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 informationCS 323: Numerical Analysis and Computing
CS 323: Numerical Analysis and Computing MIDTERM #2 Instructions: This is an open notes exam, i.e., you are allowed to consult any textbook, your class notes, homeworks, or any of the handouts from us.
More informationProgram : M.A./M.Sc. (Mathematics) M.A./M.Sc. (Final) Paper Code:MT-08 Numerical Analysis Section A (Very Short Answers Questions)
Program : M../M.Sc. (Mathematics) M../M.Sc. (Final) Paper Code:MT-08 Numerical nalysis Section (Very Short nswers Questions) 1. Write two examples of transcendental equations. (i) x 3 + sin x = 0 (ii)
More informationCOURSE Iterative methods for solving linear systems
COURSE 0 4.3. Iterative methods for solving linear systems Because of round-off errors, direct methods become less efficient than iterative methods for large systems (>00 000 variables). An iterative scheme
More informationNumerical Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by
Numerical Analysis A Comprehensive Introduction H. R. Schwarz University of Zürich Switzerland with a contribution by J. Waldvogel Swiss Federal Institute of Technology, Zürich JOHN WILEY & SONS Chichester
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationSAZ3C NUMERICAL AND STATISTICAL METHODS Unit : I -V
SAZ3C NUMERICAL AND STATISTICAL METHODS Unit : I -V UNIT-I Introduction Mathematical Preliminaries Errors: Computations, Formula Errors in a Series Approximation Roots of Equations Linear Equations Bisection
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationExamination paper for TMA4125 Matematikk 4N
Department of Mathematical Sciences Examination paper for TMA45 Matematikk 4N Academic contact during examination: Anne Kværnø a, Louis-Philippe Thibault b Phone: a 9 66 38 4, b 9 3 0 95 Examination date:
More informationHomework and Computer Problems for Math*2130 (W17).
Homework and Computer Problems for Math*2130 (W17). MARCUS R. GARVIE 1 December 21, 2016 1 Department of Mathematics & Statistics, University of Guelph NOTES: These questions are a bare minimum. You should
More informationPreface. 2 Linear Equations and Eigenvalue Problem 22
Contents Preface xv 1 Errors in Computation 1 1.1 Introduction 1 1.2 Floating Point Representation of Number 1 1.3 Binary Numbers 2 1.3.1 Binary number representation in computer 3 1.4 Significant Digits
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 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 informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationNUMERICAL ANALYSIS PROBLEMS
NUMERICAL ANALYSIS PROBLEMS JAMES KEESLING The problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.. Solving Equations Problem.
More informationFixed Points and Contractive Transformations. Ron Goldman Department of Computer Science Rice University
Fixed Points and Contractive Transformations Ron Goldman Department of Computer Science Rice University Applications Computer Graphics Fractals Bezier and B-Spline Curves and Surfaces Root Finding Newton
More informationPrevious Year Questions & Detailed Solutions
Previous Year Questions & Detailed Solutions. The rate of convergence in the Gauss-Seidal method is as fast as in Gauss Jacobi smethod ) thrice ) half-times ) twice 4) three by two times. In application
More informationComputational Methods
Numerical Computational Methods Revised Edition P. B. Patil U. P. Verma Alpha Science International Ltd. Oxford, U.K. CONTENTS Preface List ofprograms v vii 1. NUMER1CAL METHOD, ERROR AND ALGORITHM 1 1.1
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Taylor s Theorem Can often approximate a function by a polynomial The error in the approximation
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
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
More informationASSIGNMENT BOOKLET. Numerical Analysis (MTE-10) (Valid from 1 st July, 2011 to 31 st March, 2012)
ASSIGNMENT BOOKLET MTE-0 Numerical Analysis (MTE-0) (Valid from st July, 0 to st March, 0) It is compulsory to submit the assignment before filling in the exam form. School of Sciences Indira Gandhi National
More informationNumerical Programming I (for CSE)
Technische Universität München WT / Fakultät für Mathematik Prof. Dr. M. Mehl B. Gatzhammer February 7, Numerical Programming I (for CSE) Repetition ) Floating Point Numbers and Rounding a) Let f : R R
More informationNumerical Differential Equations: IVP
Chapter 11 Numerical Differential Equations: IVP **** 4/16/13 EC (Incomplete) 11.1 Initial Value Problem for Ordinary Differential Equations We consider the problem of numerically solving a differential
More informationMathematical Methods for Numerical Analysis and Optimization
Biyani's Think Tank Concept based notes Mathematical Methods for Numerical Analysis and Optimization (MCA) Varsha Gupta Poonam Fatehpuria M.Sc. (Maths) Lecturer Deptt. of Information Technology Biyani
More informationNumerical Methods. Scientists. Engineers
Third Edition Numerical Methods for Scientists and Engineers K. Sankara Rao Numerical Methods for Scientists and Engineers Numerical Methods for Scientists and Engineers Third Edition K. SANKARA RAO Formerly,
More informationEuler s Method, cont d
Jim Lambers MAT 461/561 Spring Semester 009-10 Lecture 3 Notes These notes correspond to Sections 5. and 5.4 in the text. Euler s Method, cont d We conclude our discussion of Euler s method with an example
More informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
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 informationReview of matrices. Let m, n IN. A rectangle of numbers written like A =
Review of matrices Let m, n IN. A rectangle of numbers written like a 11 a 12... a 1n a 21 a 22... a 2n A =...... a m1 a m2... a mn where each a ij IR is called a matrix with m rows and n columns or an
More informationNumerical Methods for Ordinary Differential Equations
Numerical Methods for Ordinary Differential Equations Answers of the exercises C Vuik, S van Veldhuizen and S van Loenhout 08 Delft University of Technology Faculty Electrical Engineering, Mathematics
More informationThe problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.
NUMERICAL ANALYSIS PRACTICE PROBLEMS JAMES KEESLING The problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.. Solving Equations
More informationNumerical Methods for Engineers. and Scientists. Applications using MATLAB. An Introduction with. Vish- Subramaniam. Third Edition. Amos Gilat.
Numerical Methods for Engineers An Introduction with and Scientists Applications using MATLAB Third Edition Amos Gilat Vish- Subramaniam Department of Mechanical Engineering The Ohio State University Wiley
More informationOutline. 1 Interpolation. 2 Polynomial Interpolation. 3 Piecewise Polynomial Interpolation
Outline Interpolation 1 Interpolation 2 3 Michael T. Heath Scientific Computing 2 / 56 Interpolation Motivation Choosing Interpolant Existence and Uniqueness Basic interpolation problem: for given data
More informationUSHA RAMA COLLEGE OF ENGINEERING & TECHNOLOGY
Code No: R007/R0 Set No. I B.Tech I Semester Supplementary Examinations, Feb/Mar 04 MATHEMATICAL METHODS ( Common to Civil Engineering, Electrical & Electronics Engineering, Computer Science & Engineering,
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 informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationNumerical Methods with MATLAB
Numerical Methods with MATLAB A Resource for Scientists and Engineers G. J. BÖRSE Lehigh University PWS Publishing Company I(T)P AN!NTERNATIONAL THOMSON PUBLISHING COMPANY Boston Albany Bonn Cincinnati
More informationAM205: Assignment 3 (due 5 PM, October 20)
AM25: Assignment 3 (due 5 PM, October 2) For this assignment, first complete problems 1, 2, 3, and 4, and then complete either problem 5 (on theory) or problem 6 (on an application). If you submit answers
More informationMAT 460: Numerical Analysis I. James V. Lambers
MAT 460: Numerical Analysis I James V. Lambers January 31, 2013 2 Contents 1 Mathematical Preliminaries and Error Analysis 7 1.1 Introduction............................ 7 1.1.1 Error Analysis......................
More information