Notions of. Antonio Mucherino. last update: August University of Rennes 1 Numerical Analysis. A.
|
|
- Anna Lynch
- 5 years ago
- Views:
Transcription
1 Notions of Antonio Mucherino University of Rennes 1 last update: August 2013 commom s
2 ... commom s
3 Aircraft and wind Suppose an aircraft flies from Paris to Rio and then it comes back. Suppose is constant during the whole travel and it is able to influence the speed of the aircraft. Paris - Rio, time t 1 = 5.1 hours, aircraft flying against wind Rio - Paris, time t 2 = 4.7 hours, aircraft flying with wind distance: 5700 miles... How can we find the average speed of the aircraft and the average speed of? commom s
4 How to solve this problem? Let x be the average speed of the aircraft, and let y be the average speed of : the actual aircraft speed is x y when it flies against the actual aircraft speed is x + y when it flies with the distance d for each travel can be computed as the product between the time (t 1 or t 2 ) and the actual speed We can define the following system of equations: { t1 (x y) = d t 2 (x + y) = d... commom s This is a linear system: t 1, t 2 and d are parameters (already known) x and y are variables
5 General form of a linear system with 2 equations: { a11 x + a 12 y = b 1 a 21 x + a 22 y = b 2... commom s And, in matrix form: ( a11 a 12 a 21 a 22 ) ( x y The coefficient matrix about the existence of solutions about the number of solutions (out of the scope of this course). ) = ( b1 b 2 ) ( ) a11 a 12 is able to provide information a 21 a 22
6 The solution for our example For the system { t1 x t 1 y = d t 2 x + t 2 y = d we can find a solution (x, y) analytically:... x = d t 1 + y y = d ( 1 t ) 2 2t 2 t 1 commom s
7 The solution for our example So, in our example: t 1 = 5.1 t 2 = 4.7 d = 5700 and therefore: x = miles/hours... y = 47.6 miles/hours Can we program a computer to make this work for us? commom s Note that, in this simple example, we did not consider the health rotation.
8 Back substitution... Suppose the coefficient matrix of our linear system is an upper triangular matrix: a 11 a 12 a 13 a 14 0 a 22 a 23 a a 33 a a 44 In this situation, we can compute: x 4 = b 4 a 44 x 1 x 2 x 3 x 4 = b 1 b 2 b 3 b 4 commom s
9 Back substitution... Suppose the coefficient matrix of our linear system is an upper triangular matrix: a 11 a 12 a 13 a 14 0 a 22 a 23 a a 33 a a 44 In this situation, we can compute: x 1 x 2 x 3 x 4 x 3 = b 3 a 34 x 4 a 33 = b 1 b 2 b 3 b 4 commom s
10 Back substitution... Suppose the coefficient matrix of our linear system is an upper triangular matrix: a 11 a 12 a 13 a 14 0 a 22 a 23 a a 33 a a 44 In this situation, we can compute: x 1 x 2 x 3 x 4 = x 2 = b 2 a 23 x 3 a 24 x 4 a 22 b 1 b 2 b 3 b 4 commom s
11 Back substitution... Suppose the coefficient matrix of our linear system is an upper triangular matrix: a 11 a 12 a 13 a 14 0 a 22 a 23 a a 33 a a 44 In this situation, we can compute: x 1 x 2 x 3 x 4 = x 1 = b 1 a 12 x 2 a 13 x 3 a 14 x 4 a 11 b 1 b 2 b 3 b 4 commom s
12 C function for back substitution... commom s void back(int n,double **a,double *x) { // n is the system dimension // a is the coefficient matrix // (must be upper triangular) // x is // the vector of known terms (input) // the solution (output) }; int i,j; for (i = n - 1; i >= 0; i--) { for (j = i+1; j < n; j++) { x[i] = x[i] - a[i][j]*x[j]; }; x[i] = x[i]/a[i][i]; };
13 Gaussian elimination... How to solve linear whose coefficient matrix is not in triangular form? Gaussian elimination : transform the system in an equivalent system whose coefficient matrix is in triangular form. Example: 2x +y z = 8 3x y +2z = 11 2x +y +2z = 3 = 2x +y z = y z = 1 z = 1 commom s
14 LAPACK LAPACK Algebra PACKage... free library for linear algebra (including linear ) it s a freely-available software package (library + sources) originally developed in Fortran, there are versions for C and C++ it s based on another library called BLAS, which contains for efficient matrix manipulations (sums, products,... ) commom s
15 Some references... Wikipedia page about linear, Online solution of linear, LAPACK, BLAS, commom s
16 ... commom s
17 Stationary points In many applications, stationary points of are of particular interest.... They might provide: the minimum and maximum points of the equilibrium point of a dynamic system (may be stable or not)... In order to find a stationary point, the derivative of a function must be computed, and roots of such a derivative must be identified: df(x) dx = 0 commom s
18 Given a function f : [a, b] Y, how to find its roots (zeros)? f(x) = 0 Examples: ax = b = x = b a... ax 2 + bx + c = 0 = x 1 = b b 2 4ac 2a x 2 = b + b 2 4ac 2a commom s
19 A simple Simplest for finding roots: Extract a predefined number of points from the function domain [a, b] and evaluate the function in all these points. One of these points can be a root (or be close to a root).... commom s This is not able to provide a good approximation of roots of having a more complex shape.
20 Bisection: basic idea... commom s is an iterative which defines a sequence of intervals {a k, b k } k=1,2,...,itmax converging to one function root. At the beginning, the whole function domain is considered: [a 0, b 0 ] = [a, b] At each iteration, the following two steps are performed: the average point of the current interval is computed: x k = a k + b k a k 2 the new interval is then defined as: { [ak+1, b k+1 ] = [a k, x k ] if f(a k )f(x k ) 0 [a k+1, b k+1 ] = [x k, b k ] otherwise
21 Bisection: basic idea In the bisection, intervals [a k, b k ] are reduced in size at each iteration, and they are supposed to converge to a function root.... commom s
22 Bisection: basic idea In the bisection, intervals [a k, b k ] are reduced in size at each iteration, and they are supposed to converge to a function root.... commom s
23 Bisection: basic idea In the bisection, intervals [a k, b k ] are reduced in size at each iteration, and they are supposed to converge to a function root.... commom s
24 Applicability... commom s can be applied to f : [a, b] Y : if all points in [a, b] can be evaluated if f is a continuous function Moreover, if at least one of the intervals [a k, b k ] is such that f(a k )f(b k ) < 0 then, the converges toward one of the roots contained in the interval. The can be stopped when a k b k < ε or f(a k ) f(b k ) < ε where ε is a small real number (tolerance).
25 C function for the bisection algorithm... commom s double bisection(double a,double b,double (*f)(double),double eps,int itmax) { // [a,b], function domain // double (*f)(double), pointer to a function // eps, tolerance // itmax, maximum number of iterations int it; double ca,cb,cx,fa,fb,fx; ca = a; cb = b; fa = f(a); fb = f(b); cx = (ca+cb)/2.0; fx = f(cx); it = 0; while (it <= itmax && fabs(fx) > eps && fabs(cb-ca) > eps && fabs(fb-fa) > eps) { it = it + 1; if (fa*fx < 0) { cb = cx; fb = fx; } else { ca = cx; fa = fx; }; cx = (ca+cb)/2.0; fx = f(cx); }; return cx; };
26 Pointers to A pointer to a function can make reference to any function of a predefined type: double (*f)(double)... commom s In the main function: double xcube(double x); double polynomial(double x); double bisection(double a,double b,double (*f)(double),double eps,int itmax); main() { int i,j; double a,b,root;... double *f(double);... }; f = xcube; root = bisection(a,b,f,0.001,100);... f = polynomial; root = bisection(a,b,f,0.001,100);...
27 Other s for finding roots... Newton s it is based on the computation of the tangent to the function in the current root approximation Secant similar to the Newton s, but the tangent is replaced by a secant (the function does not have to be differentiable in the whole domain) Lehmer-Schur extension of the bisection Brent s combination of different s, including the bisection, with the aim of speeding up the search commom s
28 Polynomial interpolation... commom s
29 ... The for ammonia reacting in hydrogen and nitrogen gases depends upon the hydrogen-nitrogen mole ratio, the pressure, and the temperature. For a 3-to-1 hydrogen-nitrogen mole ratio, the equilibrium constant K p for a range of pressures and temperatures is given by: 100 atm 200 atm 300 atm 400 atm 500 atm 400 C C C C C Encyclopedia of Chemical Technology, vol. 2, 2 nd edition, New York, Wiley, commom s
30 ... Suppose that values for K p related to 500 C and 300 atm are not available, and that, for same reason, we cannot perform any experiment to find them. 100 atm 200 atm 300 atm 400 atm 500 atm 400 C C C C C How can we find the needed values for the constant K p? Easiest solution: linear interpolation. commom s
31 ... Suppose that values for K p related to 500 C and 300 atm are not available, and that, for same reason, we cannot perform any experiment to find them. 100 atm 200 atm 300 atm 400 atm 500 atm 400 C C C C C How can we find the needed values for the constant K p? Easiest solution: linear interpolation. commom s
32 ... Suppose that values for K p related to 500 C and 300 atm are not available, and that, for same reason, we cannot perform any experiment to find them. 100 atm 200 atm 300 atm 400 atm 500 atm 400 C C C C C How can we find the needed values for the constant K p? Easiest solution: linear interpolation. commom s
33 ... Suppose that values for K p related to 500 C and 300 atm are not available, and that, for same reason, we cannot perform any experiment to find them. 100 atm 200 atm 300 atm 400 atm 500 atm 400 C C C C C How can we find the needed values for the constant K p? Easiest solution: linear interpolation. commom s
34 ... Suppose that values for K p related to 500 C and 300 atm are not available, and that, for same reason, we cannot perform any experiment to find them. 100 atm 200 atm 300 atm 400 atm 500 atm 400 C C C C C How can we find the needed values for the constant K p? Easiest solution: linear interpolation. commom s
35 ... Suppose that values for K p related to 500 C and 300 atm are not available, and that, for same reason, we cannot perform any experiment to find them. 100 atm 200 atm 300 atm 400 atm 500 atm 400 C C C C C How can we find the needed values for the constant K p? Easiest solution: linear interpolation. commom s
36 interpolation Let f : [a, b] Y be a function such that the pair (x 1,f(x 1 )) is known, with x 1 [a, b] the pair (x 2,f(x 2 )) is known, with x 2 [a, b] and x 2 > x 1 f(x) is not known for any x (x 1, x 2 )... commom s interpolation: assign to the interval (x 1, x 2 ) of f(x) the equation of the line between x 1 and x 2.
37 interpolation Let f : [a, b] Y be a function such that the pair (x 1,f(x 1 )) is known, with x 1 [a, b] the pair (x 2,f(x 2 )) is known, with x 2 [a, b] and x 2 > x 1 f(x) is not known for any x (x 1, x 2 )... commom s interpolation: assign to the interval (x 1, x 2 ) of f(x) the equation of the line between x 1 and x 2.
38 interpolation Let f : [a, b] Y be a function such that the pair (x 1,f(x 1 )) is known, with x 1 [a, b] the pair (x 2,f(x 2 )) is known, with x 2 [a, b] and x 2 > x 1 f(x) is not known for any x (x 1, x 2 )... commom s interpolation: assign to the interval (x 1, x 2 ) of f(x) the equation of the line between x 1 and x 2.
39 Quadratic interpolation Let f : [a, b] Y be a function such that the pair (x 1,f(x 1 )) is known, with x 1 [a, b] the pair (x 2,f(x 2 )) is known, with x 2 [a, b] and x 2 > x 1 the pair (x 3,f(x 3 )) is known, with x 3 [a, b] and x 3 > x 2 > x 1 f(x) is not known for any x [a, b] \ {x 1, x 2, x 3 }... commom s Quadratic interpolation: assign to the interval (x 1, x 3 ) of f(x) the equation of the parabola passing through x 1, x 2 and x 3.
40 Quadratic interpolation Let f : [a, b] Y be a function such that the pair (x 1,f(x 1 )) is known, with x 1 [a, b] the pair (x 2,f(x 2 )) is known, with x 2 [a, b] and x 2 > x 1 the pair (x 3,f(x 3 )) is known, with x 3 [a, b] and x 3 > x 2 > x 1 f(x) is not known for any x [a, b] \ {x 1, x 2, x 3 }... commom s Quadratic interpolation: assign to the interval (x 1, x 3 ) of f(x) the equation of the parabola passing through x 1, x 2 and x 3.
41 Quadratic interpolation Let f : [a, b] Y be a function such that the pair (x 1,f(x 1 )) is known, with x 1 [a, b] the pair (x 2,f(x 2 )) is known, with x 2 [a, b] and x 2 > x 1 the pair (x 3,f(x 3 )) is known, with x 3 [a, b] and x 3 > x 2 > x 1 f(x) is not known for any x [a, b] \ {x 1, x 2, x 3 }... commom s Quadratic interpolation: assign to the interval (x 1, x 3 ) of f(x) the equation of the parabola passing through x 1, x 2 and x 3.
42 Lagrangian interpolation... Let f : [a, b] Y be a function such that the pairs (x i,f(x i )) are known, with x i {x 1, x 2,...,x n } [a, b] f(x) is not known for any x [a, b] \ {x 1, x 2,...,x n } Lagrangian interpolation: assign to the interval (x 1, x n ) of f(x) the equation of the polynomial of degree n 1 passing through the n points x 1, x 2,... x n. commom s
43 Lagrangian interpolation A general polynomial of degree n 1 can be written as:... f(x) = a n 1 x n 1 + a n 2 x n a 2 x 2 + a 1 x + a 0 For the polynomial to pass through the n points (x i, y i ), we need to solve the following system of linear equations: y 1 = a n 1 x n a n 2 x n a 2 x1 2 + a 1x 1 + a 0 y 2 = a n 1 x n a n 2 x n a 2 x2 2 + a 1x 2 + a 0 y 3 = a n 1 x n a n 2 x n a 2 x3 2 + a 1x 3 + a 0... y n = a n 1 xn n 1 + a n 2 xn n a 2 xn 2 + a 1 x n + a 0 commom s
44 Lagrangian interpolation It can be proved that the system of linear equations has only one solution: {a 0, a 1, a 2,..., a n 1 } the polynomial of degree n 1 and having as coefficients the found a i s is such that: y i = f(x i ), i = 1, 2,..., n... commom s The general formula for Lagrangian interpolation is: f(x) = n n ( ) x xj y i i=1 j=1,j i x i x j
45 C function for interpolation... double interpol(int n,double *x,double *y,double p) { // n, number of available (x,y) // x, vector containing all x s // y, vector containing all y s // p, point where to evaluate lagrangian polynomial int i,j; double sum,prod; sum = 0.0; for (i=0; i<n; i++) { prod = 1.0; for (j=0; j<n; j++) { if (j!=i) { prod = prod * ((p-x[j])/(x[i]-x[j])); }; }; sum = sum + y[i]*prod; }; return sum; }; commom s
46 Regression Suppose that the form of f(x) is known a priori.... commom s If f(x) is linear, would the lagrangian polynomial be a good model? Solution: regression models.
47 Regression Suppose that the form of f(x) is known a priori.... commom s If f(x) is linear, would the lagrangian polynomial be a good model? No! Solution: regression models.
48 ... commom s
49 Archimedes Archimedes (287BC 212BC) was a Greek mathematician, physicist, engineer, inventor, and astronomer. He is generally considered to be the greatest mathematician of antiquity. Archimedes was able to approximate the area of a circle with polygons converging to the shape of the circle.... He was able to approximate the value of π to commom s
50 Solving the definite integral b a f(x)dx equals to finding the area under the curve y = f(x) and between x = a and x = b.... commom s We suppose: a and b are not ±, there are no points x [a, b] such that f( x) = ±.
51 Idea: approximate the area defined by f(x) with the area of the trapezoid ABCD:... commom s Area of trapezoid: 1 2 h (f(a) + f(b)). The area between y = f(x) and the segment DC corresponds to the error introduced with this approximation.
52 : divide [a, b] in n equal parts of length h and approximate each subinterval [x i, x i+1 ] with the area of the corresponding trapezoid:... commom s Trapezoidal formula: 1 n 2 h (f(x i ) + f(x i+1 )). i=1
53 : the C function... double trapez(double a,double b,int n, double (*f)(double)) { // interval [a,b] (input) // n, number of subintervals (input) // f, pointer to function (input) // returning value: approx. of the area defined by f(x) in [a,b] int i; double h,area; double ca,cb,fa,fb; h = (b - a)/n; ca = a; cb = ca + h; fa = f(ca); fb = f(cb); area = fa + fb; for (i = 1; i < n; i++) { ca = cb; fa = fb; cb = cb + h; fb = f(cb); area = area + fa + fb; }; return h*area/2.0; }; commom s
54 Other algorithms Simpson rule: instead of using trapezoids to approximate the areas, parabolas interpolating 3 consecutive points are employed. This simple modification increases the accuracy of the. Gaussian quadrature: subintervals of [a, b] do not have the same length but they are chosen so that the global accuracy increases.... W.S. Dorn, D.D. Mc Cracken, Methods with Fortran IV Case Studies, John Wiley & Sons, Inc., commom s
55 ... commom s
56 problems General form on an optimization problem:... min f(x) x A subject to a set of constraints: { x B g(x) = 0 x C h(x) 0 where f(x) is the objective function g(x) represents the equality constraints h(x) represents the inequality constraints commom s
57 Some s for optimization Deterministic s Simplex (may require some assumptions to be satisfied)... commom s Branch & Bound Branch & Prune... Heuristic s (no guarantees for optimality) Simulated Annealing Genetic Algorithms Tabu Search Variable Neighbourhood Search...
Jim 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 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 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 informationMS 2001: Test 1 B Solutions
MS 2001: Test 1 B Solutions Name: Student Number: Answer all questions. Marks may be lost if necessary work is not clearly shown. Remarks by me in italics and would not be required in a test - J.P. Question
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 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 informationRoot Finding (and Optimisation)
Root Finding (and Optimisation) M.Sc. in Mathematical Modelling & Scientific Computing, Practical Numerical Analysis Michaelmas Term 2018, Lecture 4 Root Finding The idea of root finding is simple we want
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 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 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 information1. Method 1: bisection. The bisection methods starts from two points a 0 and b 0 such that
Chapter 4 Nonlinear equations 4.1 Root finding Consider the problem of solving any nonlinear relation g(x) = h(x) in the real variable x. We rephrase this problem as one of finding the zero (root) of a
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 informationBACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES
No. of Printed Pages : 5 BCS-054 BACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 058b9 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES Time : 3 hours Maximum Marks
More informationSolving Non-Linear Equations (Root Finding)
Solving Non-Linear Equations (Root Finding) Root finding Methods What are root finding methods? Methods for determining a solution of an equation. Essentially finding a root of a function, that is, a zero
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 informationTHE SECANT METHOD. q(x) = a 0 + a 1 x. with
THE SECANT METHOD Newton s method was based on using the line tangent to the curve of y = f (x), with the point of tangency (x 0, f (x 0 )). When x 0 α, the graph of the tangent line is approximately the
More informationCHAPTER 4 ROOTS OF EQUATIONS
CHAPTER 4 ROOTS OF EQUATIONS Chapter 3 : TOPIC COVERS (ROOTS OF EQUATIONS) Definition of Root of Equations Bracketing Method Graphical Method Bisection Method False Position Method Open Method One-Point
More informationMATH 019: Final Review December 3, 2017
Name: MATH 019: Final Review December 3, 2017 1. Given f(x) = x 5, use the first or second derivative test to complete the following (a) Calculate f (x). If using the second derivative test, calculate
More informationUNCC 2001 Algebra II
UNCC 2001 Algebra II March 5, 2001 1. Compute the sum of the roots of x 2 5x + 6 = 0. (A) 3 (B) 7/2 (C) 4 (D) 9/2 (E) 5 (E) The sum of the roots of the quadratic ax 2 + bx + c = 0 is b/a which, for this
More informationEdexcel GCE Further Pure Mathematics (FP1) Required Knowledge Information Sheet. Daniel Hammocks
Edexcel GCE Further Pure Mathematics (FP1) Required Knowledge Information Sheet FP1 Formulae Given in Mathematical Formulae and Statistical Tables Booklet Summations o =1 2 = 1 + 12 + 1 6 o =1 3 = 1 64
More informationMathematics for Engineers. Numerical mathematics
Mathematics for Engineers Numerical mathematics Integers Determine the largest representable integer with the intmax command. intmax ans = int32 2147483647 2147483647+1 ans = 2.1475e+09 Remark The set
More informationRoot Finding: Close Methods. Bisection and False Position Dr. Marco A. Arocha Aug, 2014
Root Finding: Close Methods Bisection and False Position Dr. Marco A. Arocha Aug, 2014 1 Roots Given function f(x), we seek x values for which f(x)=0 Solution x is the root of the equation or zero of the
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 informationRoots of Equations. ITCS 4133/5133: Introduction to Numerical Methods 1 Roots of Equations
Roots of Equations Direct Search, Bisection Methods Regula Falsi, Secant Methods Newton-Raphson Method Zeros of Polynomials (Horner s, Muller s methods) EigenValue Analysis ITCS 4133/5133: Introduction
More informationNumerical Integration
Numerical Integration Sanzheng Qiao Department of Computing and Software McMaster University February, 2014 Outline 1 Introduction 2 Rectangle Rule 3 Trapezoid Rule 4 Error Estimates 5 Simpson s Rule 6
More informationSolution of Algebric & Transcendental Equations
Page15 Solution of Algebric & Transcendental Equations Contents: o Introduction o Evaluation of Polynomials by Horner s Method o Methods of solving non linear equations o Bracketing Methods o Bisection
More informationLab 11: Numerical Integration Techniques. Figure 1. From the Fundamental Theorem of Calculus, we know that if we want to calculate f ( x)
Lab 11: Numerical Integration Techniques Introduction The purpose of this laboratory experience is to develop fundamental methods for approximating the area under a curve for the definite integral. With
More informationNonlinear Equations and Continuous Optimization
Nonlinear Equations and Continuous Optimization Sanzheng Qiao Department of Computing and Software McMaster University March, 2014 Outline 1 Introduction 2 Bisection Method 3 Newton s Method 4 Systems
More informationDRAFT - Math 101 Lecture Note - Dr. Said Algarni
2 Limits 2.1 The Tangent Problems The word tangent is derived from the Latin word tangens, which means touching. A tangent line to a curve is a line that touches the curve and a secant line is a line that
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 informationMath Numerical Analysis Mid-Term Test Solutions
Math 400 - Numerical Analysis Mid-Term Test Solutions. Short Answers (a) A sufficient and necessary condition for the bisection method to find a root of f(x) on the interval [a,b] is f(a)f(b) < 0 or f(a)
More informationCollege Algebra. Basics to Theory of Equations. Chapter Goals and Assessment. John J. Schiller and Marie A. Wurster. Slide 1
College Algebra Basics to Theory of Equations Chapter Goals and Assessment John J. Schiller and Marie A. Wurster Slide 1 Chapter R Review of Basic Algebra The goal of this chapter is to make the transition
More informationNumerical techniques to solve equations
Programming for Applications in Geomatics, Physical Geography and Ecosystem Science (NGEN13) Numerical techniques to solve equations vaughan.phillips@nateko.lu.se Vaughan Phillips Associate Professor,
More information15 Nonlinear Equations and Zero-Finders
15 Nonlinear Equations and Zero-Finders This lecture describes several methods for the solution of nonlinear equations. In particular, we will discuss the computation of zeros of nonlinear functions f(x).
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 informationMidterm Review. Igor Yanovsky (Math 151A TA)
Midterm Review Igor Yanovsky (Math 5A TA) Root-Finding Methods Rootfinding methods are designed to find a zero of a function f, that is, to find a value of x such that f(x) =0 Bisection Method To apply
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 informationPART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435
PART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435 Professor Biswa Nath Datta Department of Mathematical Sciences Northern Illinois University DeKalb, IL. 60115 USA E mail: dattab@math.niu.edu
More informationChapter 2 Solutions of Equations of One Variable
Chapter 2 Solutions of Equations of One Variable 2.1 Bisection Method In this chapter we consider one of the most basic problems of numerical approximation, the root-finding problem. This process involves
More informationROOT FINDING REVIEW MICHELLE FENG
ROOT FINDING REVIEW MICHELLE FENG 1.1. Bisection Method. 1. Root Finding Methods (1) Very naive approach based on the Intermediate Value Theorem (2) You need to be looking in an interval with only one
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 informationp 1 p 0 (p 1, f(p 1 )) (p 0, f(p 0 )) The geometric construction of p 2 for the se- cant method.
80 CHAP. 2 SOLUTION OF NONLINEAR EQUATIONS f (x) = 0 y y = f(x) (p, 0) p 2 p 1 p 0 x (p 1, f(p 1 )) (p 0, f(p 0 )) The geometric construction of p 2 for the se- Figure 2.16 cant method. Secant Method The
More informationOrder of convergence
Order of convergence Linear and Quadratic Order of convergence Computing square root with Newton s Method Given a > 0, p def = a is positive root of equation Newton s Method p k+1 = p k p2 k a 2p k = 1
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 informationBisection and False Position Dr. Marco A. Arocha Aug, 2014
Bisection and False Position Dr. Marco A. Arocha Aug, 2014 1 Given function f, we seek x values for which f(x)=0 Solution x is the root of the equation or zero of the function f Problem is known as root
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 information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationHanoi Open Mathematical Competition 2017
Hanoi Open Mathematical Competition 2017 Junior Section Saturday, 4 March 2017 08h30-11h30 Important: Answer to all 15 questions. Write your answers on the answer sheets provided. For the multiple choice
More information3. Solve the following inequalities and express your answer in interval notation.
Youngstown State University College Algebra Final Exam Review (Math 50). Find all Real solutions for the following: a) x 2 + 5x = 6 b) 9 x2 x 8 = 0 c) (x 2) 2 = 6 d) 4x = 8 x 2 e) x 2 + 4x = 5 f) 36x 3
More informationNumerical Integration
Chapter 1 Numerical Integration In this chapter we examine a few basic numerical techniques to approximate a definite integral. You may recall some of this from Calculus I where we discussed the left,
More informationMath 473: Practice Problems for Test 1, Fall 2011, SOLUTIONS
Math 473: Practice Problems for Test 1, Fall 011, SOLUTIONS Show your work: 1. (a) Compute the Taylor polynomials P n (x) for f(x) = sin x and x 0 = 0. Solution: Compute f(x) = sin x, f (x) = cos x, f
More informationQuestion Bank (I scheme )
Question Bank (I scheme ) Name of subject: Applied Mathematics Subject code: 22206/22224/22210/22201 Course : CH/CM/CE/EJ/IF/EE/ME Semester: II UNIT-3 (CO3) Unit Test : II (APPLICATION OF INTEGRATION)
More informationDepartment of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
Department of Applied Mathematics and Theoretical Physics AMA 204 Numerical analysis Exam Winter 2004 The best six answers will be credited All questions carry equal marks Answer all parts of each question
More informationNONLINEAR EQUATIONS AND TAYLOR S THEOREM
APPENDIX C NONLINEAR EQUATIONS AND TAYLOR S THEOREM C.1 INTRODUCTION In adjustment computations it is frequently necessary to deal with nonlinear equations. For example, some observation equations relate
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationDerivatives and Shapes of Curves
MATH 1170 Section 43 Worksheet NAME Derivatives and Shapes of Curves In Section 42 we discussed how to find the extreme values of a function using the derivative These results say, In Chapter 2, we discussed
More informationMathematical Olympiad Training Polynomials
Mathematical Olympiad Training Polynomials Definition A polynomial over a ring R(Z, Q, R, C) in x is an expression of the form p(x) = a n x n + a n 1 x n 1 + + a 1 x + a 0, a i R, for 0 i n. If a n 0,
More information1.1: The bisection method. September 2017
(1/11) 1.1: The bisection method Solving nonlinear equations MA385/530 Numerical Analysis September 2017 3 2 f(x)= x 2 2 x axis 1 0 1 x [0] =a x [2] =1 x [3] =1.5 x [1] =b 2 0.5 0 0.5 1 1.5 2 2.5 1 Solving
More informationIntroduction Linear system Nonlinear equation Interpolation
Interpolation Interpolation is the process of estimating an intermediate value from a set of discrete or tabulated values. Suppose we have the following tabulated values: y y 0 y 1 y 2?? y 3 y 4 y 5 x
More informationAS PURE MATHS REVISION NOTES
AS PURE MATHS REVISION NOTES 1 SURDS A root such as 3 that cannot be written exactly as a fraction is IRRATIONAL An expression that involves irrational roots is in SURD FORM e.g. 2 3 3 + 2 and 3-2 are
More informationA nonlinear equation is any equation of the form. f(x) = 0. A nonlinear equation can have any number of solutions (finite, countable, uncountable)
Nonlinear equations Definition A nonlinear equation is any equation of the form where f is a nonlinear function. Nonlinear equations x 2 + x + 1 = 0 (f : R R) f(x) = 0 (x cos y, 2y sin x) = (0, 0) (f :
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 information2007 Hypatia Contest
Canadian Mathematics Competition An activity of the Centre for Education in Mathematics and Computing, University of Waterloo, Waterloo, Ontario 007 Hypatia Contest Wednesday, April 18, 007 Solutions c
More informationLesson 9 Exploring Graphs of Quadratic Functions
Exploring Graphs of Quadratic Functions Graph the following system of linear inequalities: { y > 1 2 x 5 3x + 2y 14 a What are three points that are solutions to the system of inequalities? b Is the point
More informationNonlinear equations. Norms for R n. Convergence orders for iterative methods
Nonlinear equations Norms for R n Assume that X is a vector space. A norm is a mapping X R with x such that for all x, y X, α R x = = x = αx = α x x + y x + y We define the following norms on the vector
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 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 informationTest 2 Review Math 1111 College Algebra
Test 2 Review Math 1111 College Algebra 1. Begin by graphing the standard quadratic function f(x) = x 2. Then use transformations of this graph to graph the given function. g(x) = x 2 + 2 *a. b. c. d.
More informationSec 4 Maths. SET A PAPER 2 Question
S4 Maths Set A Paper Question Sec 4 Maths Exam papers with worked solutions SET A PAPER Question Compiled by THE MATHS CAFE 1 P a g e Answer all the questions S4 Maths Set A Paper Question Write in dark
More informationax 2 + bx + c = 0 where
Chapter P Prerequisites Section P.1 Real Numbers Real numbers The set of numbers formed by joining the set of rational numbers and the set of irrational numbers. Real number line A line used to graphically
More information1.1 Basic Algebra. 1.2 Equations and Inequalities. 1.3 Systems of Equations
1. Algebra 1.1 Basic Algebra 1.2 Equations and Inequalities 1.3 Systems of Equations 1.1 Basic Algebra 1.1.1 Algebraic Operations 1.1.2 Factoring and Expanding Polynomials 1.1.3 Introduction to Exponentials
More informationAnnouncements. Topics: Homework: - sections , 6.1 (extreme values) * Read these sections and study solved examples in your textbook!
Announcements Topics: - sections 5.2 5.7, 6.1 (extreme values) * Read these sections and study solved examples in your textbook! Homework: - review lecture notes thoroughly - work on practice problems
More informationAPPM/MATH Problem Set 6 Solutions
APPM/MATH 460 Problem Set 6 Solutions This assignment is due by 4pm on Wednesday, November 6th You may either turn it in to me in class or in the box outside my office door (ECOT ) Minimal credit will
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 informationSolution Methods. Richard Lusby. Department of Management Engineering Technical University of Denmark
Solution Methods Richard Lusby Department of Management Engineering Technical University of Denmark Lecture Overview (jg Unconstrained Several Variables Quadratic Programming Separable Programming SUMT
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 informationNonlinear Equations. Chapter The Bisection Method
Chapter 6 Nonlinear Equations Given a nonlinear function f(), a value r such that f(r) = 0, is called a root or a zero of f() For eample, for f() = e 016064, Fig?? gives the set of points satisfying y
More informationAdditional exercises with Numerieke Analyse
Additional exercises with Numerieke Analyse March 10, 017 1. (a) Given different points x 0, x 1, x [a, b] and scalars y 0, y 1, y, z 1, show that there exists at most one polynomial p P 3 with p(x i )
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 information19.4 Spline Interpolation
c9-b.qxd 9/6/5 6:4 PM Page 8 8 CHAP. 9 Numerics in General 9.4 Spline Interpolation Given data (function values, points in the xy-plane) (x, ƒ ), (x, ƒ ),, (x n, ƒ n ) can be interpolated by a polynomial
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationAlgebra II Vocabulary Word Wall Cards
Algebra II Vocabulary Word Wall Cards Mathematics vocabulary word wall cards provide a display of mathematics content words and associated visual cues to assist in vocabulary development. The cards should
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 5. Nonlinear Equations
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T Heath Chapter 5 Nonlinear Equations Copyright c 2001 Reproduction permitted only for noncommercial, educational
More informationNumerical Methods for Engineers
Numerical Methods for Engineers SEVENTH EDITION Steven C Chopra Berger Chair in Computing and Engineering Tufts University Raymond P. Canal Professor Emeritus of Civil Engineering of Michiaan University
More informationOptimization. Charles J. Geyer School of Statistics University of Minnesota. Stat 8054 Lecture Notes
Optimization Charles J. Geyer School of Statistics University of Minnesota Stat 8054 Lecture Notes 1 One-Dimensional Optimization Look at a graph. Grid search. 2 One-Dimensional Zero Finding Zero finding
More informationOptimization. Totally not complete this is...don't use it yet...
Optimization Totally not complete this is...don't use it yet... Bisection? Doing a root method is akin to doing a optimization method, but bi-section would not be an effective method - can detect sign
More informationQR Decomposition. When solving an overdetermined system by projection (or a least squares solution) often the following method is used:
(In practice not Gram-Schmidt, but another process Householder Transformations are used.) QR Decomposition When solving an overdetermined system by projection (or a least squares solution) often the following
More informationRoundoff Error. Monday, August 29, 11
Roundoff Error A round-off error (rounding error), is the difference between the calculated approximation of a number and its exact mathematical value. Numerical analysis specifically tries to estimate
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 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 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 information8.8. Applications of Taylor Polynomials. Infinite Sequences and Series 8
8.8 Applications of Taylor Polynomials Infinite Sequences and Series 8 Applications of Taylor Polynomials In this section we explore two types of applications of Taylor polynomials. First we look at how
More information. (a) Express [ ] as a non-trivial linear combination of u = [ ], v = [ ] and w =[ ], if possible. Otherwise, give your comments. (b) Express +8x+9x a
TE Linear Algebra and Numerical Methods Tutorial Set : Two Hours. (a) Show that the product AA T is a symmetric matrix. (b) Show that any square matrix A can be written as the sum of a symmetric matrix
More informationCS412: Introduction to Numerical Methods
CS412: Introduction to Numerical Methods MIDTERM #1 2:30PM - 3:45PM, Tuesday, 03/10/2015 Instructions: This exam is a closed book and closed notes exam, i.e., you are not allowed to consult any textbook,
More informationComputational Methods CMSC/AMSC/MAPL 460. Solving nonlinear equations and zero finding. Finding zeroes of functions
Computational Methods CMSC/AMSC/MAPL 460 Solving nonlinear equations and zero finding Ramani Duraiswami, Dept. of Computer Science Where does it arise? Finding zeroes of functions Solving functional equations
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 5 Nonlinear Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction
More informationLINEAR SYSTEMS (11) Intensive Computation
LINEAR SYSTEMS () Intensive Computation 27-8 prof. Annalisa Massini Viviana Arrigoni EXACT METHODS:. GAUSSIAN ELIMINATION. 2. CHOLESKY DECOMPOSITION. ITERATIVE METHODS:. JACOBI. 2. GAUSS-SEIDEL 2 CHOLESKY
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 informationNumerical Methods. Roots of Equations
Roots of Equations by Norhayati Rosli & Nadirah Mohd Nasir Faculty of Industrial Sciences & Technology norhayati@ump.edu.my, nadirah@ump.edu.my Description AIMS This chapter is aimed to compute the root(s)
More informationCHAPTER-II ROOTS OF EQUATIONS
CHAPTER-II ROOTS OF EQUATIONS 2.1 Introduction The roots or zeros of equations can be simply defined as the values of x that makes f(x) =0. There are many ways to solve for roots of equations. For some
More information