Root-finding and optimisation

Size: px
Start display at page:

Download "Root-finding and optimisation"

Transcription

1 Root-finding and optimisation Reference tets: Chapra, S., R. Canale. Numerical Methods for Engineers, 3rd ed., McGraw-Hill, New York, Parts 2 and 4. Root-finding Simple mathematical equations have eplicit solutions, meaning that the variable we need to calculate can be isolated on the left hand side of the equation. Many problems we face as engineers involve solution of implicit equations. For eample, in the following polynomial, if is known then y can be solved eplicitly, whereas if y is known, solving for is an implicit problem. y = There are various ways of solving an implicit equation. The most common is to consider it as a root-finding problem. The roots of a function are where that function equals zero. For eample, using the polynomial above, solve for y = 5. This is equivalent to solving f() = 0, where f() is defined as, f ( ) = If you plot f() against you will see there are 2 roots at approimately = -0.3 and = 4.7. There are a range of root-finding methods available. Here we introduce two methods using the same eample f()

2 Newton-Raphson This requires an initial estimate of the root (for our eample, we ll choose 1 = 0.75), and estimates the gradient of the function f() at this value of. The gradient is then projected to the intersection with the -ais, which gives the second estimate of the root, 2. A third estimate is achieved in the same way, etc, until the estimate converges to the true root, in our case, f() The advantage of this method is that it usually converges to an accurate estimate very quickly. However, there are a number of eamples where it takes a long time to converge or finds the wrong root (see Chapra, pages ). For eample, what would have happened if we had used 1 = 1? Also, in general, our problems are more comple than this polynomial and we do not want to have to plot out f() every time. Therefore we generally do not know whether more than one root eists, and we might not get the answer we are looking for. Also, some functions we use are discontinuous i.e. there are some values of where gradients do not eist, in which case Newton-Raphson should not be used. Bisection method First, upper and lower bounds of are chosen (u 1 and l 1 ), so that we know the root lies within these bounds (in practical problems this is usually easy, because there are usually physical constraints on the value of ). f() is calculated at these two values of and also at their mid-point m 1. In our eample below, we have chosen l 1 = -0.85, u 1 = 0 and m 1 = The three values of f() are used to identify the sub-interval within which f() changes from positive to negative. This sub-interval is used to define the

3 second estimates of l and u. In our eample, l 2 = and u 2 = 0, so m 2 becomes This process continues until the interval converges on the root. f() -1 l 1 m 1 u 1 =0 m 2 3rd interval 2nd interval initial interval The advantage of this method is that it avoids the need to calculate gradients, and will always find a root which is between the initial estimates l 1 and u 1 (whereas the Newton-Raphson method may find a root which is far away from 1 ). However, in most problems this bisection method is less efficient that Newton-Raphson. Chapra s book describes improved versions of both these methods of root-finding, and gives instructions for developing computer code. Root-finding using Ecel: Manning s equation eample For steady flow in a rectangular open channel, Manning's equation can be used to determine the flow velocity as a function of channel characteristics during uniform, steady flow: 1 2 / 3 1/ 2 U = R S (1) n where U = mean velocity [m/s]; n = Manning roughness coefficient ; R = hydraulic radius [m] = A c /P; A c = cross-sectional area [m 2 ]; B = width [m]; H = mean depth [m]; P =

4 wetted perimeter [m]; and S = slope. Note that A c and P are related to the more fundamental parameters B and H by A c = BH (2) P = B + 2H (3) P A c B S H Q, U The flow is equal U multiplied by A c, Q = UA c (4) where Q = flow rate [m 3 /s]. Given the channel characteristics (H, B, S and n), these two equations can be solved directly for flow and velocity. A more interesting problem involves determine the depth and velocity if Q, B, S and n are given. This can be illustrated by eliminating U by substituting Eq. (1) into Eq. (4) Q = S n 12 / 53 / ( BH) ( B+ 2H) 23 / (5) Although we have one equation with one unknown, it is impossible to solve eplicitly for H. There is no way to algebraically manipulate Eq. (5) so that H is isolated on one side of the equation. However, the depth can be determined by root-finding. This is done by subtracting Q from each side to yield, f ( H) = 12 / 53 / S n ( BH) ( B+ 2H) 23 / Q = 0 (6) The value of H that makes this equation equal to zero is the root. For eample, if Q = 5 m 3 /s, B = 20 m, n = 0.03 and S = , the equation is

5 53 / ( 20H) f ( H) = = 0 (7) 23 / ( H) If you didn t know how to find roots using numerical methods, you could obtain a solution by trial-and-error. That is, you could guess values of H until you determined a value that drove the function to zero. Although this is an inefficient approach, eventually you would determine H = m. The result can be checked by substitution into Eq. (7) to give, ( ) f ( H ) = ( ) 5 / 3 2 / 3 5 = which is very close to zero. Our other unknown, the velocity can be determined by substitution back into Eq. (4), Q U = BH = 5 20(0.702) = 0.36 m/s An easier way to find the root is to use Ecel s Goal-Seek (from the Tools menu). Goalseek is based on the Newton-Raphson method. Try this yourself. Set cell refer to the cell in which you have calculated the value of f(h). To value enter the target value of f(h) = 0. By changing cell refer to the value of H (which will also be refered to within your f(h) formula). Eercise: Set up the polynomial equation in Ecel, and find both roots using Goal-seek.

