MA2264 -NUMERICAL METHODS UNIT V : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL. By Dr.T.Kulandaivel Department of Applied Mathematics SVCE

Size: px
Start display at page:

Download "MA2264 -NUMERICAL METHODS UNIT V : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL. By Dr.T.Kulandaivel Department of Applied Mathematics SVCE"

Transcription

1 MA64 -NUMERICAL METHODS UNIT V : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS B Dr.T.Kulandaivel Department of Applied Matematics SVCE

2

3 Numerical ordinar differential equations is te part of numerical analsis wic studies te numerical solution of ordinar differential equations (ODEs). It is ver important for industrial and practical engineering applications, Mecanics, Cemistr and Psics..

4 It is a numerical metods for solving ordinar differential equations (O.D.Es) wit a given conditions. d f (, ), ( ) d Talor series epansion of te function f() as te following form:. ( ) ( ) () ( ) () ()!!... Assume = ( ) () () () ()!!!...

5 Hint: - Wen te step size becomes small, te numerical solution becomes nearl te same as te eact solution. - computing approimate values of te solution f() at te points: = +, = +, = +,..

6 Solve te boundar value problem = - at =.,.4 given tat () = Given tat () &,.. Evaluate te derivatives at te point of epansion zero () ( ) ( ), () () ( ) ( ), () () ( ) ()

7 Talor epansion given b: ( ) () () ( ) ()!!!.... (.) (.) (.) ( )!! (.) ()! ( ) Here.,.8867 ( ) ( ) (.).8867 " ( ) ' ( ) ' ( ) " ( )

8 () () ( ) ( ) ( )!!!. (.) (.) (.4).8867 (.8867 ) (.8867 ) (.8867 )..!!!.67...

9 It is a first-order numerical procedure for solving ordinar differential equations (O.D.E) wit a given conditions. d f d (, ), ( ) Te general Euler formula for solving boundar value problem is: n n f( n, n )

10 Assume = n+ n Computing approimate values of te solution f() at te points: = +, = +, = +,. ( ), ( ), ( ),...

11 Solve te following boundar value problem at =.,.4,.6 and.8 given tat () =. d d Given tat., d d.4, ;,.6, 4,.8 General Euler formula is given b:

12 [ ], n n n n For n =...( ) [ ].( ). For n = [ ]..((.).).44

13 For n = [ ].44. ((.4).44 ).6

14 For n = [ 4 ].6.((.6).6).85 For n = [ 4 4 ].85.((.8).85).4

15 We can put te results in a table: n n n n

16 Tis metod tries to improve te error of Euler metod. Te general modified Euler formula for solving boundar value problem is:,,, n n n n n n f f

17 Using modified Euler s metod, find (.), if = /, () = wit =., tree decimal places are required Modified Eulr ' s formula : n n f n, n f,, n n

