AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E.

Size: px
Start display at page:

Download "AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E."

Transcription

1 AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 14 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering an Applie Science Gonzaga University SDC PUBLICATIONS Schroff Development Corporation

2 Chapter : Introuction to Numerical Methos 43 C H A P T E R INTRODUCTION TO NUMERICAL METHODS.1 THE USE OF NUMERICAL METHODS IN SCIENCE AND ENGINEERING Analysis of problems in engineering an the physical sciences typically involves four steps as follows. (1) Development of a suitable mathematical moel that realistically represents a given physical system. () Derivation of the system governing equations using physical laws such as Newton's laws of motion, conservation of energy, the laws governing electrical circuits etc. (3) Solution of the governing equations, an (4) Interpretation of the results. Because real worl problems are generally quite complex with the generation of close-form analytical solutions becoming impossible in many situations, there exists, most efinitely, a nee for the proper utilization of computer-base techniques in the solution of practical problems. The avancement of computer technology has mae the effective use of numerical methos an computer-base techniques very feasible, an thus, solutions can now be obtaine much faster than ever before an with much better than acceptable accuracy. However, there are avantages as well as isavantages associate with any numerical proceure that is resorte to, an these must be kept in min when using it.. COMPARISON OF NUMERICAL METHODS WITH ANALYTICAL METHODS While an analytical solution will be exact if it exists, a numerical metho, on the other han, will generally require iterations to generate a solution, which is only an approximation an which certainly cannot be consiere exact by any means. A isavantage associate with analytical solution techniques is that they are generally applicable only to very special cases of problems. Numerical solutions, on the contrary, will solve complex situations as well. to generate a solution can be looke upon as isavantages. While numerical techniques have several avantages incluing easy programming on a computer an the convenience with which they hanle complex problems, the initial estimate of the solution along with the many number of iterations that are sometimes require

3 44 AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD.3 SOURCES OF NUMERICAL ERRORS AND THEIR COMPUTATION It is inee possible for miscalculations to creep into a numerical solution because of various sources of error. These inclue inaccurate mathematical moeling, wrong programming, wrong input, rouning off of numbers an truncation of an infinite series. Roun-off error is the general name given to inaccuracies that affect the calculation scene when a finite number of igits are assigne to represent an actual number. In a long sequence of calculations, this roun-off error can accumulate, then propagate through the process of calculation an finally grow very rapily to a significant number. A truncation error results when an infinite series is approximate by a finite number of terms, an, typically, upper bouns are place on the size of this error. The true error is efine as the ifference between the compute value an the true value of a number. E True = X Comp X True (.1) while the relative true error is the error relative to the true value X Comp X True e r = X True (.) Expresse as a percentage, the relative true error is written as X Comp X True e r = X True 1 (%) (.3).4 TAYLOR SERIES EXPANSION The Taylor series is consiere as a basis of approximation in numerical analysis. If the value of a function of x is provie at " x ", then the Taylor series provies a means of evaluating the function at " x + h", where " x " is the starting value of the inepenent variable an " h " is the ifference between the starting value an the new value at which the function is to be approximate = fx fx + h h + fx x h + fx! x h fx 3! x (.4) This equation can be use for generating various orers of approximations as shown below. The orer of approximation is efine by the highest erivative inclue in the series. For example, If only terms up to the secon erivative are retaine in the series, the result is a secon orer approximation.

4 Chapter : Introuction to Numerical Methos 4 First orer approximation: = fx ( ) + h fx + h fx Secon orer approximation: x (.) = fx ( ) + h fx + h fx x x fx ( ) (.6) Thir orer approximation: = fx ( ) + h fx + h fx h h fx (! ) fx ( x 3! 3 ) +... x x It is to note that the significance of the higher orer terms in the Taylor series increases with the nonlinearity of the function involve as well as the ifference between the " starting x" value an the "x" value at which the function is to be approximate. Thus, the fewer the terms that are inclue in the series, the larger will be the error associate with the computation of the function value. If the function is linear, however, only terms up to the first erivative term nee to be inclue. Example.1 Using the Taylor series expansion for f(x) = -.1 x x 3 -. x -. x + 1. etermine the zeroth, first, secon, thir, fourth an fifth orer approximations of f(x + h ) where x = an h = 1,, 3, 4, an compare these with the exact solutions. (.7) h=1.: Put in the function an generate its erivatives: fx :=.1 x 4.17 x 3. x.x + 1. x := h:= 1. fprime( x) :=.6 x 3.1 x. x. <--Generate erivatives fprime( x) := 1.8 x 1. x. f3prime x 3.6 := x 1. f4prime( x) := 3.6 fprime( x) :=.

