MATH745 Fall MATH745 Fall
|
|
- Brian McGee
- 6 years ago
- Views:
Transcription
1 MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics bprotas@mcmasterca Office HH 36, Ext 46 Course Webpage: ttp://wwwmatmcmasterca/ bprotas/math745 Wat is NUMERICAL ANALYSIS? Development of COMPUTATIONAL ALGORITHMS for solutions of problems in algebra and analysis Use of metods of MATHEMATICAL ANALYSIS to determine a priori properties of tese algoritms suc as: ACCURACY, STABILITY, CONVERGENCE Application of tese algoritms to solve actual problems arising in practice MATH745 Fall 5 3 MATH745 Fall 5 4 BASICS a I PART I Finite Differences A Review ASSUMPTIONS : f : Ω R is a smoot function, ie is continuously differentiable sufficiently many times, te domain Ω = [a,b] is discretized wit a uniform grid {x = a,,x N = b}, suc tat x + x = = extensions to nonuniform grids are straigtforward PROBLEM given te nodal values of te function f, ie, f = f x, =,, N approximate te nodal values of te function derivative d f dx x = f x, =,,N a Details can be found in any standard textbook on elementary numerical analysis, eg, K Atkinson and W Han, Elementary Numerical Analysis, Wiley, 4
2 MATH745 Fall 5 5 MATH745 Fall 5 6 BASICS II Te simplest approac Derivation of finite difference formulae via TAYLOR SERIES EXPANSIONS f + = f + x + x f + x + x Rearrange te expansion! = f + f f = f + f f + x + x 3 f + 3! f + = f + f +O, were O α denotes te contribution from all terms wit powers of greater or equal α ere α = Neglecting O, we obtain a FIRST ORDER FORWARD DIFFERENCE FORMULA : = f + f BASICS III Backward difference formula is obtained by expanding f about x and proceeding as before: f = f f f + = = f f Neglected term wit te lowest power of is te LEADING ORDER APPROXIMATION ERROR, ie, Err = f x C α Te exponent α of in te leading order error represents te ORDER OF ACCURACY OF THE METHOD it tells ow quickly te approximation error vanises wen te resolution is refined Te actual value of te approximation error depends on te constant C caracterizing te function f In te examples above Err = f, ence te metods are FIRST ORDER ACCURATE MATH745 Fall 5 7 MATH745 Fall 5 8 HIGHER ORDER FORMULA I Consider two expansions: Subtracting te second from te first: Central Difference Formula f = f + f Te leading order error is ACCURATE f + = f + f f = f f f + f = f f + + = = f + f, tus te metod is SECOND ORDER Manipulating four different Taylor series expansions one can obtain a fourt order central difference formula : = f f + 8 f + f, Err = 4 3 f v APPROXIMATION OF THE SECOND DERIVATIVE Consider two expansions: Adding te two expansions f + = f + f f = f f f + + f = f f iv + Central difference formula for te second derivative: f = f + f + f f iv + = Te leading order error is ACCURATE f iv δ f = f + f + f, tus te metod is SECOND ORDER
3 MATH745 Fall 5 9 MATH745 Fall 5 AN ALTERNATIVE APPROACH I An alternative derivation of a finite difference sceme: Find an N t order accurate interpolating function px wic interpolates te function f x at te nodes x, =,,N, ie, suc tat px = f x, =,,N Differentiate te interpolating function px and evaluate at te nodes to obtain an approximation of te derivative p x f x, =,,N Example: for =,,N, let te interpolant ave te form of a quadratic polynomial p x on [x,x + ] Lagrange interpolating polynomial p x = x x x x + f + x x x x + f + x x x x f + p x = x x x + f + x x x + f + x x x f + Evaluating at x = x we obtain f x p x = f + f ie, second order accurate center difference formula AN ALTERNATIVE APPROACH II Generalization to iger orders straigtforward Example: for = 3,,N, one can use a fourt order polynomial as interpolant p x on [x,x + ] Differentiating wit respect to x and evaluating at x = x we arrive at te fourt order accurate finite difference formula = f f + 8 f + f, Err = 4 3 f v Order of accuracy of te finite difference formula is one less tan te order of te interpolating polynomial Te set of grid points needed to evaluate a finite difference formula is called STENCIL In general, iger order formulas ave larger stencils MATH745 Fall 5 MATH745 Fall 5 TAYLOR TABLE I A general metod for coosing te coefficients of a finite difference formula to ensure te igest possible order of accuracy Example: consider a one sided finite difference formula p= α p f +p, were te coefficients α p, p =,, are to be determined Form an expression for te approximation error f p= and expand it about x in te powers of α p f +p = ε TAYLOR TABLE II Expansions can be collected in a Taylor table f f f f f a f a a f + a a a a 3 6 a f + a a a a 3 6 te leftmost column contains te terms present in te expression for te approximation error te corresponding rows multiplied by te top row represent te terms obtained from expansions about x columns represent terms wit te same order in sums of columns are te contributions to te approximation error wit te given order in Te coefficients α p, p =,, can now be cosen to cancel te contributions to te approximation error wit te lowest powers of
4 MATH745 Fall 5 3 MATH745 Fall 5 4 TAYLOR TABLE III Setting te coefficients of te first tree terms to zero: a a a = a a = = a = 3 a a, a =, a = = Te resulting formula: = f f + 3 f Te approximation error determined te evaluating te first column wit non zero coefficient: a 3 6 a 3 f = 6 3 f Te formula is tus SECOND ORDER ACCURATE AN OPERATOR PERSPECTIVE I Quick review of FUNCTIONAL ANALYSIS background NORMED SPACES X: : X R suc tat x,y X Banac spaces x, x + y x + y, x = x vector spaces: finite dimensional R N vs infinite dimensional l p function spaces on Ω R N : Lebesgue spaces L p Ω, Sobolev spaces W p,q Ω Hilbert spaces: inner products, ortogonality & proections, bases, etc Linear Operators: operator norms, functionals, Riesz Teorem MATH745 Fall 5 5 MATH745 Fall 5 6 AN OPERATOR PERSPECTIVE II Assume tat f and f belong to a function space X; DIFFERENTIATION d dx : f f can ten be regarded as a LINEAR OPERATOR dx d : X X Wen f and f are approximated by teir nodal values as f = [ f f f N ] T and f = [ f f f N ]T, ten te differential operator dx d can be approximated by a DIFFERENTIATION MATRIX A R N N suc tat f = Af ; How can we determine tis matrix? Assume for simplicity tat te domain Ω is periodic, ie, f = f N and f = f N+ ; ten differentiation wit te second order center difference formula can be represented as te following matrix vector product f f = f N f N AN OPERATOR PERSPECTIVE III Using te fourt order center difference formula we would obtain a pentadiagonal system increased order of accuracy entails increased bandwidt of te differentiation matrix A A is a TOEPLITZ MATRIX, since is as constant entries along te te diagonals; in fact, it a also a CIRCULANT MATRIX wit entries a i depending only on i mod N Note tat te matrix A defined above is SINGULAR as a zero eigenvalue λ = Wy? Tis property is in fact inerited from te original continuous operator d dx wic is also singular and as a zero eigenvalue A singular matrix A does not ave an inverse at least, now in te classical sense; wat can we do to get around tis difficulty?
5 MATH745 Fall 5 7 MATH745 Fall 5 8 AN OPERATOR PERSPECTIVE IV Matrix singularity linearly dependent rows te LHS vector does not contain enoug information to determine UNIQUELY te RHS vector MATRIX DESINGULARIZATION incorporating additional information into te matrix, so tat its argument te RHS vector can be determined uniquely Example desingularization of te second order center difference differentiation matrix: in a center difference formula, even and odd nodes are decoupled knowing f, =,,N and f, one can recover f, = 3,5, ie, te odd nodes only f must also be provided ence, te zero eigenvalue as multiplicity two wen desingularizing te differentiation matrix one must modify at least two rows see, eg, sing_diff_mat_m AN OPERATOR PERSPECTIVE V Wat is WRONG wit te differentiation operator? Te differentiation operator dx d is UNBOUNDED! One usually cannot find a constant C R independent of f, suc tat d dx f x X C f X, f X For instance, f x = e ikx, so tat C = k for k Unfortunately, finite dimensional emulations of te differentiation operator te DIFFERENTIATION MATRICES inerit tis property OPERATOR NORM for matrices A = max x = Ax = max x Ax, Ax x,x x,a T Ax = max = λ max A T A = σ x x,x maxa Tus, te norm of a matrix is given by te square root of its largest SINGULAR VALUE σ max A MATH745 Fall 5 9 MATH745 Fall 5 AN OPERATOR PERSPECTIVE VI As can be rigorously proved in many specific cases, A grows witout bound as N or, tis is a reflection of te unbounded nature of te underlying dim operator Te loss of precision wen solving te system Ax = b is caracterized by te CONDITION NUMBER wit respect to inversion κ p A = A p A p for p =, κ A = σ maxa σ min A wen te condition number is large, te matrix is said to be ILL CONDITIONED solution of te system Ax = b is prone to round off errors if A is singular, κ p A = + SUBTRACTIVE CANCELLATION ERRORS SUBTRACTIVE CANCELLATION ERRORS wen comparing two numbers wic are almost te same using finite precision aritmetic, te relative round off error is proportional to te inverse of te difference between te two numbers Tus, if te difference between te two numbers is decreased by an order of magnitude, te relative accuracy wit wic tis difference may be calculated using finite precision aritmetic is also decreased by an order of magnitude Problems wit finite difference formulae wen loss of precision due to finite precision aritmetic SUBTRACTIVE CANCELLATION, eg, for double precision: 345 e 8% error 34 e 3 9% error
6 MATH745 Fall 5 MATH745 Fall 5 COMPLEX STEP DERIVATIVE a PADÉ APPROXIMATION I Consider te complex extension f z, were z = x + iy, of f x and compute te complex Taylor series expansion Take imaginary part and divide by f x + i = f + i f i 3 +O 4 f = I f x + i + 6 f +O 3 = = I f x + i Note tat te sceme is second order accurate were is conservation of complexity? Te metod doesn t suffer from cancellation errors, is easy to implement and quite useful a J N Lyness and C BMoler, Numerical differentiation of analytical functions, SIAM J Numer Anal 4, -, 967 GENERAL IDEA include in te finite difference formula not only te function values, but also te values of te FUNCTION DERIVATIVE at te adacent nodes, eg: b f + f + b f + p= α p f +p = ε Construct te Taylor table using te following expansions: f + = f + f + f f f + = f + f + f f iv + 5 f v + 6 f iv f v + NOTE need an expansion for te derivative and a iger order expansion for te function more coefficient to determine MATH745 Fall 5 3 MATH745 Fall 5 4 PADÉ APPROXIMATION II Te Taylor table f f f f b f b b b b 3 6 b 4 4 f b f + b b b b 3 6 b 4 4 a f a a a a 3 6 a 4 4 a 5 a f a a f + a a a a 3 6 a 4 4 a 5 Te algebraic system: / / / / 3 /6 3 /6 3 /6 3 /6 4 /4 4 /4 b b a a a f iv = = b b a a a f v /4 /4 = 3/4 3/4 PADÉ APPROXIMATION III Te Padé approximation: = 3 f + f 4 Leading order error 4 3 f v FOURTH ORDER ACCURATE Te approximation is NONLOCAL, in tat it requires derivatives at te adacent nodes wic are also unknowns; Tus all derivatives must be determined at once via te solution of te following algebraic system /4 /4 + = 3 4 f + f
7 MATH745 Fall 5 5 MATH745 Fall 5 6 PADÉ APPROXIMATION IV Closing te system at ENDPOINTS were neigbors are not available use a lower order one sided ie, forward or backward finite difference formula Te vector of derivatives can tus be obtained via solution of te following algebraic system were Bf = 3 Af = f = 3 B Af B is a tri diagonal matrix wit b i,i = and b i,i = b i,i+ = 4, i =,,N A is a second order accurate differentiation matrix MODIFIED WAVENUMBER ANALYSIS I How do finite differences perform at different WAVELENGTHS? Finite Difference formulae applied to THE FOURIER MODE f x = e ikx wit te exact derivative f x = ike ikx Central Difference formula: = f + f = eikx + e ikx were te modified wavenumber k sink = eik e ik e ikx = i sink f = ik f, Comparison of te modified wavenumber k wit te actual wavenumber k sows ow numerical differentiation errors affect different Fourier components of a given function MATH745 Fall 5 7 MODIFIED WAVENUMBER ANALYSIS II Fourt-order central difference formula = f f + 8 f + f ] [ 4 = i 3 sink 6 sink = e ik e ik f e ik e ik f 3 f = ik f were te modified wavenumber k [ 4 3 sink 6 sink] Fourt order Padé sceme: were Tus: = ik e ikx + = ik e ik f and = 3 4 f + f, = ik e ikx = ik e ik f ik 4 eik e ik f = 3 e ik e ik f 4 ik + cosk f = i 3 sink f = k 3sink + cosk
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 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 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 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 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 informationFall 2014 MAT 375 Numerical Methods. Numerical Differentiation (Chapter 9)
Fall 2014 MAT 375 Numerical Metods (Capter 9) Idea: Definition of te derivative at x Obviuos approximation: f (x) = lim 0 f (x + ) f (x) f (x) f (x + ) f (x) forward-difference formula? ow good is tis
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 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 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 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 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 informationNotes on Multigrid Methods
Notes on Multigrid Metods Qingai Zang April, 17 Motivation of multigrids. Te convergence rates of classical iterative metod depend on te grid spacing, or problem size. In contrast, convergence rates of
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 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 informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
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 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 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 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 informationThe Laplace equation, cylindrically or spherically symmetric case
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,
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 informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
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 informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
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 informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
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 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 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 informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (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 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 informationDigital Filter Structures
Digital Filter Structures Te convolution sum description of an LTI discrete-time system can, in principle, be used to implement te system For an IIR finite-dimensional system tis approac is not practical
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 information7 Semiparametric Methods and Partially Linear Regression
7 Semiparametric Metods and Partially Linear Regression 7. Overview A model is called semiparametric if it is described by and were is nite-dimensional (e.g. parametric) and is in nite-dimensional (nonparametric).
More informationA Finite Element Primer
A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.
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 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 informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
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 informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
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 informationPOLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY
APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky
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 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 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 informationNON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE
Global and Stocastic Analysis Vol. 4 No. 1, January 2017, 1-10 NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE K. PHANEENDRA AND E. SIVA PRASAD Abstract.
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 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 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 information1 (10) 2 (10) 3 (10) 4 (10) 5 (10) 6 (10) Total (60)
First Name: OSU Number: Last Name: Signature: OKLAHOMA STATE UNIVERSITY Department of Matematics MATH 2144 (Calculus I) Instructor: Dr. Matias Sculze MIDTERM 1 September 17, 2008 Duration: 50 minutes No
More informationMATH 155A FALL 13 PRACTICE MIDTERM 1 SOLUTIONS. needs to be non-zero, thus x 1. Also 1 +
MATH 55A FALL 3 PRACTICE MIDTERM SOLUTIONS Question Find te domain of te following functions (a) f(x) = x3 5 x +x 6 (b) g(x) = x+ + x+ (c) f(x) = 5 x + x 0 (a) We need x + x 6 = (x + 3)(x ) 0 Hence Dom(f)
More informationGeneric maximum nullity of a graph
Generic maximum nullity of a grap Leslie Hogben Bryan Sader Marc 5, 2008 Abstract For a grap G of order n, te maximum nullity of G is defined to be te largest possible nullity over all real symmetric n
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 informationINTRODUCTION TO CALCULUS LIMITS
Calculus can be divided into two ke areas: INTRODUCTION TO CALCULUS Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and
More informationCONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION
CONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION SÖREN BARTELS AND ANDREAS PROHL Abstract. Te Landau-Lifsitz-Gilbert equation describes dynamics of ferromagnetism,
More informationExercise 19 - OLD EXAM, FDTD
Exercise 19 - OLD EXAM, FDTD A 1D wave propagation may be considered by te coupled differential equations u x + a v t v x + b u t a) 2 points: Derive te decoupled differential equation and give c in terms
More informationMath Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim
Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z
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 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 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 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 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 informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
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 informationChapter 1. Density Estimation
Capter 1 Density Estimation Let X 1, X,..., X n be observations from a density f X x. Te aim is to use only tis data to obtain an estimate ˆf X x of f X x. Properties of f f X x x, Parametric metods f
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 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 informationEfficient algorithms for for clone items detection
Efficient algoritms for for clone items detection Raoul Medina, Caroline Noyer, and Olivier Raynaud Raoul Medina, Caroline Noyer and Olivier Raynaud LIMOS - Université Blaise Pascal, Campus universitaire
More informationSin, Cos and All That
Sin, Cos and All Tat James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 9, 2017 Outline Sin, Cos and all tat! A New Power Rule Derivatives
More information5.1 introduction problem : Given a function f(x), find a polynomial approximation p n (x).
capter 5 : polynomial approximation and interpolation 5 introduction problem : Given a function f(x), find a polynomial approximation p n (x) Z b Z application : f(x)dx b p n(x)dx, a a one solution : Te
More informationIntegral Calculus, dealing with areas and volumes, and approximate areas under and between curves.
Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
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 informationInvestigation of Tangent Polynomials with a Computer Algebra System The AMATYC Review, Vol. 14, No. 1, Fall 1992, pp
Investigation of Tangent Polynomials wit a Computer Algebra System Te AMATYC Review, Vol. 14, No. 1, Fall 199, pp. -7. Jon H. Matews California State University Fullerton Fullerton, CA 9834 By Russell
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 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 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 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 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 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 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 informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationMTH 119 Pre Calculus I Essex County College Division of Mathematics Sample Review Questions 1 Created April 17, 2007
MTH 9 Pre Calculus I Essex County College Division of Matematics Sample Review Questions Created April 7, 007 At Essex County College you sould be prepared to sow all work clearly and in order, ending
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 informationFunction Composition and Chain Rules
Function Composition an Cain Rules James K. Peterson Department of Biological Sciences an Department of Matematical Sciences Clemson University November 2, 2018 Outline Function Composition an Continuity
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 informationChapter 10. Function approximation Function approximation. The Lebesgue space L 2 (I)
Capter 1 Function approximation We ave studied metods for computing solutions to algebraic equations in te form of real numbers or finite dimensional vectors of real numbers. In contrast, solutions to
More informationThese error are made from replacing an infinite process by finite one.
Introduction :- Tis course examines problems tat can be solved by metods of approximation, tecniques we call numerical metods. We begin by considering some of te matematical and computational topics tat
More informationMath 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0
3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,
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 informationTHE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225
THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationHomework 1 Due: Wednesday, September 28, 2016
0-704 Information Processing and Learning Fall 06 Homework Due: Wednesday, September 8, 06 Notes: For positive integers k, [k] := {,..., k} denotes te set of te first k positive integers. Wen p and Y q
More informationOn convergence of the immersed boundary method for elliptic interface problems
On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many
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 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