SME 3023 Applied Numerical Methods
|
|
- Charleen Hicks
- 6 years ago
- Views:
Transcription
1 UNIVERSITI TEKNOLOGI MALAYSIA SME 3023 Applied Numerical Methods Solution of Nonlinear Equations Abu Hasan Abdullah Faculty of Mechanical Engineering Sept 2012 Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
2 Outline 1 Introduction Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
3 Outline 1 Introduction 2 Engineering Applications Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
4 Outline 1 Introduction 2 Engineering Applications 3 Methods Available Incremental Search Method Bisection Method Newton-Raphson Method Secant Method Fixed Point Iteration Method Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
5 Outline 1 Introduction 2 Engineering Applications 3 Methods Available Incremental Search Method Bisection Method Newton-Raphson Method Secant Method Fixed Point Iteration Method 4 Root Finding with Matlab Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
6 Outline 1 Introduction 2 Engineering Applications 3 Methods Available Incremental Search Method Bisection Method Newton-Raphson Method Secant Method Fixed Point Iteration Method 4 Root Finding with Matlab 5 Roots of Nonlinear Polynomials Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
7 Outline 1 Introduction 2 Engineering Applications 3 Methods Available Incremental Search Method Bisection Method Newton-Raphson Method Secant Method Fixed Point Iteration Method 4 Root Finding with Matlab 5 Roots of Nonlinear Polynomials 6 Bibliography Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
8 Introduction Many engineering analyses require determination of value(s) of variable x that satisfy a nonlinear equation f(x) = 0 (1) where x is known as roots of Eq. (1) or zeros of function f(x). Number of roots maybe finite or infinite depending on nature of problem and physical problem Examples of f(x) are x 4 80x = 0 tan x tanh x = 0 polynomial transcendental equation Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
9 Engineering Applications Example Problem 1 Figure 1 : Discharge of water from reservoir. Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
10 Engineering Applications Example Problem 1 Problem Statement: Water is discharge from a reservoir through a long pipe as shown in Figure 1.By neglecting the change in the level of the reservoir, the transient velocity of the water flowing from pipe, v(t), can be expressed as: v(t) t p p = tanh 2gh 2gh 2L where h is the height of the fluid in the reservoir, L is the length of the pipe, g is the acceleration due to gravity, and t is the time elapsed from the beginning of the flow. Find the value of h necessary for achieving a velocity of v = 5 m/s at time t = 3 s when L = 5 m and g = 9.81 m/s 2. Solution: Work through the example see Rao (2002). Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
11 Engineering Applications Example Problem 2 Figure 2 : Open belt drive. Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
12 Engineering Applications Example Problem 2 Problem Statement: The length of a belt in an open-belt drive, L, is given by L = p 4c 2 (D d) `DθD + dθ d (E1) where θ D = π + 2 sin 1 D d 2c ««θ d = π 2 sin 1 D d 2c (E2,E3) c is the centre distance, D is the diameter of the larger pulley, d is the diameter of the smaller pulley, θ D is the angle of contact of the belt with the larger pulley, and θ d is the angle of contact of the belt with the smaller pulley, Figure 2. If a belt having a length 11 m is used to connect the two pulleys with diameters 0.4 m and 0.2 m, determine the centre distance between the pulleys. Solution: Work through the example see Rao (2002). Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
13 Engineering Applications Example Problem 3 Figure 3 : Contact stress between spheres. Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
14 Engineering Applications Example Problem 3 Problem Statement: The shear stress induced along the z-axis when two spheres are in contact with each other, while carrying a load F, is given by where h(λ) = λ tan λ2 1 λ «0.65 h(λ) = τzx p max and λ = z a in which τ zx is the shear force, (E1) p max = 3F 2πa 2 is the maximum pressure deveoped at the centre of the contact area, and (E2) Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
15 Engineering Applications Example Problem 3 the radius of the contact area, Figure 3, is 8 1 >< + 1 «91/3 E 1 E >= 2 a = F 1 >: + 1 «>; d 1 d 2 (E3) where E 1 and E 2 are Young s moduli of the two spheres, and d 1 and d 2 are diameters of the two spheres. Poisson s ratios of the two spheres was assume to be 0.3 in deriving Eqs. (E1) and (E3). Determine the value of λ at which the shear stress, given by Eq. (E1), attains its maximum value. Solution: Work through the example see Rao (2002). Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
16 Methods Available Incremental Search Method Bisection Method Newton-Raphson Method Secant Method Regula Falsi Method Fixed Point Iteration Method Bairstow Method Muller s Method Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
17 Methods Available Incremental Search Method Value of x is incremented, by x, from an initial value, x 1, successively until a change in the sign of the function f(x) is observed. f(x) changes sign between x i and x i+1, if it has root in the interval [x i, x i+1 ] which implies f(x i ) f(x i+1 ) < 0 wherever a root is crossed. Plot of the function is usually very useful in guiding the task of finding the interval. A potential problem is the choice of increment length, x: too small, the search can be very time-consuming, too great, closely spaced roots might be missed. Algorithm 1 Sets an initial guess for x i=1, and a stepsize x. 2 Call the function f(x) to calculate its value at x i=1. 3 Increment x i+1 = x i + x. 4 Call the function f(x) to calculate its value at x i+1. 5 Compares the sign of the returned function value f(x i+1 ) to the previous f(x i ). * If the sign of f(x i+1 ) does NOT change, repeat from Step 3 again. * If the sign of f(x i+1 ) does change, the root lies between x i and x i+1. Reduce stepsize x, set x i = x i 1 and repeat from Step 3. Iterate until f(x i+1 ) f(x i ) ε where ε is a specified very small number called tolerance. Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
18 Methods Available Incremental Search Method Example 1 Problem Statement: Find the root of the equation f(x) = 1.5x 0.65 (1 + x 2 tan 1 ) 2 « x x 1 + x = 0 (E1) 2 using the incremental search method with x 1 = 0.0 and x (1) = 0.1. Solution: Work through the example. Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
19 Methods Available Bisection Method If f(x) is real and continuous in the interval prescribed by a lower bound, x L and upper bound, x U and f(x L) and f(x U) have opposite signs, such that f(x L) f(x U) < 0 (2) then there is at least one real root in the interval [x L, x U] Algorithm 1 Choose lower x L and upper x U guesses for the root such that the function changes sign over the interval. This can be checked by ensuring that f(x L) f(x U) < 0 2 An estimate of the root x R is determined by xl + xu x R = 2 3 Make the following evaluations to determine in which subinterval the root lies: * if f(x L) f(x R) = 0, the root equals x R. Terminate computation. * if f(x L) f(x R) < 0, the root lies in the lower subinterval. Therefore set x U = x R and return to step 2 * if f(x L) f(x R) > 0, the root lies in the upper subinterval. Therefore set x L = x R and return to step 2 Iterate until f(x R) ε where ε is a specified very small number called tolerance. Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
20 Methods Available Bisection Method Example 1 Problem Statement: Find the root of the equation f(x) = 1.5x 0.65 (1 + x 2 tan 1 ) 2 « x x 1 + x = 0 (E1) 2 using the bisection method with x (1) l = 0.0, x (1) u = 2.0 and ε = Solution: Modify the Fortran code, ܽ ½º, and Matlab code, غÑ, to solve Eq. (E1) above. Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
21 Methods Available Newton-Raphson Method Function f(x) is expressed using Taylor s series about an arbitrary point, x 1 f(x) = f(x 1) + (x x 1)f (x 1) + 1 2! (x x1)2 f (x 1) +... (3) where f, f, f,... are evaluated at x 1 Consider only the first two terms in the expansion f(x) = f(x 1) + (x x 1)f (x 1) = 0 (4) and set f(x) = 0 to give f(x 1) + (x x 1)f (x 1) = 0 (5) Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
22 Methods Available Newton-Raphson Method Since higher order derivative terms were neglected in the approximation of f(x) in Eq. (4), the solution of Eq. (5) yields the next approximation to the root (instead of the exact root) as x = x 2 = x 1 f(x1) f (x 1) (6) x 2 denotes improved approximation to the root. For the next improvement we use x 2 in place of x 1 on the RHS of Eq. (6) to obtain x 3 The iterative procedure of Newton-Raphson method is generalized as x i+1 = x i f(x i) : i = 1, 2,... (7) f (x i ) Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
23 Methods Available Newton-Raphson Method Example 1 Problem Statement: Find the root of the equation f(x) = 1.5x 0.65 (1 + x 2 tan 1 ) 2 « x x 1 + x = 0 (E1) 2 using the Newton-Raphson method with the starting point x 1 = 0.0 and the convergence criteria ε = Solution: Modify the Fortran code, ܽ ¾º, and Matlab code, Ò ÛØÓÒºÑ, to solve Eq. (E1) above. Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
24 Methods Available Secant Method Similar to Newton s method but is different in that the derivative f is approximated using two consecutive iterative values of f f (x i ) f(x i) f(x i 1 ) x i x i 1 (8) and Eq. (7) can be re-written as x i+1 = x i f(x i) f (x i ) = x i f(x i )» x i x i 1 f(x i ) f(x i 1 ) (9) Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
25 Methods Available Fixed Point Iteration Method Also known as Successive Substitution Method, in which function f(x) = 0 is re-written in the form of x = g(x) and iterative procedure is used on x i+1 = g(x i ) : i = 1, 2, 3,... (10) where new approximation to root, x i+1, is found using the previous one, x i. Iterative process is stopped with a convergence criterion 10 3 < ε < 10 6 x i+1 g(x i+1 ) ε (11) This method is simple but may NOT always converge with an arbitrarily chosen form of the function g(x) but condition to be satisfied for convergence to the correct root is given by g (x i+1 ) < 1 (12) Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
26 Methods Available Fixed Point Iteration Method Example 1 Problem Statement: Find the root of the equation f(x) = 1.5x (1 + x tan 1 ) 2 « x x 1 + x 2 = 0 (E1) using the fixed point iteration with the starting point x 1 = 0.0 and the convergence criteria ε = Solution: Re-arrange Eq. (E1) 1.5x (1 + x 2 = 0.65 tan 1 ) 2 «1 0.65x x 1 + x 2 = 0 1 «13 30 x(1 + x2 ) x = (1 + x2 ) 2 tan 1 x {z } g(x) so the RHS is g(x). Thus, from Eq. (E2), the iterative process can be expressed as x i+1 = 13 «1 30 (1 + x2 i )2 tan 1 13 x i 30 x i(1 + xi 2 ) (E2) Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
27 Root Finding with Matlab Example 1 Problem Statement: Find the root of the equation f(x) = tan 1 x = 0 using Matlab. Solution: Matlab Code Ü ¹¾ Ô ¼º¼½ ¾ Ô Ü ³ Ø Ò Üµ³ ÔÐÓØ Ü ¹¾ Ô ¾ Ô µ ÜÐ Ð ³Ü³µ ÝÐ Ð ³Ý ܵ³µ Ø ØÐ ³ÊÓÓØ Ò Ò ³µ ÖÓÓØ Þ ÖÓ Ü ½º µ Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
28 Root Finding with Matlab Example 2 Problem Statement: Find the root of the equation f(x) = sin 2x = 0 using Matlab. Solution: Matlab Code Ü ¹¾ Ô ¼º¼½ ¾ Ô Ü ³ Ò ¾ ܵ³ ÔÐÓØ Ü ¹¾ Ô ¾ Ô µ ÜÐ Ð ³Ü³µ ÝÐ Ð ³Ý ܵ³µ Ø ØÐ ³ÊÓÓØ Ò Ò ³µ ÖÓÓØ Þ ÖÓ Ü ¼º µ Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
29 Root Finding with Matlab Example 3 Problem Statement: Find the root of the equation f(x) = 1.5x 0.65 (1 + x 2 tan 1 ) 2 using Matlab. Solution: Matlab Code « x x 1 + x = 0 2 Ü ¹¾º¼ ¼º¼½ ¾ Ü ³½º ܺ» ½ ܺ ¾µº ¾¹¼º Ø Ò ½º»Üµµ ¼º ܺ» ½ ܺ ¾µ³ ÔÐÓØ Ü ¹¾ ¾ µ ÜÐ Ð ³Ü³µ ÝÐ Ð ³Ý ܵ³µ Ø ØÐ ³ÊÓÓØ Ò Ò ³µ ܾ Þ ÖÓ Ü ¹¼º½µ Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
30 Roots of Nonlinear Polynomials There is theorem which states that there is no such formula for general polynomials of degree higher than four. In practice we use numerical method to solve polynomial equations of degree higher than two. As a special case of f(x) = 0, consider a polynomial equation f(x) = a nx n + a n 1x n a 2x 2 + a 1x + a 0 = 0 (13) where n denote degree of polynomial, a 0, a 1, a 2... a n are coefficients of polynomial. Eq. (13) in general will have n roots Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
31 Roots of Nonlinear Polynomials If x 1, x 2... x n are roots of Eq. (13), they are related to the coefficients of polynomial as nx i=1 nx x i nx i=1 j=1,j 1 nx x i nx nx i=1 j=1,j i k=1,k j x i = a n 1 a n x i x j = a n 2 a n x i x j x k = a n 3 a n... x 1 x 2... x n 1 x n = ( 1) n a 0 a n (14) Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
32 Roots of Nonlinear Polynomials Muller s Method Your homework! Read Section 2.11, pp of SINGIRESU S. RAO (2002): Applied Numerical Methods for Engineers and Scientists, ISBN X, Prentice Hall Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
33 Roots of Nonlinear Polynomials Example 1 Problem Statement: Use Matlab to determine the roots of the following polynomial f(x) = x 4 12x 3 + 0x x = 0 Solution: Matlab Code Ô ½ ¹½¾ ¼ ¾ ½½ Ö ÖÓÓØ Ôµ ÔÔ ÔÓÐÝ Öµ ± Ò Ó ÒØ ØÓ Ñ ØÖ Ü Ô ± Ò ÖÓÓØ Ò Ò Ø Ñ ØÓ Ö ± Ë ÓÙÐ Ö ØÙÖÒ ÔÔ Ôº Ì ± ÓÙ Ñ Ø Û ÒØ ØÓ ÔÐÓØ Ø ÔÓÐÝÒÓÑ Ð Ö Ø Ü Ð Ò Ô ¼ ½¾µ ± Ë Ø Ö Ò ÓÖ Ü Ý ÔÓÐÝÚ Ð Ô Üµ ± Ú ÐÙ Ø ÔÓÐÝÒÓÑ Ð Ø ÐÐ Ú ÐÙ Ó Ü Ü Ü Ü ¼ ± ¹ Ü ÔÐÓØ Ü Ü Ü Ü Ýµ ± ÈÐÓØ Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
34 Roots of Nonlinear Polynomials Example 2 Problem Statement: Use Matlab to determine the roots of the following polynomial f(x) = x 5 3.5x x x x which has three real roots: 0.5, 1.0 and 2.0, and a pair of complex roots: 1 ± 0.5i. Solution: Matlab Code ±± ÒØ Ö Ó ÒØ Ó ÔÓÐÝÒÓÑ Ð ÒØÓ Ñ ØÖ Ü Ô ½ ¹ º ¾º ¾º½¾ ¹ º ½º¾ ±± ÓÙ Ñ Ý Ú ÐÙ Ø ÔÓÐÝÒÓÑ Ð Üµ Ý Ø Ü ½ ÓÖ Ö ÒØ Ø Øººº ÔÓÐÝÚ Ð Ô ½µ ÔÓÐÝ Ö Ôµ ±± Ö Ø ÕÙ Ö Ø ÕÒ ÖÓÑ ØÛÓ ÒÓÛÒ ÖÓÓØ ±± ܹ¼º µ Ü ½µ Ü ¾ ¼º Ü ¹ ¼º Ö ¼º ¹½º¼ Õ ÔÓÐÝ Öµ Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
35 Roots of Nonlinear Polynomials Example 2 Matlab Code (continued) ± Ú Ø ÓÖ Ò Ð ÔÓÐÝÒÓÑ Ð Üµ Ý Ø ÕÙ Ö Ø ± ÕÒ Ò Ò Ö ÙÐØ ÒØÓ ÕÙÓØ ÒØ Ò Ö Ñ Ò Ö ÓÒÚ Ô Õµ ÖÓÓØ µ ± ÊÓÓØ Ó ÕÙÓØ ÒØ ÔÓÐÝÒÓÑ Ð ÓÒÚ Õµ ± ÅÙÐØ ÔÐÝ Ý Õ ØÓ ÓÑ ÙÔ Û Ø Ø ÓÖ Ò Ð Üµ ± ÇÖ ÄÄ ÖÓÓØ Ñ Ý Ø ÖÑ Ò Ý Ö ÖÓÓØ Ôµ ± Ù ØÓ Û Ö ÓÑ Ó Ø ÖÓÓØ Ñ Ø ÝÓÙ Ñ Ý ± Û ÒØ ØÓ ÔÐÓØ ÔÓÐÝÒÓÑ Ð Üµ ÓÚ Ö Ö Ò Ý ¹½¼ Ü ½¼ Ü ¹½¼ ¼º½ ½¼ ± Ì Ò Ø Ý Üµ Ò Ü¹ Ü ººº Ò ÔÐÓØ Ý Üº ¹ º ܺ ¾º ܺ ¾¹ º Ü ½º¾ Ý ¼ Ü ¼ ÔÐÓØ Ü Ý ¼ Ü Ýµ Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
36 Bibliography 1 STEVEN C. CHAPRA, RAYMOND P. CANALE (2006): Applied Numerical Methods for Engineers, 5ed, ISBN , McGraw-Hill 2 SINGIRESU S. RAO (2002): Applied Numerical Methods for Engineers and Scientists, ISBN X, Prentice Hall 3 DAVID KINCAID AND WARD CHENEY (1991): Numerical Analysis: Mathematics of Scientific Computing, ISBN , Brooks/Cole Publishing Co. 4 STEVEN C. CHAPRA (2005): Applied Numerical Methods with MATLAB for Engineers and Scientists, ISBN , McGraw-Hill 5 WILLIAM J. PALM III (2005): Introduction to Matlab 7 for Engineers, ISBN , McGraw-Hill 6 JOHN D. ANDERSON, JR. (1995): Computational Fluid Dynamics The Basics with Applications, ISBN , McGraw-Hill Abu Hasan Abdullah (FME) SME 3023 Applied Numerical Methods Sept / 31
SKMM 3023 Applied Numerical Methods
SKMM 3023 Applied Numerical Methods Solution of Nonlinear Equations ibn Abdullah Faculty of Mechanical Engineering Òº ÙÐÐ ÚºÒÙÐÐ ¾¼½ SKMM 3023 Applied Numerical Methods Solution of Nonlinear Equations
More informationSME 3023 Applied Numerical Methods
UNIVERSITI TEKNOLOGI MALAYSIA SME 3023 Applied Numerical Methods Ordinary Differential Equations Abu Hasan Abdullah Faculty of Mechanical Engineering Sept 2012 Abu Hasan Abdullah (FME) SME 3023 Applied
More informationSKMM 3023 Applied Numerical Methods
UNIVERSITI TEKNOLOGI MALAYSIA SKMM 3023 Applied Numerical Methods Ordinary Differential Equations ibn Abdullah Faculty of Mechanical Engineering Òº ÙÐÐ ÚºÒÙÐÐ ¾¼½ SKMM 3023 Applied Numerical Methods Ordinary
More informationSKMM 3023 Applied Numerical Methods
UNIVERSITI TEKNOLOGI MALAYSIA SKMM 3023 Applied Numerical Methods Numerical Differentiation ibn Abdullah Faculty of Mechanical Engineering Òº ÙÐÐ ÚºÒÙÐÐ ¾¼½ SKMM 3023 Applied Numerical Methods Numerical
More informationComputational Methods for Engineers Programming in Engineering Problem Solving
Computational Methods for Engineers Programming in Engineering Problem Solving Abu Hasan Abdullah January 6, 2009 Abu Hasan Abdullah 2009 An Engineering Problem Problem Statement: The length of a belt
More informationSolution of Nonlinear Equations: Graphical and Incremental Sea
Outlines Solution of Nonlinear Equations: Graphical and s September 2, 2004 Outlines Part I: Solution of Nonlinear Equations Solution of Nonlinear Equations Introduction General Form of the Problem Types
More informationMMJ 1113 Computational Methods for Engineers
Faculty of Mechanical Engineering Engineering Computing Panel MMJ 1113 Computational Methods for Engineers Engineering Problem Solving Abu Hasan Abdullah Feb 2013 Ùº Òº ÙÐÐ ÚºÒÙÐÐ MMJ 1113 Computational
More informationF(jω) = a(jω p 1 )(jω p 2 ) Û Ö p i = b± b 2 4ac. ω c = Y X (jω) = 1. 6R 2 C 2 (jω) 2 +7RCjω+1. 1 (6jωRC+1)(jωRC+1) RC, 1. RC = p 1, p
ÓÖ Ò ÊÄ Ò Ò Û Ò Ò Ö Ý ½¾ Ù Ö ÓÖ ÖÓÑ Ö ÓÒ Ò ÄÈ ÐØ Ö ½¾ ½¾ ½» ½½ ÓÖ Ò ÊÄ Ò Ò Û Ò Ò Ö Ý ¾ Á b 2 < 4ac Û ÒÒÓØ ÓÖ Þ Û Ö Ð Ó ÒØ Ó Û Ð Ú ÕÙ Ö º ËÓÑ Ñ ÐÐ ÕÙ Ö Ö ÓÒ Ò º Ù Ö ÓÖ ½¾ ÓÖ Ù Ö ÕÙ Ö ÓÖ Ò ØÖ Ò Ö ÙÒØ ÓÒ
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 informationLund Institute of Technology Centre for Mathematical Sciences Mathematical Statistics
Lund Institute of Technology Centre for Mathematical Sciences Mathematical Statistics STATISTICAL METHODS FOR SAFETY ANALYSIS FMS065 ÓÑÔÙØ Ö Ü Ö Ì ÓÓØ ØÖ Ô Ð ÓÖ Ø Ñ Ò Ý Ò Ò ÐÝ In this exercise we will
More informationÇÙÐ Ò ½º ÅÙÐ ÔÐ ÔÓÐÝÐÓ Ö Ñ Ò Ú Ö Ð Ú Ö Ð ¾º Ä Ò Ö Ö Ù Ð Ý Ó ËÝÑ ÒÞ ÔÓÐÝÒÓÑ Ð º Ì ÛÓ¹ÐÓÓÔ ÙÒÖ Ö Ô Û Ö Ö ÖÝ Ñ ¹ ÝÓÒ ÑÙÐ ÔÐ ÔÓÐÝÐÓ Ö Ñ
ÅÙÐ ÔÐ ÔÓÐÝÐÓ Ö Ñ Ò ÝÒÑ Ò Ò Ö Ð Ö Ò Ó Ò Ö ÀÍ ÖÐ Òµ Ó Ò ÛÓÖ Û Ö Ò ÖÓÛÒ Ö Ú ½ ¼¾º ¾½ Û Åº Ä Ö Ö Ú ½ ¼¾º ¼¼ Û Äº Ñ Ò Ëº Ï ÒÞ ÖÐ Å ÒÞ ½ º¼ º¾¼½ ÇÙÐ Ò ½º ÅÙÐ ÔÐ ÔÓÐÝÐÓ Ö Ñ Ò Ú Ö Ð Ú Ö Ð ¾º Ä Ò Ö Ö Ù Ð Ý Ó ËÝÑ
More informationSeminar to the lecture Computer-based Engineering Mathematics
Seminar to the lecture Computer-based Engineering Mathematics N T S Prof. Dr.-Ing. A. Czylwik Mohammad Abdelqader, M.Sc. Room: BA 249, Tel: +49-203-379-3474 E-Mail: abdelqader@nts.uni-duisburg-essen.de
More information2 Hallén s integral equation for the thin wire dipole antenna
Ú Ð Ð ÓÒÐ Ò Ø ØØÔ»» Ѻ Ö Ùº º Ö ÁÒغ º ÁÒ Ù ØÖ Ð Å Ø Ñ Ø ÎÓк ÆÓº ¾ ¾¼½½µ ½ ¹½ ¾ ÆÙÑ Ö Ð Ñ Ø Ó ÓÖ Ò ÐÝ Ó Ö Ø ÓÒ ÖÓÑ Ø Ò Û Ö ÔÓÐ ÒØ ÒÒ Ëº À Ø ÑÞ ¹Î ÖÑ ÞÝ Ö Åº Æ Ö¹ÅÓ Êº Ë Þ ¹Ë Ò µ Ô ÖØÑ ÒØ Ó Ð ØÖ Ð Ò Ò
More informationLecture 16: Modern Classification (I) - Separating Hyperplanes
Lecture 16: Modern Classification (I) - Separating Hyperplanes Outline 1 2 Separating Hyperplane Binary SVM for Separable Case Bayes Rule for Binary Problems Consider the simplest case: two classes are
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 informationSolution of Nonlinear Equations: Graphical and Incremental Sea
Outlines Solution of Nonlinear Equations: Graphical and Incremental Search Methods September 2, 2004 Outlines Part I: Review of Previous Lecture Part II: Sample Problems Solved with Numerical Methods Part
More informationScientific Computing. Roots of Equations
ECE257 Numerical Methods and Scientific Computing Roots of Equations Today s s class: Roots of Equations Polynomials Polynomials A polynomial is of the form: ( x) = a 0 + a 1 x + a 2 x 2 +L+ a n x n f
More informationSOLVING EQUATIONS OF ONE VARIABLE
1 SOLVING EQUATIONS OF ONE VARIABLE ELM1222 Numerical Analysis Some of the contents are adopted from Laurene V. Fausett, Applied Numerical Analysis using MATLAB. Prentice Hall Inc., 1999 2 Today s lecture
More informationNUMERICAL AND STATISTICAL COMPUTING (MCA-202-CR)
NUMERICAL AND STATISTICAL COMPUTING (MCA-202-CR) Autumn Session UNIT 1 Numerical analysis is the study of algorithms that uses, creates and implements algorithms for obtaining numerical solutions to problems
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 informationBEKG 2452 NUMERICAL METHODS Solution of Nonlinear Equations
BEKG 2452 NUMERICAL METHODS Solution of Nonlinear Equations Ser Lee Loh a, Wei Sen Loi a a Fakulti Kejuruteraan Elektrik Universiti Teknikal Malaysia Melaka Lesson Outcome Upon completion of this lesson,
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 informationcha1873x_p02.qxd 3/21/05 1:01 PM Page 104 PART TWO
cha1873x_p02.qxd 3/21/05 1:01 PM Page 104 PART TWO ROOTS OF EQUATIONS PT2.1 MOTIVATION Years ago, you learned to use the quadratic formula x = b ± b 2 4ac 2a to solve f(x) = ax 2 + bx + c = 0 (PT2.1) (PT2.2)
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 informationCE 205 Numerical Methods. Some of the analysis methods you have used so far.. Algebra Calculus Differential Equations etc.
CE 205 Numerical Methods Dr. Charisma Choudhury Lecture 1 March 30, 2009 Objective Some of the analysis methods you have used so far.. Algebra Calculus Differential Equations etc. Often not possible to
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 informationZeroes of Transcendental and Polynomial Equations. Bisection method, Regula-falsi method and Newton-Raphson method
Zeroes of Transcendental and Polynomial Equations Bisection method, Regula-falsi method and Newton-Raphson method PRELIMINARIES Solution of equation f (x) = 0 A number (real or complex) is a root of the
More informationChapter 1. Root Finding Methods. 1.1 Bisection method
Chapter 1 Root Finding Methods We begin by considering numerical solutions to the problem f(x) = 0 (1.1) Although the problem above is simple to state it is not always easy to solve analytically. This
More informationChapter 4. Solution of Non-linear Equation. Module No. 1. Newton s Method to Solve Transcendental Equation
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Mathematics National Institute of Technology Durgapur Durgapur-713209 email: anita.buie@gmail.com 1 . Chapter 4 Solution of Non-linear
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 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 informationx 0, x 1,...,x n f(x) p n (x) = f[x 0, x 1,..., x n, x]w n (x),
ÛÜØ Þ ÜÒ Ô ÚÜ Ô Ü Ñ Ü Ô Ð Ñ Ü ÜØ º½ ÞÜ Ò f Ø ÚÜ ÚÛÔ Ø Ü Ö ºÞ ÜÒ Ô ÚÜ Ô Ð Ü Ð Þ Õ Ô ÞØÔ ÛÜØ Ü ÚÛÔ Ø Ü Ö L(f) = f(x)dx ÚÜ Ô Ü ÜØ Þ Ü Ô, b] Ö Û Þ Ü Ô Ñ ÒÖØ k Ü f Ñ Df(x) = f (x) ÐÖ D Ü Ü ÜØ Þ Ü Ô Ñ Ü ÜØ Ñ
More informationAn Example file... log.txt
# ' ' Start of fie & %$ " 1 - : 5? ;., B - ( * * B - ( * * F I / 0. )- +, * ( ) 8 8 7 /. 6 )- +, 5 5 3 2( 7 7 +, 6 6 9( 3 5( ) 7-0 +, => - +< ( ) )- +, 7 / +, 5 9 (. 6 )- 0 * D>. C )- +, (A :, C 0 )- +,
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 information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 4
2.29 Spring 2015 Lecture 4 Review Lecture 3 Truncation Errors, Taylor Series and Error Analysis Taylor series: 2 3 n n i1 i i i i i n f( ) f( ) f '( ) f ''( ) f '''( )... f ( ) R 2! 3! n! n1 ( n1) Rn f
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 informationWhat are Numerical Methods? (1/3)
What are Numerical Methods? (1/3) Numerical methods are techniques by which mathematical problems are formulated so that they can be solved by arithmetic and logic operations Because computers excel at
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 informationAll that begins... peace be upon you
All that begins... peace be upon you Faculty of Mechanical Engineering Department of Thermo Fluids SKMM 2323 Mechanics of Fluids 2 «An excerpt (mostly) from White (2011)» ibn Abdullah May 2017 Outline
More informationCONVEX OPTIMIZATION OVER POSITIVE POLYNOMIALS AND FILTER DESIGN. Y. Genin, Y. Hachez, Yu. Nesterov, P. Van Dooren
CONVEX OPTIMIZATION OVER POSITIVE POLYNOMIALS AND FILTER DESIGN Y. Genin, Y. Hachez, Yu. Nesterov, P. Van Dooren CESAME, Université catholique de Louvain Bâtiment Euler, Avenue G. Lemaître 4-6 B-1348 Louvain-la-Neuve,
More informationA Language for Task Orchestration and its Semantic Properties
DEPARTMENT OF COMPUTER SCIENCES A Language for Task Orchestration and its Semantic Properties David Kitchin, William Cook and Jayadev Misra Department of Computer Science University of Texas at Austin
More informationGreen s function, wavefunction and Wigner function of the MIC-Kepler problem
Green s function, wavefunction and Wigner function of the MIC-Kepler problem Tokyo University of Science Graduate School of Science, Department of Mathematics, The Akira Yoshioka Laboratory Tomoyo Kanazawa
More informationNON-LINEAR ALGEBRAIC EQUATIONS Lec. 5.1: Nonlinear Equation in Single Variable
NON-LINEAR ALGEBRAIC EQUATIONS Lec. 5.1: Nonlinear Equation in Single Variable Dr. Niket Kaisare Department of Chemical Engineering IIT Madras NPTEL Course: MATLAB Programming for Numerical Computations
More informationNumerical Methods. Root Finding
Numerical Methods Solving Non Linear 1-Dimensional Equations Root Finding Given a real valued function f of one variable (say ), the idea is to find an such that: f() 0 1 Root Finding Eamples Find real
More informationu x + u y = x u . u(x, 0) = e x2 The characteristics satisfy dx dt = 1, dy dt = 1
Õ 83-25 Þ ÛÐ Þ Ð ÚÔÜØ Þ ÝÒ Þ Ô ÜÞØ ¹ 3 Ñ Ð ÜÞ u x + u y = x u u(x, 0) = e x2 ÝÒ Þ Ü ÞØ º½ dt =, dt = x = t + c, y = t + c 2 We can choose c to be zero without loss of generality Note that each characteristic
More informationµ(, y) Computing the Möbius fun tion µ(x, x) = 1 The Möbius fun tion is de ned b y and X µ(x, t) = 0 x < y if x6t6y 3
ÈÖÑÙØØÓÒ ÔØØÖÒ Ò Ø ÅÙ ÙÒØÓÒ ÙÖ ØÒ ÎØ ÂÐÒ Ú ÂÐÒÓÚ Ò ÐÜ ËØÒÖÑ ÓÒ ÒÖ Ì ØÛÓµ 2314 ½¾ ½ ¾ ¾½ ¾ ½ ½¾ ¾½ ½¾ ¾½ ½ Ì ÔÓ Ø Ó ÔÖÑÙØØÓÒ ÛºÖºØº ÔØØÖÒ ÓÒØÒÑÒØ ½ 2314 ½¾ ½ ¾ ¾½ ¾ ½ ½¾ ¾½ ½¾ ¾½ Ì ÒØÖÚÐ [12,2314] ½ ¾ ÓÑÔÙØÒ
More informationSolution of Nonlinear Equations
Solution of Nonlinear Equations In many engineering applications, there are cases when one needs to solve nonlinear algebraic or trigonometric equations or set of equations. These are also common in Civil
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2018 3 Lecture 3 3.1 General remarks March 4, 2018 This
More informationNumerical Study of Some Iterative Methods for Solving Nonlinear Equations
International Journal of Engineering Science Invention ISSN (Online): 2319 6734, ISSN (Print): 2319 6726 Volume 5 Issue 2 February 2016 PP.0110 Numerical Study of Some Iterative Methods for Solving Nonlinear
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 informationby Martin Mendez, UASLP Copyright 2006 The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
Chapter 5 by Martin Mendez, 1 Roots of Equations Part Why? b m b a + b + c = 0 = aa 4ac But a 5 4 3 + b + c + d + e + f = 0 sin + = 0 =? =? by Martin Mendez, Nonlinear Equation Solvers Bracketing Graphical
More informationINTRODUCTION TO NUMERICAL ANALYSIS
INTRODUCTION TO NUMERICAL ANALYSIS Cho, Hyoung Kyu Department of Nuclear Engineering Seoul National University 3. SOLVING NONLINEAR EQUATIONS 3.1 Background 3.2 Estimation of errors in numerical solutions
More informationA Derivative Free Hybrid Equation Solver by Alloying of the Conventional Methods
DOI: 1.15415/mjis.213.129 A Derivative Free Hybrid Equation Solver by Alloying of the Conventional Methods Amit Kumar Maheshwari Advanced Materials and Processes Research Institute (CSIR, Bhopal, India
More informationProving observational equivalence with ProVerif
Proving observational equivalence with ProVerif Bruno Blanchet INRIA Paris-Rocquencourt Bruno.Blanchet@inria.fr based on joint work with Martín Abadi and Cédric Fournet and with Vincent Cheval June 2015
More informationFinding the Roots of f(x) = 0. Gerald W. Recktenwald Department of Mechanical Engineering Portland State University
Finding the Roots of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab:
More informationFinding the Roots of f(x) = 0
Finding the Roots of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab:
More informationNumerical solutions of nonlinear systems of equations
Numerical solutions of nonlinear systems of equations Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan E-mail: min@math.ntnu.edu.tw August 28, 2011 Outline 1 Fixed points
More informationRadu Alexandru GHERGHESCU, Dorin POENARU and Walter GREINER
È Ö Ò Ò Ù Ò Ò Ò ÖÝ ÒÙÐ Ö Ý Ø Ñ Radu Alexandru GHERGHESCU, Dorin POENARU and Walter GREINER Radu.Gherghescu@nipne.ro IFIN-HH, Bucharest-Magurele, Romania and Frankfurt Institute for Advanced Studies, J
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 informationNumerical Solution of f(x) = 0
Numerical Solution of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@pdx.edu ME 350: Finding roots of f(x) = 0 Overview Topics covered in these slides
More informationProblem 1 (From the reservoir to the grid)
ÈÖÓ º ĺ ÙÞÞ ÐÐ ÈÖÓ º ʺ ³ Ò Ö ½ ½¹¼ ¼¹¼¼ ËÝ Ø Ñ ÅÓ Ð Ò ÀË ¾¼½ µ Ü Ö ÌÓÔ ÀÝ ÖÓ Ð ØÖ ÔÓÛ Ö ÔÐ ÒØ À Èȵ ¹ È ÖØ ÁÁ Ð ÖÒ Ø Þº ÇØÓ Ö ½ ¾¼½ Problem (From the reservoir to the grid) The causality diagram of the
More informationSolutions of Equations in One Variable. Newton s Method
Solutions of Equations in One Variable Newton s Method Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University c 2011 Brooks/Cole,
More informationNumerical Analysis: Solving Nonlinear Equations
Numerical Analysis: Solving Nonlinear Equations Mirko Navara http://cmp.felk.cvut.cz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office 104a
More informationComparative Analysis of Convergence of Various Numerical Methods
Journal of Computer and Mathematical Sciences, Vol.6(6),290-297, June 2015 (An International Research Journal), www.compmath-journal.org ISSN 0976-5727 (Print) ISSN 2319-8133 (Online) Comparative Analysis
More informationA three point formula for finding roots of equations by the method of least squares
Journal of Applied Mathematics and Bioinformatics, vol.2, no. 3, 2012, 213-233 ISSN: 1792-6602(print), 1792-6939(online) Scienpress Ltd, 2012 A three point formula for finding roots of equations by the
More informationJuan Juan Salon. EH National Bank. Sandwich Shop Nail Design. OSKA Beverly. Chase Bank. Marina Rinaldi. Orogold. Mariposa.
( ) X é X é Q Ó / 8 ( ) Q / ( ) ( ) : ( ) : 44-3-8999 433 4 z 78-19 941, #115 Z 385-194 77-51 76-51 74-7777, 75-5 47-55 74-8141 74-5115 78-3344 73-3 14 81-4 86-784 78-33 551-888 j 48-4 61-35 z/ zz / 138
More informationNotes for Numerical Analysis Math 5465 by S. Adjerid Virginia Polytechnic Institute and State University. (A Rough Draft)
Notes for Numerical Analysis Math 5465 by S. Adjerid Virginia Polytechnic Institute and State University (A Rough Draft) 1 2 Contents 1 Error Analysis 5 2 Nonlinear Algebraic Equations 7 2.1 Convergence
More informationAppendix 8 Numerical Methods for Solving Nonlinear Equations 1
Appendix 8 Numerical Methods for Solving Nonlinear Equations 1 An equation is said to be nonlinear when it involves terms of degree higher than 1 in the unknown quantity. These terms may be polynomial
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More informationA New Modification of Newton s Method
A New Modification of Newton s Method Vejdi I. Hasanov, Ivan G. Ivanov, Gurhan Nedjibov Laboratory of Mathematical Modelling, Shoumen University, Shoumen 971, Bulgaria e-mail: v.hasanov@@fmi.shu-bg.net
More informationChapter 5 Roots of Equations: Bracketing Models. Gab-Byung Chae
Chapter 5 Roots of Equations: Bracketing Models Gab-Byung Chae 2008 4 17 2 Chapter Objectives Studying Bracketing Methods Understanding what roots problems are and where they occur in engineering and science.
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 informationAn Improved Hybrid Algorithm to Bisection Method and Newton-Raphson Method
Applied Mathematical Sciences, Vol. 11, 2017, no. 56, 2789-2797 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.710302 An Improved Hybrid Algorithm to Bisection Method and Newton-Raphson
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 informationNonlinearity Root-finding Bisection Fixed Point Iteration Newton s Method Secant Method Conclusion. Nonlinear Systems
Nonlinear Systems CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Nonlinear Systems 1 / 24 Part III: Nonlinear Problems Not all numerical problems
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 informationComputational Methods. Solving Equations
Computational Methods Solving Equations Manfred Huber 2010 1 Solving Equations Solving scalar equations is an elemental task that arises in a wide range of applications Corresponds to finding parameters
More informationApplications of Discrete Mathematics to the Analysis of Algorithms
Applications of Discrete Mathematics to the Analysis of Algorithms Conrado Martínez Univ. Politècnica de Catalunya, Spain May 2007 Goal Given some algorithm taking inputs from some set Á, we would like
More informationArbeitstagung: Gruppen und Topologische Gruppen Vienna July 6 July 7, Abstracts
Arbeitstagung: Gruppen und Topologische Gruppen Vienna July 6 July 7, 202 Abstracts ÁÒÚ Ö Ð Ñ Ø Ó Ø¹Ú ÐÙ ÙÒØ ÓÒ ÁÞØÓ Ò ÞØÓ º Ò ÙÒ ¹Ñ º ÙÐØÝ Ó Æ ØÙÖ Ð Ë Ò Ò Å Ø Ñ Ø ÍÒ Ú Ö ØÝ Ó Å Ö ÓÖ ÃÓÖÓ ½ ¼ Å Ö ÓÖ ¾¼¼¼
More informationProblem 1 (From the reservoir to the grid)
ÈÖÓ º ĺ ÙÞÞ ÐÐ ÈÖÓ º ʺ ³ Ò Ö ½ ½¹¼ ¹¼¼ ËÝ Ø Ñ ÅÓ Ð Ò ÀË ¾¼½ µ Ü Ö ËÓÐÙØ ÓÒ ÌÓÔ ÀÝ ÖÓ Ð ØÖ ÔÓÛ Ö ÔÐ ÒØ À Èȵ ¹ È ÖØ ÁÁ Ð ÖÒ Ø Þº ÇØÓ Ö ¾ ¾¼½ Problem 1 (From the reservoir to the grid) The causality diagram
More informationNumerical Methods Lecture 3
Numerical Methods Lecture 3 Nonlinear Equations by Pavel Ludvík Introduction Definition (Root or zero of a function) A root (or a zero) of a function f is a solution of an equation f (x) = 0. We learn
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 informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
3 Numerical Solution of Nonlinear Equations and Systems 3.1 Fixed point iteration Reamrk 3.1 Problem Given a function F : lr n lr n, compute x lr n such that ( ) F(x ) = 0. In this chapter, we consider
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 informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationCourse in. Geometric nonlinearity. Nonlinear FEM. Computational Mechanics, AAU, Esbjerg
Course in Nonlinear FEM Geometric nonlinearity Nonlinear FEM Outline Lecture 1 Introduction Lecture 2 Geometric nonlinearity Lecture 3 Material nonlinearity Lecture 4 Material nonlinearity it continued
More informationINRIA Sophia Antipolis France. TEITP p.1
ÌÖÙ Ø ÜØ Ò ÓÒ Ò ÓÕ Ä ÙÖ ÒØ Ì ÖÝ INRIA Sophia Antipolis France TEITP p.1 ÅÓØ Ú Ø ÓÒ Ï Ý ØÖÙ Ø Ó ÑÔÓÖØ ÒØ Å ÒÐÝ ÈÖÓÚ Ò ÌÖÙØ Ø Ò ÑÔÐ Ã Ô ÈÖÓÚ Ò ÌÖÙ Ø È Ó ÖÖÝ Ò ÈÖÓÓ µ Ò Ö ØÝ ÓÑ Ò ËÔ ÔÔÐ Ø ÓÒ TEITP p.2 ÇÙØÐ
More informationMATH 3795 Lecture 12. Numerical Solution of Nonlinear Equations.
MATH 3795 Lecture 12. Numerical Solution of Nonlinear Equations. Dmitriy Leykekhman Fall 2008 Goals Learn about different methods for the solution of f(x) = 0, their advantages and disadvantages. Convergence
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 informationarxiv:hep-ph/ v1 10 May 2001
New data and the hard pomeron A Donnachie Department of Physics, Manchester University P V Landshoff DAMTP, Cambridge University DAMTP-200-38 M/C-TH-0/03 arxiv:hep-ph/005088v 0 May 200 Abstract New structure-function
More informationLA PRISE DE CALAIS. çoys, çoys, har - dis. çoys, dis. tons, mantz, tons, Gas. c est. à ce. C est à ce. coup, c est à ce
> ƒ? @ Z [ \ _ ' µ `. l 1 2 3 z Æ Ñ 6 = Ð l sl (~131 1606) rn % & +, l r s s, r 7 nr ss r r s s s, r s, r! " # $ s s ( ) r * s, / 0 s, r 4 r r 9;: < 10 r mnz, rz, r ns, 1 s ; j;k ns, q r s { } ~ l r mnz,
More informationPH Nuclear Physics Laboratory Gamma spectroscopy (NP3)
Physics Department Royal Holloway University of London PH2510 - Nuclear Physics Laboratory Gamma spectroscopy (NP3) 1 Objectives The aim of this experiment is to demonstrate how γ-ray energy spectra may
More informationHence a root lies between 1 and 2. Since f a is negative and f(x 0 ) is positive The root lies between a and x 0 i.e. 1 and 1.
The Bisection method or BOLZANO s method or Interval halving method: Find the positive root of x 3 x = 1 correct to four decimal places by bisection method Let f x = x 3 x 1 Here f 0 = 1 = ve, f 1 = ve,
More information! " # $! % & '! , ) ( + - (. ) ( ) * + / 0 1 2 3 0 / 4 5 / 6 0 ; 8 7 < = 7 > 8 7 8 9 : Œ Š ž P P h ˆ Š ˆ Œ ˆ Š ˆ Ž Ž Ý Ü Ý Ü Ý Ž Ý ê ç è ± ¹ ¼ ¹ ä ± ¹ w ç ¹ è ¼ è Œ ¹ ± ¹ è ¹ è ä ç w ¹ ã ¼ ¹ ä ¹ ¼ ¹ ±
More informationPlanning for Reactive Behaviors in Hide and Seek
University of Pennsylvania ScholarlyCommons Center for Human Modeling and Simulation Department of Computer & Information Science May 1995 Planning for Reactive Behaviors in Hide and Seek Michael B. Moore
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 informationApplied Mathematics Letters. Combined bracketing methods for solving nonlinear equations
Applied Mathematics Letters 5 (01) 1755 1760 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Combined bracketing methods for
More informationScientific Computing. Roots of Equations
ECE257 Numerical Methods and Scientific Computing Roots of Equations Today s s class: Roots of Equations Bracketing Methods Roots of Equations Given a function f(x), the roots are those values of x that
More informationThe iteration formula for to find the root of the equation
SRI RAMAKRISHNA INSTITUTE OF TECHNOLOGY, COIMBATORE- 10 DEPARTMENT OF SCIENCE AND HUMANITIES SUBJECT: NUMERICAL METHODS & LINEAR PROGRAMMING UNIT II SOLUTIONS OF EQUATION 1. If is continuous in then under
More informationCNED - Numerical Calculus. Differential Equations
Coordinating unit: Teaching unit: Academic year: Degree: ECTS credits: 2017 295 - EEBE - Barcelona East School of Engineering 749 - MAT - Department of Mathematics BACHELOR'S DEGREE IN ELECTRICAL ENGINEERING
More information