5 46 AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD term1 h := f( x) term := h fprime( x) term3 := fprime x h 3 term4 := 6 term h 4 := f4prime x 4 term6 h := 1 f3prime x ftaylor := term1 <---- one-term or zero-orer approximation ftaylor1 := term1 + term <---- first orer approximation with two terms fprime ( x ) ftaylor := term1 + term + term3 <---secon orer approximation with 3 terms ftaylor3 term1 term := + + term3 + term4 <---thir orer approximation with 4 terms ftaylor4 := term1 + term + term3 + term4 + term <---fourth orer approximation with terms fifth orer approximation with 6 terms ftaylor := term1 + term + term3 + term4 + term + term6 <----- x := x + h x = 1 ftaylor = 1. ftaylor1 = 1 ftaylor =.7 ftaylor3 =.8 ftaylor4 =.43 ftaylor =.43 <-- These are the zero- fifth orer approximations of the given function f(x) using the Taylor series. f1 =.43 <---EXACT ANSWER USING FUNCTION GIVEN. err := f( x) ftaylor err =.8 ftaylor1 err1 := f x err1 =.7 err := f( x) ftaylor err =.3 These are errors ( ifferences between exact <-- values an approximations ) for the above zero - fifth orer approximations. err3 := f( x) ftaylor3 err3 =.1 err4 := f( x) ftaylor4 err4 = err := f( x) ftaylor err = Similarly, by using h=, 3, 4,, the zeroth- fifth orer approximations for f(), f(3), f(4), f() an the associate errors can be etermine. These are given in Tables.1 an. Plots of the various Taylor series approximations of the given function an associate errors are generate below an are presente in Figs..1 an.

6 Chapter : Introuction to Numerical Methos 47 x fx := x :=,.1.. ftaylor x := <-- zeroth-orer approximation ftaylor1( x) := ftaylor( x) + ( x x) fprime( x) <---first-orer approximation ( x x) ftaylor( x) := ftaylor1( x) + fprime( x) <--secon-orer approximation ( x x) 3 ftaylor3( x) := ftaylor( x) + f3prime( x) <---thir-orer approximation 6 ftaylor4 x ftaylor x ( x x) 4 := ftaylor3( x) + f4prime( x) <---fourth-orer approximation 4 ftaylor4 x ( x x) := + fprime( x) <---fifth- orer approximation 1 Errors generate with the various approximations are as follows Zero orer approximation: err( x) := fx ftaylor( x) First orer approximation: Secon orer approximation: err1( x) := fx ftaylor1( x) err( x) := fx ftaylor( x) Thir orer approximation: Fourth orer approximation: err3( x) := fx ftaylor3( x) err4( x) := fx ftaylor4( x) Fifth orer approximation: err ( x) := fx ftaylor ( x) The various approximations generate by the above calculations an the associate errors are compare in the Table.1.

7 48 AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Taylor series approx of given function Function approximate by Taylor series ftaylor( x) ftaylor3( x) ftaylor( x) 1 Zeroth orer approx Thir orer approx Fifth orer approx/ given function x x- value Figure.1. Taylor series approximation of given function x.1 := errors as function of x,.. err( x) err3( x) err( x) 1 Errors gen ue to Taylor-series approx Zero orer approx Thir orer approx Fifth orer approx x- value Figure.. Errors associate with the various Taylor series approximations x

8 Chapter : Introuction to Numerical Methos 49 The various approximations generate by the above calculations an the associate errors are compare in the following tables. h := 1.. x :=, 1.. Table.1 Various orers of approximation generate by Taylor series approach versus true values of given function zeroth orer first orer secon orer thir orer fourth orer fifth orer True Value h = x = ftaylor( x) ftaylor1 = ( x) =ftaylor( x) Table. ftaylor3 = ( x) =ftaylor4( x) Errors associate with the ifferent orers of approximation =ftaylor ( x) =f( x) = h = zeroth orer first orer secon orer thir orer fourth orer fifth orer x = err( x) = err1( x) = err( x) = err3( x) = err4( x) = err ( x) = PROBLEMS.1. Using the Taylor series expansion for cos x, which is given as f(x) = cos x = 1- x / + x 4 / 4, etermine the one-term, two-term an three-term approximations of f(x + h ), where x =