18 ]] (.5 )[ [.5, (.) ] / [.,..,, f f f f.9548 ] (.)[.9548 /.5 ] (.5 ) (.)[.5 [.5,.5 ] (.) ]] (.5 )[ [.5, (.) f f

19 Te accurac of tis metod is equivalent to five terms of Talor s series. Te general Rung- kutta formula for solving boundar value problem is: n k (k k ) k 4 n. 6 Net Te Main Menu ا Previous

20 were k f( n, n ) k f ( ),( ) n n k ) k f ( ),( ) n n k k f ( ),( k 4 n n )

21 Find (.) and (.), if \ =, ()=, take =., four decimal places are required To get = (.) k f(, ).( ) n n. k f ( ),( ).(..5) n n k.5 k f ( n ),(n ).(.5.5) k.5

22 (, ),(,k ).(.5.). 5 k f 4 n n k (k k ) k 4 6. (.5.5 ) (.).5

23 To get = (.) k f(, ).(.5.) n n.5 k f ( n ),(n ).(.5.5) k.6 k f ( n ),(n ).(..5) k.6

24 (,),(,k ).(.45.). k f 4 n n k (k k ) k (.6 6.6). (.) =.45

25 Milne s Metod It is a multi-step metod were we assume tat te solution to te given IVP is known at te past four equall spaced point t, t, t and t.

26 In general, Milne s predictorpredictorcorrector pair can be written as 4 P : n n ( n n n ) C : n n ( n 4 n n )

27 Using Milne s Metod, find (4.4) given 5 + -=, Given (4)=, (4.) =.49 and (4.) =.97, (4.) =.4, find (.4).

28 .4..97,.49,, , 4., 4., 4,, 5 : 4 Given Solution ) 49. ( 45. ). 4 5 ( ) 4. ( ). 4 5 ( ) 97. ( ). 4 5 ( ) 49. ( 5 ' ' '

29 B Milne s predictor formula 4 ' ' ' 4, P 4(.) [(.49).467 (.45].897 (.897) (4.4) ' 4 4

30 4, C ' 4 ' ' [.467 4(.45).47]

31 Eample: Using Milne s metod, find (.4) given = +, ()=. Use Talor series to get te values of (.), (.), (.). Given ' Here,.,.,., 4.4,.

32 ' " ' " "' ' ' " ' ) ( : ') ( '' ' '' "' : ' ' " : '.67 (.) () 6. ()..... ) (! ) (! ) (! ) ( ) (

33 )(.587 ( ) (.)( ) (.67 ) (.)(.67 ' ' '' ' 6.4 ) ( 4.87 ' " ' " '''

34 !!! (.).767. (4.87 ) ( 6.4) 6

35 (.767)(.885).767 (.)(.885).885 (.767) (.)(.767) ' ' '' ' ) ( ' " ' " '''

36 .5 (.) ).6875 (8 6. ) ( !!! ] [ 4, ln '.,,,,, 4, 4 formula espredictor BMi wewillfind Nowknowing p

37 ,.776 (.5) (.)(.5).885 (.767) (.)(.767).587 (.67) (.)(.67) ' ' ' Now Now,.897 (.776)].885 [(.587) 4(.) 4, p

38 , ' ln, sin (.897) (.4)(.897) , formula s corrector e Mi g U p.8698 (.4) 4.996] 4(.776) [ ] 4 [ 4, 4, p c

39 Adam s predictor and corrector mtod mtod:: It is anoter predictor predictor--corrector metod, were we use te fact tat te solution to te given initial value problem is known at past four equall spaced points n, n, n, n Te task is to compute te values of at n

40 Let us consider te differential equation d d f (, ) In general, AdamAdam-Moulton predictor corrector pair can be written as P : n n 55 n 59 n 7 n 9 n 4 C : n n 9 n 9 n 5 n n 4

41 In general, AdamAdam-Basfort predictor -corrector formula can be written as P : n n 55 n 59 n 7 n 9 n 4 C : n n 9 n 9 n 5 n n 4

42 Problem: Using Adam s Basfort metod find (4.4) given 5 + =,Y94)= (4.)=.49,(4.)=.97, (4.)=.4. Solution: Given 4, 5,, Let 4.,..49, 4.,.97, 4.,.4.

43 Adam' s n, p 4, p putting n Now, 5 predictor n [55 n 4, we ave [ formula , is n 7 7, 5 5 n n ] ] ().48,

44 U sin g tese 4, p.4 values in (), we get. [55(.45) 59(.467) 4 7(.48) 9(.5)].4 (4.4)

45 Adam' s corrector n, c 4, c putting n Now, 4 n [9 n 4, we ave [ formula is n 5 5 n n ] ] ()

46 Eqn() becomes. 4, c.4 4 [9(.47) 9(.45) 5(.467).48]..4 (.669).87 4 (4.4).87

47 Using AdamAdam-Moulton predictor predictor-corrector metod, find te solution of te initial value problem d t, dt () at t =., taking =.. Compare it wit te analtical solution.

48 Solution In order to use Adam s P-C metod, we require te solution of te given differential equation at te past four equall spaced points, use R-K metod of 4 t points, for wic we order wic is self starting.

49 Tus taking t =, =, =., we compute k =., k =.8, 8, k =.98, 98, k4 =.596, 596, and get (k k k k4 ).859 6

50 Now, we take t =.4, =.468, 468, =., and compute k =.66, 66, k =.697, 697,k =.76, 76, k4 =.757 to get (.6) (k k k k4 )

51 Now, we use Adam s PP-C pair to calculate ( (.8) and (.) as follows: P : n n 55 n 59 n 7 n 9 n 4 C : n n 9 n 9 n 5 n n 4 Tus p 4 ()

52 From te given differential equation, we ave t. Terefore, t. t.7859 t.8 t.7779

53 Hence, from Eq. ( (), we get. (.8) p 4.45 Now to obtain te corrector value of at t =.8, we use (.8) c 4 c ()

54 But, p 9 4 9( 4 t4 ) 9[.45 (.8) ].759 Terefore,. 4 (.8) c.44 () Proceeding similarl, we get (.) p 5 p

55 Noting tat 4 4 t.74, 4 we calculate p Now, te corrector formula for computing 5 is given b (.) c 5 c (4)

56 From te given differential equation, we ave t. Terefore, t. t.7859 t.8 t.7779

57 Hence, from Eq. ( (), we get. (.8) p 4.45 Now to obtain te corrector value of at t =.8, we use (.8) c 4 c ()

58 But, p 9 4 9( 4 t4 ) 9[.45 (.8) ].759 Terefore,. 4 (.8) c.44 () Proceeding similarl, we get (.) p 5 p

59 Noting tat 4 4 t.74, 4 we calculate p Now, te corrector formula for computing 5 is given b (.) c 5 c (4)

60 9 9 t 5 But, p Tus, finall we get. 5 (.)

61

62 Te numerical solution of a differential equation can be sown to converge to its eact solution, if te step size is ver small.

63 Te numerical solution of a differential equation is said to be stable if te error do not grow eponentiall as we compute from one step to anoter.

64 Stabilit consideration are ver important in finding te numerical solutions of te differential equations eiter b single-stepstep metods or b using multi-step metods.

Chapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods

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

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1

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

Initial Value Problems for. Ordinary Differential Equations

Initial Value Problems for. Ordinary Differential Equations Initial Value Problems for Ordinar Differential Equations INTRODUCTION Equations which are composed of an unnown function and its derivatives are called differential equations. It becomes an initial value

More information

158 Calculus and Structures

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

(5.5) Multistep Methods

(5.5) Multistep Methods (5.5) Mulstep Metods Consider te inial-value problem for te ordinary differenal equaon: y t f t, y, a t b, y a. Let y t be te unique soluon. In Secons 5., 5. and 5.4, one-step numerical metods: Euler Metod,

More information

Ordinary Differential Equations n

Ordinary Differential Equations n Numerical Analsis MTH63 Ordinar Differential Equations Introduction Talor Series Euler Method Runge-Kutta Method Predictor Corrector Method Introduction Man problems in science and engineering when formulated

More information

3.1 Extreme Values of a Function

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

232 Calculus and Structures

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

Solving Second Order Linear Dirichlet and Neumann Boundary Value Problems by Block Method

Solving Second Order Linear Dirichlet and Neumann Boundary Value Problems by Block Method IAENG International Journal o Applied Matematics 43: IJAM_43 04 Solving Second Order Linear Diriclet and Neumann Boundar Value Problems b Block Metod Zanaria Abdul Majid Mod Mugti Hasni and Norazak Senu

More information

University Mathematics 2

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

Computational Methods CMSC/AMSC/MAPL 460. Ordinary differential equations

Computational Methods CMSC/AMSC/MAPL 460. Ordinary differential equations Computational Methods CMSC/AMSC/MAPL 460 Ordinar differential equations Ramani Duraiswami, Dept. of Computer Science Several slides adapted from Prof. ERIC SANDT, TAMU ODE: Previous class Standard form

More information

Chapter XI. Solution of Ordinary Differential Equations

Chapter XI. Solution of Ordinary Differential Equations dy d 7 Capter XI Solution of Ordinary Differential Equations Several approaces are described for te solution of ODE (ordinary differential equations). Tese include: Te Taylor series Metod Te Euler Metod

More information

Runge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t))

Runge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t)) Runge-Kutta metods Wit orders of Taylor metods yet witout derivatives of f (t, y(t)) First order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous

