ETNA Kent State University
|
|
- Matthew Mitchell
- 5 years ago
- Views:
Transcription
1 Electronic Transactions on Numerical Analysis. Volume 34, pp , Copyrigt 2008,. ISSN ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS CARMEN ARÉVALO Dedicated to Víctor Pereyra on te occasion of is 70t birtday Abstract. Wen differential-algebraic equations of index 3 or iger are solved wit backward differentiation formulas, te solution can ave gross errors in te first few steps, even if te initial values are equal to te exact solution and even if te stepsize is kept constant. Tis raises te question of wat are consistent initial values for te difference equations. Here we study ow to cange te exact initial values into wat we call numerically consistent initial values for te implicit Euler metod. Key words. ig index differential-algebraic equations, consistent initial values AMS subject classifications. 65L05 1. Introduction. A differential-algebraic equation (DAE) as te form F(t, x, ẋ) = 0, were te matrix F/ ẋ is singular. Here we sall consider te differential-algebraic equation of te form ṗ = U(t, v) v = F(t, p, v) + G(t, p, v)λ (1.1a) (1.1b) 0 = R(t, p), (1.1c) were t R, p R n, v R m, λ R s, and U : R R m R n, F : R R n R m R m, G : R R n R m R m s and R : R R n R s. Assume tat s min(m, n) in order to avoid an over-determined system. Te system (1.1) as index 3 if R/ p U/ v G is nonsingular for all t [t 0, T]. For simplicity we take t 0 = 0. Te algebraic variables appear linearly as in (1.1b) in important classes of pysical problems. Tis condition is fulfilled, for example, by te Euler-Lagrange equations of multibody mecanics, wic ave applications in biomecanics, te dynamics of macinery, robotics, and veicle design. To solve tis type of DAE, several tecniques ave been considered. One proposition as been to solve te system in its original formulation using a backward differentiation formula (BDF), as implemented in te DAE solver DASSL [4]; but suc a variable-step-size, variableorder code based on BDF metods presents some essential difficulties wen solving iger index DAEs, especially in te accuracy of te algebraic variables [2]. Te initial values (p 0, v 0, λ 0 ) are said to be consistent if te DAE as a differentiable solution (p(t), v(t), λ(t)) in te interval [0, T] suc tat (p(0), v(0), λ(0)) = (p 0, v 0, λ 0 ). Brenan and Engquist [3] defined numerically consistent starting values to order k + 1 for te k-step BDF applied to (1.1) as starting values [p k 1,..., p 0 ], [v k 1,..., v 0 ], [λ k 1,..., λ 0 ] Received Marc 31, Accepted August 6, Publised online on December 13, Recommended by Godela Scerer. Numerical Analysis, Centre for Matematical Sciences, Lund University, Box 118, SE Lund, Sweden. (Carmen.Arevalo@na.lu.se) 14
2 A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES 15 suc tat p j p(t j ) K 1 k+1 v j v(t j ) K 2 k+1 R(t j, p j ) K 3 k+2 (1.2a) (1.2b) (1.2c) for some constants K 1, K 2, K 3 and j = 0, 1,...,k 1. Tey proved tat te k-step BDF (k = 1,...,6) wit constant step size converges globally wit O( k ) accuracy to te solutions, but only after after k+1 steps, and provided te metod uses numerically consistent starting values to order k + 1. In particular, for te implicit Euler metod te state variables p and v ave O() accuracy after te first step, but te error in te algebraic variable λ is O(1), even wen te initial values are exact. By approximating te errors in te algebraic variable, a correction mecanism was devised in order to obtain O() accuracy even in te initial k steps [2]. Tese formulas correct te errors locally and produce O() accurate algebraic variables. Our view is tat tese O(1) errors are caused by initial values tat are inconsistent wit te difference equations. We attempt to redefine numerically consistent initial values as initial values tat are consistent wit te difference equations, as opposed to tose consistent wit te differential equation, as in (1.2). We explain te O(1) errors in te algebraic variables as a result of starting te BDF solver wit initial values tat are consistent wit te differential equation, but not wit te difference equations. Once tis is establised, we develop a sceme to construct numerically consistent initial values. We illustrate our results by solving te same problems tat appear in [1]. 2. Numerically consistent initial values. Initial values tat are numerically consistent wit te difference equation generated by an order p metod sould produce O( p ) accurate solutions. DEFINITION 2.1. Te set (x 0, x 1,..., x k 1 ) is a numerically consistent initial value at t = t 0 for a differential equation f(t, x, ẋ) = 0, solved wit a k-step multistep metod of order p, if te metod initiated wit (x 0, x 1,...,x k 1 ) generates te approximation x k suc tat x(t k ) x k = O( p ). (2.1) For differential-algebraic equations (DAEs) in semi-explicit form ẋ = f(t, x, λ) 0 = g(t, x, λ) of index less tan 2, consistent initial conditions are also numerically consistent; but for iger index DAEs, consistent initial conditions do not necessarily produce an approximation satisfying (2.1). For index 3 DAEs, in fact, BDF metods approximate te algebraic variables wit O(1) errors after te first step, even if te initial values are exact. 3. Euler-Lagrange equations and te implicit Euler formula. Te index 3 formulation of constrained multibody systems is ṗ = v M v = F(p, v) + G(p) T λ (3.1a) (3.1b) 0 = R(p), (3.1c)
3 16 C. ARÉVALO were M is te positive definite mass matrix, G = R/ p, and G(p)M 1 G(p) T is nonsingular. To simplify notation, we will denote G(p 1 ) by G 1 and F(p 1, v 1 ) by F 1. Wen te implicit Euler formula is applied to (3.1) using exact initial values, te algebraic variable λ 1 as an O(1) error ǫ λ = λ(t 1 ) λ 1 tat can be approximated wit O() accuracy by [2] δλ = (G 1 M 1 G T 1 ) 1 G 1 M 1 F 1 + λ 1. Te implicit Euler metod requires initial values only for state variables p and v. Taking te initial condition p 0 = p(0) + O( 3 ) v 0 = v(0) M 1 G T 1 δλ + O(2 ), (3.2a) (3.2b) te first step is p 1 = p 0 + v 1 v 1 = v 0 + (I B)M 1 F 1 + M 1 G T 1 δλ 0 = R(p 1 ), were te matrix B is dependent on p 1 and is defined by B = M 1 G T 1 (G 1 M 1 G T 1 ) 1 G 1. Te matrix B is a projector and B M 1 G T 1 = M 1 G T 1. We denote te global errors at t 1 by ǫ p = p(t 1 ) p 1, ǫ v = v(t 1 ) v 1 and ǫ λ = λ(t 1 ) λ 1. Te O() accurate correction term δλ may be expressed in terms of te exact value of te algebraic variable as δλ = 1 2 ((G 1M 1 G T 1 ) 1 G 1 M 1 F 1 + λ(t 1 )). By expanding p(t) and v(t) about t 1 and evaluating at t = 0 we obtain ǫ v = M 1 G T 1 (δλ + ǫλ ) + O( 2 ) ǫ p = ǫ v 2 2 M 1 (F(p(t 1 ), v(t 1 )) + G(p(t 1 )) T λ(t 1 )) + O( 3 ). Tese expressions lead to ǫ p = 2 2 M 1 (G T 1 (G 1 M 1 G T 1 ) 1 G 1 M 1 F 1 F(t 1 )) + 2 M 1 G T 1 ǫ λ + O( 3 ), and multiplying by G 1 we can conclude tat G 1 ǫ p = 2 G 1 M 1 G T 1 ǫλ + O( 3 ). If we expand R(p) about p 1 and use te fact tat ǫ p = O( 2 ) togeter wit R(p(t 1 )) = 0 and R(p 1 ) = 0, we finally obtain G 1 ǫ p = O( 4 ), tus implying tat ǫ λ = O(). In oter words, te initial condition (3.2) is numerically consistent for (3.1) solved wit te implicit Euler metod.
4 A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES 17 Te initial value (3.2b) may be written as Equation (3.3) in te form v 0 = v(0) + B(v(0) v 1 ) + O( 2 ). (3.3) v 0 v(0) + O( 2 ) = B (v 1 v(0)) (3.4) can be interpreted as an indication tat te initial values for v consistent wit te differential equation are also consistent wit te difference equation only in te case tat v(0) is in te null space of G 1. In practical terms, tis would require tat λ(t 0 ) = (G 0 M 1 G T 0 ) 1 G 0 M 1 F 0 + O(). As te cange in initial values for te state variable v is O(), it will not affect te accuracy of te solution in any of te state variables. 4. Te more general case. Te more general problem we consider ere is ṗ = U(t, q) q = F(t, p, q) + G(t, p, q)λ (4.1a) (4.1b) 0 = R(t, p). (4.1c) Tis is an index 3 DAE if R p U q G is non-singular. According to [2], te algebraic variable Λ is estimated by te first step of te implicit Euler metod wit an O(1) error ǫ Λ = Λ(t 1 ) Λ 1, wic is approximated wit O() accuracy by δλ = (R p U q G) 1 R p [U q (F + GΛ)U t ], were all functions are evaluated at t = t 1. Te first implicit Euler step is p 1 = p 0 + U(t 1, q 1 ) q 1 = q 0 + (F(t 1, p 1, q 1 ) + G(t 1, p 1, q 1 )Λ 1 ) 0 = R(t 1, p 1 ) wit numerically consistent initial condition p 0 = p(0) + O( 3 ) q 0 = q(0) GδΛ + O( 2 ). (4.2a) (4.2b) We may write (4.2b) as q 0 = q(0) AU q (q 1 q(0)) AU t + O( 2 ) wit A = G(R p U q G) 1 R p. Te matrix AU q is a projector wit AU q G = G. Te last equation can be written as q 0 q(0) = AU q q 1 q(0) AU t + O(), mirroring equation (3.4) wen U is independent of t. In tis case, te numerically consistent initial values for q are consistent wit te differential equation if AU q q = 0, tat is, only if q is in te null space of AU q. In practical terms, tis would require tat Λ(0) = (R p U q G) 1 R p U q F t=0.
5 18 C. ARÉVALO 5. Computational results. We illustrate our results wit te following two index 3 DAEs from [1]. PROBLEM 5.1. Te system of Euler-Lagrange equations as te solution ẍ = 2y + xλ ÿ = 2x + yλ (5.1a) (5.1b) 0 = x 2 + y 2 1 (5.1c) x = sin(1 + t) 2 y = cos(1 + t) 2 u = 2(1 + t)cos(1 + t) 2 v = 2(1 + t)sin(1 + t) 2 λ = 4(1 + t) 2. PROBLEM 5.2. Tis problem is defined in te form (4.1) wit p = [ x y z ] T q = [ u v w ] T Λ = [ λ β ] T U = [ 2u v w 1 ] T F = [ y 2x + y sin t 2 4yt 4zt sint 2] T [ ] T x 0 2z G = 0 2y 1 R = [ x 2 + y 2 + z 2 1 z 0.5 ] T. Its solution is 3 x = 2 cost2 y = 3 2 sin t2 z = u = 2 t sin t2 v = 3t cost 2 w = 1 λ = 2t 2 β = 0.5 sint 2. Problems 5.1 and 5.2 were solved wit implicit Euler using = and = After te first step te algebraic variable λ sows O(1) errors. Problem 5.1 was solved wit t 0 = 0, so G(P 0 ) T λ 0 = [sin(1) cos(1)] λ 0. Terefore, λ(0) would ave to be close to zero (but is instead equal to -4) for te exact initial values to be numerically consistent. Te numerically consistent initial values were calculated from (3.3) and te step was redone. Te exact initial values [ u v ] = [ ] were canged to [ ] for = and to [ ] for = Tis resulted in algebraic variables of O() accuracy (5.1). For Problem 5.2 wit t 0 = 1 and = 0.001, te initial values for q were canged from [ ] to [ ], producing algebraic variables of O() accuracy, as sown in Table 5.2. Te results presented in Figure 5.1 sow te O(1) errors wen Problem 5.1 was solved using exact initial values, and te O() errors in te algebraic variable obtained by using numerically consistent initial values.
6 A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES 19 TABLE 5.1 Absolute errors in te algebraic variable λ for Problem 5.1 using implicit Euler wit step sizes and O(1) errors and teir corrected values are displayed in bold face. = = t j numerically numerically exact IV consistent IV exact IV consistent IV TABLE 5.2 Absolute errors in te algebraic variable λ for Problem 5.2 using implicit Euler wit step sizes and O(1) errors and teir corrected values are displayed in bold face. = = t j numerically numerically exact IV consistent IV exact IV consistent IV Implicit Euler solution wit exact initial values log global error = = = time Implicit Euler solution wit numerically consistent initial values log global error = = = time FIG Implicit Euler used to solve an index 3 DAE test problem (see text) using numerically consistent initial values. Global errors wit exact initial values (top) and numerically consistent initial values (bottom). REFERENCES [1] C. ARÉVALO, Matcing te structure of DAEs and multistep metods, P. D. tesis, Department of Computer Science, Lund University, Lund, Sweden, [2] C. ARÉVALO AND P. LÖDSTEDT, Improving te accuracy of BDF metods for index 3 differential-algebraic equations, BIT, 35 (1995), pp [3] K. BRENAN AND B. ENGQUIST, Backward differentiation approximations of nonlinear differential/algebraic systems, Mat. Comp., 51 (1988), pp [4] L. PETZOLD, Differential/algebraic equations are not ODEs, SIAM J. Sci. Statist. Comput., 3 (1982), pp
Numerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationUniversity Mathematics 2
University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More informationAMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.
Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error
More informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow
1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation
More informationReflection Symmetries of q-bernoulli Polynomials
Journal of Nonlinear Matematical Pysics Volume 1, Supplement 1 005, 41 4 Birtday Issue Reflection Symmetries of q-bernoulli Polynomials Boris A KUPERSHMIDT Te University of Tennessee Space Institute Tullaoma,
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationSection 3.1: Derivatives of Polynomials and Exponential Functions
Section 3.1: Derivatives of Polynomials and Exponential Functions In previous sections we developed te concept of te derivative and derivative function. Te only issue wit our definition owever is tat it
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Matematics National Institute of Tecnology Durgapur Durgapur-7109 email: anita.buie@gmail.com 1 . Capter 8 Numerical Solution of Ordinary
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More informationBlanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS
Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationRecent Progress in the Integration of Poisson Systems via the Mid Point Rule and Runge Kutta Algorithm
Recent Progress in te Integration of Poisson Systems via te Mid Point Rule and Runge Kutta Algoritm Klaus Bucner, Mircea Craioveanu and Mircea Puta Abstract Some recent progress in te integration of Poisson
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationComputer Derivations of Numerical Differentiation Formulae. Int. J. of Math. Education in Sci. and Tech., V 34, No 2 (March-April 2003), pp
Computer Derivations o Numerical Dierentiation Formulae By Jon H. Matews Department o Matematics Caliornia State University Fullerton USA Int. J. o Mat. Education in Sci. and Tec. V No (Marc-April ) pp.8-87.
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationNUMERICAL DIFFERENTIATION
NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationInfluence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations
International Journal of Difference Equations ISSN 0973-6069, Volume 12, Number 2, pp. 281 302 (2017) ttp://campus.mst.edu/ijde Influence of te Stepsize on Hyers Ulam Stability of First-Order Homogeneous
More informationMATHEMATICS FOR ENGINEERING DIFFERENTIATION TUTORIAL 1 - BASIC DIFFERENTIATION
MATHEMATICS FOR ENGINEERING DIFFERENTIATION TUTORIAL 1 - BASIC DIFFERENTIATION Tis tutorial is essential pre-requisite material for anyone stuing mecanical engineering. Tis tutorial uses te principle of
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More information1 Solutions to the in class part
NAME: Solutions to te in class part. Te grap of a function f is given. Calculus wit Analytic Geometry I Exam, Friday, August 30, 0 SOLUTIONS (a) State te value of f(). (b) Estimate te value of f( ). (c)
More informationName: Answer Key No calculators. Show your work! 1. (21 points) All answers should either be,, a (finite) real number, or DNE ( does not exist ).
Mat - Final Exam August 3 rd, Name: Answer Key No calculators. Sow your work!. points) All answers sould eiter be,, a finite) real number, or DNE does not exist ). a) Use te grap of te function to evaluate
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationHOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS
HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3
More informationQuaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016.
Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim1887@aol.com rev 1 Aug 8, 216 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial
More informationBob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk
Bob Brown Mat 251 Calculus 1 Capter 3, Section 1 Completed 1 Te Tangent Line Problem Te idea of a tangent line first arises in geometry in te context of a circle. But before we jump into a discussion of
More informationLogarithmic functions
Roberto s Notes on Differential Calculus Capter 5: Derivatives of transcendental functions Section Derivatives of Logaritmic functions Wat ou need to know alread: Definition of derivative and all basic
More informationLab 6 Derivatives and Mutant Bacteria
Lab 6 Derivatives and Mutant Bacteria Date: September 27, 20 Assignment Due Date: October 4, 20 Goal: In tis lab you will furter explore te concept of a derivative using R. You will use your knowledge
More informationAPPENDIXES. Let the following constants be established for those using the active Mathcad
3 APPENDIXES Let te following constants be establised for tose using te active Matcad form of tis book: m.. e 9.09389700 0 3 kg Electron rest mass. q.. o.6077330 0 9 coul Electron quantum carge. µ... o.5663706
More informationMODIFIED DIFFERENTIAL EQUATIONS. Dedicated to Prof. Michel Crouzeix
ESAIM: PROCEEDINGS, September 2007, Vol.21, 16-20 Gabriel Caloz & Monique Dauge, Editors MODIFIED DIFFERENTIAL EQUATIONS Pilippe Cartier 1, Ernst Hairer 2 and Gilles Vilmart 1,2 Dedicated to Prof. Micel
More informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationPractice Problem Solutions: Exam 1
Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More information2.1 THE DEFINITION OF DERIVATIVE
2.1 Te Derivative Contemporary Calculus 2.1 THE DEFINITION OF DERIVATIVE 1 Te grapical idea of a slope of a tangent line is very useful, but for some uses we need a more algebraic definition of te derivative
More informationSome Review Problems for First Midterm Mathematics 1300, Calculus 1
Some Review Problems for First Midterm Matematics 00, Calculus. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd,
More information(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?
Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More information2.3 Product and Quotient Rules
.3. PRODUCT AND QUOTIENT RULES 75.3 Product and Quotient Rules.3.1 Product rule Suppose tat f and g are two di erentiable functions. Ten ( g (x)) 0 = f 0 (x) g (x) + g 0 (x) See.3.5 on page 77 for a proof.
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationChapter 2. Limits and Continuity 16( ) 16( 9) = = 001. Section 2.1 Rates of Change and Limits (pp ) Quick Review 2.1
Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 969) Quick Review..... f ( ) ( ) ( ) 0 ( ) f ( ) f ( ) sin π sin π 0 f ( ). < < < 6. < c c < < c 7. < < < < < 8. 9. 0. c < d d < c
More informationParametric Spline Method for Solving Bratu s Problem
ISSN 749-3889 print, 749-3897 online International Journal of Nonlinear Science Vol4202 No,pp3-0 Parametric Spline Metod for Solving Bratu s Problem M Zarebnia, Z Sarvari 2,2 Department of Matematics,
More informationDIRECTLY SOLVING SECOND ORDER LINEAR BOUNDARY VALUE PROBLEMS OF ORDINARY DIFFERENTIAL EQUATIONS. Ra ft Abdelrahim 1, Z. Omar 2
International Journal of Pure and Applied Matematics Volume 6 No. 6, -9 ISSN: - (printed version); ISSN: -95 (on-line version) url: ttp://www.ijpam.eu doi:.7/ijpam.v6i. PAijpam.eu DIRECTLY SOLVING SECOND
More informationAnalytic Functions. Differentiable Functions of a Complex Variable
Analytic Functions Differentiable Functions of a Complex Variable In tis capter, we sall generalize te ideas for polynomials power series of a complex variable we developed in te previous capter to general
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationTeaching Differentiation: A Rare Case for the Problem of the Slope of the Tangent Line
Teacing Differentiation: A Rare Case for te Problem of te Slope of te Tangent Line arxiv:1805.00343v1 [mat.ho] 29 Apr 2018 Roman Kvasov Department of Matematics University of Puerto Rico at Aguadilla Aguadilla,
More information3.1 Extreme Values of a Function
.1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find
More informationMath 312 Lecture Notes Modeling
Mat 3 Lecture Notes Modeling Warren Weckesser Department of Matematics Colgate University 5 7 January 006 Classifying Matematical Models An Example We consider te following scenario. During a storm, a
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More informationChapter 2 Limits and Continuity
4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationOn One Justification on the Use of Hybrids for the Solution of First Order Initial Value Problems of Ordinary Differential Equations
Pure and Applied Matematics Journal 7; 6(5: 74 ttp://wwwsciencepublisinggroupcom/j/pamj doi: 648/jpamj765 ISSN: 6979 (Print; ISSN: 698 (Online On One Justiication on te Use o Hybrids or te Solution o First
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips
More informationParabolic PDEs: time approximation Implicit Euler
Part IX, Capter 53 Parabolic PDEs: time approximation We are concerned in tis capter wit bot te time and te space approximation of te model problem (52.4). We adopt te metod of line introduced in 52.2.
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationQuantum Mechanics Chapter 1.5: An illustration using measurements of particle spin.
I Introduction. Quantum Mecanics Capter.5: An illustration using measurements of particle spin. Quantum mecanics is a teory of pysics tat as been very successful in explaining and predicting many pysical
More informationGELFAND S PROOF OF WIENER S THEOREM
GELFAND S PROOF OF WIENER S THEOREM S. H. KULKARNI 1. Introduction Te following teorem was proved by te famous matematician Norbert Wiener. Wiener s proof can be found in is book [5]. Teorem 1.1. (Wiener
More informationf a h f a h h lim lim
Te Derivative Te derivative of a function f at a (denoted f a) is f a if tis it exists. An alternative way of defining f a is f a x a fa fa fx fa x a Note tat te tangent line to te grap of f at te point
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationTwo Step Hybrid Block Method with Two Generalized Off-step Points for Solving Second Ordinary Order Differential Equations Directly
Global Journal of Pure and Applied Matematics. ISSN 0973-768 Volume 2, Number 2 (206), pp. 59-535 Researc India Publications ttp://www.ripublication.com/gjpam.tm Two Step Hybrid Block Metod wit Two Generalized
More informationPreface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.
Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationLecture 21. Numerical differentiation. f ( x+h) f ( x) h h
Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More information2.3 More Differentiation Patterns
144 te derivative 2.3 More Differentiation Patterns Polynomials are very useful, but tey are not te only functions we need. Tis section uses te ideas of te two previous sections to develop tecniques for
More informationSolving Continuous Linear Least-Squares Problems by Iterated Projection
Solving Continuous Linear Least-Squares Problems by Iterated Projection by Ral Juengling Department o Computer Science, Portland State University PO Box 75 Portland, OR 977 USA Email: juenglin@cs.pdx.edu
More information7.1 Using Antiderivatives to find Area
7.1 Using Antiderivatives to find Area Introduction finding te area under te grap of a nonnegative, continuous function f In tis section a formula is obtained for finding te area of te region bounded between
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
More information