9 AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD ra an h =. 1, ra, an compare these with the exact solution. Using Mathca, the inepenent variable x. generate plots of the various Taylor series approximations an associate errors as functions of. Develop a Taylor series expansion of the following function: f(x) = x - 6 x 4 + 3x + 9. Use x =3 as the base an h as the increment. Using Mathca, evaluate the series for h=.1,...1., aing terms incrementally as in Problem.1. Compare the various Taylor series approximations obtaine with true values in a table. Generate plots of the approximations an associate errors as functions of x..3 Given the following function: f(x) = x 3-3 x + x + 1, etermine f ( x + h ) with the help of a Taylor series expansion, where x = an h =.4. Compare the true value of f (.4 ) with estimates obtaine by resorting to (a) one term only (b) two terms (c) three terms an () four terms of the series..4 Given the following function fx = 3x 3 6x + 1 x + use a Taylor series expansion to etermine the zeroth, first, secon an thir orer approximations of f(x +h) where x = an h =.. solution. Compare these with the exact. By eveloping a Taylor series expansion for f(x) = e x about x =, etermine the fourth-orer approximation of e. an compare it with the exact solution..6. By eveloping a Taylor series expansion for f(x) = ln(-x) about x =, etermine the fourth-orer approximation of ln (.) an compare it with the exact solution

10 Chapter : Introuction to Numerical Methos 1.7. By eveloping a Taylor series expansion for f(x) = x 3 e - x about x = 1, etermine the thir-orer approximation of f(1.) an compare it with the exact solution..8. By eveloping a Taylor series expansion for f(x) = e cos x about x =, etermine the fourth- orer approximation of f (π) an compare it with the exact solution...9. By eveloping a Taylor series expansion for f(x) = (x - ) 1/ about x =, etermine the thir-orer approximation of f (.), that is, (.) 1/, an compare it with the exact solution..1. Given the function f(x) = x - x. + 6, use a Taylor series expansion to etermine the first, secon, thir an fourth orer approximations of f (. ) by resorting to x = an h =.. Compare these with the exact solution..11. Given the function f(x) = 6 x 3-9 x + x + 4, use a Taylor series expansion to etermine the zeroth, first, secon an thir orer approximations of f (x + h ) where x = 3 an h = 1. Compare these with the exact solution..1. Given the function f(x) = 4 x 4-7 x 3 + x - 6 x + 9 use a Taylor series expansion to etermine the zeroth, first, secon, thir an fourth orer

11 AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD approximations of f (x + h ) ) where x = 3 an h =.. Compare these with the exact solution. Calculate errors abn generate calculations to three ecimal places..13. Given the function f(x) = 8 x 3-1 x + x + 4, use a Taylor series expansion to etermine the zeroth, first, secon an thir orer approximations of f (x + h ) where x = an h = 1. Compare these with the exact solution..14. Given the function f(x) = 1 + x + x /! + x 3 / 3! +x 4 / 4! use a Taylor series expansion to etermine the zeroth, first, secon, thir an fourth orer approximations of f (x + h ) where x = an h =.. Compare these with the exact solution. Generate answers correct to four ecimal places..1. Given the function f(x) = x + x 3 / 3 + x / 1 use a Taylor series expansion to etermine the zeroth, first, secon, thir an fourth orer approximations of f (x + h ) where x = an h =.8. Compare these with the exact solution by computing percentage errors. Generate answers correct to four ecimal places..16. Given the function f(x) = sin (x) use a Taylor series expansion to etermine the fifth orer approximation of f (x + h ) where x = an h =. raians. Compare your answer with the true value. Generate answers correct to four ecimal places..17. Given the function f(x) = 3 x - 6 x. + 9, solution by computing percentage errors. Generate answers correct to four ecimal places. use a Taylor series expansion to etermine the zeroth, first, secon, thir an fourth orer approximations of f (x + h ) whereo x = 3 an h = 1. Compare these with the exact

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 13. Khyruddin Akbar Ansari, Ph.D., P.E.

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 13. Khyruddin Akbar Ansari, Ph.D., P.E. AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 13 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering Gonzaga University SDC PUBLICATIONS Schroff