More information

Computational Methods CMSC/AMSC/MAPL 460. Ordinary differential equations

Computational Methods CMSC/AMSC/MAPL 460. Ordinary differential equations Computational Methods CMSC/AMSC/MAPL 460 Ordinar differential equations Ramani Duraiswami, Dept. of Computer Science Several slides adapted from Prof. ERIC SANDT, TAMU ODE: Previous class Applications

More information

NUMERICAL DIFFERENTIAL 1

NUMERICAL DIFFERENTIAL 1 NUMERICAL DIFFERENTIAL Ruge-Kutta Metods Ruge-Kutta metods are ver popular ecause o teir good eiciec; ad are used i most computer programs or dieretial equatios. Te are sigle-step metods as te Euler metods.

More information

Two Step Hybrid Block Method with Two Generalized Off-step Points for Solving Second Ordinary Order Differential Equations Directly

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

SKP Engineering College

SKP Engineering College S.K.P. Engineering College, Tiruvannamalai SKP Engineering College Tiruvannamalai 666 A Course Material on By S.Danasekar Assistant Professor Department of Matematics S.K.P. Engineering College, Tiruvannamalai

More information

1 ode.mcd. Find solution to ODE dy/dx=f(x,y). Instructor: Nam Sun Wang

1 ode.mcd. Find solution to ODE dy/dx=f(x,y). Instructor: Nam Sun Wang Fin solution to ODE /=f(). Instructor: Nam Sun Wang oe.mc Backgroun. Wen a sstem canges wit time or wit location, a set of ifferential equations tat contains erivative terms "/" escribe suc a namic sstem.

