2.29 Numerical Fluid Mechanics Fall 2011 Lecture 3
|
|
- Christina Floyd
- 6 years ago
- Views:
Transcription
1 Numerical Fluid Mechaics Fall 2011 Lecture 3 REVIEW Lectures 1-2 Approximatio ad roud-off errors ˆx a xˆ Absolute ad relative errors: E a xˆ a ˆx, a ˆx a xˆ Iterative schemes ad stop criterio: ˆx 1 a ˆx 1 For digits: 10 s 2 Number represetatios Floatig-Poit represetatio: x = m b e Quatizig errors ad their cosequeces Arithmetic operatios (e.g. radd.m) Errors of arithmetic/umerical operatios Large umber of additio/subtractios, large & small umbers, ier products, etc Recursio algorithms (Hero, Horer s scheme, Bessel fuctios, etc): Order of computatios matter Roud-off error growth ad (i)-stability s Numerical Fluid Mechaics PFJL Lecture 3, 1 1
2 Numerical Fluid Mechaics Fall 2011 Lecture 3: Outlie Trucatio Errors, Taylor Series ad Error Aalysis Taylor series: f ( x i 1 ) f ( x i ) x f '( x i ) x 1 ( 1 R f ) ( ) 1! x 2 x 3 x f ''( x i ) f '''(x i )... 2! 3!! f (x i ) R Use of Taylor Series to derive fiite differece schemes (first-order Euler scheme ad forward, backward ad cetered differeces) Geeral error propagatio formulas ad error estimatio, with examples Cosider y f ( x 1, x 2, x 3,..., x ). I f i' s are magitudes of errors o x i 's, what is the error o y? The Differetial Formula: f( x 1,..., x ) y i i 1 x i f 2 The Stadard Error (statistical formula): 2 E( s y) i i 1 x i Error cacellatio (e.g. subtractio of errors of the same sig) Coditio umber: x f '( x ) K p f ( x ) Well-coditioed problems vs. well-coditioed algorithms Referece: Chapra ad Caale, Numerical stability Chaps 3, 4 ad 5 Numerical Fluid Mechaics PFJL Lecture 3, 2 2
3 Numerical Fluid Mechaics: Lectures 3 ad 4 Outlie Coditio Numbers Roots of oliear equatios Chaps 3, 4 ad 5 Bracketig Methods Example: Hero s formula Bisectio False Positio Ope Methods Ope-poit Iteratio (Geeral method or Picard Iteratio) Examples Covergece Criteria Order of Covergece Newto-Raphso Covergece speed ad examples Secat Method Examples Covergece ad efficiecy Extesio of Newto-Raphso to systems of oliear equatios Roots of Polyomial (all real/complex roots) Ope methods (applicatios of the above for complex umbers) Special Methods (e.g. Muller s ad Bairstow s methods) Referece: Chapra ad Caale, Numerical Fluid Mechaics PFJL Lecture 3, 3 3
4 Trucatio Errors, Taylor Series ad Error Aalysis Taylor Series: Provides a mea to predict a fuctio at oe poit i terms of its values ad derivatives at aother poit (i the form of a polyomial) Hece, ay smooth fuctios ca be approximated by a polyomial Taylor Series (Mea for itegrals theorems): 2 3 x x x i 1 i x f '(x i ) f ''(x i ) f '''(x i )... f (x i ) R f (x ) f (x ) 1 x ( 1) R f ( ) 1! = costat + lie + parabola + etc 2! 3!! Numerical Fluid Mechaics PFJL Lecture 3, 4 4
5 Derivatio of Geeral Differetial Error Propagatio Formula Uivariate Case Recall: Hece, f i 1 ( x ) f ( x ) x f R 1 x 1! i f ( 1) '( x i ) x 2! 1 ( ) O ( x ) y f f( x i 1) f( x i ) For xf'( x ) x 1, f x i x 2! f 2 '( x i ) 2 f ''( x i ) f ''( x i ) x 3! 3 x 3! 2 O ( x ) x f '( x ) 3 f '''( x i )... f '''( x i )... i x! x! f ( x ) R i f ( x i ) R Thus, for a error o x equal to Δx such that x 1, we have a error o y equal to : y y f x f '(x i ) x f '( x i ) f '( x i ) Multivariate case y f (x, x, x,..., x ) f( x 1,..., x ) For x i 1, y i 1 x i i Derivatio doe i class o the board Numerical Fluid Mechaics PFJL Lecture 3, 5 5
6 Geeral Error Propagatio Formula (The Differetial Formula) For this large plotted x, secod derivative ot egligible Absolute Errors? Fuctio of oe variable x = x -x y~f (x) x Geeral Error Propagatio Formula x x Numerical Fluid Mechaics PFJL Lecture 3, 6 6
7 Error Propagatio Example with Differetial Approach: Multiplicatios Multiplicatio Error Propagatio Formula => => => ε r y ε r i => Relative Errors Add for Multiplicatio Aother example, more geeral case: ε r y ε r i Numerical Fluid Mechaics PFJL Lecture 3, 7 7
8 Statistical Error Propagatio Estimate: Error Expectatios ad the Stadard Error Example: Additio a. Choppig Error Expectatio (o average, half the iterval is chopped) b. Roudig Numerical Fluid Mechaics PFJL Lecture 3, 8 8
9 Statistical Error Propagatio Estimate: Error Expectatios ad the Stadard Error y y( x 1, x 2, x 3,..., x ) Stadard Error (assumig Gaussia errors o x i, compute error stadard deviatio o y) 2 Example: Additio Stadard Error is a more accurate measure of expected errors: It is likely that the error will be of that magitude. However, the geeral error propagatio formula is a upper boud Numerical Fluid Mechaics PFJL Lecture 3, 9 9
10 Error Propagatio: Error Cacellatio (occurs whe errors of the same sig are substracted) Example: Cosider this fuctio of oe variable Defie: For 2 terms (z 1 ad z 2 ): Ge. error (max): Stad. Error : Direct error (sigle y term): Error cacellatio 2 2 error if 4 digits Numerical Fluid Mechaics PFJL Lecture 3, 10 10
11 The Coditio Number The coditio of a mathematical problem relates to its sesitivity to chages i its iput values A computatio is umerically ustable if the ucertaity of the iput values are magified by the umerical method Cosiderig x ad f(x), the coditio umber is the ratio of the relative error i f(x) to that i x. Usig first-order Taylor series f ( x ) f ( x ) f '( x)( x x ) Relative error i f(x): Relative error i x: (x x ) x f ( x ) f ( x ) f '( x)( x x ) f( x) f ( x ) Coditio Nb = Ratio of relative errors: K p x f '( x) f ( x ) Numerical Fluid Mechaics PFJL Lecture 3, 11 11
12 Error Propagatio Coditio Number: aother derivatio ad a example y = f(x) x x y y x x (1 ) y y (1 ) f ( x ) f( x)(1 ) Problem Coditio Number Problem ill-coditioed Error cacellatio example With first-order Taylor series Well-coditioed problem Numerical Fluid Mechaics PFJL Lecture 3, 12 12
13 Error Propagatio Problem vs. Algorithm Coditio Numbers a. Problem Coditio Number b. Algorithm (4 Sigificat Digits) Coditio Number 2 K A is the algorithm coditio umber, which may be much larger tha K P due to digitized umber represetatio. Solutio: Higher precisio Rewrite algorithm Well-coditioed Algorithm (base 10) Algorithm Coditio Number Numerical Fluid Mechaics PFJL Lecture 3, 13 13
14 Roots of Noliear Equatios Fidig roots of equatios (x f (x) = 0) ubiquitous i egieerig problems, icludig umerical fluid dyamics applicatios I geeral: fidig roots of fuctios (implicit i ukow parameters/roots) Fluid Applicatio Examples Implicit itegratio schemes Statioary poits i fluid equatios, stability characterizatios, etc Fluid dyamical system properties Fluid egieerig system desig Most realistic root fidig problems require umerical methods Two types of root problems/methods: 1.Determie the real roots of algebraic/trascedetal equatios (these methods usually eed some sort of prior kowledge o a root locatio) 2.Determie all real ad complex roots of polyomials (all roots at oce) Numerical Fluid Mechaics PFJL Lecture 3, 14 14
15 Roots of Noliear Equatios Bracketig Methods Bracket root, systematically reduce width of bracket, track error for covergece Methods: Bisectio ad False-Positio (Regula Falsi) Ope Methods Systematic Trial ad Error schemes, do t require a bracket (ca start from a sigle value) Computatioally efficiet, but do t always coverge Methods: Ope-poit Iteratio (Geeral method or Picard Iteratio), Newto-Raphso, Secat Method Extesio of Newto-Raphso to systems of oliear equatios Roots of Polyomials (all real/complex roots) Ope methods (applicatios of the above for complex umbers) Special Methods (e.g. Muller s ad Bairstow s methods) Numerical Fluid Mechaics PFJL Lecture 3, 15 15
16 Bracketig Methods: Graphical Approach A fuctio typically chages sig at the viciity of the root bracket the root f(x) f(x) f(x) f(x) x x x x x t x a x t x a x t x a xt xa (i) (ii) (i) (ii) f(x) f(x) x x x t x a x t x a (iii) (iv) These graphs show geeral ways i which roots may occur o a bouded iterval. I (i) ad (ii), if the fuctio is positive at both boudaries, there will be either a eve umber of roots or o roots; i the other two graphs, if the fuctio has opposite sigs at the boudary poits, there will be a odd umber of roots withi the iterval. Image by MIT OpeCourseWare. These graphs show some exceptios to the geeral cases illustrated o the left. I (i), a multiple root occurs where the fuctio is tagetial to the x axis. There are a eve umber of itersectios for the iterval. I (ii) is depicted a discotiuous fuctio with a eve umber of roots ad ed poits of opposite sigs. Fidig the roots i such cases requires special strategies. Numerical Fluid Mechaics PFJL Lecture 3, 16 16
17 Roots of Noliear Equatios: Bracketig Methods If Example Square root Hero s Priciple Guess root a=2; =6; hero.m g=2; % Number of Digits dig=5; sq(1)=g; for i=2: sq(i)= 0.5*radd(sq(i-1),a/sq(i-1),dig); ed ' i value ' [ [1:]' sq'] hold off plot([0 ],[sqrt(a) sqrt(a)],'b') hold o plot(sq,'r') plot(a./sq,'r-.') plot((sq-sqrt(a))/sqrt(a),'g') leged('sqrt','x','s/x','relative Err') grid o If i value Mea is better guess ( ) / 2 Iteratio Formula ( ) / Numerical Fluid Mechaics PFJL Lecture 3, 17 17
18 Roots of Noliear Equatios: Bracketig Methods Bisectio Icremetal search method defied as: Keep dividig i two the iterval withi which at least oe root lies Algorithm = +1 f(x) yes x o Numerical Fluid Mechaics PFJL Lecture 3, 18 18
19 Roots of Noliear Equatios: Bracketig Methods Bisectio Algorithm = +1 % Root fidig by bi-sectio f=ilie(' a*x -1','x','a'); a=2 bisect.m figure(1); clf; hold o x=[0 1.5]; eps=1e-3; err=max(abs(x(1)-x(2)),abs(f(x(1),a)-f(x(2),a))); while (err>eps & f(x(1),a)*f(x(2),a) <= 0) xo=x; x=[xo(1) 0.5*(xo(1)+xo(2))]; if ( f(x(1),a)*f(x(2),a) > 0 ) x=[0.5*(xo(1)+xo(2)) xo(2)] ed x err=max(abs(x(1)-x(2)),abs(f(x(1),a)-f(x(2),a))); b=plot(x,f(x,a),'.b'); set(b,'markersize',20); grid o; ed Stop criterio above based o x~0 ad f(x)~0 yes o Numerical Fluid Mechaics PFJL Lecture 3, 19 19
20 Roots of Noliear Equatios: Bracketig Methods Bisectio Stop criteria: whe relative estimated error or relative added value is egligible: a 1 xˆ r xˆ r s xˆ r Each successive iteratio halves the maximum error, a geeral formula thus relates the error ad the umber of iteratios: x 0 log 2 E ad, where x 0 is a iitial error ad E ad the desired error, Numerical Fluid Mechaics PFJL Lecture 3, 20 20
21 Aother Bracketig Method: The False Positio Method (Regula Falsi) Bisectio is a brute force scheme Somewhat iefficiet approach Does t accout for magitudes of f(x L ) ad f(x U ) f( xu ) x x x ( L x U ) r U f( xl ) f( x U ) False Positio f( x ) f( x ) f( x ) f( x ) x x ) 0 U x A graphical U x L depictio of the false positio method, Regula U Falsi. Shaded i gree are the r L r U triagles used to derive the method. U L r U ( r f( xl ) f( x ) (x x ) (x x ) f(x) x r f(x u ) x L x u x f( xu ) x x x ( L x U ) r U f( x ) f( x ) L U f(x L ) Image by MIT OpeCourseWare. Numerical Fluid Mechaics PFJL Lecture 3, 21 21
22 Commets o Bracketig Methods (Bisectio ad False Postio) Bisectio ad False Positio: Coverget methods, but oe eeds to miimize fuctio evaluatios Error of False Positio ca decrease much faster tha that of Bisectio, but if the fuctio is far from a straight lie, False Positio will be slow (oe of the ed poit remais fixed) Remedy: Modified False Positio e.g. if ed poit remais fixed for a few iteratios, reduce ed-poit fuctio value at ed poit, e.g. divide iterval by 2 Oe ca ot really geeralize for Root Fidig methods (results are fuctio deped) Always substitute estimate of root x r at the ed to check how close estimate is to f(x r )=0 Numerical Fluid Mechaics PFJL Lecture 3, 22 22
23 Roots of Noliear Equatios: Ope Methods Fixed Poit Iteratio (Geeral Method or Picard Iteratio) No-liear Equatio Goal: Covergig series Rewrite Problem Example Example: Cube root % f(x) = x^3 - a = 0 % g(x) = x + C*(x^3 - a) a=2; cube.m =10; g=1.0; C=-0.1; sq(1)=g; for i=2: sq(i)= sq(i-1) + C*(sq(i-1)^3 -a); ed hold off plot([0 ],[a^(1./3.) a^(1/3.)],'b') hold o plot(sq,'r') plot( (sq-a^(1./3.))/(a^(1./3.)),'g') grid o Iteratio Numerical Fluid Mechaics PFJL Lecture 3, 23 23
24 Roots of Noliear Equatios: Bracketig vs. Ope Methods i Graphics f(x) f(x) x i xi+1 x x l x u x f(x) x l x u x l x u x i xi+1 x x l x u x l x u Image by MIT OpeCourseWare. A graphical depictio of the fudametal differeces betwee the bracketig method (show o the right) ad the ope methods (show o the left). O the right, the root is costraied withi the iterval; o the left, a formula is used to project iteratively betwee x i ad x i+1 ad ca diverge or coverge rapidly. Numerical Fluid Mechaics PFJL Lecture 3, 24 24
25 Ope Methods (Fixed Poit Iteratio) Stop-criteria Urealistic stop-criteria: cotiue while Realistic stop-criteria Machie Epsilo Use a or combiatio of the two criteria f(x) flat f(x) f(x) steep f(x) x x Caot require Caot require Numerical Fluid Mechaics PFJL Lecture 3, 25 25
26 MIT OpeCourseWare Numerical Fluid Mechaics Fall 2011 For iformatio about citig these materials or our Terms of Use, visit:
THE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.
THE SOLUTION OF NONLINEAR EQUATIONS f( ) = 0. Noliear Equatio Solvers Bracketig. Graphical. Aalytical Ope Methods Bisectio False Positio (Regula-Falsi) Fied poit iteratio Newto Raphso Secat The root of
More informationCS537. Numerical Analysis and Computing
CS57 Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 Jauary 9 9 What is the Root May physical system ca be
More informationME NUMERICAL METHODS Fall 2007
ME - 310 NUMERICAL METHODS Fall 2007 Group 02 Istructor: Prof. Dr. Eres Söylemez (Rm C205, email:eres@metu.edu.tr ) Class Hours ad Room: Moday 13:40-15:30 Rm: B101 Wedesday 12:40-13:30 Rm: B103 Course
More informationRoot Finding COS 323
Root Fidig COS 323 Remider Sig up for Piazza Assigmet 0 is posted, due Tue 9/25 Last time.. Floatig poit umbers ad precisio Machie epsilo Sources of error Sesitivity ad coditioig Stability ad accuracy
More informationPC5215 Numerical Recipes with Applications - Review Problems
PC55 Numerical Recipes with Applicatios - Review Problems Give the IEEE 754 sigle precisio bit patter (biary or he format) of the followig umbers: 0 0 05 00 0 00 Note that it has 8 bits for the epoet,
More informationChapter 2 The Solution of Numerical Algebraic and Transcendental Equations
Chapter The Solutio of Numerical Algebraic ad Trascedetal Equatios Itroductio I this chapter we shall discuss some umerical methods for solvig algebraic ad trascedetal equatios. The equatio f( is said
More informationSection A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics
X0/70 NATIONAL QUALIFICATIONS 005 MONDAY, MAY.00 PM 4.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
More informationCS321. Numerical Analysis and Computing
CS Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 September 8 5 What is the Root May physical system ca
More information*X203/701* X203/701. APPLIED MATHEMATICS ADVANCED HIGHER Numerical Analysis. Read carefully
X0/70 NATIONAL QUALIFICATIONS 006 MONDAY, MAY.00 PM.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More informationNotes on iteration and Newton s method. Iteration
Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationMath 312 Lecture Notes One Dimensional Maps
Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 21-23 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More information(a) (b) All real numbers. (c) All real numbers. (d) None. to show the. (a) 3. (b) [ 7, 1) (c) ( 7, 1) (d) At x = 7. (a) (b)
Chapter 0 Review 597. E; a ( + )( + ) + + S S + S + + + + + + S lim + l. D; a diverges by the Itegral l k Test sice d lim [(l ) ], so k l ( ) does ot coverge absolutely. But it coverges by the Alteratig
More informationMIDTERM 3 CALCULUS 2. Monday, December 3, :15 PM to 6:45 PM. Name PRACTICE EXAM SOLUTIONS
MIDTERM 3 CALCULUS MATH 300 FALL 08 Moday, December 3, 08 5:5 PM to 6:45 PM Name PRACTICE EXAM S Please aswer all of the questios, ad show your work. You must explai your aswers to get credit. You will
More informationWe are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n
Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationMath 257: Finite difference methods
Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...
More informationENGI Series Page 6-01
ENGI 3425 6 Series Page 6-01 6. Series Cotets: 6.01 Sequeces; geeral term, limits, covergece 6.02 Series; summatio otatio, covergece, divergece test 6.03 Stadard Series; telescopig series, geometric series,
More informationMath 113 Exam 4 Practice
Math Exam 4 Practice Exam 4 will cover.-.. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationChapter 6 Infinite Series
Chapter 6 Ifiite Series I the previous chapter we cosidered itegrals which were improper i the sese that the iterval of itegratio was ubouded. I this chapter we are goig to discuss a topic which is somewhat
More informationLesson 10: Limits and Continuity
www.scimsacademy.com Lesso 10: Limits ad Cotiuity SCIMS Academy 1 Limit of a fuctio The cocept of limit of a fuctio is cetral to all other cocepts i calculus (like cotiuity, derivative, defiite itegrals
More informationME 501A Seminar in Engineering Analysis Page 1
Accurac, Stabilit ad Sstems of Equatios November 0, 07 Numerical Solutios of Ordiar Differetial Equatios Accurac, Stabilit ad Sstems of Equatios Larr Caretto Mecaical Egieerig 0AB Semiar i Egieerig Aalsis
More informationANALYSIS OF EXPERIMENTAL ERRORS
ANALYSIS OF EXPERIMENTAL ERRORS All physical measuremets ecoutered i the verificatio of physics theories ad cocepts are subject to ucertaities that deped o the measurig istrumets used ad the coditios uder
More informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
More informationComputational Fluid Dynamics. Lecture 3
Computatioal Fluid Dyamics Lecture 3 Discretizatio Cotiued. A fourth order approximatio to f x ca be foud usig Taylor Series. ( + ) + ( + ) + + ( ) + ( ) = a f x x b f x x c f x d f x x e f x x f x 0 0
More informationAssignment 1 : Real Numbers, Sequences. for n 1. Show that (x n ) converges. Further, by observing that x n+2 + x n+1
Assigmet : Real Numbers, Sequeces. Let A be a o-empty subset of R ad α R. Show that α = supa if ad oly if α is ot a upper boud of A but α + is a upper boud of A for every N. 2. Let y (, ) ad x (, ). Evaluate
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More informationResearch Article A New Second-Order Iteration Method for Solving Nonlinear Equations
Abstract ad Applied Aalysis Volume 2013, Article ID 487062, 4 pages http://dx.doi.org/10.1155/2013/487062 Research Article A New Secod-Order Iteratio Method for Solvig Noliear Equatios Shi Mi Kag, 1 Arif
More informationTwo-step Extrapolated Newton s Method with High Efficiency Index
Jour of Adv Research i Damical & Cotrol Systems Vol. 9 No. 017 Two-step Etrapolated Newto s Method with High Efficiecy Ide V.B. Kumar Vatti Dept. of Egieerig Mathematics Adhra Uiversity Visakhapatam Idia.
More information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
More informationRoot Finding COS 323
Root Fidig COS 33 Why Root Fidig? Solve or i ay equatio: b where? id root o g b 0 Might ot be able to solve or directly e.g., e -0. si3-0.5 Evaluatig might itsel require solvig a dieretial equatio, ruig
More informationQuiz. Use either the RATIO or ROOT TEST to determine whether the series is convergent or not.
Quiz. Use either the RATIO or ROOT TEST to determie whether the series is coverget or ot. e .6 POWER SERIES Defiitio. A power series i about is a series of the form c 0 c a c a... c a... a 0 c a where
More informationZeros of Polynomials
Math 160 www.timetodare.com 4.5 4.6 Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered with fidig the solutios of polyomial equatios of ay degree
More information6.003 Homework #3 Solutions
6.00 Homework # Solutios Problems. Complex umbers a. Evaluate the real ad imagiary parts of j j. π/ Real part = Imagiary part = 0 e Euler s formula says that j = e jπ/, so jπ/ j π/ j j = e = e. Thus the
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More informationNonlinear regression
oliear regressio How to aalyse data? How to aalyse data? Plot! How to aalyse data? Plot! Huma brai is oe the most powerfull computatioall tools Works differetly tha a computer What if data have o liear
More informationThe axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.
5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =
More information-ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION
NEW NEWTON-TYPE METHOD WITH k -ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION R. Thukral Padé Research Cetre, 39 Deaswood Hill, Leeds West Yorkshire, LS7 JS, ENGLAND ABSTRACT The objective
More informationMath 113, Calculus II Winter 2007 Final Exam Solutions
Math, Calculus II Witer 7 Fial Exam Solutios (5 poits) Use the limit defiitio of the defiite itegral ad the sum formulas to compute x x + dx The check your aswer usig the Evaluatio Theorem Solutio: I this
More informationThe Growth of Functions. Theoretical Supplement
The Growth of Fuctios Theoretical Supplemet The Triagle Iequality The triagle iequality is a algebraic tool that is ofte useful i maipulatig absolute values of fuctios. The triagle iequality says that
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 informationTime-Domain Representations of LTI Systems
2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costat-coefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. State-variable
More informationMa 530 Introduction to Power Series
Ma 530 Itroductio to Power Series Please ote that there is material o power series at Visual Calculus. Some of this material was used as part of the presetatio of the topics that follow. What is a Power
More informationx x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula
NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS. If g( is cotiuous i [a,b], te uder wat coditio te iterative (or iteratio metod = g( as a uique solutio i [a,b]? '( i [a,b].. Wat
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam will cover.-.9. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for you
More informationNumerical Solution of Non-linear Equations
Numerical Solutio of Noliear Equatios. INTRODUCTION The most commo reallife problems are oliear ad are ot ameable to be hadled by aalytical methods to obtai solutios of a variety of mathematical problems.
More informationREGRESSION (Physics 1210 Notes, Partial Modified Appendix A)
REGRESSION (Physics 0 Notes, Partial Modified Appedix A) HOW TO PERFORM A LINEAR REGRESSION Cosider the followig data poits ad their graph (Table I ad Figure ): X Y 0 3 5 3 7 4 9 5 Table : Example Data
More informationOptimization Methods MIT 2.098/6.255/ Final exam
Optimizatio Methods MIT 2.098/6.255/15.093 Fial exam Date Give: December 19th, 2006 P1. [30 pts] Classify the followig statemets as true or false. All aswers must be well-justified, either through a short
More informationSome New Iterative Methods for Solving Nonlinear Equations
World Applied Scieces Joural 0 (6): 870-874, 01 ISSN 1818-495 IDOSI Publicatios, 01 DOI: 10.589/idosi.wasj.01.0.06.830 Some New Iterative Methods for Solvig Noliear Equatios Muhammad Aslam Noor, Khalida
More informationMATH301 Real Analysis (2008 Fall) Tutorial Note #7. k=1 f k (x) converges pointwise to S(x) on E if and
MATH01 Real Aalysis (2008 Fall) Tutorial Note #7 Sequece ad Series of fuctio 1: Poitwise Covergece ad Uiform Covergece Part I: Poitwise Covergece Defiitio of poitwise covergece: A sequece of fuctios f
More informationECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015
ECE 8527: Itroductio to Machie Learig ad Patter Recogitio Midterm # 1 Vaishali Ami Fall, 2015 tue39624@temple.edu Problem No. 1: Cosider a two-class discrete distributio problem: ω 1 :{[0,0], [2,0], [2,2],
More informationCarleton College, Winter 2017 Math 121, Practice Final Prof. Jones. Note: the exam will have a section of true-false questions, like the one below.
Carleto College, Witer 207 Math 2, Practice Fial Prof. Joes Note: the exam will have a sectio of true-false questios, like the oe below.. True or False. Briefly explai your aswer. A icorrectly justified
More informationFinally, we show how to determine the moments of an impulse response based on the example of the dispersion model.
5.3 Determiatio of Momets Fially, we show how to determie the momets of a impulse respose based o the example of the dispersio model. For the dispersio model we have that E θ (θ ) curve is give by eq (4).
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationFINALTERM EXAMINATION Fall 9 Calculus & Aalytical Geometry-I Questio No: ( Mars: ) - Please choose oe Let f ( x) is a fuctio such that as x approaches a real umber a, either from left or right-had-side,
More informationThe Method of Least Squares. To understand least squares fitting of data.
The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve
More informationZ ß cos x + si x R du We start with the substitutio u = si(x), so du = cos(x). The itegral becomes but +u we should chage the limits to go with the ew
Problem ( poits) Evaluate the itegrals Z p x 9 x We ca draw a right triagle labeled this way x p x 9 From this we ca read off x = sec, so = sec ta, ad p x 9 = R ta. Puttig those pieces ito the itegralrwe
More informationNumerical Methods in Geophysics: Implicit Methods
Numerical Methods i Geophysics: What is a implicit scheme? Explicit vs. implicit scheme for Newtoia oolig rak-nicholso Scheme (mixed explicit-implicit Explicit vs. implicit for the diffusio equatio Relaxatio
More information1.010 Uncertainty in Engineering Fall 2008
MIT OpeCourseWare http://ocw.mit.edu.00 Ucertaity i Egieerig Fall 2008 For iformatio about citig these materials or our Terms of Use, visit: http://ocw.mit.edu.terms. .00 - Brief Notes # 9 Poit ad Iterval
More informationMTH 142 Exam 3 Spr 2011 Practice Problem Solutions 1
MTH 42 Exam 3 Spr 20 Practice Problem Solutios No calculators will be permitted at the exam. 3. A pig-pog ball is lauched straight up, rises to a height of 5 feet, the falls back to the lauch poit ad bouces
More informationMAT1026 Calculus II Basic Convergence Tests for Series
MAT026 Calculus II Basic Covergece Tests for Series Egi MERMUT 202.03.08 Dokuz Eylül Uiversity Faculty of Sciece Departmet of Mathematics İzmir/TURKEY Cotets Mootoe Covergece Theorem 2 2 Series of Real
More information3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials
Math 60 www.timetodare.com 3. Properties of Divisio 3.3 Zeros of Polyomials 3.4 Complex ad Ratioal Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered
More informationAP Calculus BC Review Applications of Derivatives (Chapter 4) and f,
AP alculus B Review Applicatios of Derivatives (hapter ) Thigs to Kow ad Be Able to Do Defiitios of the followig i terms of derivatives, ad how to fid them: critical poit, global miima/maima, local (relative)
More informationMTH Assignment 1 : Real Numbers, Sequences
MTH -26 Assigmet : Real Numbers, Sequeces. Fid the supremum of the set { m m+ : N, m Z}. 2. Let A be a o-empty subset of R ad α R. Show that α = supa if ad oly if α is ot a upper boud of A but α + is a
More informationSOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.
SOLUTIONS TO EXAM 3 Problem Fid the sum of the followig series 2 + ( ) 5 5 2 5 3 25 2 2 This series diverges Solutio: Note that this defies two coverget geometric series with respective radii r 2/5 < ad
More informationx x x Using a second Taylor polynomial with remainder, find the best constant C so that for x 0,
Math Activity 9( Due with Fial Eam) Usig first ad secod Taylor polyomials with remaider, show that for, 8 Usig a secod Taylor polyomial with remaider, fid the best costat C so that for, C 9 The th Derivative
More informationMATH 31B: MIDTERM 2 REVIEW
MATH 3B: MIDTERM REVIEW JOE HUGHES. Evaluate x (x ) (x 3).. Partial Fractios Solutio: The umerator has degree less tha the deomiator, so we ca use partial fractios. Write x (x ) (x 3) = A x + A (x ) +
More informationSection 5.5. Infinite Series: The Ratio Test
Differece Equatios to Differetial Equatios Sectio 5.5 Ifiite Series: The Ratio Test I the last sectio we saw that we could demostrate the covergece of a series a, where a 0 for all, by showig that a approaches
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This
More information15.081J/6.251J Introduction to Mathematical Programming. Lecture 21: Primal Barrier Interior Point Algorithm
508J/65J Itroductio to Mathematical Programmig Lecture : Primal Barrier Iterior Poit Algorithm Outlie Barrier Methods Slide The Cetral Path 3 Approximatig the Cetral Path 4 The Primal Barrier Algorithm
More informationAdvanced Analysis. Min Yan Department of Mathematics Hong Kong University of Science and Technology
Advaced Aalysis Mi Ya Departmet of Mathematics Hog Kog Uiversity of Sciece ad Techology September 3, 009 Cotets Limit ad Cotiuity 7 Limit of Sequece 8 Defiitio 8 Property 3 3 Ifiity ad Ifiitesimal 8 4
More informationAdditional Notes on Power Series
Additioal Notes o Power Series Mauela Girotti MATH 37-0 Advaced Calculus of oe variable Cotets Quick recall 2 Abel s Theorem 2 3 Differetiatio ad Itegratio of Power series 4 Quick recall We recall here
More informationWe will conclude the chapter with the study a few methods and techniques which are useful
Chapter : Coordiate geometry: I this chapter we will lear about the mai priciples of graphig i a dimesioal (D) Cartesia system of coordiates. We will focus o drawig lies ad the characteristics of the graphs
More informationPRELIM PROBLEM SOLUTIONS
PRELIM PROBLEM SOLUTIONS THE GRAD STUDENTS + KEN Cotets. Complex Aalysis Practice Problems 2. 2. Real Aalysis Practice Problems 2. 4 3. Algebra Practice Problems 2. 8. Complex Aalysis Practice Problems
More informationDECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan
Mathematical ad Computatioal Applicatios, Vol. 9, No. 3, pp. 30-40, 04 DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS Muhammad Aslam Noor, Khalida Iayat Noor ad Asif Waheed
More informationSIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road QUESTION BANK (DESCRIPTIVE)
QUESTION BANK 8 SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayaavaam Road 5758 QUESTION BANK (DESCRIPTIVE Subject with Code : (6HS6 Course & Brach: B.Tech AG Year & Sem: II-B.Tech& I-Sem
More informationMath 128A: Homework 1 Solutions
Math 8A: Homework Solutios Due: Jue. Determie the limits of the followig sequeces as. a) a = +. lim a + = lim =. b) a = + ). c) a = si4 +6) +. lim a = lim = lim + ) [ + ) ] = [ e ] = e 6. Observe that
More informationPAPER : IIT-JAM 2010
MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure
More information4x 2. (n+1) x 3 n+1. = lim. 4x 2 n+1 n3 n. n 4x 2 = lim = 3
Exam Problems (x. Give the series (, fid the values of x for which this power series coverges. Also =0 state clearly what the radius of covergece is. We start by settig up the Ratio Test: x ( x x ( x x
More informationPolynomial Functions. New Section 1 Page 1. A Polynomial function of degree n is written is the form:
New Sectio 1 Page 1 A Polyomial fuctio of degree is writte is the form: 1 P x a x a x a x a x a x a 1 1 0 where is a o-egative iteger expoet ad a 0 ca oly take o values where a, a, a,..., a, a. 0 1 1.
More informationCalculus 2 - D. Yuen Final Exam Review (Version 11/22/2017. Please report any possible typos.)
Calculus - D Yue Fial Eam Review (Versio //7 Please report ay possible typos) NOTE: The review otes are oly o topics ot covered o previous eams See previous review sheets for summary of previous topics
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationChapter 7: The z-transform. Chih-Wei Liu
Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability
More information1. (25 points) Use the limit definition of the definite integral and the sum formulas 1 to compute
Math, Calculus II Fial Eam Solutios. 5 poits) Use the limit defiitio of the defiite itegral ad the sum formulas to compute 4 d. The check your aswer usig the Evaluatio Theorem. ) ) Solutio: I this itegral,
More informationMath 341 Lecture #31 6.5: Power Series
Math 341 Lecture #31 6.5: Power Series We ow tur our attetio to a particular kid of series of fuctios, amely, power series, f(x = a x = a 0 + a 1 x + a 2 x 2 + where a R for all N. I terms of a series
More informationx c the remainder is Pc ().
Algebra, Polyomial ad Ratioal Fuctios Page 1 K.Paulk Notes Chapter 3, Sectio 3.1 to 3.4 Summary Sectio Theorem Notes 3.1 Zeros of a Fuctio Set the fuctio to zero ad solve for x. The fuctio is zero at these
More informationIt is often useful to approximate complicated functions using simpler ones. We consider the task of approximating a function by a polynomial.
Taylor Polyomials ad Taylor Series It is ofte useful to approximate complicated fuctios usig simpler oes We cosider the task of approximatig a fuctio by a polyomial If f is at least -times differetiable
More informationStudy the bias (due to the nite dimensional approximation) and variance of the estimators
2 Series Methods 2. Geeral Approach A model has parameters (; ) where is ite-dimesioal ad is oparametric. (Sometimes, there is o :) We will focus o regressio. The fuctio is approximated by a series a ite
More informationProblem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =
Mathematics Summer Wilso Fial Exam August 8, ANSWERS Problem 1 (a) Fid the solutio to y +x y = e x x that satisfies y() = 5 : This is already i the form we used for a first order liear differetial equatio,
More informationIntroduction to Signals and Systems, Part V: Lecture Summary
EEL33: Discrete-Time Sigals ad Systems Itroductio to Sigals ad Systems, Part V: Lecture Summary Itroductio to Sigals ad Systems, Part V: Lecture Summary So far we have oly looked at examples of o-recursive
More informationSeptember 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1
September 0 s (Edecel) Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright
More informationNumerical Methods for Ordinary Differential Equations
Numerical Methods for Ordiary Differetial Equatios Braislav K. Nikolić Departmet of Physics ad Astroomy, Uiversity of Delaware, U.S.A. PHYS 460/660: Computatioal Methods of Physics http://www.physics.udel.edu/~bikolic/teachig/phys660/phys660.html
More informationHigher-order iterative methods by using Householder's method for solving certain nonlinear equations
Math Sci Lett, No, 7- ( 7 Mathematical Sciece Letters A Iteratioal Joural http://dxdoiorg/785/msl/5 Higher-order iterative methods by usig Householder's method for solvig certai oliear equatios Waseem
More informationMA131 - Analysis 1. Workbook 3 Sequences II
MA3 - Aalysis Workbook 3 Sequeces II Autum 2004 Cotets 2.8 Coverget Sequeces........................ 2.9 Algebra of Limits......................... 2 2.0 Further Useful Results........................
More information