More information

TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS

TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS MISN-0-4 TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS f(x ± ) = f(x) ± f ' (x) + f '' (x) 2 ±... 1! 2! = 1.000 ± 0.100 + 0.005 ±... TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS by Peter Signell 1.

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

Linear and quadratic approximation

Linear and quadratic approximation Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

Differentiation ( , 9.5)

Differentiation ( , 9.5) Chapter 2 Differentiation (8.1 8.3, 9.5) 2.1 Rate of Change (8.2.1 5) Recall that the equation of a straight line can be written as y = mx + c, where m is the slope or graient of the line, an c is the

More information

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E.

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E. AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathcad Release 14 Khyruddin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering and Applied Science Gonzaga University

More information

ECE 422 Power System Operations & Planning 7 Transient Stability

ECE 422 Power System Operations & Planning 7 Transient Stability ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of

More information

Solutions to MATH 271 Test #3H

Solutions to MATH 271 Test #3H Solutions to MATH 71 Test #3H This is the :4 class s version of the test. See pages 4 7 for the 4:4 class s. (1) (5 points) Let a k = ( 1)k. Is a k increasing? Decreasing? Boune above? Boune k below? Convergant

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

2Algebraic ONLINE PAGE PROOFS. foundations

2Algebraic ONLINE PAGE PROOFS. foundations Algebraic founations. Kick off with CAS. Algebraic skills.3 Pascal s triangle an binomial expansions.4 The binomial theorem.5 Sets of real numbers.6 Surs.7 Review . Kick off with CAS Playing lotto Using

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

MATH , 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208

MATH , 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208 MATH 321-03, 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208 Instructor: Brent Deschamp Email: brent.eschamp@ssmt.eu Office: McLaury 316B Phone:

More information

Single Variable Calculus Warnings

Single Variable Calculus Warnings Single Variable Calculus Warnings These notes highlight number of common, but serious, first year calculus errors. Warning. The formula g(x) = g(x) is vali only uner the hypothesis g(x). Discussion. In

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

More from Lesson 6 The Limit Definition of the Derivative and Rules for Finding Derivatives.