More information

1 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

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

1.5 Functions and Their Rates of Change

1.5 Functions and Their Rates of Change 66_cpp-75.qd /6/8 4:8 PM Page 56 56 CHAPTER Introduction to Functions and Graps.5 Functions and Teir Rates of Cange Identif were a function is increasing or decreasing Use interval notation Use and interpret

More information

Logarithmic functions

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

Continuity and Differentiability Worksheet

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

The Verlet Algorithm for Molecular Dynamics Simulations

The Verlet Algorithm for Molecular Dynamics Simulations Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical

More information

Applied Linear Statistical Models. Simultaneous Inference Topics. Simultaneous Estimation of β 0 and β 1 Issues. Simultaneous Inference. Dr.

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

Higher Derivatives. Differentiable Functions

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

Sixth and Fourth Order Compact Finite Difference Schemes for Two and Three Dimension Poisson Equation with Two Methods to derive These Schemes

Sixth and Fourth Order Compact Finite Difference Schemes for Two and Three Dimension Poisson Equation with Two Methods to derive These Schemes Basra Journal of Scienec (A) Vol.(),-0, 00 Sit and Fourt Order Compact Finite Difference Scemes for Two and Tree Dimension Poisson Equation wit Two Metods to derive Tese Scemes Akil J. Harfas Huda A. Jalob

More information

Section 3.1: Derivatives of Polynomials and Exponential Functions

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

Exponentials and Logarithms Review Part 2: Exponentials

Exponentials and Logarithms Review Part 2: Exponentials Eponentials and Logaritms Review Part : Eponentials Notice te difference etween te functions: g( ) and f ( ) In te function g( ), te variale is te ase and te eponent is a constant. Tis is called a power

More information

ORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD

ORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD Numercal Analss or Engneers German Jordanan Unverst ORDINARY DIFFERENTIAL EQUATIONS We wll eplore several metods o solvng rst order ordnar derental equatons (ODEs and we wll sow ow tese metods can be appled

More information

Fourth Order RK-Method

Fourth Order RK-Method Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)

More information

d y f f dy Numerical Solution of Ordinary Differential Equations Consider the 1 st order ordinary differential equation (ODE) . dx

d y f f dy Numerical Solution of Ordinary Differential Equations Consider the 1 st order ordinary differential equation (ODE) . dx umerical Solutio o Ordiar Dieretial Equatios Cosider te st order ordiar dieretial equatio ODE d. d Te iitial coditio ca be tae as. Te we could use a Talor series about ad obtai te complete solutio or......!!!

More information

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.

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

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT

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

WYSE Academic Challenge 2004 Sectional Mathematics Solution Set

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

1 2 x Solution. The function f x is only defined when x 0, so we will assume that x 0 for the remainder of the solution. f x. f x h f x.

1 2 x Solution. The function f x is only defined when x 0, so we will assume that x 0 for the remainder of the solution. f x. f x h f x. Problem. Let f x x. Using te definition of te derivative prove tat f x x Solution. Te function f x is only defined wen x 0, so we will assume tat x 0 for te remainder of te solution. By te definition of