6 Optimisation Optimisation is one of the most important challenges in hydrology and environmental engineering. The term optimisation algorithm covers a wide range of numerical methods of finding an optimal (i.e. maimised or minimised) solution. A optimisation algorithm consists of at least 2 components: 1. a series of formulae that produce an output given inputs such as data and parameters; 2. a measure of how good the output is (called one of the following: objective function; criteria; goodness of fit; cost function). The aim is usually to identify the input values which produce the optimum (maimum or minimum) objective function value. There are many different methods of optimisation that range from simple automated trial and error to very comple algorithms. In this workshop we introduce 2 simple types and use the optimisation tool provided in Ecel. Before we start, some of you will realise that you have already done some optimisation - in the 2 nd Ecel workbook. There, you used regression to optimise the parameters of a pipe orifice equation. In that case, the objective function which you minimised was the errors between the observed and calculated orifice flow. Regression, however, is limited to fitting linear models to data we need to look at more general methods. Newton s method This is almost the same as the Newton-Raphson method which we used for root-finding. When a continuous function reaches a maimum or minimum, its first derivative becomes zero. For a univariate function (i.e. a function with only one variable input,, and objective function y), y = f () dy = f '( ) = 0 d So all we have to do in order to optimise y is locate the root of f (). To maimise y from the previous polynomial eample, y = dy d 3 2 = f '( ) =

7 f'() The major problem in this eample is that there are three roots and two maima! This is a general problem with the Newton method, and so it should only be used if you know your starting point 1 is close to the optimum. Another disadvantage is the need to calculate derivatives, so the function needs to be continuous in the vicinity of the root. Ecel includes an optimisation tool called Solver 1. Solver uses a comple version of Newton s method, which combines a quasi-newton algorithm with a conjugate gradient search (see Chapra, Chapters 14 and 15). Solver does not have problems with distinguishing between maima and minima, and allows multivariate optimisation (i.e. of functions with more than one variable input). However, it still has problems distinguishing between two maima, as in our polynomial eample. Make a copy of your spreadsheet which includes the polynomial equation. Use Solver to maimise the value of y. Go to Tools, Solver: 1 You might have to install Solver. Go to Tools on the Ecel menubar, and choose Add-ins, and tick Solver add-in, then OK.

8 Your target cell is the one you want to be maimised (i.e. the one with the polynomial formula in it), and the cell you want to change will have the value of in it. There are no constraints in this problem. Note that if you use a starting point of less than about 1.85, then Solver will not find the true maimum. (Note - Solver can also be used to find roots by specifying a value ( Value of: ) rather than a maimum or minimum). Optimisation using random sampling The first step in this method is to define upper and lower limits of all the input variables, between which you know the optimum values lie. The second step is to generate a large number of random samples from within these limits and run the formulue using each sample. Finally, locate the optimum value of the objective function and corresponding input values (e.g. by sorting). The sampled optimum will be an approimation of the true optimum. A smaller number of samples will tend to produce a poorer approimation, and so will a wide difference between the upper and lower limits. Therefore, after applying the method initially, it is common to then go back and refine the estimates of the upper and lower limits and repeat. This method is very robust, because it does not rely on gradients and has no problems when there are multiple maima and minima it will (approimately) find them all. In the polynomial eample, lets assume upper and lower limits, u = 5 and l = 0, and take 8 random samples of : f() l u I have been lucky and one of my samples gave a good approimation of the maimum y (note - in practical problems we would not usually know the true maimum, so we would not know how good our approimation was). However, it can be seen that many more samples would be required to be sure of a good approimation of the maimum y. In

9 comple practical problems, hundreds of thousands of samples and function evaluations might be needed, in which case this method may be too computationally inefficient. Constrained optimisation: eample using Ecel s Solver In this eample, you will optimise a simplified water abstraction and wastewater treatment system to minimise the associated costs. The system is illustrated below, Drinking water abstraction = A Sewage treatment plant: pollution load = W flow = Q flow = Q-A Distance downstream of STP = The input variables which you will optimise are: A percentage reduction in the abstraction rate A A percentage reduction in the pollution load W The associated costs are: The cost of other sources of drinking water to replace the reduction in A The cost of installing a new STP to reduce W The cost due to loss of amenity because of pollution in the river Cost functions Assume that the cost of replacing the abstraction is, 2 XA = 2000R where R is the percentage reduction in A. Assume that the cost of the new STP is, XW = S 0.75

