Computational Methods for Engineers Programming in Engineering Problem Solving
|
|
- Leslie Cox
- 6 years ago
- Views:
Transcription
1 Computational Methods for Engineers Programming in Engineering Problem Solving Abu Hasan Abdullah January 6, 2009 Abu Hasan Abdullah 2009
2 An Engineering Problem Problem Statement: The length of a belt in an open-belt drive, L, is given by L = q 4c 2 (D d) `DθD + dθ d 2 (1) where «D d θ D = π + 2 sin 1 2c «D d θ d = π 2 sin 1 2c 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, see Figure-2.8 of Rao (2002). 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. Abu Hasan Abdullah
3 An Engineering Problem Figure 1: Open belt drive. Abu Hasan Abdullah
4 Analysis of Engineering Problem 1. Problem Statement: Recognise and understand the problem (what is it that needed to be solved?). 2. Governing Equations or Mathematical Models: Identify parameters affecting the problem, make the necessary assumptions, develop mathematical model or governing equations (based on theories from Engineering Mathematics and other Engineering Subjects). 3. Solution: Solution of the governing equations may make use of the computer programming (why?). 4. Verification: Verify and interpret the solution (right/wrong?). Abu Hasan Abdullah
5 Analysis of Engineering Problem Problem Statement The length of a belt in an open-belt drive, L, is given by q L = 4c 2 (D d) `DθD + dθ d 2 (2) where «D d θ D = π + 2 sin 1 2c «D d θ d = π 2 sin 1 2c 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, see Figure-2.8 of Rao (2002). 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. Abu Hasan Abdullah
6 Analysis of Engineering Problem Problem Statement Figure 2: Open belt drive. Abu Hasan Abdullah
7 Analysis of Engineering Problem Mathematical Model Defined as a formulation or equation that expresses the essential features of a physical system or process in mathematical terms. Its simplest form can be represented as a functional relationship thus where Dependent variable = f(independent variables, parameters, forcing functions) dependent variable: a characteristic that reflects the behaviour/state of system independent variables: dimensions (time, space, mass) along which the system s behaviour that is being determined parameters: reflective of system s properties or composition forcing functions: external influences acting on the system Abu Hasan Abdullah
8 Analysis of Engineering Problem Mathematical Model Mathematical model ranges from a simple algebraic relationship to large complicated set of DE. Mathematical models (a.k.a. governing equations) are derived by applying physical laws such as Equilibrium Equation Newton s Law of Motion Conservation Laws: Mass, Momentum, Energy Equation of State Abu Hasan Abdullah
9 Analysis of Engineering Problem Solution Solution of the governing equation or mathematical model may appear as Transcendental Functions Linear or Nonlinear Algebraic Equations Homogeneous Equations leading to an Eigenvalue Problem Ordinary or Partial Differential Equations Equations involving Integrals or Derivatives which are either closed-form or open-ended. Abu Hasan Abdullah
10 Analysis of Engineering Problem Solution: Closed-form Closed-form mathematical expression I 1 = b a xe x2 dx = = 1 2 e b e a2 = 1 2 [ 12 e x2] b a (e a2 e b2) leads to analytical solution Abu Hasan Abdullah
11 Analysis of Engineering Problem Solution: Open-ended Open-ended mathematical expressions, e.g. I 1 = Z b a need to be approximated numerically f(x)dx = Z b a e x2 dx Nowadays, approximated numerical solutions are done by developing a computer program. Because numerical methods deal extensively with approximations connected with the manipulation of numbers, accuracy, precision and error feature prominently in programming the solution. We shall cover these later! Abu Hasan Abdullah
12 Analysis of Engineering Problem Solution: Computer Program Steps in computer program development: Algorithm Design: Listing down of the sequence of steps to define the problem at hand. Techniques available: algorithm, flowchart, pseudocode Program Coding: Writing these steps in a computer language. Debugging: Testing the program to ensure that it is error-free and reliable. Documentation: Making the program easy to understand and use through manual or guide. See SME 1013 Programming for Engineers for details. Abu Hasan Abdullah
13 Analysis of Engineering Problem Solution: Computer Program Algorithm Algorithm: A general sequence of the logical steps in solving a specific problem. Flowchart: A graphical representation of the algorithm. Better suited for visualizing complex algorithms. Pseudocode: Uses code-like statements in place of the graphical symbols of flowchart. Easier to develop a program with it than with a flowchart. Elements of good algorithm Each step must be deterministic i.e. not ambiguous. The process must end after a finite number of steps. The algorithm must be general enough to deal with any contingency. Abu Hasan Abdullah
14 Analysis of Engineering Problem Solution: Computer Program Algorithm Find roots of equation ax 2 + bx + c = 0 using the quadratic formula x = b ± b 2 4ac 2a 1. Start 2. Read coefficients a, b and c 3. Implement quadratic formula (Avoid division by zero, allow for complex roots) 4. Display solution i.e. values of x 5. Stop Abu Hasan Abdullah
15 Analysis of Engineering Problem Solution: Computer Program Pseudocode Find roots of equation ax 2 + bx + c = 0 using the quadratic formula x = b ± b 2 4ac 2a DO READ a, b, c root1 = (-b + SQRT(b^2-4ac)/(2a) root2 = (-b - SQRT(b^2-4ac)/(2a) PRINT root1, root2 PRINT Try again? Answer yes or no READ response IF response = no EXIT ENDDO Abu Hasan Abdullah
16 Analysis of Engineering Problem Solution: Computer Program Flowchart Start/Stop Process Input/Output Refers to separate flowchart Decision Connector Off-page Connector Preparation (for Loop, etc.) Abu Hasan Abdullah
17 Analysis of Engineering Problem Solution: Computer Program Coding A program is a sequence of instructions to the computer for it to solve a particular problem. A set of programs is called code. Programs are written in some programming language, e.g. Fortran, Matlab, Basic, C++, Pascal, Java. Programs are stored in files which are a sequence of bytes which is given a name and stored on a disk. Abu Hasan Abdullah
18 Analysis of Engineering Problem Solution: Computer Program Coding A program is a file containing a sequence of statements, each of which tells the computer to do a specific action. Once a program is run or executed the commands are followed and actions occur in a sequential manner. If the program is designed to interact with the outside world, then it must have input and output. Abu Hasan Abdullah
19 Analysis of Engineering Problem Solution: Computer Program Coding A program is said to have a bug if it contains a mistake or it does not function in the way it is intended to. Bugs can happen both in the logic of the program, as well as in the commands. In order that the program perform the exact actions it is intended to do, before the actual program is written an algorithm for solving the problem must first be outlined. Abu Hasan Abdullah
20 Example Problem 1 Problem Statement: Assuming that the thrust T of a screw propeller is dependent upon diameter D, speed of advance v, fluid density ρ, rotational speed of propeller N and coefficient of viscosity µ, derive and expression that relates all the parameters involved and solve for T. Governing Equation: Through dimensional analysis ( µ T = ρv 2 D 2 f ρvd, ND ) v Abu Hasan Abdullah
21 Example Problem 2 Problem Statement: Given temperature in degrees Fahrenheit, the temperature in degrees Kelvin is to be computed and shown. Governing Equation: From Physics, these two temperature scales are related through ( ) TF 32 T k = and the parameters involved in this problem are T K and T F Abu Hasan Abdullah
22 Example Problem 2 Algorithm: 1. Start 2. Get the temperature in Fahrenheit 3. Compute the temperature in Kelvin using the formula: T k = ( ) TF Show the temperature in Kelvin 5. Stop Abu Hasan Abdullah
23 Example Problem 2 Pseudocode: Start Read TF TK = (TF-32)/ Print TK Stop Abu Hasan Abdullah
24 Example Problem 3 Problem Statement: Determine the mass of the bungee jumper with a drag coefficient of 0.25 kg/m to have a velocity of 36 m/s after 4 s of free fall. Note: The acceleration of gravity is 9.81 m/s. Governing Equation: From Physics and Mechanics the fall velocity as function of time is ( ) gm gcd v(t) = tanh c d m t or, on re-arranging in terms of a function of mass, f(m), f(m) = ( ) gm gcd tanh c d m t v(t) = 0 Abu Hasan Abdullah
25 Example Problem 3 Algorithm: The governing equation indicates a root-finding problem. 1. Choose lower m_l and upper m_u guesses for the root such that the function changes sign over the interval. This can be checked by ensuring that f(m_l) x f(m_u) < An estimate of the root m_r is determined by m R = m L + m U 2 3. Make the following evaluations to determine in which subinterval the root lies: If f(m_l) x f(m_u) < 0, root lies in lower subinterval; repeat step 2. If f(m_l) x f(m_u) > 0, root lies in upper subinterval; repeat step 2. If f(m_l) x f(m_u) = 0, root equals m_r; terminate computation. Abu Hasan Abdullah
26 Example Problem 3 Computer Program: Fortran code for root finding C*********************************************************************** C BISECTION METHOD * C*********************************************************************** C TO FIND A SOLUTION TO F(X)=0 GIVEN THE CONTINOUS FUNCTION C F ON THE INTERVAL <A,B>, WHERE F(A) AND F(B) HAVE C OPPOSITE SIGNS: C C INPUT: ENDPOINTS A,B; TOLERANCE TOL; C MAXIMUM INTERATIONS N0. C OUTPUT: APPROXIMATE SOLUTION P OR A C MESSAGE THAT THE ALGORITHM FAILS. C CHARACTER NAME*14,NAME1*14,AA*1 INTEGER INP,OUP,FLAG LOGICAL OK REAL A,B,FA,FB,X,TOL INTEGER N0 Abu Hasan Abdullah
27 C*********************************************************************** C DEFINE YOUR FUNCTION HERE * C*********************************************************************** C DEFINE F C! F(X)=X**3+4*X**2-10 c c OPEN(UNIT=5,FILE= CON,ACCESS= SEQUENTIAL ) OPEN(UNIT=6,FILE= CON,ACCESS= SEQUENTIAL ) WRITE(6,*) This is the Bisection Method. WRITE(6,*) Has the function F been created in the program? WRITE(6,*) Enter Y or N WRITE(6,*) READ(5,*) AA IF(( AA.EQ. Y ).OR. ( AA.EQ. y )) THEN OK =.FALSE. 10 IF (OK) GOTO 11 WRITE(6,*) Input endpoints A < B separated by blank WRITE(6,*) READ(5,*) A, B IF (A.GT.B) THEN X = A A = B Abu Hasan Abdullah
28 B = X ENDIF IF (A.EQ.B) THEN WRITE(6,*) A cannot equal B WRITE(6,*) ELSE FA = F( A ) FB = F( B ) IF ( FA * FB.GT. 0.0 ) THEN WRITE(6,*) F(A) and F(B) have same sign WRITE(6,*) ELSE OK =.TRUE. ENDIF ENDIF GOTO OK =.FALSE. 12 IF (OK) GOTO 13 WRITE(6,*) Input tolerance WRITE(6,*) READ(5,*) TOL Abu Hasan Abdullah
29 IF (TOL.LE.0.0) THEN WRITE(6,*) Tolerance must be positive WRITE(6,*) ELSE OK =.TRUE. ENDIF GOTO OK =.FALSE. 14 IF (OK) GOTO 15 WRITE(6,*) Input maximum number of iterations WRITE(6,*) - no decimal point WRITE(6,*) READ(5,*) N0 IF ( N0.LE. 0 ) THEN WRITE(6,*) Must be positive integer WRITE(6,*) ELSE OK =.TRUE. ENDIF GOTO CONTINUE Abu Hasan Abdullah
30 ELSE WRITE(6,*) The program will end so that the function F WRITE(6,*) can be created OK =.FALSE. ENDIF IF (.NOT.OK) GOTO 40 WRITE(6,*) Select output destinations: WRITE(6,*) 1. Screen WRITE(6,*) 2. Text file WRITE(6,*) Enter 1 or 2 WRITE(6,*) READ(5,*) FLAG IF ( FLAG.EQ. 2 ) THEN WRITE(6,*) Input the file name in the form - WRITE(6,*) drive:name.ext WRITE(6,*) with the name contained within quotes WRITE(6,*) as example: A:OUTPUT.DTA WRITE(6,*) READ(5,*) NAME1 OUP = 3 OPEN(UNIT=OUP,FILE=NAME1,STATUS= NEW ) Abu Hasan Abdullah
31 ELSE OUP = 6 ENDIF WRITE(6,*) Select amount of output WRITE(6,*) 1. Answer only WRITE(6,*) 2. All intermediate approximations WRITE(6,*) Enter 1 or 2 WRITE(6,*) READ(5,*) FLAG WRITE(OUP,*) BISECTION METHOD IF (FLAG.EQ.2) THEN WRITE(OUP,004) 004 FORMAT(3X, I,15X, P,12X, F(P) ) ENDIF C STEP 1 I=1 C STEP IF (I.GT.N0) GOTO 020 C STEP 3 C COMPUTE P(I) P=A+(B-A)/2 Abu Hasan Abdullah
32 FP=F(P) IF (FLAG.EQ.2) THEN WRITE(OUP,005) I,P,FP 005 FORMAT(1X,I3,2X,E15.8,2X,E15.8) ENDIF C STEP 4 C C STEP 5 I=I+1 C STEP 6 C IF( ABS(FP).LE.1.0E-20.OR. (B-A)/2.LT. TOL) THEN PROCEDURE COMPLETED SUCCESSFULLY WRITE(OUP,002) P, I, TOL GOTO 040 ENDIF COMPUTE A(I) AND B(I) IF( F(A)*FP.GT. 0) THEN A=P ELSE B=P ENDIF GOTO 016 Abu Hasan Abdullah
33 020 CONTINUE C STEP 7 C PROCEDURE COMPLETED UNSUCCESSFULLY WRITE(OUP,003) N0 IF(OUP.NE.6) WRITE(6,003) N0 040 CLOSE(UNIT=5) CLOSE(UNIT=OUP) IF (OUP.NE.6) CLOSE(UNIT=6) STOP 002 FORMAT(1X, THE APPROXIMATE SOLUTION IS,1X *,E15.8,1X, AFTER,1X,I2,1X, ITERATIONS, WITH TOLERANCE *,1X,E15.8) 003 FORMAT(1X, THE METHOD FAILS AFTER,1X,I2,1X, *ITERATIONS ) END Abu Hasan Abdullah
34 References 1. Singiresu S. Rao (2002): Applied Numerical Methods for Engineers and Scientists, ISBN X, Prentice Hall 2. Steven C. Chapra, Raymond P. Canale (2006): Applied Numerical Methods for Engineers and Scientists, 10ed, ISBN , McGraw-Hill 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 Abu Hasan Abdullah
SME 3023 Applied Numerical Methods
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
More informationSKMM 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 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 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 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 informationLecture 44. Better and successive approximations x2, x3,, xn to the root are obtained from
Lecture 44 Solution of Non-Linear Equations Regula-Falsi Method Method of iteration Newton - Raphson Method Muller s Method Graeffe s Root Squaring Method Newton -Raphson Method An approximation to the
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 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 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 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 informationVariable. Peter W. White Fall 2018 / Numerical Analysis. Department of Mathematics Tarleton State University
Newton s Iterative s Peter W. White white@tarleton.edu Department of Mathematics Tarleton State University Fall 2018 / Numerical Analysis Overview Newton s Iterative s Newton s Iterative s Newton s Iterative
More informationGoals for This Lecture:
Goals for This Lecture: Learn the Newton-Raphson method for finding real roots of real functions Learn the Bisection method for finding real roots of a real function Look at efficient implementations of
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 informationNumerical Methods School of Mechanical Engineering Chung-Ang University
Part 2 Chapter 5 Roots: Bracketing Methods Prof. Hae-Jin Choi hjchoi@cau.ac.kr 1 Overview of Part 2 l To find the roots of general second order polynomial, the quadratic formula is used b b 2 2 - ± - 4ac
More informationIntro to Scientific Computing: How long does it take to find a needle in a haystack?
Intro to Scientific Computing: How long does it take to find a needle in a haystack? Dr. David M. Goulet Intro Binary Sorting Suppose that you have a detector that can tell you if a needle is in a haystack,
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 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 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 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 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 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 informationCS 221 Lecture 9. Tuesday, 1 November 2011
CS 221 Lecture 9 Tuesday, 1 November 2011 Some slides in this lecture are from the publisher s slides for Engineering Computation: An Introduction Using MATLAB and Excel 2009 McGraw-Hill Today s Agenda
More informationCOMPUTER PROGRAMMING
Prof. Dr. Namık Kemal ÖZTORUN Lecture Notes for Computer Programming Course Page 1 / 36 COMPUTER PROGRAMMING References: WITH FORTRAN Lecture Notes prepared by Prof. Dr. Namık Kemal ÖZTORUN Dr. Faruk TOKDEMİR,
More informationChapter 2. Review of Mathematics. 2.1 Exponents
Chapter 2 Review of Mathematics In this chapter, we will briefly review some of the mathematical concepts used in this textbook. Knowing these concepts will make it much easier to understand the mathematical
More informationbc7f2306 Page 1 Name:
Name: Questions 1 through 4 refer to the following: Solve the given inequality and represent the solution set using set notation: 1) 3x 1 < 2(x + 4) or 7x 3 2(x + 1) Questions 5 and 6 refer to the following:
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 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 informationChapter 8. Exploring Polynomial Functions. Jennifer Huss
Chapter 8 Exploring Polynomial Functions Jennifer Huss 8-1 Polynomial Functions The degree of a polynomial is determined by the greatest exponent when there is only one variable (x) in the polynomial Polynomial
More informationNumerical Methods Dr. Sanjeev Kumar Department of Mathematics Indian Institute of Technology Roorkee Lecture No 7 Regula Falsi and Secant Methods
Numerical Methods Dr. Sanjeev Kumar Department of Mathematics Indian Institute of Technology Roorkee Lecture No 7 Regula Falsi and Secant Methods So welcome to the next lecture of the 2 nd unit of this
More informationINTRODUCTION, FOUNDATIONS
1 INTRODUCTION, FOUNDATIONS 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 Information
More informationBisected Direct Quadratic Regula Falsi
Applied Mathematical Sciences, Vol. 4, 2010, no. 15, 709-718 Bisected Direct Quadratic Regula Falsi Robert G. Gottlieb 1 and Blair F. Thompson Odyssey Space Research 1120 NASA Parkway, Suite 505 Houston,
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 informationMath Introduction to Numerical Methods - Winter 2011 Homework 2 Assigned: Friday, January 14, Due: Thursday, January 27,
Math 371 - Introduction to Numerical Methods - Winter 2011 Homework 2 Assigned: Friday, January 14, 2011. Due: Thursday, January 27, 2011.. Include a cover page. You do not need to hand in a problem sheet.
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 informationLearning Objectives for Math 165
Learning Objectives for Math 165 Chapter 2 Limits Section 2.1: Average Rate of Change. State the definition of average rate of change Describe what the rate of change does and does not tell us in a given
More informationCHAPTER 2 EXTRACTION OF THE QUADRATICS FROM REAL ALGEBRAIC POLYNOMIAL
24 CHAPTER 2 EXTRACTION OF THE QUADRATICS FROM REAL ALGEBRAIC POLYNOMIAL 2.1 INTRODUCTION Polynomial factorization is a mathematical problem, which is often encountered in applied sciences and many of
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 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 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 informationChapter 8 Linear Algebraic Equations
PowerPoint to accompany Introduction to MATLAB for Engineers, Third Edition William J. Palm III Chapter 8 Linear Algebraic Equations Copyright 2010. The McGraw-Hill Companies, Inc. This work is only for
More informationCLASS NOTES Computational Methods for Engineering Applications I Spring 2015
CLASS NOTES Computational Methods for Engineering Applications I Spring 2015 Petros Koumoutsakos Gerardo Tauriello (Last update: July 2, 2015) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material
More informationFurther Pure Mathematics FP1
Write your name here Surname Other names Pearson Edexcel International Advanced Level Centre Number Further Pure Mathematics FP1 Advanced/Advanced Subsidiary Wednesday 29 January 2014 Morning Time: 1 hour
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 informationMath 2142 Homework 5 Part 1 Solutions
Math 2142 Homework 5 Part 1 Solutions Problem 1. For the following homogeneous second order differential equations, give the general solution and the particular solution satisfying the given initial conditions.
More informationOutline. Math Numerical Analysis. Intermediate Value Theorem. Lecture Notes Zeros and Roots. Joseph M. Mahaffy,
Outline Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research
More informationMath Numerical Analysis
Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationOutline. Additional Nonlinear Systems. Abstract. Finding Equilibrium Points Numerically. Newton s Method
Outline Finding Equilibrium Points Numerically Additional Nonlinear Systems James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University June 13, 2017
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... Introduce the topic of numerical methods Consider the Error analysis and sources of errors Introduction A numerical method which
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 informationSTOP, a i+ 1 is the desired root. )f(a i) > 0. Else If f(a i+ 1. Set a i+1 = a i+ 1 and b i+1 = b Else Set a i+1 = a i and b i+1 = a i+ 1
53 17. Lecture 17 Nonlinear Equations Essentially, the only way that one can solve nonlinear equations is by iteration. The quadratic formula enables one to compute the roots of p(x) = 0 when p P. Formulas
More informationNumerical Methods in Informatics
Numerical Methods in Informatics Lecture 2, 30.09.2016: Nonlinear Equations in One Variable http://www.math.uzh.ch/binf4232 Tulin Kaman Institute of Mathematics, University of Zurich E-mail: tulin.kaman@math.uzh.ch
More informationCompanion. Jeffrey E. Jones
MATLAB7 Companion 1O11OO1O1O1OOOO1O1OO1111O1O1OO 1O1O1OO1OO1O11OOO1O111O1O1O1O1 O11O1O1O11O1O1O1O1OO1O11O1O1O1 O1O1O1111O11O1O1OO1O1O1O1OOOOO O1111O1O1O1O1O1O1OO1OO1OO1OOO1 O1O11111O1O1O1O1O Jeffrey E.
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 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 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 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 informationCore Mathematics C4. You must have: Mathematical Formulae and Statistical Tables (Pink)
Write your name here Surname Other names Pearson Edexcel GCE Centre Number Core Mathematics C4 Advanced Candidate Number Friday 23 June 2017 Morning Time: 1 hour 30 minutes Paper Reference 6666/01 You
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 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 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 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 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 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 informationGoal: Approximate the area under a curve using the Rectangular Approximation Method (RAM) RECTANGULAR APPROXIMATION METHODS
AP Calculus 5. Areas and Distances Goal: Approximate the area under a curve using the Rectangular Approximation Method (RAM) Exercise : Calculate the area between the x-axis and the graph of y = 3 2x.
More informationFP1 practice papers A to G
FP1 practice papers A to G Paper Reference(s) 6667/01 Edexcel GCE Further Pure Mathematics FP1 Advanced Subsidiary Practice Paper A Time: 1 hour 30 minutes Materials required for examination Mathematical
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 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 informationTop Down Design. Gunnar Gotshalks 03-1
Top Down Design 03-1 Top Down Description Top down is also known as step wise refinement and functional decomposition Given an operation, there are only the following three choices for refinement» Sequence
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationhexadecimal-to-decimal conversion
OTHER NUMBER SYSTEMS: octal (digits 0 to 7) group three binary numbers together and represent as base 8 3564 10 = 110 111 101 100 2 = (6X8 3 ) + (7X8 2 ) + (5X8 1 ) + (4X8 0 ) = 6754 8 hexadecimal (digits
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 informationExaminers Report/ Principal Examiner Feedback. June GCE Core Mathematics C2 (6664) Paper 1
Examiners Report/ Principal Examiner Feedback June 011 GCE Core Mathematics C (6664) Paper 1 Edexcel is one of the leading examining and awarding bodies in the UK and throughout the world. We provide a
More information2 Discrete Dynamical Systems (DDS)
2 Discrete Dynamical Systems (DDS) 2.1 Basics A Discrete Dynamical System (DDS) models the change (or dynamics) of single or multiple populations or quantities in which the change occurs deterministically
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 informationMath 471. Numerical methods Root-finding algorithms for nonlinear equations
Math 471. Numerical methods Root-finding algorithms for nonlinear equations overlap Section.1.5 of Bradie Our goal in this chapter is to find the root(s) for f(x) = 0..1 Bisection Method Intermediate value
More informationNumerical Methods - Preliminaries
Numerical Methods - Preliminaries Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Preliminaries 2013 1 / 58 Table of Contents 1 Introduction to Numerical Methods Numerical
More informationCHAPTER 10 Zeros of Functions
CHAPTER 10 Zeros of Functions An important part of the maths syllabus in secondary school is equation solving. This is important for the simple reason that equations are important a wide range of problems
More informationRegister machines L2 18
Register machines L2 18 Algorithms, informally L2 19 No precise definition of algorithm at the time Hilbert posed the Entscheidungsproblem, just examples. Common features of the examples: finite description
More informationphysicsandmathstutor.com Paper Reference Core Mathematics C3 Advanced Thursday 15 January 2009 Morning Time: 1 hour 30 minutes
Centre No. Candidate No. Paper Reference(s) 6665/01 Edexcel GCE Core Mathematics C3 Advanced Thursday 15 January 2009 Morning Time: 1 hour 30 minutes Materials required for examination Mathematical Formulae
More informationZeros of Functions. Chapter 10
Chapter 10 Zeros of Functions An important part of the mathematics syllabus in secondary school is equation solving. This is important for the simple reason that equations are important a wide range of
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 informationChapter 1A -- Real Numbers. iff. Math Symbols: Sets of Numbers
Fry Texas A&M University! Fall 2016! Math 150 Notes! Section 1A! Page 1 Chapter 1A -- Real Numbers Math Symbols: iff or Example: Let A = {2, 4, 6, 8, 10, 12, 14, 16,...} and let B = {3, 6, 9, 12, 15, 18,
More informationMCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory
MCE/EEC 647/747: Robot Dynamics and Control Lecture 8: Basic Lyapunov Stability Theory Reading: SHV Appendix Mechanical Engineering Hanz Richter, PhD MCE503 p.1/17 Stability in the sense of Lyapunov A
More informationTHE TEACHER UNDERSTANDS THE REAL NUMBER SYSTEM AND ITS STRUCTURE, OPERATIONS, ALGORITHMS, AND REPRESENTATIONS
The real number SySTeM C O M P E T E N C Y 1 THE TEACHER UNDERSTANDS THE REAL NUMBER SYSTEM AND ITS STRUCTURE, OPERATIONS, ALGORITHMS, AND REPRESENTATIONS This competency section reviews some of the fundamental
More informationPLC Papers. Created For:
PLC Papers Created For: Algebra and proof 2 Grade 8 Objective: Use algebra to construct proofs Question 1 a) If n is a positive integer explain why the expression 2n + 1 is always an odd number. b) Use
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 informationAn Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations
An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical
More informationOutline. Scientific Computing: An Introductory Survey. Nonlinear Equations. Nonlinear Equations. Examples: Nonlinear Equations
Methods for Systems of Methods for Systems of Outline Scientific Computing: An Introductory Survey Chapter 5 1 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationAlgorithms and Flowcharts
1 CHAPTER ONE Algorithms and Flowcharts Introduction The program is a set of instruction gave to the computer to execute successive operations leads to solve specific problem. In general to solve any problem
More informationTHE TEACHER UNDERSTANDS THE REAL NUMBER SYSTEM AND ITS STRUCTURE, OPERATIONS, ALGORITHMS, AND REPRESENTATIONS
THE REAL NUMBER SYSTEM C O M P E T E N C Y 1 THE TEACHER UNDERSTANDS THE REAL NUMBER SYSTEM AND ITS STRUCTURE, OPERATIONS, ALGORITHMS, AND REPRESENTATIONS This competency section reviews some of the fundamental
More informationAccel Alg E. L. E. Notes Solving Quadratic Equations. Warm-up
Accel Alg E. L. E. Notes Solving Quadratic Equations Warm-up Solve for x. Factor. 1. 12x 36 = 0 2. x 2 8x Factor. Factor. 3. 2x 2 + 5x 7 4. x 2 121 Solving Quadratic Equations Methods: (1. By Inspection)
More informationAdvanced Mathematics Unit 2 Limits and Continuity
Advanced Mathematics 3208 Unit 2 Limits and Continuity NEED TO KNOW Expanding Expanding Expand the following: A) (a + b) 2 B) (a + b) 3 C) (a + b)4 Pascals Triangle: D) (x + 2) 4 E) (2x -3) 5 Random Factoring
More informationAdvanced Mathematics Unit 2 Limits and Continuity
Advanced Mathematics 3208 Unit 2 Limits and Continuity NEED TO KNOW Expanding Expanding Expand the following: A) (a + b) 2 B) (a + b) 3 C) (a + b)4 Pascals Triangle: D) (x + 2) 4 E) (2x -3) 5 Random Factoring
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 informationLecture 5. September 4, 2018 Math/CS 471: Introduction to Scientific Computing University of New Mexico
Lecture 5 September 4, 2018 Math/CS 471: Introduction to Scientific Computing University of New Mexico 1 Review: Office hours at regularly scheduled times this week Tuesday: 9:30am-11am Wed: 2:30pm-4:00pm
More information5 Finding roots of equations
Lecture notes for Numerical Analysis 5 Finding roots of equations Topics:. Problem statement. Bisection Method 3. Newton s Method 4. Fixed Point Iterations 5. Systems of equations 6. Notes and further
More informationNov : Lecture 22: Differential Operators, Harmonic Oscillators
14 Nov. 3 005: Lecture : Differential Operators, Harmonic Oscillators Reading: Kreyszig Sections:.4 (pp:81 83),.5 (pp:83 89),.8 (pp:101 03) Differential Operators The idea of a function as something that
More informationGetting to the Roots of Quadratics
NAME BACKGROUND Graphically: The real roots of a function are the x-coordinates of the points at which the graph of the function intercepts/crosses the x-axis. For a quadratic function, whose graph is
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 information