More from Lesson 6 The Limit Definition of the Derivative and Rules for Finding Derivatives. Math 1314 ONLINE More from Lesson 6 The Limit Definition of the Derivative an Rules for Fining Derivatives Eample 4: Use the Four-Step Process for fining the erivative of the function Then fin f (1) f(

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1 We consier ifferential equations of the form y + a()y = b(), (1) y( 0 ) = y 0, where a() an b() are functions. Observe that this class of equations inclues equations

More information

A Weak First Digit Law for a Class of Sequences

A Weak First Digit Law for a Class of Sequences International Mathematical Forum, Vol. 11, 2016, no. 15, 67-702 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1288/imf.2016.6562 A Weak First Digit Law for a Class of Sequences M. A. Nyblom School of

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

3.2 Differentiability

3.2 Differentiability Section 3 Differentiability 09 3 Differentiability What you will learn about How f (a) Might Fail to Eist Differentiability Implies Local Linearity Numerical Derivatives on a Calculator Differentiability

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

YORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics. MATH A Test #2. June 25, 2014 SOLUTIONS

YORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics. MATH A Test #2. June 25, 2014 SOLUTIONS YORK UNIVERSITY Faculty of Science Department of Mathematics an Statistics MATH 505 6.00 A Test # June 5, 04 SOLUTIONS Family Name (print): Given Name: Stuent No: Signature: INSTRUCTIONS:. Please write

More information

TIME-DELAY ESTIMATION USING FARROW-BASED FRACTIONAL-DELAY FIR FILTERS: FILTER APPROXIMATION VS. ESTIMATION ERRORS

TIME-DELAY ESTIMATION USING FARROW-BASED FRACTIONAL-DELAY FIR FILTERS: FILTER APPROXIMATION VS. ESTIMATION ERRORS TIME-DEAY ESTIMATION USING FARROW-BASED FRACTIONA-DEAY FIR FITERS: FITER APPROXIMATION VS. ESTIMATION ERRORS Mattias Olsson, Håkan Johansson, an Per öwenborg Div. of Electronic Systems, Dept. of Electrical

More information

ECEn 621 Computer Arithmetic. Method Comparisons

ECEn 621 Computer Arithmetic. Method Comparisons ECEn 61 Computer Arithmetic Chapter 7: Iterative Approximation Slie #1 Metho Comparisons Digit Recurrence Metho Computes Multiplication Division/Square Root Log/Exp Operations Use Shifting Aition Simple

More information

An inductance lookup table application for analysis of reluctance stepper motor model

An inductance lookup table application for analysis of reluctance stepper motor model ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR

More information

Section 7.1: Integration by Parts

Section 7.1: Integration by Parts Section 7.1: Integration by Parts 1. Introuction to Integration Techniques Unlike ifferentiation where there are a large number of rules which allow you (in principle) to ifferentiate any function, the

More information

6 Wave equation in spherical polar coordinates

6 Wave equation in spherical polar coordinates 6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.

More information

Mathcad Lecture #5 In-class Worksheet Plotting and Calculus

Mathcad Lecture #5 In-class Worksheet Plotting and Calculus Mathca Lecture #5 In-class Worksheet Plotting an Calculus At the en of this lecture, you shoul be able to: graph expressions, functions, an matrices of ata format graphs with titles, legens, log scales,

More information

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

Fall 2016: Calculus I Final

Fall 2016: Calculus I Final Answer the questions in the spaces provie on the question sheets. If you run out of room for an answer, continue on the back of the page. NO calculators or other electronic evices, books or notes are allowe

More information

11.7. Implicit Differentiation. Introduction. Prerequisites. Learning Outcomes

11.7. Implicit Differentiation. Introduction. Prerequisites. Learning Outcomes Implicit Differentiation 11.7 Introuction This Section introuces implicit ifferentiation which is use to ifferentiate functions expresse in implicit form (where the variables are foun together). Examples

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Final Exam Study Guide and Practice Problems Solutions

Final Exam Study Guide and Practice Problems Solutions Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

We G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies

We G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll

More information

Average value of position for the anharmonic oscillator: Classical versus quantum results

Average value of position for the anharmonic oscillator: Classical versus quantum results verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive

More information

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0. Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?

More information

Polynomial Inclusion Functions

Polynomial Inclusion Functions Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl

More information

New Statistical Test for Quality Control in High Dimension Data Set

New Statistical Test for Quality Control in High Dimension Data Set International Journal of Applie Engineering Research ISSN 973-456 Volume, Number 6 (7) pp. 64-649 New Statistical Test for Quality Control in High Dimension Data Set Shamshuritawati Sharif, Suzilah Ismail

More information

First Order Linear Differential Equations

First Order Linear Differential Equations LECTURE 6 First Orer Linear Differential Equations A linear first orer orinary ifferential equation is a ifferential equation of the form ( a(xy + b(xy = c(x. Here y represents the unknown function, y

More information

1 The Derivative of ln(x)

1 The Derivative of ln(x) Monay, December 3, 2007 The Derivative of ln() 1 The Derivative of ln() The first term or semester of most calculus courses will inclue the it efinition of the erivative an will work out, long han, a number

More information

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu

More information

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0. Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?

More information

Thermal runaway during blocking

Thermal runaway during blocking Thermal runaway uring blocking CES_stable CES ICES_stable ICES k 6.5 ma 13 6. 12 5.5 11 5. 1 4.5 9 4. 8 3.5 7 3. 6 2.5 5 2. 4 1.5 3 1. 2.5 1. 6 12 18 24 3 36 s Thermal runaway uring blocking Application

More information

MATH2231-Differentiation (2)

MATH2231-Differentiation (2) -Differentiation () The Beginnings of Calculus The prime occasion from which arose my iscovery of the metho of the Characteristic Triangle, an other things of the same sort, happene at a time when I ha

More information

Derivative of a Constant Multiple of a Function Theorem: If f is a differentiable function and if c is a constant, then

Derivative of a Constant Multiple of a Function Theorem: If f is a differentiable function and if c is a constant, then Bob Brown Math 51 Calculus 1 Chapter 3, Section Complete 1 Review of the Limit Definition of the Derivative Write the it efinition of the erivative function: f () Derivative of a Constant Multiple of a

More information

Module 2. DC Circuit. Version 2 EE IIT, Kharagpur

Module 2. DC Circuit. Version 2 EE IIT, Kharagpur Moule 2 DC Circuit Lesson 9 Analysis of c resistive network in presence of one non-linear element Objectives To unerstan the volt (V ) ampere ( A ) characteristics of linear an nonlinear elements. Concept

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Math Chapter 2 Essentials of Calculus by James Stewart Prepared by Jason Gaddis

Math Chapter 2 Essentials of Calculus by James Stewart Prepared by Jason Gaddis Math 231 - Chapter 2 Essentials of Calculus by James Stewart Prepare by Jason Gais Chapter 2 - Derivatives 21 - Derivatives an Rates of Change Definition A tangent to a curve is a line that intersects

More information

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (II)

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (II) Laplace s Equation in Cylinrical Coorinates an Bessel s Equation (II Qualitative properties of Bessel functions of first an secon kin In the last lecture we foun the expression for the general solution

More information

Calculus Math Fall 2012 (Cohen) Lecture Notes

Calculus Math Fall 2012 (Cohen) Lecture Notes Calculus Math 70.200 Fall 202 (Cohen) Lecture Notes For the purposes of this class, we will regar calculus as the stuy of limits an limit processes. Without yet formally recalling the efinition of a limit,

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

d dx [xn ] = nx n 1. (1) dy dx = 4x4 1 = 4x 3. Theorem 1.3 (Derivative of a constant function). If f(x) = k and k is a constant, then f (x) = 0.

d dx [xn ] = nx n 1. (1) dy dx = 4x4 1 = 4x 3. Theorem 1.3 (Derivative of a constant function). If f(x) = k and k is a constant, then f (x) = 0. Calculus refresher Disclaimer: I claim no original content on this ocument, which is mostly a summary-rewrite of what any stanar college calculus book offers. (Here I ve use Calculus by Dennis Zill.) I

More information

18 EVEN MORE CALCULUS

18 EVEN MORE CALCULUS 8 EVEN MORE CALCULUS Chapter 8 Even More Calculus Objectives After stuing this chapter you shoul be able to ifferentiate an integrate basic trigonometric functions; unerstan how to calculate rates of change;

More information

Solutions to the Exercises of Chapter 9

Solutions to the Exercises of Chapter 9 9A. Vectors an Forces Solutions to the Exercises of Chapter 9. F = 5 sin 5.9 an F = 5 cos 5 4.8.. a. By the Pythagorean theorem, the length of the vector from to (, ) is + = 5. So the magnitue of the force

More information

Chapter 2 Derivatives

Chapter 2 Derivatives Chapter Derivatives Section. An Intuitive Introuction to Derivatives Consier a function: Slope function: Derivative, f ' For each, the slope of f is the height of f ' Where f has a horizontal tangent line,

More information

The Sokhotski-Plemelj Formula

The Sokhotski-Plemelj Formula hysics 24 Winter 207 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),

More information

Module FP2. Further Pure 2. Cambridge University Press Further Pure 2 and 3 Hugh Neill and Douglas Quadling Excerpt More information

Module FP2. Further Pure 2. Cambridge University Press Further Pure 2 and 3 Hugh Neill and Douglas Quadling Excerpt More information 5548993 - Further Pure an 3 Moule FP Further Pure 5548993 - Further Pure an 3 Differentiating inverse trigonometric functions Throughout the course you have graually been increasing the number of functions

More information

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique J. Basic. Appl. Sci. Res., ()9-99,, TextRoa Publication ISSN 9-434 Journal of Basic an Applie Scientific Research www.textroa.com A Comparison between a Conventional Power System Stabilizer (PSS) an Novel

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

Outline. MS121: IT Mathematics. Differentiation Rules for Differentiation: Part 1. Outline. Dublin City University 4 The Quotient Rule

Outline. MS121: IT Mathematics. Differentiation Rules for Differentiation: Part 1. Outline. Dublin City University 4 The Quotient Rule MS2: IT Mathematics Differentiation Rules for Differentiation: Part John Carroll School of Mathematical Sciences Dublin City University Pattern Observe You may have notice the following pattern when we

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Scott E. Grasman 1, Zaki Sari 2 and Tewfik Sari 3

Scott E. Grasman 1, Zaki Sari 2 and Tewfik Sari 3 RAIRO Operations Research RAIRO Oper. Res. 41 (27) 455 464 DOI: 1.151/ro:2731 NEWSVENDOR SOLUTIONS WITH GENERAL RANDOM YIELD DISTRIBUTIONS Scott E. Grasman 1, Zaki Sari 2 an Tewfik Sari 3 Abstract. Most

More information

February 21 Math 1190 sec. 63 Spring 2017

February 21 Math 1190 sec. 63 Spring 2017 February 21 Math 1190 sec. 63 Spring 2017 Chapter 2: Derivatives Let s recall the efinitions an erivative rules we have so far: Let s assume that y = f (x) is a function with c in it s omain. The erivative

More information

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements Aaptive Gain-Scheule H Control of Linear Parameter-Varying Systems with ime-delaye Elements Yoshihiko Miyasato he Institute of Statistical Mathematics 4-6-7 Minami-Azabu, Minato-ku, okyo 6-8569, Japan

More information

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable

More information

The Press-Schechter mass function

The Press-Schechter mass function The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Math 2163, Practice Exam II, Solution

Math 2163, Practice Exam II, Solution Math 63, Practice Exam II, Solution. (a) f =< f s, f t >=< s e t, s e t >, an v v = , so D v f(, ) =< ()e, e > =< 4, 4 > = 4. (b) f =< xy 3, 3x y 4y 3 > an v =< cos π, sin π >=, so

More information

Define each term or concept.

Define each term or concept. Chapter Differentiation Course Number Section.1 The Derivative an the Tangent Line Problem Objective: In this lesson you learne how to fin the erivative of a function using the limit efinition an unerstan

More information

Adjustable Fractional-Delay Filters Utilizing the Farrow Structure and Multirate Techniques

Adjustable Fractional-Delay Filters Utilizing the Farrow Structure and Multirate Techniques Ajustable Fractional-Delay Filters Utilizing the Farrow Structure an Multirate Techniques Håkan Johansson (hakanj@isy.liu.se)* an Ewa Hermanowicz (hewa@eti.pg.ga.pl)** *Division of Electronics Systems,

More information

Taylor Expansions in 2d

Taylor Expansions in 2d Taylor Expansions in 2 In your irst year Calculus course you evelope a amily o ormulae or approximating a unction F(t) or t near any ixe point t 0. The cruest approximation was just a constant. F(t 0 +

More information

Modeling time-varying storage components in PSpice

Modeling time-varying storage components in PSpice Moeling time-varying storage components in PSpice Dalibor Biolek, Zenek Kolka, Viera Biolkova Dept. of EE, FMT, University of Defence Brno, Czech Republic Dept. of Microelectronics/Raioelectronics, FEEC,

More information

1 Lecture 20: Implicit differentiation

1 Lecture 20: Implicit differentiation Lecture 20: Implicit ifferentiation. Outline The technique of implicit ifferentiation Tangent lines to a circle Derivatives of inverse functions by implicit ifferentiation Examples.2 Implicit ifferentiation

More information

SYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d

SYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d SYNTHESIS ON THE ASSESSMENT OF HIPS ONTRATION OEFFIIENT Marius MILEA Technical University Gheorghe Asachi of Iași, Machine Manufacturing an Inustrial Management Faculty, B-ul D Mangeron, nr 59, A, 700050,

More information

A note on the Mooney-Rivlin material model

A note on the Mooney-Rivlin material model A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel

More information

The Chain Rule. d dx x(t) dx. dt (t)

The Chain Rule. d dx x(t) dx. dt (t) The Chain Rule The Problem You alreay routinely use the one imensional chain rule t f xt = f x xt x t t in oing computations like t sint2 = cost 2 2t In this example, fx = sinx an xt = t 2. We now generalize

More information

Quantum mechanical approaches to the virial

Quantum mechanical approaches to the virial Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from

More information

Section The Chain Rule and Implicit Differentiation with Application on Derivative of Logarithm Functions

Section The Chain Rule and Implicit Differentiation with Application on Derivative of Logarithm Functions Section 3.4-3.6 The Chain Rule an Implicit Differentiation with Application on Derivative of Logarithm Functions Ruipeng Shen September 3r, 5th Ruipeng Shen MATH 1ZA3 September 3r, 5th 1 / 3 The Chain

More information