10 where S is the percentage reduction in W. Assume that the cost associated with river pollution is, XC = 1000D where D is the length of river in which the dissolved oygen concentration is less than 5g/m 3. Mathematical model We need a mathematical model in order to estimate D. The model to use is called the Streeter-Phelps model of river dissolved oygen, = c k L k d 0 c s ka d e kd v e ka v where c (g/m 3 ) is the concentration of dissolved oygen at distance downstream of the STP, c s (g/m 3 ) is the natural concentration of dissolved oygen, k d is the pollution degradation rate (s -1 ), k a is the natural rearation rate (s -1 ), L 0 is the concentration of pollution in the river at the outfall of the STP (i.e. at = 0), and v (m/s) is the river velocity. Using a simple mass balance equation where the STP outfall meets the river, L 0 can be calculated as, W ' L 0 = Q' where W is the pollution load after the reduction, and Q is the flow after the abstraction, W ' = W (1 S /100) Q' = Q A ( 1 R /100) Assume v is the following function of flow, v = 0.2 Q [ '] Use the following inputs and parameter values,

11 c s = 10 g/m 3 kd = 2 day -1 [= s -1 ] ka = 5 day -1 [= s -1 ] W = kg/day Q = 5 m 3 /s A = 2 m 3 /s Ecel spreadsheet First you should set up a spreadsheet which contains entries for all the parameters and inputs. Also include initial estimates for S and R, and the equations for XA, XW, W, Q, L 0 and v. For eample: You then need to calculate D. To do so, calculate c for a column of values, with ranging from 0 up to, say, 50000m at intervals of 500m. You should find that c(=0) = 10g/m 3 and then c decreases to a minimum then increases again back up to 10g/m 3. Using Ecel s IF function, calculate D and hence XC. Hence calculate the total cost = XC + XA + XW. You can now try to use Solver to minimise this cost by changing S and R. You should use the following constraints in Solver: 0 R S Min ( c ) 0 [My answer for the optimum cost was just under 5m, with S = 52% and R = 0]

Appendix F. + 1 Ma 1. 2 Ma Ma Ma ln + K = 0 (4-173)

Appendix F. + 1 Ma 1. 2 Ma Ma Ma ln + K = 0 (4-173) 5:39p.m. Page:949 Trimsize:8.5in 11in Appendix F F.1 MICROSOFT EXCEL SOLVER FOR NON-LINEAR EQUATIONS The Solver is an optimization package that finds a maximum, minimum, or specified value of a target

More information

2.29 Numerical Fluid Mechanics Spring 2015 Lecture 4

2.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 information

Numerical Methods. Root Finding

Numerical 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 information

Optimisation : Maxima & Minima

Optimisation : Maxima & Minima Optimisation : Maima & Minima Problem Solving Differentiation can be used to solve problems which require maimum or minimum values. Problems typically cover topics such as areas, volumes and rates of change.

More information

In order to master the techniques explained here it is vital that you undertake plenty of practice exercises so that they become second nature.

In order to master the techniques explained here it is vital that you undertake plenty of practice exercises so that they become second nature. Maima and minima In this unit we show how differentiation can be used to find the maimum and minimum values of a function. Because the derivative provides information about the gradient or slope of the

More information

Review: Properties of Exponents (Allow students to come up with these on their own.) m n m n. a a a. n n n m. a a a. a b a

Review: Properties of Exponents (Allow students to come up with these on their own.) m n m n. a a a. n n n m. a a a. a b a Algebra II Notes Unit Si: Polynomials Syllabus Objectives: 6. The student will simplify polynomial epressions. Review: Properties of Eponents (Allow students to come up with these on their own.) Let a

More information

Fundamentals of Algebra, Geometry, and Trigonometry. (Self-Study Course)

Fundamentals of Algebra, Geometry, and Trigonometry. (Self-Study Course) Fundamentals of Algebra, Geometry, and Trigonometry (Self-Study Course) This training is offered eclusively through the Pennsylvania Department of Transportation, Business Leadership Office, Technical

More information

Math 2414 Activity 1 (Due by end of class July 23) Precalculus Problems: 3,0 and are tangent to the parabola axis. Find the other line.

Math 2414 Activity 1 (Due by end of class July 23) Precalculus Problems: 3,0 and are tangent to the parabola axis. Find the other line. Math 44 Activity (Due by end of class July 3) Precalculus Problems: 3, and are tangent to the parabola ais. Find the other line.. One of the two lines that pass through y is the - {Hint: For a line through

More information

Basic methods to solve equations

Basic methods to solve equations Roberto s Notes on Prerequisites for Calculus Chapter 1: Algebra Section 1 Basic methods to solve equations What you need to know already: How to factor an algebraic epression. What you can learn here:

More information

Math 140 Final Sample A Solutions. Tyrone Crisp

Math 140 Final Sample A Solutions. Tyrone Crisp Math 4 Final Sample A Solutions Tyrone Crisp (B) Direct substitution gives, so the limit is infinite. When is close to, but greater than,, the numerator is negative while the denominator is positive. So

More information

Chapter 6 Overview: Applications of Derivatives

Chapter 6 Overview: Applications of Derivatives Chapter 6 Overview: Applications of Derivatives There are two main contets for derivatives: graphing and motion. In this chapter, we will consider the graphical applications of the derivative. Much of

More information

Lesson #33 Solving Incomplete Quadratics

Lesson #33 Solving Incomplete Quadratics Lesson # Solving Incomplete Quadratics A.A.4 Know and apply the technique of completing the square ~ 1 ~ We can also set up any quadratic to solve it in this way by completing the square, the technique

More information

M151B Practice Problems for Exam 1

M151B Practice Problems for Exam 1 M151B Practice Problems for Eam 1 Calculators will not be allowed on the eam. Unjustified answers will not receive credit. 1. Compute each of the following its: 1a. 1b. 1c. 1d. 1e. 1 3 4. 3. sin 7 0. +

More information

In #1-5, find the indicated limits. For each one, if it does not exist, tell why not. Show all necessary work.

In #1-5, find the indicated limits. For each one, if it does not exist, tell why not. Show all necessary work. Calculus I Eam File Fall 7 Test # In #-5, find the indicated limits. For each one, if it does not eist, tell why not. Show all necessary work. lim sin.) lim.) 3.) lim 3 3-5 4 cos 4.) lim 5.) lim sin 6.)

More information

Methods for Advanced Mathematics (C3) Coursework Numerical Methods

Methods for Advanced Mathematics (C3) Coursework Numerical Methods Woodhouse College 0 Page Introduction... 3 Terminolog... 3 Activit... 4 Wh use numerical methods?... Change of sign... Activit... 6 Interval Bisection... 7 Decimal Search... 8 Coursework Requirements on

More information

A. Incorrect! Apply the rational root test to determine if any rational roots exist.

A. Incorrect! Apply the rational root test to determine if any rational roots exist. College Algebra - Problem Drill 13: Zeros of Polynomial Functions No. 1 of 10 1. Determine which statement is true given f() = 3 + 4. A. f() is irreducible. B. f() has no real roots. C. There is a root

More information

Definition: Quadratic equation: A quadratic equation is an equation that could be written in the form ax 2 + bx + c = 0 where a is not zero.

Definition: Quadratic equation: A quadratic equation is an equation that could be written in the form ax 2 + bx + c = 0 where a is not zero. We will see many ways to solve these familiar equations. College algebra Class notes Solving Quadratic Equations: Factoring, Square Root Method, Completing the Square, and the Quadratic Formula (section

More information

Algebra Concepts Equation Solving Flow Chart Page 1 of 6. How Do I Solve This Equation?

Algebra Concepts Equation Solving Flow Chart Page 1 of 6. How Do I Solve This Equation? Algebra Concepts Equation Solving Flow Chart Page of 6 How Do I Solve This Equation? First, simplify both sides of the equation as much as possible by: combining like terms, removing parentheses using

More information

Math Methods 12. Portfolio Assignment 4 Type II. 1. The general formula for the binomial expansion of ( x + y) n is given by 2! 3!

Math Methods 12. Portfolio Assignment 4 Type II. 1. The general formula for the binomial expansion of ( x + y) n is given by 2! 3! Math Methods Portfolio Assignment 4 Type II Name: Rajesh Swaminathan Block: D Date: 4-Mar-5 ANALYSIS OF A QUARTIC FUNCTION. The general formula for the binomial epansion of ( + y) n is given by nn ( )

More information

STATIC LECTURE 4: CONSTRAINED OPTIMIZATION II - KUHN TUCKER THEORY

STATIC LECTURE 4: CONSTRAINED OPTIMIZATION II - KUHN TUCKER THEORY STATIC LECTURE 4: CONSTRAINED OPTIMIZATION II - KUHN TUCKER THEORY UNIVERSITY OF MARYLAND: ECON 600 1. Some Eamples 1 A general problem that arises countless times in economics takes the form: (Verbally):

More information

Linear Equations in Linear Algebra

Linear Equations in Linear Algebra 1 Linear Equations in Linear Algebra 1.1 SYSTEMS OF LINEAR EQUATIONS LINEAR EQUATION,, 1 n A linear equation in the variables equation that can be written in the form a a a b 1 1 2 2 n n a a is an where

More information

MEI Core 2. Sequences and series. Section 1: Definitions and Notation

MEI Core 2. Sequences and series. Section 1: Definitions and Notation Notes and Eamples MEI Core Sequences and series Section : Definitions and Notation In this section you will learn definitions and notation involving sequences and series, and some different ways in which

More information

Math 2414 Activity 1 (Due by end of class Jan. 26) Precalculus Problems: 3,0 and are tangent to the parabola axis. Find the other line.

Math 2414 Activity 1 (Due by end of class Jan. 26) Precalculus Problems: 3,0 and are tangent to the parabola axis. Find the other line. Math Activity (Due by end of class Jan. 6) Precalculus Problems: 3, and are tangent to the parabola ais. Find the other line.. One of the two lines that pass through y is the - {Hint: For a line through

More information

Chapter 9 Regression. 9.1 Simple linear regression Linear models Least squares Predictions and residuals.

Chapter 9 Regression. 9.1 Simple linear regression Linear models Least squares Predictions and residuals. 9.1 Simple linear regression 9.1.1 Linear models Response and eplanatory variables Chapter 9 Regression With bivariate data, it is often useful to predict the value of one variable (the response variable,

More information

Engineering Hydrology (ECIV 4323) CHAPTER FOUR. Stream flow measurement. Instructors: Dr. Yunes Mogheir Dr. Ramadan Al Khatib

Engineering Hydrology (ECIV 4323) CHAPTER FOUR. Stream flow measurement. Instructors: Dr. Yunes Mogheir Dr. Ramadan Al Khatib Engineering Hydrology (ECIV 4323) CHAPTER FOUR Stream flow measurement Instructors: Dr. Yunes Mogheir Dr. Ramadan Al Khatib -١ 4.1 Introduction - Surface water hydrology deals with the movement of water

More information

Math 180, Final Exam, Spring 2008 Problem 1 Solution. 1. For each of the following limits, determine whether the limit exists and, if so, evaluate it.

Math 180, Final Exam, Spring 2008 Problem 1 Solution. 1. For each of the following limits, determine whether the limit exists and, if so, evaluate it. Math 80, Final Eam, Spring 008 Problem Solution. For each of the following limits, determine whether the limit eists and, if so, evaluate it. + (a) lim 0 (b) lim ( ) 3 (c) lim Solution: (a) Upon substituting

More information

Lesson Goals. Unit 2 Functions Analyzing Graphs of Functions (Unit 2.2) Graph of a Function. Lesson Goals

Lesson Goals. Unit 2 Functions Analyzing Graphs of Functions (Unit 2.2) Graph of a Function. Lesson Goals Unit Functions Analzing Graphs of Functions (Unit.) William (Bill) Finch Mathematics Department Denton High School Lesson Goals When ou have completed this lesson ou will: Find the domain and range of

More information

APPLICATIONS OF DIFFERENTIATION

APPLICATIONS OF DIFFERENTIATION 4 APPLICATIONS OF DIFFERENTIATION APPLICATIONS OF DIFFERENTIATION 4.8 Newton s Method In this section, we will learn: How to solve high degree equations using Newton s method. INTRODUCTION Suppose that

More information

Finding Roots of Equations

Finding Roots of Equations Finding Roots of Equations Solution Methods Overview Bisection/Half-interval Search Method of false position/regula Falsi Secant Method Newton Raphson Iteration Method Many more. Open Methods Bracketing

More information

Scientific Computing. Roots of Equations

Scientific 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 information

Presented by: Civil Engineering Academy

Presented by: Civil Engineering Academy Presented by: Civil Engineering Academy Open-Channel Flow Uniform Flow (See CERM Ch. 19) Characterized by constant depth volume, and cross section. It can be steady or unsteady Non-uniform Flow *Not on

More information

dx dx x sec tan d 1 4 tan 2 2 csc d 2 ln 2 x 2 5x 6 C 2 ln 2 ln x ln x 3 x 2 C Now, suppose you had observed that x 3

dx dx x sec tan d 1 4 tan 2 2 csc d 2 ln 2 x 2 5x 6 C 2 ln 2 ln x ln x 3 x 2 C Now, suppose you had observed that x 3 CHAPTER 8 Integration Techniques, L Hôpital s Rule, and Improper Integrals Section 8. Partial Fractions Understand the concept of a partial fraction decomposition. Use partial fraction decomposition with

More information

Further algebra. polynomial identities

Further algebra. polynomial identities 8 8A Polynomial identities 8B Partial fractions 8C Simultaneous equations areas of study The solution of simultaneous equations arising from the intersection of a line with a parabola, circle or rectangular

More information

Exact and Approximate Numbers:

Exact and Approximate Numbers: Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.

More information

171, Calculus 1. Summer 1, CRN 50248, Section 001. Time: MTWR, 6:30 p.m. 8:30 p.m. Room: BR-43. CRN 50248, Section 002

171, Calculus 1. Summer 1, CRN 50248, Section 001. Time: MTWR, 6:30 p.m. 8:30 p.m. Room: BR-43. CRN 50248, Section 002 171, Calculus 1 Summer 1, 018 CRN 5048, Section 001 Time: MTWR, 6:0 p.m. 8:0 p.m. Room: BR-4 CRN 5048, Section 00 Time: MTWR, 11:0 a.m. 1:0 p.m. Room: BR-4 CONTENTS Syllabus Reviews for tests 1 Review

More information

Today 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. 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 information

12.10 Lagrange Multipliers

12.10 Lagrange Multipliers .0 Lagrange Multipliers In the last two sections we were often solving problems involving maimizing or minimizing a function f subject to a 'constraint' equation g. For eample, we minimized the cost of

More information

3.6 Determinants. 3.6 Determinants 1

3.6 Determinants. 3.6 Determinants 1 3.6 Determinants 1 3.6 Determinants We said in Section 3.3 that a 2 2 matri a b c d is invertible if and only if its erminant, ad - bc, is nonzero, and we saw the erminant used in the formula for the inverse

More information

Equations and Inequalities

Equations and Inequalities Equations and Inequalities Figure 1 CHAPTER OUTLINE.1 The Rectangular Coordinate Systems and Graphs. Linear Equations in One Variable.3 Models and Applications. Comple Numbers.5 Quadratic Equations.6 Other

More information

Chapter III. Unconstrained Univariate Optimization

Chapter III. Unconstrained Univariate Optimization 1 Chapter III Unconstrained Univariate Optimization Introduction Interval Elimination Methods Polynomial Approximation Methods Newton s Method Quasi-Newton Methods 1 INTRODUCTION 2 1 Introduction Univariate

More information

CHAPTER 1-2: SHADOW PRICES

CHAPTER 1-2: SHADOW PRICES Essential Microeconomics -- CHAPTER -: SHADOW PRICES An intuitive approach: profit maimizing firm with a fied supply of an input Shadow prices 5 Concave maimization problem 7 Constraint qualifications

More information

Polynomial Functions of Higher Degree

Polynomial Functions of Higher Degree SAMPLE CHAPTER. NOT FOR DISTRIBUTION. 4 Polynomial Functions of Higher Degree Polynomial functions of degree greater than 2 can be used to model data such as the annual temperature fluctuations in Daytona

More information

Core Connections Algebra 2 Checkpoint Materials

Core Connections Algebra 2 Checkpoint Materials Core Connections Algebra 2 Note to Students (and their Teachers) Students master different skills at different speeds. No two students learn eactly the same way at the same time. At some point you will

More information

Intro to Nonlinear Optimization

Intro to Nonlinear Optimization Intro to Nonlinear Optimization We now rela the proportionality and additivity assumptions of LP What are the challenges of nonlinear programs NLP s? Objectives and constraints can use any function: ma

More information

ACCUPLACER MATH 0311 OR MATH 0120

ACCUPLACER MATH 0311 OR MATH 0120 The University of Teas at El Paso Tutoring and Learning Center ACCUPLACER MATH 0 OR MATH 00 http://www.academics.utep.edu/tlc MATH 0 OR MATH 00 Page Factoring Factoring Eercises 8 Factoring Answer to Eercises

More information

Calculus of Variation An Introduction To Isoperimetric Problems

Calculus of Variation An Introduction To Isoperimetric Problems Calculus of Variation An Introduction To Isoperimetric Problems Kevin Wang The University of Sydney SSP Working Seminars, MATH2916 May 4, 2013 Contents I Lagrange Multipliers 2 1 Single Constraint Lagrange

More information

Math 1314 Lesson 4 Limits

Math 1314 Lesson 4 Limits Math 1314 Lesson 4 Limits What is calculus? Calculus is the study of change, particularly, how things change over time. It gives us a framework for measuring change using some fairly simple models. In

More information

Closed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis.

Closed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis. OPEN CHANNEL FLOW Open channel flow is a flow of liquid, basically water in a conduit with a free surface. The open channel flows are driven by gravity alone, and the pressure gradient at the atmospheric

More information

SAMPLE. Applications of differentiation

SAMPLE. Applications of differentiation Objectives C H A P T E R 10 Applications of differentiation To be able to find the equation of the tangent and the normal at a given point on a curve. To be able to calculate the angles between straight

More information

Equations and Inequalities

Equations and Inequalities Equations and Inequalities Figure 1 CHAPTER OUTLINE 1 The Rectangular Coordinate Systems and Graphs Linear Equations in One Variable Models and Applications Comple Numbers Quadratic Equations 6 Other Types

More information

Comparison of Average Energy Slope Estimation Formulas for One-dimensional Steady Gradually Varied Flow

Comparison of Average Energy Slope Estimation Formulas for One-dimensional Steady Gradually Varied Flow Archives of Hydro-Engineering and Environmental Mechanics Vol. 61 (2014), No. 3 4, pp. 89 109 DOI: 10.1515/heem-2015-0006 IBW PAN, ISSN 1231 3726 Comparison of Average Energy Slope Estimation Formulas

More information

PACKET Unit 4 Honors ICM Functions and Limits 1

PACKET Unit 4 Honors ICM Functions and Limits 1 PACKET Unit 4 Honors ICM Functions and Limits 1 Day 1 Homework For each of the rational functions find: a. domain b. -intercept(s) c. y-intercept Graph #8 and #10 with at least 5 EXACT points. 1. f 6.

More information

Iteration & Fixed Point

Iteration & Fixed Point Iteration & Fied Point As a method for finding the root of f this method is difficult, but it illustrates some important features of iterstion. We could write f as f g and solve g. Definition.1 (Fied Point)

More information

9.1. Solving Quadratic Equations. Investigation: Rocket Science CONDENSED LESSON

9.1. Solving Quadratic Equations. Investigation: Rocket Science CONDENSED LESSON CONDENSED LESSON 9.1 Solving Quadratic Equations In this lesson you will look at quadratic functions that model projectile motion use tables and graphs to approimate solutions to quadratic equations solve

More information

Prep for College Algebra

Prep for College Algebra Prep for College Algebra This course covers the topics outlined below. You can customize the scope and sequence of this course to meet your curricular needs. Curriculum (219 topics + 85 additional topics)

More information

Open Channel Flow I - The Manning Equation and Uniform Flow COURSE CONTENT

Open Channel Flow I - The Manning Equation and Uniform Flow COURSE CONTENT Open Channel Flow I - The Manning Equation and Uniform Flow Harlan H. Bengtson, PhD, P.E. COURSE CONTENT 1. Introduction Flow of a liquid may take place either as open channel flow or pressure flow. Pressure

More information

Find the following limits. For each one, if it does not exist, tell why not. Show all necessary work.

Find the following limits. For each one, if it does not exist, tell why not. Show all necessary work. Calculus I Eam File Spring 008 Test #1 Find the following its. For each one, if it does not eist, tell why not. Show all necessary work. 1.) 4.) + 4 0 1.) 0 tan 5.) 1 1 1 1 cos 0 sin 3.) 4 16 3 1 6.) For

More information

Prep for College Algebra with Trigonometry

Prep for College Algebra with Trigonometry Prep for College Algebra with Trigonometry This course covers the topics outlined below. You can customize the scope and sequence of this course to meet your curricular needs. Curriculum (246 topics +

More information

Unconstrained Multivariate Optimization

Unconstrained Multivariate Optimization Unconstrained Multivariate Optimization Multivariate optimization means optimization of a scalar function of a several variables: and has the general form: y = () min ( ) where () is a nonlinear scalar-valued

More information

Unit 11 - Solving Quadratic Functions PART ONE

Unit 11 - Solving Quadratic Functions PART ONE Unit 11 - Solving Quadratic Functions PART ONE PREREQUISITE SKILLS: students should be able to add, subtract and multiply polynomials students should be able to factor polynomials students should be able

More information

Part Two. Diagnostic Test

Part Two. Diagnostic Test Part Two Diagnostic Test AP Calculus AB and BC Diagnostic Tests Take a moment to gauge your readiness for the AP Calculus eam by taking either the AB diagnostic test or the BC diagnostic test, depending

More information

Summer Packet Honors PreCalculus

Summer Packet Honors PreCalculus Summer Packet Honors PreCalculus Honors Pre-Calculus is a demanding course that relies heavily upon a student s algebra, geometry, and trigonometry skills. You are epected to know these topics before entering

More information

CHAPTER 8 Quadratic Equations, Functions, and Inequalities

CHAPTER 8 Quadratic Equations, Functions, and Inequalities CHAPTER Quadratic Equations, Functions, and Inequalities Section. Solving Quadratic Equations: Factoring and Special Forms..................... 7 Section. Completing the Square................... 9 Section.

More information

Lecture 10: River Channels

Lecture 10: River Channels GEOG415 Lecture 10: River Channels 10-1 Importance of channel characteristics Prediction of flow was the sole purpose of hydrology, and still is a very important aspect of hydrology. - Water balance gives

More information

Chapter 4. Solution of a Single Nonlinear Algebraic Equation

Chapter 4. Solution of a Single Nonlinear Algebraic Equation Single Nonlinear Algebraic Equation - 56 Chapter 4. Solution of a Single Nonlinear Algebraic Equation 4.1. Introduction Life, my fris, is nonlinear. As such, in our roles as problem-solvers, we will be

More information

Fluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur

Fluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Fluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture - 42 Flows with a Free Surface Part II Good morning. I welcome you to this session

More information

Pure Core 2. Revision Notes

Pure Core 2. Revision Notes Pure Core Revision Notes June 06 Pure Core Algebra... Polynomials... Factorising... Standard results... Long division... Remainder theorem... 4 Factor theorem... 5 Choosing a suitable factor... 6 Cubic

More information

y intercept Gradient Facts Lines that have the same gradient are PARALLEL

y intercept Gradient Facts Lines that have the same gradient are PARALLEL CORE Summar Notes Linear Graphs and Equations = m + c gradient = increase in increase in intercept Gradient Facts Lines that have the same gradient are PARALLEL If lines are PERPENDICULAR then m m = or

More information

Section 6.2 Long Division of Polynomials

Section 6.2 Long Division of Polynomials Section 6. Long Division of Polynomials INTRODUCTION In Section 6.1 we learned to simplify a rational epression by factoring. For eample, + 3 10 = ( + 5)( ) ( ) = ( + 5) 1 = + 5. However, if we try to

More information

Troy High School AP Calculus Summer Packet

Troy High School AP Calculus Summer Packet Troy High School AP Calculus Summer Packet As instructors of AP Calculus, we have etremely high epectations of students taking our courses. We epect a certain level of independence to be demonstrated by

More information

NUMERICAL METHODS FOR SOLVING EQUATIONS

NUMERICAL METHODS FOR SOLVING EQUATIONS Mathematics Revision Guides Numerical Methods for Solving Equations Page of M.K. HOME TUITION Mathematics Revision Guides Level: AS / A Level AQA : C3 Edecel: C3 OCR: C3 NUMERICAL METHODS FOR SOLVING EQUATIONS

More information

Chapter 27 AB Calculus Practice Test

Chapter 27 AB Calculus Practice Test Chapter 7 AB Calculus Practice Test The Eam AP Calculus AB Eam SECTION I: Multiple-Choice Questions DO NOT OPEN THIS BOOKLET UNTIL YOU ARE TOLD TO DO SO. At a Glance Total Time 1 hour and 45 minutes Number

More information

is the intuition: the derivative tells us the change in output y (from f(b)) in response to a change of input x at x = b.

is the intuition: the derivative tells us the change in output y (from f(b)) in response to a change of input x at x = b. Uses of differentials to estimate errors. Recall the derivative notation df d is the intuition: the derivative tells us the change in output y (from f(b)) in response to a change of input at = b. Eamples.

More information

Partial Fractions. dx dx x sec tan d 1 4 tan 2. 2 csc d. csc cot C. 2x 5. 2 ln. 2 x 2 5x 6 C. 2 ln. 2 ln x

Partial Fractions. dx dx x sec tan d 1 4 tan 2. 2 csc d. csc cot C. 2x 5. 2 ln. 2 x 2 5x 6 C. 2 ln. 2 ln x 460_080.qd //04 :08 PM Page CHAPTER 8 Integration Techniques, L Hôpital s Rule, and Improper Integrals Section 8. Partial Fractions Understand the concept of a partial fraction decomposition. Use partial

More information

Computer Aided Design of Thermal Systems (ME648)

Computer Aided Design of Thermal Systems (ME648) Computer Aided Design of Thermal Systems (ME648) PG/Open Elective Credits: 3-0-0-9 Updated Syallabus: Introduction. Basic Considerations in Design. Modelling of Thermal Systems. Numerical Modelling and

More information

f x, y x 2 y 2 2x 6y 14. Then

f x, y x 2 y 2 2x 6y 14. Then SECTION 11.7 MAXIMUM AND MINIMUM VALUES 645 absolute minimum FIGURE 1 local maimum local minimum absolute maimum Look at the hills and valles in the graph of f shown in Figure 1. There are two points a,

More information

Solutions to the Exercises of Chapter 8

Solutions to the Exercises of Chapter 8 8A Domains of Functions Solutions to the Eercises of Chapter 8 1 For 7 to make sense, we need 7 0or7 So the domain of f() is{ 7} For + 5 to make sense, +5 0 So the domain of g() is{ 5} For h() to make

More information

BC Calculus Diagnostic Test

BC Calculus Diagnostic Test BC Calculus Diagnostic Test The Eam AP Calculus BC Eam SECTION I: Multiple-Choice Questions DO NOT OPEN THIS BOOKLET UNTIL YOU ARE TOLD TO DO SO. At a Glance Total Time hour and 5 minutes Number of Questions

More information

MATH 1325 Business Calculus Guided Notes

MATH 1325 Business Calculus Guided Notes MATH 135 Business Calculus Guided Notes LSC North Harris By Isabella Fisher Section.1 Functions and Theirs Graphs A is a rule that assigns to each element in one and only one element in. Set A Set B Set

More information

Polynomials and Polynomial Functions

Polynomials and Polynomial Functions Unit 5: Polynomials and Polynomial Functions Evaluating Polynomial Functions Objectives: SWBAT identify polynomial functions SWBAT evaluate polynomial functions. SWBAT find the end behaviors of polynomial

More information

18 19 Find the extreme values of f on the region described by the inequality. 20. Consider the problem of maximizing the function

18 19 Find the extreme values of f on the region described by the inequality. 20. Consider the problem of maximizing the function 940 CHAPTER 14 PARTIAL DERIVATIVES 14.8 EXERCISES 1. Pictured are a contour map of f and a curve with equation t, y 8. Estimate the maimum and minimum values of f subject to the constraint that t, y 8.

More information

Learning Targets: Standard Form: Quadratic Function. Parabola. Vertex Max/Min. x-coordinate of vertex Axis of symmetry. y-intercept.

Learning Targets: Standard Form: Quadratic Function. Parabola. Vertex Max/Min. x-coordinate of vertex Axis of symmetry. y-intercept. Name: Hour: Algebra A Lesson:.1 Graphing Quadratic Functions Learning Targets: Term Picture/Formula In your own words: Quadratic Function Standard Form: Parabola Verte Ma/Min -coordinate of verte Ais of

More information

SNAP Centre Workshop. Solving Systems of Equations

SNAP Centre Workshop. Solving Systems of Equations SNAP Centre Workshop Solving Systems of Equations 35 Introduction When presented with an equation containing one variable, finding a solution is usually done using basic algebraic manipulation. Example

More information

I have not checked this review sheet for errors, hence there maybe errors in this document. thank you.

I have not checked this review sheet for errors, hence there maybe errors in this document. thank you. I have not checked this review sheet for errors, hence there maybe errors in this document. thank you. Class test II Review sections 3.7(differentials)-5.5(logarithmic differentiation) ecluding section

More information

Introduction. So, why did I even bother to write this?

Introduction. So, why did I even bother to write this? Introduction This review was originally written for my Calculus I class, but it should be accessible to anyone needing a review in some basic algebra and trig topics. The review contains the occasional

More information

Functions and Their Graphs

Functions and Their Graphs Functions and Their Graphs 015 College Board. All rights reserved. Unit Overview In this unit you will study polynomial and rational functions, their graphs, and their zeros. You will also learn several

More information

Example - Newton-Raphson Method

Example - Newton-Raphson Method Eample - Newton-Raphson Method We now consider the following eample: minimize f( 3 3 + -- 4 4 Since f ( 3 2 + 3 3 and f ( 6 + 9 2 we form the following iteration: + n 3 ( n 3 3( n 2 ------------------------------------

More information

Objectives. By the time the student is finished with this section of the workbook, he/she should be able

Objectives. By the time the student is finished with this section of the workbook, he/she should be able FUNCTIONS Quadratic Functions......8 Absolute Value Functions.....48 Translations o Functions..57 Radical Functions...61 Eponential Functions...7 Logarithmic Functions......8 Cubic Functions......91 Piece-Wise

More information

Unit 2. Quadratic Functions and Modeling. 24 Jordan School District

Unit 2. Quadratic Functions and Modeling. 24 Jordan School District Unit Quadratic Functions and Modeling 4 Unit Cluster (F.F.4, F.F.5, F.F.6) Unit Cluster (F.F.7, F.F.9) Interpret functions that arise in applications in terms of a contet Analyzing functions using different

More information

Review: A Cross Section of the Midterm. MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.

Review: A Cross Section of the Midterm. MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Review: A Cross Section of the Midterm Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Find the limit, if it eists. 4 + ) lim - - ) A) - B) -

More information

Definition 8.1 Two inequalities are equivalent if they have the same solution set. Add or Subtract the same value on both sides of the inequality.

Definition 8.1 Two inequalities are equivalent if they have the same solution set. Add or Subtract the same value on both sides of the inequality. 8 Inequalities Concepts: Equivalent Inequalities Linear and Nonlinear Inequalities Absolute Value Inequalities (Sections.6 and.) 8. Equivalent Inequalities Definition 8. Two inequalities are equivalent

More information

4.1 & 4.2 Student Notes Using the First and Second Derivatives. for all x in D, where D is the domain of f. The number f()

4.1 & 4.2 Student Notes Using the First and Second Derivatives. for all x in D, where D is the domain of f. The number f() 4.1 & 4. Student Notes Using the First and Second Derivatives Deinition A unction has an absolute maimum (or global maimum) at c i ( c) ( ) or all in D, where D is the domain o. The number () c is called

More information

4.3 How derivatives affect the shape of a graph. The first derivative test and the second derivative test.

4.3 How derivatives affect the shape of a graph. The first derivative test and the second derivative test. Chapter 4: Applications of Differentiation In this chapter we will cover: 41 Maimum and minimum values The critical points method for finding etrema 43 How derivatives affect the shape of a graph The first

More information

Solutions to Math 41 Final Exam December 9, 2013

Solutions to Math 41 Final Exam December 9, 2013 Solutions to Math 4 Final Eam December 9,. points In each part below, use the method of your choice, but show the steps in your computations. a Find f if: f = arctane csc 5 + log 5 points Using the Chain

More information

Closed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis.

Closed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis. OPEN CHANNEL FLOW Open channel flow is a flow of liquid, basically water in a conduit with a free surface. The open channel flows are driven by gravity alone, and the pressure gradient at the atmospheric

More information

Algebra 8.6 Simple Equations

Algebra 8.6 Simple Equations Algebra 8.6 Simple Equations 1. Introduction Let s talk about the truth: 2 = 2 This is a true statement What else can we say about 2 that is true? Eample 1 2 = 2 1+ 1= 2 2 1= 2 4 1 = 2 2 4 2 = 2 4 = 4

More information

It is convenient to introduce some notation for this type of problems. I will write this as. max u (x 1 ; x 2 ) subj. to. p 1 x 1 + p 2 x 2 m ;

It is convenient to introduce some notation for this type of problems. I will write this as. max u (x 1 ; x 2 ) subj. to. p 1 x 1 + p 2 x 2 m ; 4 Calculus Review 4.1 The Utility Maimization Problem As a motivating eample, consider the problem facing a consumer that needs to allocate a given budget over two commodities sold at (linear) prices p

More information

Identifying end behavior of the graph of a polynomial function

Identifying end behavior of the graph of a polynomial function 56 Polnomial Functions 3.1.1 Eercises For a link to all of the additional resources available for this section, click OSttS Chapter 3 materials. In Eercises 1-10, find the degree, the leading term, the

More information

Optimal scaling of the random walk Metropolis on elliptically symmetric unimodal targets

Optimal scaling of the random walk Metropolis on elliptically symmetric unimodal targets Bernoulli 15(3), 2009, 774 798 DOI: 10.3150/08-BEJ176 Optimal scaling of the random walk Metropolis on elliptically symmetric unimodal targets CHRIS SHERLOCK 1 and GARETH ROBERTS 2 1 Department of Mathematics

More information