The Laplace equation, cylindrically or spherically symmetric case
|
|
- Alexis Small
- 5 years ago
- Views:
Transcription
1 Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential, waves or vibrations are just a few typical applications were u = (Laplace equation) u = q (Poisson equation) occur. In tese equations, denotes te Laplace-Operator, wic in treedimensional cartesian xyz-coordinates @y Te cylindrical coordinate system uses polar coordinates r and to locate a point in te x y plane and te coordinate z for te eigt of te point above te plane. Te Laplace operator in tis coordinate system @r + + For functions depending on r (distance from z axis) only (cylindrical symmetry), tis reduces : () Similarly, in sperical symmetry (were r measures te distance from te : () Tese differential operators (,) occur in te following examples, wic terefore in fact are just special cases of partial differential equations, reduced to ordinary differential equations by symmetry. Many features of our simple examples are also typical for real Laplace- and Poisson equations. We use examples in cylindrical symmetry ere just because one-dimensional versions of te Laplace operator () in cartesian geometry would not include first derivatives. Solutions in tat case are so simple tat tere is not muc to learn from. (On te oter and: wat s wrong wit simple examples? We sall discuss also te simple cases.)
2 Steady-state diffusion or eat conduction A cylindrically symmetric Laplace equation, cf. () governs stationary eat conduction or diffusion wit axial symmetry. An example would be a cylindrical pipe, interior radius r, exterior radius r. Te temperature = (r) in te wall of te pipe follows te equation rr + r r = for r int < r < r ext : (4) Diriclet boundary conditions specify te temperature at te interior and exterior walls, (r int ) = int (r ext ) = ext ( inner wall temperature ), ( outer wall temperature ), (Note: tere are several ways to denote derivatives of a function like (r). d (r) dr ; ; rr all stand for te same differential operator. Depending on context, one way or anoter may be more convenient to work wit.) For tis problem, te exact solution is known. We will use it to ceck te accuracy of our numerical approximations. (r) = int log(r=r ext ) + ext log(r int =r) : log(r int =r ext ) Alternatively, we migt prescribe te eat flow rate q at te inner boundary (A Diriclet condition at te rigt and a Neumann condition at te left end). In tat case, Te exact solution ere is r (r int ) = ; (r ext ) = ext : (r) = ext + r int log r r ext : For tis differential equation, Neumann boundary conditions on bot sides are no good coice. If te corresponding rates are equal, te solution is unique only up to an arbitrary additive constant. If te rates are not equal, no steady state can develop; tus, no solution exists Tis simple example demonstrates tat boundary value problems do not necessarily ave a solution. (in contrast to initial value problems, were existence of a solution is guaranteed in te vicinity of te starting point for a wide class of functions) Diriclet and Neumann conditions are also called boundary conditions of te first and second kind, respectively. Tere are tird-kind boundary
3 conditions as well (some attac to tem te names Robin or Sturm, or call tem radiation conditions), of te form c p(r) + c p (r) =. Also nonlinear boundary conditions may occur, like c p(r) 4 + c p (r) =, wic would model termal radiation according to Stefan-Boltzmann s law. Cooling of cylindrical objects To find solutions of te time-dependent eat equation in cylindrical symmetry, t = a rr + r r we set (r; t) = e at u(r), wic reduces te problem to an ordinary differential equation for u(r), u rr + r u r + u = ; < r < R (5) wic is (after multiplying by r and a transformation of variables r r) a form of Bessel s differential equation. Symmetry dictates a Neumann boundary condition u r = at te center r =. For an outer Diriclet boundary condition u = at r = R, it turns out tat nontrivial solutions exist for certain values of only (te trivial solution would be u(r) = for all r). Te solution involves J, te Bessel function of order zero of te first kind, u(r) = AJ (r) were A is an arbitrary constant (remember: te problem is a omogeneous one; any multiple of a solution is also a solution). Te outer boundary condition ten requires to be a root of J (R) = Steady one-dimensional convection-diffusion equation Let c = c(x) denote te concentration of some quantity in a flow wit velocity w and a coefficient of diffusion D. Ten Dc wc = for < x < L; c = c at x = ; c = c L at x = L;
4 models a steady-state situation were convective transport is balanced by diffusion. Tis problem as te exact solution x c(r) = c + epe L e Pe (c L c ): Here Pe is te Peclet number, Pe = wl=d. Ferzinger and Perić comment on tis problem: Because it is so simple, tis problem is often used as a test of numerical metods, including bot discretization and solution scemes. [... ] Tere are few actual flows in wic tis balance plays an important role. Normally, convection is balanced by eiter a pressure gradient or diffusion in te direction normal to te flow. [... ] Indeed, use of tis problem as a test case as probably produced more poor coices of metod tan any oter in te field. Discretization on an Equidistant Grid, Diriclet Conditions We discuss te example + r = for = (r); < r int < r < r ext (r int ) = int ; (r ext ) = ext : Te first step in te application of finite difference metods is te selection of a set of mes points in te interval (r int ; r ext ). Let us simply take n equidistant radius points r j, r j = r int + j wit mes size = r ext r int ; j = ; : : : ; n n + Boundary points are r = r int and r n+ = r ext. Next, we represent te differential equation by difference equations involving values of te unknown pressure at te mes points. A straigtforward way to do so is to replace te derivatives by appropriate difference quotients. Standard approximations are te difference formulae f (x) f (x) f(x + ) f(x ) f(x + ) f(x) + f(x ) Notation: For te value (r) of te unknown temperature at te mes point r j we write j. At mes point r j+ = r j +, it is j+, similarly j at r j. Tis way we may write j ( j+ j ) 4
5 j ( j+ j + j ) Substitution leads to difference equations, one per interior mes point. Te temperature values j at all interior mes points fulfill approximately r j j + j + r j j+ = for j = ; : : : ; n: (6) Tis is a system of n linear algebraic equations. In matrix form, A x = b wit (7) A = 6 4 ( + r ) ( r ) ( + r ) ( r ) ( + r ) ( rn ) ( + rn ) ( rn ) 7 5 ; ~x = 6 4. n ; ~ b = ( r ) int. ( + r n ) ext Note tat tis matrix is not symmetric, but diagonally dominant as long as < r. : 7 5 Exercises Exercise Solve te diriclet problem (4) for steady-state radial flow for two cases.. Strong influence of radial geometry: r int = :; r ext = ; int = ; ext = : 5
6 . moderate influence of radial geometry: To evaluate te results, plot r int = ; r ext = ; int = ; ext = :. Te error as a function of r for some values of n, say n = ; ;.. Te maximum error as a function of n for a range of n as large as computer power allows. Sow te error in a log-log-plot and estimate te asymptotic beavior as n. An important quantity is te slope of te curve. Wat we sall find out: For large n te maximum error decreases proportional to n. Tis quadratic convergence, owever, does not become apparent unless te grid size is comparable to te interior radius r int. Te memory required is proportional to n. Can you specify te memory requirements of your program more precisely? Exercise Derive a difference approximation for te convection-diffusion equation and set up te system of linear equations. Is te matrix symmetric? diagonally dominant? Try to solve te simple case D = ; L = ; w = ; c = ; c L = for various, and compare results wit te exact solution. Cange w = and try again. Evaluate your results wit similar plots as in Exercise. Exercise Te matrix (7) in our radial flow example is not symmetric. Sow tat by multiplying te i-t equation by r i a symmetric matrix results. Systems wit symmetric matrices are teoretically and numerically in many ways easier to andle tan non-symmetric ones. Implement te symmetric version and ceck wit te original version. Exercise 4 Test te sensitivity of te system (7) to small perturbations in te rigt-and side. In our MATLAB code, insert % o r i g i n a l s o l u t i o n p = A\b ; %s o l u t i o n wit p e r t u r b e d r i g t and s i d e % t r y d i f f e r e n t p e r t u r b a t i o n s e r r=r a n d ( n, ). ; %e r r =(rand (n,).5).; b = b+e r r ; p e r r = A\b ; 6
7 Te relative error in b is norm(err)/norm(b), wile te relative error in te solution is norm(p-perr)/norm(p). How muc larger is te relative error in te solution in relation to te relative error of te rigt-and side? Calculate tis quotient (a measure of te sensitivity of te linear system) for different values of n and estimate te beavior for large n. Sow a loglog-plot. In MATLAB, c = cond(a) returns te -norm condition number, te ratio of te largest singular value of A to te smallest. Te condition number of a matrix measures te sensitivity of te solution of a system of linear equations to errors in te data. It gives an indication of te accuracy of te results from matrix inversion and te linear equation solution. It also is an upper bound on te sensitivity, as calculated before. Values of cond(a) near indicate a well-conditioned matrix: errors in te data generate errors of te same magnitude in te solution. Add a plot of te condition number depending on n to your loglog-plot. Is te condition number a quite sarp or a rater pessimistic estimate on te sensitivity in tis example? Discretization on an Equidistant Grid, Mixed Boundary Conditions Just a few modifications in our first program are necessary to implement a Neumann condition at te inner radius and a Diriclet condition at te outer boundary. We are now going to solve our second example of flow in porous media, te problem p + r p = for p = p(r); r w < r < r e p (r w ) = ; (8) p(r e ) = p e : We may use te same grid and te same difference formula as before, but tere is one problem: For te first interior gridpoint r, te difference approximation is r p + p + r p = : (9) However, te boundary condition does not give us any value for p, it just specifies te derivative tere. We can exploit tis information by using a one-sided approximation of te first derivative, f (x) = f(x) + f(x + ) 7 f (x) + O( ): ()
8 In our case, we set p p = ; or p = p : () Inserting for p in Equation (9), we obtain a difference equation for te first interior gridpoint at a Neumann boundary, + p + p = r r r : () Tere is an inconsistency wit using tis approximation along wit te rest of te difference equations. Because of te one-sided approximation (), is a fist-order approximation only. At all oter points, te difference equations are second-order approximations to te original differential equation. Te first-order approximation of te Neumann condition does indeed contaminate our final results in te wole region. In general, te order of approximation of a problem is te lowest order used in any parts of te problem. Sometimes, owever, a bad approximation at just a few points will not affect te overall accuracy. To get a second-order accurate approximation at r, we can replace te second derivative in using te differential equation. In our case, p (r ) = p p p (r ) = r p (r ) = r p (r ) + O( ) from te differential equation from te boundary condition From tese equations follows p = p + r () More commonly, owever, Neumann boundaries are treated by modifying te grid. Instead of r = r w ; r = r w + ; te modified grid at a Neumann boundary as All gridpoints are given by r = r w ; r = r w + r i = r w + (i =) wit = r e r w n + = ; i = ; : : : ; n 8
9 Te point r is now outside of te region were our problem is defined. (It is called a gost point.) As a result of tis modification, te approximation p = p is of second order. However, tis sceme does not compute te pressure p wf directly, since tere is no gridpoint at r w. We may set p wf = (p + p )= + O( ). Exercises Exercise 5 Solve te mixed-type problem (8) and test te tree ways to implement Neumann boundary conditions. Data: r w = :; r e = ; p (r w ) = = ; p e = : To evaluate te results, plot te maximum error as a function of n in a log-log-plot and estimate te asymptotic beavior as n. Te difference between first-order and second-order accuracy sould become clearly apparent. Exercise 6 Solve for r w = ; and r e = te inomogeneos differential equation wit omogeneous mixed-type boundary conditions (quasi-steadystate flow), p + r p = for r w < r < r e ; p (r e ) = ; p(r w ) = : Evaluate te accuracy of your solution as in Exercise. (If you ave coded tis exercise, te modifications in te code sould be rater marginal) Exercise 7 Discretize te differential equation (5). Assume simple values for a and. Wat appens to te rigt-and side of te system? Wat solution gives MATLAB s A\b? Any ideas ow to proceed in tis case? Deriving Difference Approximations for Derivatives Consider a smoot function f(x). ( Smoot ere means tat as many derivatives of f exixt as we may need.) We will illustrate, by way of example, a metod to find approximations to derivatives. It is sometimes called te metod of undetermined coefficients. 9
10 Example: Derive an approximation for f at gridpoint x using values at gridpoints x ; x and x of an equidistant grid wit spacing. Expand f and f in Taylor series around x f = f f + f 6 f f iv + O(5 ) f = f + f + f + 6 f f iv + O(5 ) We now seek a linear combination of te tree values f ; f and f approximating f in te form f af + bf + cf : To determine te coefficients a; b and c, we insert te Taylor series and collect terms in f ; f ; : : :. We find f (a + b + c)f + ( a + c)f + ( a + c)f (a + c)f iv + O(5 ) (a + c)f + Since we ave tree degrees of freedom (tree undetermined coefficients), we can plan on setting te coefficients of f and f equal to zero and te coefficient of f equal to one. We are left wit te following system of equations. a + b + c = a + c = a + c = Solving te above system of equations gives a = ; b = ; c = : As an additional bonus, tis solution makes te coefficient of f equal to zero too. Te coefficient of f iv turns out to be =. Hence, we get te approximation f f + f = f + f iv + O( ): Tis is te well-known centered difference approximation tat we ave used already.
11 As an additional exercise, you may derive a five-point centered approximation to f. Te system of equations for five coefficients a; b; c; d and e in tis case will be a + b + c + d + e = a b + d + e = 4a + b + d + 4e = 8a b + d + 8e = 6= 6a + b + d + 6e = (Using an algebraic manipulator suc as Maple or Matematica makes solving suc systems muc easier) Te solution is a = ; b = ; c = ; d = ; e = : An alternative metod: find an interpolating polynomial p(x) for te data points (x i ; f i ). Differentiate te polynomial wit respect to x and evaluate at x. If you do tis by and, you sould use te Lagrange interpolation formula: Te polynom p(x) tat interpolates te n + pairs of values (x ; f ); (x ; f ); : : : ; (x n ; f n ) is given by p(x) = f L (x) + f L (x) + + f n L n (x) ; were L k (x) = (x x )(x x ) (x x k )(x x k+ ) (x x n ) (x k x )(x k x ) (x k x k )(x k x k+ ) (x k x n ) for eac k = ; ; : : : ; n. In Matematica, tis metod conveniently finds difference approximations for iger derivatives. Te commands data={{-,f[-]},{-,f[-]},{,f[]},{,f[]},{,f[]}}; poly=interpolatingpolynomial[data,x]; appr=d[poly,{x,}]/.x-> //Togeter
12 find an interpolating polynomial and evaluate its tird derivative at x =. Matematica gives appr = f( ) + f( ) f() + f( ) : We may find te truncation error by te command Series[appr,{,,}] wic produces te output f () () + f (5) () 4 + O() 4 : Since many tables provide difference approximations for f ; f ; : : :, it is normally not necessary to derive basic formulae. However, for irregular grids, complicated equations or some specific boundary conditions, tese metods are of great value. Non-Equidistant Grid Write f = f(x ); f = f(x); f + = f(x + + ); Ten f (x) = = r i r i ; + = r i+ r i : f (x) = (f + f ) + + (f f ) f + f f f f f ( + ) + + ( ) f iv Note tat in te second approximation, for + 6= now tere is a loworder truncation error term involving f. Equidistant grids usually provide approximations wit iger-order truncation errors. So you ave to balance te gain by finer messize against possibly iger truncation error. Linear solvers Te Tomas Algoritm for tridiagonal Systems Te so-called Tomas Algoritm is just a form of elimination for solving tridiagonal systems of linear equations. Llewellyn H. Tomas used it around 95
13 to solve elliptic partial differential equations. Te attribution to Tomas seems to be more common in some engineering disciplines tan it is in numerical analysis. Let A be a tridiagonal matrix and b te rigt-and side of a system Ax = b. A = 6 4 d e c d e c d e c n d n e n ; b = b b b. b n : 7 5 c n d n b n Transform it to upper triangular form, by Gaussian elimination. Use backsubstitution to solve te resulting equivalent system Ax = b e b e b e b A = ; b = :. 6 4 e 7 6 n 5 4 b 7 n 5 Te expressions for te entries in A and b follow from a step-by-step application of te usual elimination procedure. for i = ; : : : ; n: e i = e = e d ; b = b d ; e i d i c i e i ; b = b i c i b i i ; d i c i e i Backsubstitution: x n = b n, und f"ur i = n ; : : : ; : Tis algoritm is stable if 8 >< >: d i > ; i = ; ; : : : ; n; d > je j; x i = b i e i x i+: d i jc i j + je i j und c i 6= ; e i 6= ; i = ; : : : ; n ; d n jc n j: b n
14 Essentially tese conditions require A to be (weakly) diagonally dominant. Tey are sufficient but not necessary. 4
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 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 informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
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 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 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 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 informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
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 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 informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
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 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 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 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 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 informationMath 34A Practice Final Solutions Fall 2007
Mat 34A Practice Final Solutions Fall 007 Problem Find te derivatives of te following functions:. f(x) = 3x + e 3x. f(x) = x + x 3. f(x) = (x + a) 4. Is te function 3t 4t t 3 increasing or decreasing wen
More informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
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 informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
More informationCalculus I Practice Exam 1A
Calculus I Practice Exam A Calculus I Practice Exam A Tis practice exam empasizes conceptual connections and understanding to a greater degree tan te exams tat are usually administered in introductory
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 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 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 informationHigher Derivatives. Differentiable Functions
Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.
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 informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More information1. Which one of the following expressions is not equal to all the others? 1 C. 1 D. 25x. 2. Simplify this expression as much as possible.
004 Algebra Pretest answers and scoring Part A. Multiple coice questions. Directions: Circle te letter ( A, B, C, D, or E ) net to te correct answer. points eac, no partial credit. Wic one of te following
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 informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
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 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 informationPolynomials 3: Powers of x 0 + h
near small binomial Capter 17 Polynomials 3: Powers of + Wile it is easy to compute wit powers of a counting-numerator, it is a lot more difficult to compute wit powers of a decimal-numerator. EXAMPLE
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 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 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 informationarxiv: v1 [math.na] 3 Nov 2011
arxiv:.983v [mat.na] 3 Nov 2 A Finite Difference Gost-cell Multigrid approac for Poisson Equation wit mixed Boundary Conditions in Arbitrary Domain Armando Coco, Giovanni Russo November 7, 2 Abstract In
More informationMAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points
MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim
More informationTime (hours) Morphine sulfate (mg)
Mat Xa Fall 2002 Review Notes Limits and Definition of Derivative Important Information: 1 According to te most recent information from te Registrar, te Xa final exam will be eld from 9:15 am to 12:15
More information. If lim. x 2 x 1. f(x+h) f(x)
Review of Differential Calculus Wen te value of one variable y is uniquely determined by te value of anoter variable x, ten te relationsip between x and y is described by a function f tat assigns a value
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 informationAverage Rate of Change
Te Derivative Tis can be tougt of as an attempt to draw a parallel (pysically and metaporically) between a line and a curve, applying te concept of slope to someting tat isn't actually straigt. Te slope
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 informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
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 informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More information5.1 We will begin this section with the definition of a rational expression. We
Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go
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 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 informationPrecalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!
Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More information1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist
Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter
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 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 informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
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 informationALGEBRA AND TRIGONOMETRY REVIEW by Dr TEBOU, FIU. A. Fundamental identities Throughout this section, a and b denotes arbitrary real numbers.
ALGEBRA AND TRIGONOMETRY REVIEW by Dr TEBOU, FIU A. Fundamental identities Trougout tis section, a and b denotes arbitrary real numbers. i) Square of a sum: (a+b) =a +ab+b ii) Square of a difference: (a-b)
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 information2.3 Algebraic approach to limits
CHAPTER 2. LIMITS 32 2.3 Algebraic approac to its Now we start to learn ow to find its algebraically. Tis starts wit te simplest possible its, and ten builds tese up to more complicated examples. Fact.
More informationMA455 Manifolds Solutions 1 May 2008
MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using
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 informationNotes on wavefunctions II: momentum wavefunctions
Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles
More informationDerivatives. By: OpenStaxCollege
By: OpenStaxCollege Te average teen in te United States opens a refrigerator door an estimated 25 times per day. Supposedly, tis average is up from 10 years ago wen te average teenager opened a refrigerator
More informationSECTION 1.10: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES
(Section.0: Difference Quotients).0. SECTION.0: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES Define average rate of cange (and average velocity) algebraically and grapically. Be able to identify, construct,
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
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 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 information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
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 informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
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 information, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim
Mat 132 Differentiation Formulas Stewart 2.3 So far, we ave seen ow various real-world problems rate of cange and geometric problems tangent lines lead to derivatives. In tis section, we will see ow to
More informationSection 15.6 Directional Derivatives and the Gradient Vector
Section 15.6 Directional Derivatives and te Gradient Vector Finding rates of cange in different directions Recall tat wen we first started considering derivatives of functions of more tan one variable,
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 informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
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 informationRecall from our discussion of continuity in lecture a function is continuous at a point x = a if and only if
Computational Aspects of its. Keeping te simple simple. Recall by elementary functions we mean :Polynomials (including linear and quadratic equations) Eponentials Logaritms Trig Functions Rational Functions
More informationChapter 1D - Rational Expressions
- Capter 1D Capter 1D - Rational Expressions Definition of a Rational Expression A rational expression is te quotient of two polynomials. (Recall: A function px is a polynomial in x of degree n, if tere
More informationMATH1131/1141 Calculus Test S1 v8a
MATH/ Calculus Test 8 S v8a October, 7 Tese solutions were written by Joann Blanco, typed by Brendan Trin and edited by Mattew Yan and Henderson Ko Please be etical wit tis resource It is for te use of
More informationApplied Linear Statistical Models. Simultaneous Inference Topics. Simultaneous Estimation of β 0 and β 1 Issues. Simultaneous Inference. Dr.
Applied Linear Statistical Models Simultaneous Inference Dr. DH Jones Simultaneous Inference Topics Simultaneous estimation of β 0 and β 1 Bonferroni Metod Simultaneous estimation of several mean responses
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More information1 Upwind scheme for advection equation with variable. 2 Modified equations: numerical dissipation and dispersion
1 Upwind sceme for advection equation wit variable coefficient Consider te equation u t + a(x)u x Applying te upwind sceme, we ave u n 1 = a (un u n 1), a 0 u n 1 = a (un +1 u n ) a < 0. CFL condition
More information1watt=1W=1kg m 2 /s 3
Appendix A Matematics Appendix A.1 Units To measure a pysical quantity, you need a standard. Eac pysical quantity as certain units. A unit is just a standard we use to compare, e.g. a ruler. In tis laboratory
More informationChapter Seven The Quantum Mechanical Simple Harmonic Oscillator
Capter Seven Te Quantum Mecanical Simple Harmonic Oscillator Introduction Te potential energy function for a classical, simple armonic oscillator is given by ZÐBÑ œ 5B were 5 is te spring constant. Suc
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 informationFinancial Econometrics Prof. Massimo Guidolin
CLEFIN A.A. 2010/2011 Financial Econometrics Prof. Massimo Guidolin A Quick Review of Basic Estimation Metods 1. Were te OLS World Ends... Consider two time series 1: = { 1 2 } and 1: = { 1 2 }. At tis
More informationMath 1241 Calculus Test 1
February 4, 2004 Name Te first nine problems count 6 points eac and te final seven count as marked. Tere are 120 points available on tis test. Multiple coice section. Circle te correct coice(s). You do
More informationThe Derivative The rate of change
Calculus Lia Vas Te Derivative Te rate of cange Knowing and understanding te concept of derivative will enable you to answer te following questions. Let us consider a quantity wose size is described by
More informationCOMPUTATIONAL COMPARISON BETWEEN THE TAYLOR HOOD AND THE CONFORMING CROUZEIX RAVIART ELEMENT
Proceedings of ALGORITMY 5 pp. 369 379 COMPUTATIONAL COMPARISON BETWEEN E TAYLOR HOOD AND E CONFORMING OUZEIX RAVIART ELEMENT ROLF KRAHL AND EBERHARD BÄNSCH Abstract. Tis paper is concerned wit te computational
More informationMATH 1020 Answer Key TEST 2 VERSION B Fall Printed Name: Section #: Instructor:
Printed Name: Section #: Instructor: Please do not ask questions during tis exam. If you consider a question to be ambiguous, state your assumptions in te margin and do te best you can to provide te correct
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationIntroduction to Machine Learning. Recitation 8. w 2, b 2. w 1, b 1. z 0 z 1. The function we want to minimize is the loss over all examples: f =
Introduction to Macine Learning Lecturer: Regev Scweiger Recitation 8 Fall Semester Scribe: Regev Scweiger 8.1 Backpropagation We will develop and review te backpropagation algoritm for neural networks.
More informationPolynomial Functions. Linear Functions. Precalculus: Linear and Quadratic Functions
Concepts: definition of polynomial functions, linear functions tree representations), transformation of y = x to get y = mx + b, quadratic functions axis of symmetry, vertex, x-intercepts), transformations
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 6. Differential Calculus 6.. Differentiation from First Principles. In tis capter, we will introduce
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
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 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 information