More information

Pre-Calculus Review Preemptive Strike

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

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019

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

2.3. Applying Newton s Laws of Motion. Objects in Equilibrium

2.3. Applying Newton s Laws of Motion. Objects in Equilibrium Appling Newton s Laws of Motion As ou read in Section 2.2, Newton s laws of motion describe ow objects move as a result of different forces. In tis section, ou will appl Newton s laws to objects subjected

More information

Recall from our discussion of continuity in lecture a function is continuous at a point x = a if and only if

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

Numerical Solution of the Initial Value Problem of the Ordinary Differential Equation with Singular Point by Multi-Step Integration Method

Numerical Solution of the Initial Value Problem of the Ordinary Differential Equation with Singular Point by Multi-Step Integration Method Nuerical o te Initial Value Proble o te Ordinar Dierential Equation wit Singular Point b MultiStep Integration Metod KANOKNAPA ERAWUN Kase Bundit Universit Suanluang Bangkok Tailand Abstract. In [] M.

More information

5.1 We will begin this section with the definition of a rational expression. We

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

AMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.

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

lim 1 lim 4 Precalculus Notes: Unit 10 Concepts of Calculus

lim 1 lim 4 Precalculus Notes: Unit 10 Concepts of Calculus Syllabus Objectives: 1.1 Te student will understand and apply te concept of te limit of a function at given values of te domain. 1. Te student will find te limit of a function at given values of te domain.

More information

Chapter 4: Numerical Methods for Common Mathematical Problems

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

1 + t5 dt with respect to x. du = 2. dg du = f(u). du dx. dg dx = dg. du du. dg du. dx = 4x3. - page 1 -

