Department of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
|
|
- Joanna Bailey
- 5 years ago
- Views:
Transcription
1 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 24 marks are available for each question For two-part questions, part (a) : bookwork part (b) : unseen problem For question 3, part (a) : bookwork parts (b) and (c) : unseen problems For questions 5 and 8 parts (a) and (b) : bookwork part (c) : an unseen problem
2 Page 2 of 5 110AMA (a) Given an n th 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 describing any features whose meaning is not conveyed directly by the notation used. How might you be able to obtain an upper bound value to the magnitude of the error? (b) Given f(x) = x ln(x) x, obtain the cubic Taylor polynomial P 3 (x) about x 0 = 1. Show that f (4) (x) < 16 if x 1 < 0.5 and hence obtain an upper bound to the error, as a function of x, if x 1 < 0.5. Use this information to estimate an upper bound to 1.5 (f(x) P 3 (x)) dx. SOLUTION (a) The Lagrange formula for the error is given by 0.5 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) 2 f (x 0 ) + (x x 0) 3 f (x 0 ). 2 6 Since our function f(x) = x ln(x) x, it follows that f (x) = ln(x), f (x) = 1, x f (x) = 1. x 2 Substituting our choice for x 0 = 1, we thus find that f(1) = 1 f (1) = 0 f (1) = 1 f (1) = 1. The Taylor polynomial P 3 (x) is thus given by P 3 (x) = 1 + (x 1)2 2 (x 1)3. 6
3 In order to find a bound for the error, we need to calculate the fourth derivative of f(x) f (4) (x) = 2 x 3. We need to find a bound for the case that x 1 < 1/2. The maximum value for the derivative is obtained when x is set to 1/2. We then obtain a maximum value for f (4) (x) on the interval 1/2 x 3/2 of 16. The upper bound for the error in our Taylor polynomial is thus given by E(x) < 16 (x 1)4. 24 The upper bound to the error in the integral is then found as follows: (f(x) P 3 (x)) dx < (x 1)4 [ 16 dx = (x 1) (f(x) P 3 (x)) dx < ] = 1 120
4 2. (a) Derive the Newton-Raphson method for solving equations in one variable. Compared to the false-position method, name one advantage and two disadvantages of the Newton-Raphson method. (b) Use three iterations of the Newton-Raphson method to approximate the solution to f(x) = e x tan x = 0 in the interval ( π/2, π/2). A good initial guess to this problem would be x 4/3. How accurate do you estimate the final answer to be? 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, P 1 (x) = 0, 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 convergence rate is quadratic rather than linear. 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) = e x 1 cos 2 x Our initial guess is p 0 = 4/3 = The next guess thus becomes p 1 = p 0 f(p 0) f (p 0 ) = =
5 The next guess becomes p 2 = p 1 f(p 1) f (p 1 ) The next guess becomes = = p 3 = p 2 f(p 2) f (p 2 ) = = One way for estimating the accuracy in the guess is to examine the difference between consecutive guesses. This suggests that the results is accurate to within (Alternatively, one can estimate from the function value and the derivative an accuracy within , although this corresponds to an additional Newton-Raphson step).
6 3. (a) Why is pivoting important in finding the solution to a system of linear equations? What is meant by LU decomposition? How can you use LU decomposition to determine the determinant of a matrix? (b) Use LU decomposition to solve the following system of linear equations, x 1 x 2 x 3 = (c) Determine A, L and U, where A = LU are the matrix in part(b) and its decomposition. Are these norms consistent with each other? Provide reasons for your answer SOLUTION (a) The essential matrix element in step i is element a ii, the element of equation i on the diagonal, also known as the pivot. In pivoting, one aims to make the pivot the largest element in the column (excluding the part of the column above the diagonal). Pivoting is important, since the accuracy of the calculations is improved when this element is as large as possible. A small pivot may lead to a significant round-off error. In LU decomposition, the matrix A is split into a lower triangular matrix L and an upper triangular matrix U, such that LU = A. Since the matrices L and U have n 2 +n unknown entries and the matrix A has n 2 known entries, we need to choose n entries. A common choice is to set the diagonal entries of either L or U equal to 1. LU decomposition is useful in calculating the determinant of matrix A, since you can calculate the determinant of A by multiplying all the elements on the diagonals of L and U. If you made the choice of all diagonal elements of U are equal to 1, then you simply need to multiply only the diagonal elements of L. (b) If we want to solve the set of equations Ax = b using LU decomposition, then we first decompose the matrix A into two matrices L and U, so that LU = A. We then first solve Ly = b using forward-substitution and then solve Ux = y using backsubstitution. So, first we need to decompose A = l l 21 l 22 0 l 31 l 32 l 33 1 u 12 u u , where we have chosen the diagonal elements of U to be equal to 1.
7 The first row of A gives the following set of equations: l 11 = 2 l 11 u 12 = 4 u 12 = 2 l 11 u 13 = 4 u 13 = 2 The second row of A gives the following set of equations: l 21 = 1 l 21 u 12 + l 22 = 1 l 22 = 1 l 21 u 13 + l 22 u 23 = 1 u 23 = 1 The third row of A gives the following set of equations: l 31 = 2 l 31 u 12 + l 32 = 4 l 32 = 0 l 31 u 13 + l 32 u 23 + l 33 = 2 l 33 = 2 Now we have to solve Ly = b y 1 y 2 y 3 = Hence y 1 = 2/2 = 1 y 2 = 0 y 1 = 1 y 3 = (4 + 2y 1 )/2 = 1 Now we have to solve Ux = y x 1 x 2 x 3 = Hence x 3 = 1 x 2 = 1 + x 3 = 2
8 x 1 = 1 2x 3 + 2x 2 = 1 Substituting the answer back into the original problem gives and hence the answer is correct = (c) The norm is given by the maximum sum of absolute values along a row, A = max A ji j Hence A = 10. Similarly L = 4, and U = 5. These norms are consistent with each other, since A L U. i 2 0 4,
9 Page 3 of 5 110AMA (a) Briefly describe Neville s algorithm to perform an interpolation of a tabulated function. What advantage does Neville s method have over computing the Lagrange polynomial directly? (b) In inverse interpolation, one asks the following question: given a function y = f(x) tabulated at a discrete set of points, at what value of x does y assume a specified value? By reversing the roles of the dependent variable y and independent variable x, use inverse interpolation with Neville s algorithm to find the value of x for which y(x) = 0, given the data points in the table below. Use the 4 most appropriate data points. x k y k SOLUTION (a) Neville s algorithm constructs the Lagrange interpolating polynomial through iterated interpolation. Suppose P ms and P mt are two polynomials. P ms agrees with a function f(x) at a set of points x m and a point x s. P mt agrees with a function f(x) at the identical set of points x m and a point x t. We can then construct a polynomial that agrees with f(x) at the set x m, x s and x t by applying the formula P mst = (x x s)p mt (x) (x x t )P ms (x) x t x s The advantage of Neville s algorithm when we want to know the function at one particular point, is that if we want to improve the approximation through adding data points, these additional data points can be added significantly easier. (b) Let s start by rewriting the table by renaming x k as Y k and y k as X k. X k Y k Now we have a table that we can apply Neville s method to in the regular fashion. The four most appropriate points are the closest points to 0 : , , and Hence P 1 = 1, P 2 = 1, P 3 = 3, and P 4 = 5
10 We first determine P 12, P 23 and P 34 using Neville s algorithm and X = 0, P stu = (X X s)p tu (X X t )P su X t X s. P 12 = (0 ( ))1 (0 ( ))( 1) ( ) = P 23 = P 34 = P 123 and P 234 are then given by (0 ( ))3 ( ) ( ) ( )5 (0 5.75) = , = , P 123 = (0 ( )) ( )3.0090) ( ( ) = P 234 = Finally, P 1234 is given by (0 ( )) (0 5.75)2.3700) (5.75 ( ) = P 1234 = (0 ( )) (0 5.75)2.4897) (5.75 ( )) =
11 5. In the fixed-point method, one tries to find the root of a given equation by rewriting the given equation into the form x = f(x). Under certain conditions for f (x), one can get an improved estimate x 1 from an initial guess x 0 by calculating x 1 = f(x 0 ). In these circumstances, repeated calculation of x n+1 = f(x n ) gives the root to the equation. The fixed-point method converges linearly. (a) Give the definition of the convergence rate. (b) Linearly convergent methods can be accelerated using Aitken s 2 method. Derive Aitken s 2 method from the definition of the convergence rate. (c) Using an initial guess of x 0 = 0 for the problem x = cos(x), gives an improved guess of x 1 = cos(x 0 ) = cos(0) = 1. Apply the fixed-point method for a further three steps. Estimate from these three steps how many iterations you would need to achieve an accuracy of (d) Apply Aitken s 2 algorithm to the three new estimates from part (c) to find an improved estimate to the solution of x = cos(x). SOLUTION (a) The convergence rate α is the particular value of α for which p p n+1 lim n p p n C, α where p is the exact answer, p n a series of guesses converging to p and C is a constant. (b) A linear convergence means α = 1. If the signs of p p n either alternate or remain the same, we can write C p p n+1 p p n p p n+2 p p n+1. If we now assume the difference between the two approximations to be negligible, we get (p p n+1 ) 2 = (p p n )(p p n+2 ) Hence or p 2 2pp n+1 + p 2 n+1 = p 2 p(p n + p n+2 ) + p n p n+2 p(p n 2p n+1 +p n+2 ) = p n p n+2 p 2 n+1 = p n (p n 2p n+1 +p n+2 ) p 2 n+2p n p n+1 p 2 n+1 Thus p = p n (p n p n+1 ) 2 (p n 2p n+1 + p n+2 )
12 (c) We have x 1 = 1, and we need three additional iterations. Following the scheme outlined in the question, x 2 = cos x 1 = cos 1 = x 3 = cos x 2 = cos = x 4 = cos x 3 = cos = To estimate the number of iterations needed for a precision of 10 6, we consider the difference between consecutive guesses. x 3 x 2 = , and : x 4 x 3 = We want this difference to become Assuming linear convergence, each iteration reduces the difference between consecutive guesses by a factor In order to reach an accuracy of 10 6, we would need a further N guesses, where N is given by N < 10 6, :N 28 (d) The acceleration using Aitken s alogrithm then gives us x = x 2 (x 2 x 3 ) 2 (x 2 2x 3 + x 4 ) =
13 Page 4 of 5 110AMA (a) Define the finite difference operators E, and D. Derive the Newton forward difference formula for interpolation. Derive the relationship between D and, with h the stepsize. hd = ln(1 + ), (b) Given the tabulated function f(x) below, construct a difference table up to 3 rd differences. Use this difference table to find an approximation to f(0.12) and f(0.85) employing the most appropriate finite-difference formula for interpolation. Justify your choice of finite-difference formula. x k f k Also use the forward-difference formula for the derivative at a point p, f p = 1 ( ) p 2 p h p 3..., (1) to estimate the derivative of f(x) at x = 0.2. SOLUTION (a) The finite difference operators are defined as follows Ef(x) = f(x + h) f(x) = f(x + h) f(x) = (E 1)f(x) Df(x) = f (x) This leads to the following expressions for E : E = 1 + If we want to know the value of f(x) at a point x = x 0 + ph, then we have = ( f(x) = f(x 0 + ph) = E p f(x 0 ) = (1 + ) p f(x 0 ) 1 + p + = f(x 0 ) + p 0 + p(p 1) p(p 1) ) p(p 1)(p 2) f(x 0 ) 6 p(p 1)(p 2)
14 In order to relate and D, we can use Taylor s polynomial: f(x + h) = Ef(x) = ( + 1)f(x) Hence or f(x + h) = f(x) + hf (x) + h2 2! f (x) +... = f(x) + hdf(x) + h2 D 2 f(x) +... = e hd f(x) 2! + 1 = e hd hd = ln( + 1) (b) First, we have to make a difference table In order to determine f(0.15), we need to use the forward difference formula, starting from 0.1 with p = 0.5, since this formula uses the most accurate data near the top of the table. f(0.15) = f(0.1) f(0.1) f(0.1) f(0.1) = = In order to determine f(0.88), we need to use the backward difference formula, f(x 0 + ph) = f(x 0 ) + p 0 + p(p + 1) p(p + 1)(p + 2)
15 starting from 0.9 with p = 0.2, since this formula uses the most accurate data at the end of the table. f(0.88) = f(0.9) 0.2 f(0.9) f(0.9) f(0.9) = = To evaluate the derivative at x = 0.2, we use the equation given f p = 1 h ( ) p 2 p The step-size used in the table h = 0.1. For x = 0.2, the finite-difference table shows that 2 = , 2 2 = and 3 2 = Hence the derivative at x = 0.2 is given by f 2 = 10( ) = 2.370
16 7. (a) What is meant by the least-squares polynomial approximation to a function f(x) and what is meant by the minimax approximation to a function f(x)? Explain why the Legendre polynomials are useful in obtaining least-squares polynomial approximations. (b) Obtain the quadratic least-squares approximation to the function f(x) = e x in the interval [0,2] using the Legendre polynomials. SOLUTION (a) The least-squares polynomial approximation to a function f(x) is the polynomial P n (x) of degree n, which minimizes the L 2 norm of (f(x) P n (x)) (over a given range ([a, b]). The minimax polynomial approximation to a function f(x) is the polynomial P n (x) of degree n, which minimizes the maximum value of (f(x) P n (x)) (over a given range ([a, b]). The Legendre polynomials are well suited for least-squares fitting, since they are orthogonal functions over [-1,1] and a weight function w(x) = 1, (P i, P j ) = 0 if i and j are not equal. In least-squares polynomial fitting using functions φ i (x), one needs to solve a set of equations for a : Ua = b, where U = (φ i, φ j ) and b = (φ i, f). The matrix U becomes a diagonal matrix when orthogonal polynomials are used. The coefficients a can then be found directly by a i = (φ i, f) (φ i, φ i ). (b) In order to obtain the least-squares approximation, we use the Legendre polynomials up to order 2. These polynomials give the least-squares approximation on the interval [-1,1] following a k = 2k f(x)p k (x)dx. First we need to transform the interval [0,2] over which we want to obtain the leas-squares approximation to the interval [-1,1]. This can be achieved by the transformation t = x 1, or x = t + 1. The transformed function F(t) is thus F (t) = e t 1.
17 Now we can apply the formula for the coefficients and the Legendre polynomials: P 0 (t) = 1, P 1 (t) = t, P 2 (t) = 1 2 (3t2 1). So a 1 = a 0 = e t 1 dt = 1 2 [ e t 1 ] 1 1 = 1 e 2 2 te t 1 dt = 3 2 = 3 2 ( 1 ) [ te t 1 ] e t 1 dt = 1 ( e e 2) = 3e 2 a 2 = (3t 2 1)e t 1 dt 3t 2 e t 1 dt = [ 3t 2 e t 1 ] = 3e e te t 1 dt = a 2 = 5 4 (3 15e e 2 ) = 5 2 (1 7e 2 ) Thus the least-squares approximation is given by Q(t) = 1 e 2 1 3e 2 t (1 7e 2 ) 3t2 1 2 Transforming this back into a least-squares fit for the original interval [0,2], we use t = x 1 and obtain q(x) = 1 e 2 2 3e 2 (x 1) (1 7e 2 )(3(x 1) 2 1) = = 3 15e 2 + ( e 2 )x + ( e 2 )x 2.
18 8. (a) What is meant by adaptive quadrature? (b) Suppose that S 0, S 1, S 2 are the approximations to the integral of a function f(x) using Simpson s method over the intervals [a, b], [a, c] and [c, b] respectively with c the midpoint of the interval [a, b]. Since the error in Simpson s method is given by E 0 = show that b a if, approximately, f(x)dx S 0 = h5 90 f (4) (µ), with h = (b a)/2 b f(x)dx S 1 S 2 < ɛ a S 0 S 1 S 2 < 15ɛ. (c) Use adaptive quadrature and Simpson s rule to calculate π/2 with an (absolute) precision of SOLUTION 0 sin xdx (a) In adaptive quadrature, the step-size used in the quadrature is varied throughout the interval of integration, so that the approximation error is evenly distributed over the interval. Generally this means a smaller step size where the function varies rapidly and a larger step size where the function varies slower. (b) We use two methods to evaluate the integral. Let h = (b a)/2. Then and b a b a f(x)dx = S 1 We can find a µ so that f(x)dx = S 0 h5 90 f (4) (µ), h f (4) (µ ) + S 2 h f (4) (µ ). h f (4) (µ ) + h f (4) (µ ) = h5 + h f (4) (µ). We now make the approximation that µ = µ, and use the fact that c is the midpoint of [a, b]. b a f(x)dx S 1 + S 2 2 h f (4) (µ) (2)
19 If we now set the two expressions for the integral equal to each other, then we get S 0 h5 90 f (4) (µ) = S 1 + S f (4) (µ) or S 0 S 1 S 2 = 15 h f (4) (µ). If we substitute this into equation 2, we obtain b a h 5 f(x)dx S 1 + S (S 0 S 1 S 2 ), which demonstrates what we needed to show. (c) First we use Simpson s rule, b a f(x)dx b a 6 over the entire interval : π/2 0 ( f(a) + 4f( a + b ) 2 ) + f(b) sin xdx = π (0 + 4 sin(π/4) + sin(π/2)) = Then we use Simpson s rule over the subintervals [0, π/4] and [π/4, π/2] separately. π/2 π/4 π/4 0 sin xdx = π (0 + 4 sin(π/8) + sin(π/4)) = sin xdx = π (sin(π/4) + 4 sin(3π/8) + sin(π/2)) = Adding the integrals together gives We now have to check whether the difference of the two calculations for the integral from 0 to π/2 is less than 10 times the specified tolerance, (We use 10 rather than 15 as a safety margin: µ = µ is only an assumption.) The difference is Therefore, within , π/2 0 sin xdx = ,
20 9. (a) Explain how you can use the (Runge-Kutta) midpoint approach to solve the second-order differential equation d 2 y dt 2 = dy dt + y + 2t, when you know y and y at t=0. Give the equation(s) to which you apply the midpoint method, and indicate the order in which calculations need to be performed. Note: you do not need to solve this equation. (b) Given dy dt = 2y + t and y(1) = 2, use the mid-point method to obtain a value for y(2) using a stepsize of 0.25 and four significant figures. SOLUTION (a) In order to apply the midpoint approach to a second-order differential equation, you transform the differential equation into a set of two firstorder equations. Define u 1 = y Then, obviously, and, less obviously, u 2 = y u 1 = f 1 = u 2 u 2 = f 2 = u 2 + u 1 + 2t. When we apply the midpoint method to this approach, we first calculate the derivative of both u 1 and u 2 at t = t 0, multiplied by the step-size k 1u1 = hf 1 (u 1, u 2, t 0 ) k 1u2 = hf 2 (u 1, u 2, t 0 ). Using these derivatives, we can then estimate the derivative at the midpoint k 2u1 = hf 1 (u 1 + k 1u1 /2, u 2 + k 1u2 /2, t + h/2) k 2u2 = hf 2 (u 1 + k 1u1 /2, u 2 + k 1u2 /2, t + h/2)
21 (b) It is essential that this second step is only performed after we have completed the first step. Once we have calculated these derivatives, we can now propagate the solution: u 11 = u 10 + k 2u1 u 21 = u 20 + k 2u2 Now we just repeat the midpoint method until we have reached the end of the propagation interval. dy dt = 2y + t, y(1) = 2, a stepsize of 0.25 and we need to use the midpoint method : y(t + h) = y(t) + hf(t + h/2, y + hf(t, y)/2). So, we start with (x, y) = (1, 2) y(1.25) = y(1) f(1.125, f(1, 2)) f(1, 2) = 3 y(1.25) = f(1.125, 1.625) = = y(1.5) = y(1.25) f(1.375, f(1.25, )) f(1.25, ) = y(1.5) = f(1.375, ) = = y(1.75) = y(1.5) f(1.625, f(1.5, )) f(1.5, ) = y(1.75) = f(1.625, ) = = y(2) = y(1.75) f(1.875, f(1.75, )) f(1.75, ) = y(2) = f(1.875, ) = =
(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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationToday s class. Linear Algebraic Equations LU Decomposition. Numerical Methods, Fall 2011 Lecture 8. Prof. Jinbo Bi CSE, UConn
Today s class Linear Algebraic Equations LU Decomposition 1 Linear Algebraic Equations Gaussian Elimination works well for solving linear systems of the form: AX = B What if you have to solve the linear
More informationMATH 1242 FINAL EXAM Spring,
MATH 242 FINAL EXAM Spring, 200 Part I (MULTIPLE CHOICE, NO CALCULATORS).. Find 2 4x3 dx. (a) 28 (b) 5 (c) 0 (d) 36 (e) 7 2. Find 2 cos t dt. (a) 2 sin t + C (b) 2 sin t + C (c) 2 cos t + C (d) 2 cos t
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 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 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 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 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 informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 20, 2014 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
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 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 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 informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 17, 2017 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
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 information(x x 0 )(x x 1 )... (x x n ) (x x 0 ) + y 0.
> 5. Numerical Integration Review of Interpolation Find p n (x) with p n (x j ) = y j, j = 0, 1,,..., n. Solution: p n (x) = y 0 l 0 (x) + y 1 l 1 (x) +... + y n l n (x), l k (x) = n j=1,j k Theorem Let
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 informationChapter 3: Root Finding. September 26, 2005
Chapter 3: Root Finding September 26, 2005 Outline 1 Root Finding 2 3.1 The Bisection Method 3 3.2 Newton s Method: Derivation and Examples 4 3.3 How To Stop Newton s Method 5 3.4 Application: Division
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 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 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 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 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 informationSolution: In standard form (i.e. y + P (t)y = Q(t)) we have y t y = cos(t)
Math 380 Practice Final Solutions This is longer than the actual exam, which will be 8 to 0 questions (some might be multiple choice). You are allowed up to two sheets of notes (both sides) and a calculator,
More informationMatrix Factorization Reading: Lay 2.5
Matrix Factorization Reading: Lay 2.5 October, 20 You have seen that if we know the inverse A of a matrix A, we can easily solve the equation Ax = b. Solving a large number of equations Ax = b, Ax 2 =
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 informationComputational Methods. Systems of Linear Equations
Computational Methods Systems of Linear Equations Manfred Huber 2010 1 Systems of Equations Often a system model contains multiple variables (parameters) and contains multiple equations Multiple equations
More information3.1 Introduction. Solve non-linear real equation f(x) = 0 for real root or zero x. E.g. x x 1.5 =0, tan x x =0.
3.1 Introduction Solve non-linear real equation f(x) = 0 for real root or zero x. E.g. x 3 +1.5x 1.5 =0, tan x x =0. Practical existence test for roots: by intermediate value theorem, f C[a, b] & f(a)f(b)
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 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 informationExam 3 MATH Calculus I
Trinity College December 03, 2015 MATH 131-01 Calculus I By signing below, you attest that you have neither given nor received help of any kind on this exam. Signature: Printed Name: Instructions: Show
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 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 informationMath 113 Winter 2005 Key
Name Student Number Section Number Instructor Math Winter 005 Key Departmental Final Exam Instructions: The time limit is hours. Problem consists of short answer questions. Problems through are multiple
More informationJune 2011 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations
June 20 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations The topics covered in this exam can be found in An introduction to differential equations
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 informationSpring 2017 Midterm 1 04/26/2017
Math 2B Spring 2017 Midterm 1 04/26/2017 Time Limit: 50 Minutes Name (Print): Student ID This exam contains 10 pages (including this cover page) and 5 problems. Check to see if any pages are missing. Enter
More informationNumerical and Statistical Methods
F.Y. B.Sc.(IT) : Sem. II Numerical and Statistical Methods Time : ½ Hrs.] Prelim Question Paper Solution [Marks : 75 Q. Attempt any THREE of the following : [5] Q.(a) What is a mathematical model? With
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 informationSolutions to Math 41 Final Exam December 10, 2012
Solutions to Math 4 Final Exam December,. ( points) Find each of the following limits, with justification. If there is an infinite limit, then explain whether it is or. x ln(t + ) dt (a) lim x x (5 points)
More informationInterpolating Accuracy without underlying f (x)
Example: Tabulated Data The following table x 1.0 1.3 1.6 1.9 2.2 f (x) 0.7651977 0.6200860 0.4554022 0.2818186 0.1103623 lists values of a function f at various points. The approximations to f (1.5) obtained
More informationNumerical and Statistical Methods
F.Y. B.Sc.(IT) : Sem. II Numerical and Statistical Methods Time : ½ Hrs.] Prelim Question Paper Solution [Marks : 75 Q. Attempt any THREE of the following : [5] Q.(a) What is a mathematical model? With
More informationUniversity of Connecticut Department of Mathematics
University of Connecticut Department of Mathematics Math 1131 Sample Exam 2 Fall 2015 Name: Instructor Name: Section: TA Name: Discussion Section: This sample exam is just a guide to prepare for the actual
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 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 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 informationLinear Algebra Section 2.6 : LU Decomposition Section 2.7 : Permutations and transposes Wednesday, February 13th Math 301 Week #4
Linear Algebra Section. : LU Decomposition Section. : Permutations and transposes Wednesday, February 1th Math 01 Week # 1 The LU Decomposition We learned last time that we can factor a invertible matrix
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 informationInterpolation Theory
Numerical Analysis Massoud Malek Interpolation Theory The concept of interpolation is to select a function P (x) from a given class of functions in such a way that the graph of y P (x) passes through the
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 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 informationMath 112 (Calculus I) Final Exam
Name: Student ID: Section: Instructor: Math 112 (Calculus I) Final Exam Dec 18, 7:00 p.m. Instructions: Work on scratch paper will not be graded. For questions 11 to 19, show all your work in the space
More informationToday s class. Numerical differentiation Roots of equation Bracketing methods. Numerical Methods, Fall 2011 Lecture 4. Prof. Jinbo Bi CSE, UConn
Today s class Numerical differentiation Roots of equation Bracketing methods 1 Numerical Differentiation Finite divided difference First forward difference First backward difference Lecture 3 2 Numerical
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 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 informationExample 1 Which of these functions are polynomials in x? In the case(s) where f is a polynomial,
1. Polynomials A polynomial in x is a function of the form p(x) = a 0 + a 1 x + a 2 x 2 +... a n x n (a n 0, n a non-negative integer) where a 0, a 1, a 2,..., a n are constants. We say that this polynomial
More informationSECTION A. f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes.
SECTION A 1. State the maximal domain and range of the function f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes. 2. By evaluating f(0),
More informationSimpson s 1/3 Rule Simpson s 1/3 rule assumes 3 equispaced data/interpolation/integration points
CE 05 - Lecture 5 LECTURE 5 UMERICAL ITEGRATIO COTIUED Simpson s / Rule Simpson s / rule assumes equispaced data/interpolation/integration points Te integration rule is based on approximating fx using
More informationNumerical Methods of Approximation
Contents 31 Numerical Methods of Approximation 31.1 Polynomial Approximations 2 31.2 Numerical Integration 28 31.3 Numerical Differentiation 58 31.4 Nonlinear Equations 67 Learning outcomes In this Workbook
More informationFinal Examination. CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Fall 2013), Stanford University
Final Examination CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Fall 2013), Stanford University The exam runs for 3 hours. The exam contains eight problems. You must complete the first
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 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 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 informationCOURSE Numerical integration of functions
COURSE 6 3. Numerical integration of functions The need: for evaluating definite integrals of functions that has no explicit antiderivatives or whose antiderivatives are not easy to obtain. Let f : [a,
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 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 information18.01 EXERCISES. Unit 3. Integration. 3A. Differentials, indefinite integration. 3A-1 Compute the differentials df(x) of the following functions.
8. EXERCISES Unit 3. Integration 3A. Differentials, indefinite integration 3A- Compute the differentials df(x) of the following functions. a) d(x 7 + sin ) b) d x c) d(x 8x + 6) d) d(e 3x sin x) e) Express
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 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 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 information2t t dt.. So the distance is (t2 +6) 3/2
Math 8, Solutions to Review for the Final Exam Question : The distance is 5 t t + dt To work that out, integrate by parts with u t +, so that t dt du The integral is t t + dt u du u 3/ (t +) 3/ So the
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 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 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 information8.1 Introduction. Consider the initial value problem (IVP):
8.1 Introduction Consider the initial value problem (IVP): y dy dt = f(t, y), y(t 0)=y 0, t 0 t T. Geometrically: solutions are a one parameter family of curves y = y(t) in(t, y)-plane. Assume solution
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 informationProblems on Numerical Methods for Engineers. Pedro Fortuny Ayuso. EPIG, Gijón. Universidad de Oviedo address:
Problems on Numerical Methods for Engineers Pedro Fortuny Ayuso EPIG, Gijón. Universidad de Oviedo E-mail address: fortunypedro@uniovi.es CC BY: Copyright c 2011 2015 Pedro Fortuny Ayuso This work is licensed
More informationX. Numerical Methods
X. Numerical Methods. Taylor Approximation Suppose that f is a function defined in a neighborhood of a point c, and suppose that f has derivatives of all orders near c. In section 5 of chapter 9 we introduced
More information