Chapter XI. Solution of Ordinary Differential Equations
|
|
- Egbert McGee
- 5 years ago
- Views:
Transcription
1 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 Te Improved Euler Metod Te Runge-Kutta Metods. We will begin our presentation on te numerical solution of a first-order ODE and later on etend te concepts to any order ODE.. Problem Statement Given = f (, y); y( ) = y 0 0 (.) were y( 0 ) = 0 is te initial condition needed to solve te problem. Tat is, y = y 0 at = 0. Determine y as a function of in eiter tabulated form or grap. Table 0. sows te tabulated solution required. Table. Tabulated results n- n- y y 0 y... y n- y n- Were i = i + and ) is given. 0. Te Taylor Series Metod Epand y() about 0 using Taylor series:
2 74 y ( 0) y ( ) = y ( 0) + y ( 0)( 0) + ( 0) +! ( 4) y ( 0) y ( ) +!! ( ) ! 4 Let 0 =, terefore: 4 y ( ) = y0 + f ( 0, y0 ) + f ( 0, y 0) + f ( 0, y 0) + f ( 0, y 0) +!! 4! Note tat f ( 0, y0) f ( 0, y0) df ( 0, y0) = d + dy y We will drop te ( 0,y 0 ) from ere on. Dividing by d we get! (.) Similarly one can derive: df d = f f + y f = F (.) d f d d f d " i d f i d " F F = + y f = F F F = + y f = F Fi Fi = + y f = F i (.4) Eample. Epress te solution of
3 dy d 75 = y wit y( 0) =. Solution. From Equations (.) and (.4) we can write: f = 0 y0 = F f = ( 0 y0) = y0 = F f = + ( 0 y0) = 0 y0 = F f = ( y ) = + + y = F " Substituting 0 = 0 and y 0 = - we get: Tat is: f =, f =, f =, f =,! 4 y ( ) = + +!!! 4! For small we can neglect terms iger tan 4 and ence given one can calculate y(). Setting 0 = and y 0 = y() we can repeat te above steps to calculate y() and so on.. Euler Metod In te Euler metod te first two terms of te Taylor series are used. y ( ) = y ( ) + f(, y) or y ( ) = y ( ) + f(, y) i+ i i i were i+ = i +, i = 0,,,! Euler metod can be grapically described as sown in Fig... Te slope at ( 0, y 0 ) is used to calculate y( ). Using and y( ) calculate f(, y ) and from tere calculate y( ) and so on.
4 76 Figure. Euler Metod..4 Improved Euler Metod In te improved Euler metod we use tree terms of te Taylor series: y ( i+ ) = y ( i) + f( i, yi) + f ( i, y i)! i = 0!,,,, Te derivative of f( i, y i ) can be calculated using forward difference as follows: were f (, y ) * y i+ is calculated using Euler metod as follows: i i * f ( i+, yi+ ) f ( i, yi) Hence: * yi+ = yi + f ( i, yi) y ( ) = y ( ) + f(, y) + i+ i i i ( ) f i+, yi + f ( i, yi) f ( i, yi)
5 77 wic can be placed in te form: [ ( )] y ( i+ ) = y ( i) + f ( i, y i) + f i+, y i + f ( i, y i) were = y ( i ) + ( k + k ) k = f (, y ) i (, ) k = f + y + k i i i Procedure: 5. Calculate k = f ( i, yi) 6. Calculate y * i = yi + k * 7. Calculate k = f ( i +, yi ) 8. Calculate yi+ = yi + k + k ( ) Te improved Euler metod can be described grapically as sown Figure.. Figure. Improved Euler Metod.
6 78 From Fig.. te average slope f dy d = * f (, y ) + f (, y ) 0 0 = y wit y( 0) =. is used in Euler metod to calculate te net value yi+ = yi + f, wic is te same approac derived in tis section for te Improved Euler metod. Eample. Calculate y(0.) and y(0.) using =0. for te following ODE Solution. Using f (, y) = y and te initial conditions 0 = 0, y0 = we follow te procedure described above.. Calculate k = f (, 0 ) = () 0 ( ) = *. Calculate y = y + k 0 = + 0.* = 09.. * Calculate k = f ( , y ) = f ( 0., 09. ) = Calculate y y 0 y 0. 0 k k 0. ( ) ( ) = Repeating te procedure to obtain y(0.).. Calculate k = f (., ) = * Calculate y = y + k = ( 0. 75) = * Calculate k = f ( + 0., y ) = f ( 0., ) = Calculate 0. y = y(.) 0 = y(.) 0 + ( k + k ) =. + ( ) = Hence:
7 79 y( 0) = y(.) 0 = y(.) 0 = Runge Kutta Metods Te main reason for te development of alternative metods to te Taylor series approac is to avoid te differentiations required by te series. Neiter Euler nor Improved metods require tat we differentiate any function. Te Runge Kutta Metods are basically generalizations of te Improved Euler Metod. For second order Runge-Kutta metod te following is assumed: y ( + ) y ( ) + ak + ak were k = f (, y) k = f ( + m, y + mk ) (.5) Te problem is now stated as follows: obtain te unknowns a, a, m suc tat y( + ) y ( + ) = y ( ) + y ( ) + y ( )! } obtained using te above equation matces te Taylor Series epansion of y( + ) truncated after te first derivative. Te Taylor Series epansion wit te terms truncated after te second derivative is: f y y f y f (, y ) (, ) (, ) y f (, y ) + + = + + y f (, y) f (, y) f (, y) + ( ) + y + ( y)! y y +! (.6) Since te derivatives in te Taylor series epansion are in terms of f(,y) we need to epand f ( + m, y + mk) as a series in terms of f(,y). Tis calls for te use of Taylor series epansion in two variables wic is given by: Hence we can write:
8 80 k = [ f (, y) + mf (, y) + mk f (, y) +!] (.7) y f (, y) Were f(, y) = and fy(, y) = [ y] [ y] y( + ) = y( ) + a f (, y) + a f (, y) + mf (, y) + mk f = y + a f + a f + mf + mff = y + ( a + a ) f + a m ( f + ff ) y ( ) = f df (, y) f f y ( ) = = + = + d y f f ff y f (, y) y Replacing k wit f(,y) and k wit te rigt-and side of Eq.(.7) in Equation (.5) we get: Notice tat we dropped te (,y) just for compactness sake. In te Taylor series epansion (Eq..6) replace: Hence y (.8) ( ff y) y ( + ) = y ( ) + f + f + (.9) Comparing (.8) and (.9) we get: a + a = am= Tese are two equations in tree unknowns and ence we ave te coice to select one of te variables. Coosing m=, terefore Hence a second order Runge-Kutta metod is: a = a =
9 8 y ( + ) = y ( ) + ( k + k) wic turns out to be te Improved Euler Metod. Te Improved Euler metod is actually a special case of te Runge-Kutta metods. We will now tackle te fourt order Runge-Kutta metod. In tis case we assume te Runge-Kutta formulas are: Wat is required is to find te coefficients y ( + ) y ( ) + ak + bk + ck + dk k = f (, y) k = f ( + m, y+ mk) k = f ( + n, y+ mk) k = f ( + p, y+ pk ) 4 4 abcdmnp,,,,,, } duplicate te Taylor series troug te term in 4. Tat is: Since dy d suc tat te above formulas 4 y ( ) y ( ) y ( ) y ( ) y ( ) y ( iv) + = ( )!! 4! = y ( ) = f (, y) ten y ( ) = f f f ( ) = + = + y f f f f y Let F = f + fyf y ( ) = F Similarly and y ( ) = f + ff + f f + f ( f + ff ) y yy y y Let F = f + ff + f f y ( ) = F + f F y yy y
10 8 (iv) y ( ) = f + ff + f f + f f + f ( f + ff + f f ) Let y yy yyy y y yy + ( f + ff )( f + ff ) + f ( f + ff ) y y yy y y F = f + ff + f f + f f y yy yyy (iv) y ( ) = F + f F + F( f + ff ) + f F y y yy y wic allows te Taylor series to be written as: y ( + ) = y ( ) + F + ( F + ff y ) [ y ( y yy) y ] 4 F + f F + f + ff F + f F +! Turning now to te various k values, similar computations produce k = f k = f + mf + m F + m F +! 6 k = f + nf + n F + mnf yf + ( ) n F + m nfyf+ 6mn ( fy + ffyy) F } +! 6 k4 = f + pf + p F + npf yf + ( ) p F + n pf yf + 6np ( f y + ff yy) F + 6mnpf yf } +! 6 Combining tese equations as suggested by te Runge-Kutta formula,
11 8 y( + ) = y( ) + ( a + b+ c+ d) f + ( bm+ cn+ dp) F 4 + ( bm + cn + dp ) F + ( bm + cn + dp ) F ( cmn + dnp) f yf + ( cm n + dn p) f yf ( cmn + dnp ) ( f + ff ) F + dmnp f F +! y yy y Comparison wit te Taylor series results in te eigt conditions a + b+ c+ d = cmn+ dnp = bm + cn + dp = cmn + dnp = bm + cn + dp = cm n + dn p = bm + cn + dp = dmnp = Since we ave eigt equations in seven unknowns we can select one of te unknowns and solve for te rest. A classical solution set is m = n = p = a = d = b = c = 6 leading to te Runge-Kutta formulas y ( + ) = y ( ) + ( k + k + k + k4) 6 k = f (, y) k = f ( +, y + k ) k = f ( +, y + k ) k = f ( +, y + k ) 4
12 84.7 Solving Higer Order ODE An n t order ODE can be broken into n simultaneous first order ODE using a state variable approac wic will be eplained troug an eample. Eample. Using =0. determine y at t=0. and 0. for te following ODE: d y dy + + y = 0. dt dt given y() 0 = y () 0 = 00. Solve using (a) Euler. (b) Improved Euler. (c) Fourt order Runge-Kutta. Solution. Let = y = y = y = = y = 0. y y = 0. or in matri form, d 0. dt = 0. dx or = G (, t X ) dt X( 0) = = t = 0 y() 0 0 () 0 = y 0 Tis approac applies to any order ODE as follows. For an n t order ODE n d y a d n y a d n y n + n + n +! + ay n = f(, ty) dt dt dt ( n ) given y() 0 = α, y () 0 = α,! y () 0 = α n
13 85 To place in te state variable form let n = y = y = y " = y ( n ) (a) Euler We can terefore write: = = = 4 " = f (, t ) a a! a n n n n Te coefficients of te ODE can be eiter constants or functions of t and/or y. Procedure At t= t i. Calculate K = G( t i, X i ). Calculate Xi+ = Xi + G( ti, Xi) Hence at t=0. K = 0. = t = = At t= (.) 0. X = X 0 + K =. 0 + = = y 0. = y (.) 0 y(.) 0 = 00.
14 86. K = = t = = X = X + K = 0. + y(.) 0 = (.) = = y 09. = y (.) 0 (b) Improved Euler Metod Procedure At t= t i. Calculate K = G( t i, X i ). Calculate XEST = Xi + K. Calculate K G(, XEST) = t i + 4. Calculate Xi+ = Xi + [ K + K] Hence at t=0. K = G( t0, X 0 ) = = t = = XEST = X 0 + =. 0 + = Κ. K = G( t = t0 +, XEST) = #$% = 09. At t= Xi = X0 + ( K + K) = 0 + y(.) 0 = =
15 K = G( t, X) = G(., ) = = XEST = X + K = = 085. & K = G( t +, XEST) = G(., 085. = = X = X + ( K + K) = = y(.) 0 = (c) Runge-Kutta 4 t order Metod Procedure At t= t i. Calculate K = G( t i, X i ). Calculate XEST = Xi + K. Calculate K = G( ti +, XEST) 4. Calculate XEST = Xi + K 5. Calculate K = G( ti +, XEST)
16 88 6. Calculate XEST = Xi + K 7. Calculate K4 = G( ti +, XEST) 6 8. Calculate Xi+ = Xi + [ K + K + K + K4] Te calculations follow a similar procedure to te Improved Euler metod and are left to te reader as an eercise. Eercise. Use te Runge-Kutta metod to solve te following ODE for te points y(0.), y(0.) and y(0.). Use =0.. d y ty dy t + + e y = 4y e dt dt y( 0) = 05. y ( 0) = 0. t Now we return back to program development. Te following is a C++ program for solving ODE using eiter Euler or Fourt order Runge-Kutta metod. To verify tat te program works it is tested against an ODE wit known analytical solution. Te ODE used is: Its analytical solution is given by: d y dy + + y = 00., y( 0) = y ( 0) = 0. dt dt 05. yt () = e cos( t) + sin( t) #include <iostream.> #include <iomanip.> //input/output manipulation for formatting using setw, etc. #include <conio.> #include <mat.> class ODE
17 89 // Te default member classification in classes is private. //Private variables double *; //state space vector [0],[],... //Te k's in te Runge-Kutta metod. // eac k is a vector of dimension = order of ODE double *k; double *k; double *k; double *k4; double *est; void(*g)(double tp,double *p,double *dot); //function for calculating te k's double t; //independent variable int neqns; //number of equations=order of ODE double ; //step lengt public: ODE(void (*GP)(double, double *,double *), int n) G=GP; neqns=n; =new double[neqns]; k=new double[neqns]; k=new double[neqns]; k=new double[neqns]; k4=new double[neqns]; est=new double[neqns]; } ~ODE() delete []; delete []k; delete []k; delete []k; delete []k4; delete []est; } void Init(double *p, double tp, double dt); double get_t()
18 90 return t; } double* get() return ; } void Euler(); void RK4(); // 4t order Runge-Kutta }; void ODE::Init(double *p, double tp, double dt) for(int i=0; i<neqns; i++) [i]=p[i]; t=tp; =dt; } void ODE::Euler() (*G)(t,,k); for(int i=0; i<neqns; i++) [i]=[i]+*k[i]; t=t+; } //Fourt order Runge-Kutta void ODE::RK4() int i; (*G)(t,,k); for(i=0; i<neqns; i++) est[i]=[i]+0.5**k[i]; (*G)(t+0.5*,est,k); for(i=0; i<neqns; i++) est[i]=[i]+ 0.5**k[i];
19 9 (*G)(t+0.5*,est,k); for(i=0; i<neqns; i++) est[i]=[i]+*k[i]; (*G)(t+,est,k4); for(i=0; i<neqns; i++) [i]=[i]+ *(k[i]+.0*k[i]+.0*k[i]+k4[i])/6.0; t=t+; } //User supplied routine // ODE y''(t)+y'(t)+y = 0, intitial conditions y'(0)=y(0)= // Analytical solution y=ep(-0.5t)*(cos(sqrt()*t/)+sqrt()*sin(sqrt()*t/)) void Vector( double t, double *, double *dot) dot[0]=[]; dot[]=0.0-[0]-[]; } // int main() ODE eqn(vector, ); //Replace (*G) wit Vector and set neqns= double *, t,, yp[0], tp[0]; =new double[]; //initial conditions [0]=.0; []=.0; t=0.0; =0.; eqn.init(,t,); //set initial conditions tp[0]=0.0; yp[0]=[0]; int i; for(i=; i<0; i++)
20 9 eqn.rk4(); =eqn.get(); yp[i]=[0]; tp[i]=eqn.get_t(); } // print solution cout << "From Runge Kutta" << " " << "t" << " " << "From Analytical solution" << endl; for(i=0; i<0; i++) cout << setprecision(7) << setw(0) << yp[i] << setprecision() << setw(0) << tp[i] << setprecision(7) << setw(0) << (ep(-0.5*tp[i])*(cos(0.866*tp[i])+.7*sin(0.866*tp[i])))<< endl; } getc(); return ; } Print-out From Runge Kutta t From Analytical solution Problems. Solve te following ODE s in 0 # t # 5 using te improved Euler metod wit =0. (write your own program) d y a. + 8y = 0, y( 0) =, y'( 0) = 0 dt
21 b. 9 d y dy y = sin( t), y( 0) = 0, y'( 0) = dt dt c. ( e y) d y t + = t, y( 0) =, y'( 0) = 0 dt. Te number of bacterial cells (P) in a given reactor is related to time in days (t) as described by te following matematical model dp dt = 0. P P If at time t=0, P=0 6 determine te number of cells wen t=0 days. Use te fourt-order Runge-Kutta metod and a time increment of day.. Wen a veicle travels over a roug road, a vertical displacement will result. Consider te following igly idealized model: Figure. Problem. Suppose tat a manufacturer is interested in determining te largest vertical displacement for a veicle traveling at a velocity of 55 miles per our over a sinusoidal surface, given
22 94 m=0 lb-s /in. k=0 lb/in. c=88 lb-s/in. L=50 ft. )= in. Te matematical model describing tis beavior is given by: m d dt c d + + k = mϖ sinϖt dt were ϖ = πvlwere v is te velocity of te veicle. Determine te maimum displacement using te Runge-Kutta 4 t order metod and )t=0.0 s. Assume tat te veicle was initially at rest ( () 0 = '( 0) = 0). 4. Use Taylor series to solve te following second-order ordinary differential equation at t=0.: d d = 0 dt dt wit te initial conditions ( 0) =, '( 0) =.
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 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 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 informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More informationUniversity Mathematics 2
University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at
More 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 informationRunge-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 informationlim 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 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 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 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 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 informationDifferentiation. 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 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 informationThe 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 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 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 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 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 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 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 informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More 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 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 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 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 informationMA2264 -NUMERICAL METHODS UNIT V : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL. By Dr.T.Kulandaivel Department of Applied Mathematics SVCE
MA64 -NUMERICAL METHODS UNIT V : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS B Dr.T.Kulandaivel Department of Applied Matematics SVCE Numerical ordinar differential equations is te part
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 informationBob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk
Bob Brown Mat 251 Calculus 1 Capter 3, Section 1 Completed 1 Te Tangent Line Problem Te idea of a tangent line first arises in geometry in te context of a circle. But before we jump into a discussion of
More 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 information1.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 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. AB Calculus Introduction
1. AB Calculus Introduction Before we get into wat calculus is, ere are several eamples of wat you could do BC (before calculus) and wat you will be able to do at te end of tis course. Eample 1: On April
More informationCHAPTER (A) When x = 2, y = 6, so f( 2) = 6. (B) When y = 4, x can equal 6, 2, or 4.
SECTION 3-1 101 CHAPTER 3 Section 3-1 1. No. A correspondence between two sets is a function only if eactly one element of te second set corresponds to eac element of te first set. 3. Te domain of a function
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 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 informationTest 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 informationSection 3: The Derivative Definition of the Derivative
Capter 2 Te Derivative Business Calculus 85 Section 3: Te Derivative Definition of te Derivative Returning to te tangent slope problem from te first section, let's look at te problem of finding te slope
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 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 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 informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationNUMERICAL DIFFERENTIATION
NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat
More informationThe 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 informationNumerical 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 information10 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 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 information1 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 informationSection 2.1 The Definition of the Derivative. We are interested in finding the slope of the tangent line at a specific point.
Popper 6: Review of skills: Find tis difference quotient. f ( x ) f ( x) if f ( x) x Answer coices given in audio on te video. Section.1 Te Definition of te Derivative We are interested in finding te slope
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
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 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 information11.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(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 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 informationKEY CONCEPT: THE DERIVATIVE
Capter Two KEY CONCEPT: THE DERIVATIVE We begin tis capter by investigating te problem of speed: How can we measure te speed of a moving object at a given instant in time? Or, more fundamentally, wat do
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 informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
More 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 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 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 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 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 information11-19 PROGRESSION. A level Mathematics. Pure Mathematics
SSaa m m pplle e UCa ni p t ter DD iff if erfe enren tiatia tiotio nn - 9 RGRSSIN decel Slevel andmatematics level Matematics ure Matematics NW FR 07 Year/S Year decel S and level Matematics Sample material
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 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 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 informationDerivatives of Exponentials
mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function
More informationCalculus 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 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 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 informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More information0.1 The Slope of Nonlinear Function
WEEK Reading [SB],.3-.7, pp. -38 0. Te Slope of Nonlinear Function If we want approimate a nonlinear function y = f() by a linear one around some point 0, te best approimation is te line tangent to te
More informationSFU 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 informationLines, Conics, Tangents, Limits and the Derivative
Lines, Conics, Tangents, Limits and te Derivative Te Straigt Line An two points on te (,) plane wen joined form a line segment. If te line segment is etended beond te two points ten it is called a straigt
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 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 informationChapter 2 Limits and Continuity. Section 2.1 Rates of Change and Limits (pp ) Section Quick Review 2.1
Section. 6. (a) N(t) t (b) days: 6 guppies week: 7 guppies (c) Nt () t t t ln ln t ln ln ln t 8. 968 Tere will be guppies ater ln 8.968 days, or ater nearly 9 days. (d) Because it suggests te number o
More informationExponentials 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 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 information(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?
Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is
More informationMTH-112 Quiz 1 Name: # :
MTH- Quiz Name: # : Please write our name in te provided space. Simplif our answers. Sow our work.. Determine weter te given relation is a function. Give te domain and range of te relation.. Does te equation
More information= h. Geometrically this quantity represents the slope of the secant line connecting the points
Section 3.7: Rates of Cange in te Natural and Social Sciences Recall: Average rate of cange: y y y ) ) ), ere Geometrically tis quantity represents te slope of te secant line connecting te points, f (
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 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 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 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 informationSection 2: The Derivative Definition of the Derivative
Capter 2 Te Derivative Applied Calculus 80 Section 2: Te Derivative Definition of te Derivative Suppose we drop a tomato from te top of a 00 foot building and time its fall. Time (sec) Heigt (ft) 0.0 00
More informationf a h f a h h lim lim
Te Derivative Te derivative of a function f at a (denoted f a) is f a if tis it exists. An alternative way of defining f a is f a x a fa fa fx fa x a Note tat te tangent line to te grap of f at te point
More 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 informationDifferentiation. introduction to limits
9 9A Introduction to limits 9B Limits o discontinuous, rational and brid unctions 9C Dierentiation using i rst principles 9D Finding derivatives b rule 9E Antidierentiation 9F Deriving te original unction
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 informationSection 2.4: Definition of Function
Section.4: Definition of Function Objectives Upon completion of tis lesson, you will be able to: Given a function, find and simplify a difference quotient: f ( + ) f ( ), 0 for: o Polynomial functions
More informationChapter 1 Functions and Graphs. Section 1.5 = = = 4. Check Point Exercises The slope of the line y = 3x+ 1 is 3.
Capter Functions and Graps Section. Ceck Point Exercises. Te slope of te line y x+ is. y y m( x x y ( x ( y ( x+ point-slope y x+ 6 y x+ slope-intercept. a. Write te equation in slope-intercept form: x+
More informationDEFINITION OF A DERIVATIVE
DEFINITION OF A DERIVATIVE Section 2.1 Calculus AP/Dual, Revised 2017 viet.dang@umbleisd.net 2.1: Definition of a Derivative 1 DEFINITION A. Te derivative of a function allows you to find te SLOPE OF THE
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More informationLimits and an Introduction to Calculus
Limits and an Introduction to Calculus. Introduction to Limits. Tecniques for Evaluating Limits. Te Tangent Line Problem. Limits at Infinit and Limits of Sequences.5 Te Area Problem In Matematics If a
More informationMATH 111 CHAPTER 2 (sec )
MATH CHAPTER (sec -0) Terms to know: function, te domain and range of te function, vertical line test, even and odd functions, rational power function, vertical and orizontal sifts of a function, reflection
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 informationLogarithmic functions
Roberto s Notes on Differential Calculus Capter 5: Derivatives of transcendental functions Section Derivatives of Logaritmic functions Wat ou need to know alread: Definition of derivative and all basic
More information