1 + t5 dt with respect to x. du = 2. dg du = f(u). du dx. dg dx = dg. du du. dg du. dx = 4x3. - page 1 - Eercise. Find te derivative of g( 3 + t5 dt wit respect to. Solution: Te integrand is f(t + t 5. By FTC, f( + 5. Eercise. Find te derivative of e t2 dt wit respect to. Solution: Te integrand is f(t e t2.

More information

Differential equations. Differential equations

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

Seepage Analysis through Earth Dam Based on Finite Difference Method

Seepage Analysis through Earth Dam Based on Finite Difference Method J. Basic. Appl. Sci. Res., (11)111-1, 1 1, TetRoad Publication ISSN -44 Journal of Basic and Applied Scientific Researc www.tetroad.com Seepage Analysis troug Eart Dam Based on Finite Difference Metod

More information

Click here to see an animation of the derivative

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

An Accurate Self-Starting Initial Value Solvers for. Second Order Ordinary Differential Equations

An Accurate Self-Starting Initial Value Solvers for. Second Order Ordinary Differential Equations International Journal of Contemporar Matematical Sciences Vol. 9, 04, no. 5, 77-76 HIKARI Ltd, www.m-iari.com ttp://dx.doi.org/0.988/icms.04.4554 An Accurate Self-Starting Initial Value Solvers for Second

More information

SFU UBC UNBC Uvic Calculus Challenge Examination June 5, 2008, 12:00 15:00

SFU UBC UNBC Uvic Calculus Challenge Examination June 5, 2008, 12:00 15:00 SFU UBC UNBC Uvic Calculus Callenge Eamination June 5, 008, :00 5:00 Host: SIMON FRASER UNIVERSITY First Name: Last Name: Scool: Student signature INSTRUCTIONS Sow all your work Full marks are given only

More information

Derivation Of The Schwarzschild Radius Without General Relativity

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

MAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016

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

(4.2) -Richardson Extrapolation

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

ETNA Kent State University

ETNA 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

Chapter 2. Limits and Continuity 16( ) 16( 9) = = 001. Section 2.1 Rates of Change and Limits (pp ) Quick Review 2.1

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

11.6 DIRECTIONAL DERIVATIVES AND THE GRADIENT VECTOR

11.6 DIRECTIONAL DERIVATIVES AND THE GRADIENT VECTOR SECTION 11.6 DIRECTIONAL DERIVATIVES AND THE GRADIENT VECTOR 633 wit speed v o along te same line from te opposite direction toward te source, ten te frequenc of te sound eard b te observer is were c is

More information

MATH1151 Calculus Test S1 v2a

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

136 Calculus and Structures

136 Calculus and Structures 6 Calculus and Structures CHAPTER 9 THE DERIVATIVE Calculus and Structures 7 Coprigt Capter 9 THE DERIVATIVE 9. THE ORIGINS OF CALCULUS All great discoveries ave teir origins wit an important problem tat

More information

(a 1 m. a n m = < a 1/N n

(a 1 m. a n m = < a 1/N n Notes on a an log a Mat 9 Fall 2004 Here is an approac to te eponential an logaritmic functions wic avois any use of integral calculus We use witout proof te eistence of certain limits an assume tat certain

More information

10 Derivatives ( )

10 Derivatives ( ) Instructor: Micael Medvinsky 0 Derivatives (.6-.8) Te tangent line to te curve yf() at te point (a,f(a)) is te line l m + b troug tis point wit slope Alternatively one can epress te slope as f f a m lim

More information

Improved Rotated Finite Difference Method for Solving Fractional Elliptic Partial Differential Equations

Improved Rotated Finite Difference Method for Solving Fractional Elliptic Partial Differential Equations American Scientific Researc Journal for Engineering, Tecnolog, and Sciences (ASRJETS) ISSN (Print) 33-44, ISSN (Online) 33-44 Global Societ of Scientific Researc and Researcers ttp://asrjetsjournal.org/

More information

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

Integral Calculus, dealing with areas and volumes, and approximate areas under and between curves.

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

Polynomial Interpolation

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

Order of Accuracy. ũ h u Ch p, (1)

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

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001

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

Section 15.6 Directional Derivatives and the Gradient Vector

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

Printed Name: Section #: Instructor:

Printed Name: Section #: Instructor: Printed Name: Section #: Instructor: Please do not ask questions during tis eam. 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 information

Some Review Problems for First Midterm Mathematics 1300, Calculus 1

Some 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

Differentiation. Area of study Unit 2 Calculus

Differentiation. Area of study Unit 2 Calculus Differentiation 8VCE VCEco Area of stud Unit Calculus coverage In tis ca 8A 8B 8C 8D 8E 8F capter Introduction to limits Limits of discontinuous, rational and brid functions Differentiation using first

More information

Math 242: Principles of Analysis Fall 2016 Homework 7 Part B Solutions

Math 242: Principles of Analysis Fall 2016 Homework 7 Part B Solutions Mat 22: Principles of Analysis Fall 206 Homework 7 Part B Solutions. Sow tat f(x) = x 2 is not uniformly continuous on R. Solution. Te equation is equivalent to f(x) = 0 were f(x) = x 2 sin(x) 3. Since

More information

Test 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 =

Test 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 = Test Review Find te determinant of te matrix below using (a cofactor expansion and (b row reduction Answer: (a det + = (b Observe R R R R R R R R R Ten det B = (((det Hence det Use Cramer s rule to solve:

More information

How to Find the Derivative of a Function: Calculus 1

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

Continuity. Example 1

Continuity. Example 1 Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)

More information

Polynomial Interpolation

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

Numerical Differentiation

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 information

LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION

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

CHAPTER 2 MODELING OF THREE-TANK SYSTEM

CHAPTER 2 MODELING OF THREE-TANK SYSTEM 7 CHAPTER MODELING OF THREE-TANK SYSTEM. INTRODUCTION Te interacting tree-tank system is a typical example of a nonlinear MIMO system. Heiming and Lunze (999) ave regarded treetank system as a bencmark

More information

Math Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim

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

A = h w (1) Error Analysis Physics 141

A = 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 information

Implicit-explicit variational integration of highly oscillatory problems

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

Chapter 2 Limits and Continuity

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

= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)

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

The Derivative as a Function

The Derivative as a Function Section 2.2 Te Derivative as a Function 200 Kiryl Tsiscanka Te Derivative as a Function DEFINITION: Te derivative of a function f at a number a, denoted by f (a), is if tis limit exists. f (a) f(a + )

More information

Physically Based Modeling: Principles and Practice Implicit Methods for Differential Equations

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

An Approximation to the Solution of the Brusselator System by Adomian Decomposition Method and Comparing the Results with Runge-Kutta Method

An Approximation to the Solution of the Brusselator System by Adomian Decomposition Method and Comparing the Results with Runge-Kutta Method Int. J. Contemp. Mat. Sciences, Vol. 2, 27, no. 2, 983-989 An Approximation to te Solution of te Brusselator System by Adomian Decomposition Metod and Comparing te Results wit Runge-Kutta Metod J. Biazar

More information

Taylor Series and the Mean Value Theorem of Derivatives

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

4.2 - Richardson Extrapolation

4.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

Solutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014

Solutions 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

multistep methods Last modified: November 28, 2017 Recall that we are interested in the numerical solution of the initial value problem (IVP):

multistep methods Last modified: November 28, 2017 Recall that we are interested in the numerical solution of the initial value problem (IVP): MATH 351 Fall 217 multistep methods http://www.phys.uconn.edu/ rozman/courses/m351_17f/ Last modified: November 28, 217 Recall that we are interested in the numerical solution of the initial value problem

More information

Copyright c 2008 Kevin Long

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

Math 128A Spring 2003 Week 11 Solutions Burden & Faires 5.6: 1b, 3b, 7, 9, 12 Burden & Faires 5.7: 1b, 3b, 5 Burden & Faires 5.

Math 128A Spring 2003 Week 11 Solutions Burden & Faires 5.6: 1b, 3b, 7, 9, 12 Burden & Faires 5.7: 1b, 3b, 5 Burden & Faires 5. Math 128A Spring 2003 Week 11 Solutions Burden & Faires 5.6: 1b, 3b, 7, 9, 12 Burden & Faires 5.7: 1b, 3b, 5 Burden & Faires 5.8: 1b, 3b, 4 Burden & Faires 5.6. Multistep Methods 1. Use all the Adams-Bashforth

More information

MATH 1020 TEST 2 VERSION A FALL 2014 ANSWER KEY. Printed Name: Section #: Instructor:

MATH 1020 TEST 2 VERSION A FALL 2014 ANSWER KEY. Printed Name: Section #: Instructor: ANSWER KEY Printed Name: Section #: Instructor: Please do not ask questions during tis eam. If you consider a question to be ambiguous, state your assumptions in te margin and do te best you can to provide

More information

Combining functions: algebraic methods

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

ESTIMATION OF TIME-DOMAIN FREQUENCY STABILITY BASED ON PHASE NOISE MEASUREMENT

ESTIMATION OF TIME-DOMAIN FREQUENCY STABILITY BASED ON PHASE NOISE MEASUREMENT ESTIMATION OF TIME-DOMAIN FREQUENCY STABILITY BASED ON PHASE NOISE MEASUREMENT P. C. Cang, H. M. Peng, and S. Y. Lin National Standard Time & Frequenc Laborator, TL, Taiwan, Lane 55, Min-Tsu Road, Sec.

More information

A STATIC PDE APPROACH FOR MULTI-DIMENSIONAL EXTRAPOLATION USING FAST SWEEPING METHODS

A STATIC PDE APPROACH FOR MULTI-DIMENSIONAL EXTRAPOLATION USING FAST SWEEPING METHODS A STATIC PDE APPROACH FOR MULTI-DIMENSIONAL EXTRAPOLATION USING FAST SWEEPING METHODS TARIQ ASLAM, SONGTING LUO, AND HONGKAI ZHAO Abstract. A static Partial Differential Equation (PDE) approac is presented

More information

Control Errors in CFD!

Control Errors in CFD! Control Errors in CFD! Bernard Müller Department of Scientic Computing, Information Tecnolog, Uppsala Universit, Box 337, S-751 5 Uppsala, Sweden, E-mail bernd@tdb.uu.se, URL ttp://www.tdb.uu.se/ bernd/

More information

Calculus I - Spring 2014

Calculus I - Spring 2014 NAME: Calculus I - Spring 04 Midterm Exam I, Marc 5, 04 In all non-multiple coice problems you are required to sow all your work and provide te necessary explanations everywere to get full credit. In all

More information

Function Composition and Chain Rules

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

Lecture 10: Carnot theorem

Lecture 10: Carnot theorem ecture 0: Carnot teorem Feb 7, 005 Equivalence of Kelvin and Clausius formulations ast time we learned tat te Second aw can be formulated in two ways. e Kelvin formulation: No process is possible wose

More information