Exponentially-Fitted Runge-Kutta Nystrom Method of Order Three for Solving Oscillatory Problems ABSTRACT 1. INTRODUCTION
|
|
- Abel Patrick
- 5 years ago
- Views:
Transcription
1 Malaysian Journal of Mathematical Sciences 8(S): 7-4 (04) Special Issue: International Conference on Mathematical Sciences and Statistics 0 (ICMSS0) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: Exponentially-Fitted Runge-Kutta Nystrom Method of Order Three for Solving Oscillatory Problems * Faranak Rabiei, Fudziah Ismail and Norazak Senu Department of Mathematics, Faculty of Sciences and Institute for Mathematical Research, Universiti Putra Malaysia, 4400 UPM Serdang, Selangor, Malaysia faranak_rabiei@upm.edu.my *Corresponding author ABSTRACT In this paper the exponentially fitted explicit Runge-Kutta Nystrom method is proposed for solving special second-order ordinary differential equations where the solution is oscillatory. The exponentially fitting is based on given Runge-Kutta Nystrom (RKN) method of order three at a cost of three function evaluations per step. Here, we also developed the trigonometrically-fitted RKN method for solving initial value problems with oscillating solutions. The numerical results compared with the existing explicit RKN method of order three which indicates that the exponentially fitted explicit Runge-Kutta Nystrom method is more efficient than the classical RKN method. Keywords: Exponentially-fitting, Runge-Kutta Nystrom method, oscillating solutions.. INTRODUCTION Consider the second order ordinary differential equation (ODEs) y( x) f ( x, y( x)), x[ x0 X ], y( x ) y, y( x ) y For approximating solution of problem in equation (E) many different methods have been studied. One of the most common method for solving numerically () is Runge-Kutta Nystrom (RKN) methods (see Dormand (996)). Exponentially fitting is a procedure for an efficient of exponential, ()
2 Faranak Rabiei, Fudziah Ismail & Norazak Senu trigonometric or hyperbolic function for solving ODEs whose solutions are oscillatory. Simos (00) and Franco (004) proposed the different approaches for exponentially fitted on explicit Runge-Kutta Nystrom methods. In this paper we propose the exponentially and trigonometrically fitting for third order, -stages Runge-Kutta Nystrom method with minimum phase-lag derived by Norazak et al. (00) also the exponentially fitting for RKN method of order three with normal set of coefficients. In section, for exponentially fitting the RKN method integrates exactly the functions exp( iwx) where w R. In Section, the construction of trigonometrically fitting RKN methods by integrating exactly the functions exp( iwx) where wrandi, is proposed. The number of test problems are experienced and the numerical results compared with existing RKN method of order derived by Dormand (996), are given in last section. The explicit RKN method is given by s s n n n i i n n i i i i y y hy h b k, y y h b k, i n n n i n i n ij j j k f ( x, y ), k f ( x c h, y hc y h a k ), i,, s. () c [ c, c,, c ], b [ b, b,, b ], b [ b, b,, b ], and ( s s) T T T where s s s matrix A [ a ij ]. The characteristic equation for RKN formulas () is given by Franco (004), tr( D( z )) det( D( z )) 0 A( z ) B( z ) where D and A, A, B, Bare polynomial in A'( z ) B'( z ) method with characteristic equation where R z ( ) tr( D) Rz ( ) ( z) z cos ( ), and Sz ( ) ( z) S( z ), z. For RKN S( z ) det( D), ( z) is phase-lag and ( z) is q dissipation. RKN method has phase-lag order q when ( z) O( z ) and r RKN method has dissipation order when ( z) O( z ). 8 Malaysian Journal of Mathematical Sciences
3 Exponentially-Fitted Runge-Kutta Nystrom Method of Order Three for Solving Oscillatory Problems Norazak et al. (00) developed the explicit RKN method of algebraic order p and phase-lag order q 6 at a cost of three function evaluations per step of integration. Here, we chose the values of free parameters from his method.. EXPONENTIALLY-FITTED THIRD ORDER RKN METHOD To construct the exponentially-fitted RKN method, we require that all stages of third order RKN method integrate exactly to the function exp( wx). We have e c v v a 0, e c v v a 0, cv cv cv cv ( ) 0, v cv cv v cv cv v cv cv ( ) 0, cv cv e c v v a a e e c v v ( a a e ) 0, e v v ( b b e b e ) 0, e v v ( b b e b e ) 0, v cv cv e v b b e b e e v( b b e b e ) 0, where v wh, w R. Using relations: v v v v e e e e cosh( v), sinh( v). the following order conditions are obtained: cosh( c v) v a, cosh( c v) v ( a a cosh( c v)), sinh( cv) cv v a sinh( cv)), cosh( v) v ( b b cosh( cv) b cosh( cv)), sinh( v) v v ( b sinh( cv) b sinh( cv)), cosh( v) v( b sinh( c v) b sinh( c v)), sinh( v) v( b b cosh( c v) b cosh( c v)). We set the values of parameters c, c and a from RKN 6 methods with minimum phase-lag derived by Norazak et al. (00), into above order conditions, so we have (EF-N method): Malaysian Journal of Mathematical Sciences 9
4 Faranak Rabiei, Fudziah Ismail & Norazak Senu a , a 0.805, b , b , b , b , b , b b. Also by choosing c, c and a 0 as free parameters, we have (EF method): a , a , b , b , b , b , b , b b.. TRIGNOMETRICALLY-FITTED THIRD ORDER RKN METHOD To construct the Trigonometrically-fitted RKN method, we require that all stages of third order RKN method integrate exactly to the function exp( iwx) where i, so we have e ic v v a 0, e ic v v a 0, icv ic v icv icv icv icv ( ) 0, e icv v a ae icv icv e icv v ( a ae ) 0, iv icv icv e iv v ( b be be ) 0, iv ( ic v ic v iv icv ic v ) 0, ie i v b be be iv icv icv ie i v( b be be ) 0, e ic v v a a e ( ) 0, e iv v b b e b e ( ) 0, iv iv iv iv e e e e where v wh, w R. Using cos( v), sin( v), the i following order conditions are obtained cos( c v) v a, cos( c v) v ( a a cos( c v)), sin( c v) c v v a sin( c v), cos( v) v ( b b cos( c v) b cos( c v)), sin( v) v v ( b sin( c v) b sin( c v)), cos( v) v( b sin( cv) b sin( cv)), sin( v) v( b b cos( c v) b cos( c v)). 0 Malaysian Journal of Mathematical Sciences
5 Exponentially-Fitted Runge-Kutta Nystrom Method of Order Three for Solving Oscillatory Problems We set the values of parameters c, c and a from RKN 6 methods with minimum phase-lag derived by Norazak et al. (00) into above order coditions, so we have (TF-N Method): a , a , b , b , b , b , b , b b. 4. NUMERICAL RESULTS To illustrate the efficiency of new methods the following problems are solved for x [0,50] and h,,, and numerical results 4 0 are compared with existing third order RKN method derived by Dormand (996). In addition, the following abbreviations are used in Figures -. EF-N: exponentially-fitted Runge-Kutta Nystrom method of order three with minimum phase-lag. EF: exponentially-fitted Runge-Kutta Nystrom method of order three when a 0. TF-N: trigonometrically-fitted Runge-Kutta Nystrom method of order three with minimum phase-lag. RKND: classical Runge-Kutta Nystrom method of order three derived by Dormand (996). Problem. (Norazak et al (00)). y y 0.00cos( x), y(0), y(0) 0. Exact solution: y( x) cos( x) xsin( x). Problem. (Franco (006)). y y 0.00cos( x), y (0), y(0) 0, y y 0.00sin( x), y (0) 0, y (0) Malaysian Journal of Mathematical Sciences
6 Faranak Rabiei, Fudziah Ismail & Norazak Senu Exact solutions: y ( x ) cos( x ) x sin( x ), y ( x ) sin( x ) x cos( x ). Figure : Logarithm of maximum global error versusnumber of function evaluations for solving problem Figure : Logarithm of maximum global error versus number of function evaluations for solving problem Figure : Maximum global error of EF- N method (dash line) and RKND (solid line) in solving problem Figure 4: Maximum global error of EF method (dash line) and RKND (solid line) in solving problem Figure 5: Maximum global error of F-N method (dashline) and RKND (solid line) in solving problem Figure 6: Maximum global error of EF- N method (dash line) and RKND (solid line) for y in solving problem Malaysian Journal of Mathematical Sciences
7 Exponentially-Fitted Runge-Kutta Nystrom Method of Order Three for Solving Oscillatory Problems Figure 7: Maximum global error of EF- N method (dash line) and RKND (solid line) for y in solving problem Figure 8: Maximum global error of EF method (dash line) and RKND (solid line) for y in solving problem Figure 9: Maximum global error of EF method (dash line) and RKND (solid line) for y in solving problem Figure 0: Maximum global error of TF- N method (dash line) and RKND (solid line) for y in solving problem Figure : Maximum global error of TF-N method (dash line) and RKND (solid line) for y in solving problem Malaysian Journal of Mathematical Sciences
8 Faranak Rabiei, Fudziah Ismail & Norazak Senu 5. CONCLUSION In this study we constructed the exponentially- fitted and trigonometrically- fitted Runge-Kutta Nystrom method of order three with - stages. The new methods give the accurate results for problems with oscillation solutions while the step size h was not chosen as small value. Numerical examples indicate that the exponentially fitted RKN methods are more efficient than given classical RKN method. REFERENCES Dormand, J. R. (996). Numerical Method for Differential Equations (A Computational Approach). CRC Press. Inc. Franco, J. M. (004). Exponentially-Fitted Explicit Runge-Kutta Nystrom Methods. Journal of Computational and Applied Mathematics. 67:-9. Franco, J. M. (006). A Class of Explicit Two Step Hybrid Methods for Second Order IVP's. Journal of Computation and Applied Mathematics. 87: Norazak, B. S., Suleiman, M, Isamil, F. and Othman, M. (00). A Zero- Dissipative Runge-Kutta Nyströ m Methods with Minimum Pahse- Lag. Mathematical Problem in Engineering. DOI:0.555/00/ 594. Simos, T. E. (00). Exponentially-Fitted Runge-Kutta Nystrom Method for the Numerical Solution of initial value problems with oscillating solutions. Applied Mathematics Letter. 5: Malaysian Journal of Mathematical Sciences
A New Embedded Phase-Fitted Modified Runge-Kutta Method for the Numerical Solution of Oscillatory Problems
Applied Mathematical Sciences, Vol. 1, 16, no. 44, 157-178 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/ams.16.64146 A New Embedded Phase-Fitted Modified Runge-Kutta Method for the Numerical Solution
More informationFifth-Order Improved Runge-Kutta Method With Reduced Number of Function Evaluations
Australian Journal of Basic and Applied Sciences, 6(3): 9-5, 22 ISSN 99-88 Fifth-Order Improved Runge-Kutta Method With Reduced Number of Function Evaluations Faranak Rabiei, Fudziah Ismail Department
More informationImproved Extended Runge-Kutta-like Method for Solving First Order IVPs
Indian Journal of Science and Technology, Vol 9(48), DOI: 0.7485/ijst/06/v9i48/0938, December 06 ISSN (Print) : 0974-6846 ISSN (Online) : 0974-5645 Improved Extended Runge-Kutta-like Method for Solving
More informationSains Malaysiana 47(9)(2018):
Sains Malaysiana 47(9)(2018): 2223 2230 http://dx.doi.org/10.17576/jsm-2018-4709-33 Block Hybrid Method with Trigonometric-Fitting for Solving Oscillatory Problems (Kaedah Blok Hibrid dengan Penyuaian-Trigonometri
More informationEmbedded 5(4) Pair Trigonometrically-Fitted Two Derivative Runge- Kutta Method with FSAL Property for Numerical Solution of Oscillatory Problems
Embedded 5(4) Pair Trigonometrically-Fitted Two Derivative Runge- Kutta Method with FSAL Property for Numerical Solution of Oscillatory Problems N. SENU, N. A. AHMAD, F. ISMAIL & N. BACHOK Institute for
More informationA Modified Runge-Kutta-Nyström Method by using Phase Lag Properties for the Numerical Solution of Orbital Problems
Appl. Math. Inf. Sci., No., - () Applied Mathematics & Information Sciences An International Journal A Modified Runge-Kutta-Nyström Method by using Phase Lag Properties for the Numerical Solution of Orbital
More informationA NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN FILM FLOW PROBLEM
Indian J. Pure Appl. Math., 491): 151-167, March 218 c Indian National Science Academy DOI: 1.17/s13226-18-259-6 A NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN
More informationResearch Article Diagonally Implicit Block Backward Differentiation Formulas for Solving Ordinary Differential Equations
International Mathematics and Mathematical Sciences Volume 212, Article ID 767328, 8 pages doi:1.1155/212/767328 Research Article Diagonally Implicit Block Backward Differentiation Formulas for Solving
More informationarxiv: v1 [math.na] 15 Nov 2008
arxiv:0811.2481v1 [math.na] 15 Nov 2008 A Phase-Fitted Runge-Kutta-Nyström method for the Numerical Solution of Initial Value Problems with Oscillating Solutions Abstract D. F. Papadopoulos a, Z. A. Anastassi
More informationarxiv: v1 [math.na] 31 Oct 2016
RKFD Methods - a short review Maciej Jaromin November, 206 arxiv:60.09739v [math.na] 3 Oct 206 Abstract In this paper, a recently published method [Hussain, Ismail, Senua, Solving directly special fourthorder
More informationGeneralized RK Integrators for Solving Ordinary Differential Equations: A Survey & Comparison Study
Global Journal of Pure and Applied Mathematics. ISSN 7-78 Volume, Number 7 (7), pp. Research India Publications http://www.ripublication.com/gjpam.htm Generalized RK Integrators for Solving Ordinary Differential
More informationExponentially Fitted Symplectic Runge-Kutta-Nyström methods
Appl. Math. Inf. Sci. 7, No. 1, 81-85 (2013) 81 Applied Mathematics & Information Sciences An International Journal Exponentially Fitted Symplectic Runge-Kutta-Nyström methods Th. Monovasilis 1, Z. Kalogiratou
More informationEmbedded Explicit Two-Step Runge-Kutta-Nyström for Solving Special Second-Order IVPs
Applied Mathematical Sciences, Vol. 8, 14, no. 94, 4655 4673 HIKARI Ltd, www.mhikari.com http://dx.doi.org/1.1988/ams.14.4643 Embedded Explicit TwoStep RungeKuttaNyström for Solving Special SecondOrder
More informationFourth-order symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: a review
Fourth-order symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: a review G. Vanden Berghe and M. Van Daele Vakgroep Toegepaste Wiskunde en Informatica, Universiteit Gent, Krijgslaan
More informationResearch Article An Accurate Block Hybrid Collocation Method for Third Order Ordinary Differential Equations
Applied Mathematics, Article ID 549597, 9 pages http://dx.doi.org/1.1155/14/549597 Research Article An Accurate Bloc Hybrid Collocation Method for Third Order Ordinary Differential Equations Lee Ken Yap,
More informationTwo Optimized Runge-Kutta Methods for the Solution of the Schrödinger Equation
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. (8) 7-77 ISSN - Two Optimized Runge-Kutta Methods for the Solution of the Schrödinger Equation T.V. Triantafyllidis,
More informationAdebayo O. Adeniran1, Saheed O. Akindeinde2 and Babatunde S. Ogundare2 * Contents. 1. Introduction
Malaya Journal of Matematik Vol. No. 73-73 8 https://doi.org/.37/mjm/ An accurate five-step trigonometrically-fitted numerical scheme for approximating solutions of second order ordinary differential equations
More informationEfficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators
MATEMATIKA, 8, Volume 3, Number, c Penerbit UTM Press. All rights reserved Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators Annie Gorgey and Nor Azian Aini Mat Department of Mathematics,
More informationSymplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: revisited
Symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: revisited G. Vanden Berghe and M. Van Daele Vakgroep Toegepaste Wiskunde en Informatica, Universiteit Gent, Krijgslaan 8
More informationDirect Explicit Integrators of RK Type for Solving Special Fourth-Order Ordinary Differential Equations With An Application
Global Journal of Pure and Applied Mathematics. ISSN 0973-768 Volume, Number 6 (06), pp. 4687 475 Research India Publications http://www.ripublication.com/gjpam.htm Direct Explicit Integrators of RK Type
More informationSolving Delay Differential Equations By Elzaki Transform Method
SCITECH Volume 3, Issue 1 RESEARCH ORGANISATION February 22, 2017 Boson Journal of Modern Physics www.scitecresearch.com Solving Delay Differential Equations By Elzaki Transform Method Ebimene, J. Mamadu
More informationAn Optimized Symmetric 8-Step Semi-Embedded Predictor-Corrector Method for IVPs with Oscillating Solutions
Appl. Math. Inf. Sci., No., - () Applied Mathematics & Information Sciences An International Journal An Optimized Symmetric -Step Semi-Embedded Predictor-Corrector Method for IVPs with Oscillating Solutions
More informationSymplectic Runge-Kutta-Nyström Methods with Phase-Lag Oder 8 and Infinity
Appl. Math. Inf. Sci. 9, No. 3, 1105-1112 (2015) 1105 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.12785/amis/090302 Symplectic Runge-Kutta-Nyström Methods with
More informationResearch Article P-Stable Higher Derivative Methods with Minimal Phase-Lag for Solving Second Order Differential Equations
Hindawi Publishing Corporation Journal of Applied Mathematics Volume 2011, Article ID 407151, 15 pages doi:10.1155/2011/407151 Research Article P-Stable Higher Derivative Methods with Minimal Phase-Lag
More informationSolving Zhou Chaotic System Using Fourth-Order Runge-Kutta Method
World Applied Sciences Journal 21 (6): 939-944, 2013 ISSN 11-4952 IDOSI Publications, 2013 DOI: 10.529/idosi.wasj.2013.21.6.2915 Solving Zhou Chaotic System Using Fourth-Order Runge-Kutta Method 1 1 3
More informationResearch Article Solving Delay Differential Equations of Small and Vanishing Lag Using Multistep Block Method
Applied Mathematics, Article ID 348912, 10 pages http://dx.doi.org/10.1155/2014/348912 Research Article Solving Delay Differential Equations of Small and Vanishing Lag Using Multistep Bloc Method Nurul
More informationOn the reliability and stability of direct explicit Runge-Kutta integrators
Global Journal of Pure and Applied Mathematics. ISSN 973-78 Volume, Number 4 (), pp. 3959-3975 Research India Publications http://www.ripublication.com/gjpam.htm On the reliability and stability of direct
More informationOrdinary differential equation II
Ordinary Differential Equations ISC-5315 1 Ordinary differential equation II 1 Some Basic Methods 1.1 Backward Euler method (implicit method) The algorithm looks like this: y n = y n 1 + hf n (1) In contrast
More informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More information11.4. Differentiating ProductsandQuotients. Introduction. Prerequisites. Learning Outcomes
Differentiating ProductsandQuotients 11.4 Introduction We have seen, in the first three Sections, how standard functions like n, e a, sin a, cos a, ln a may be differentiated. In this Section we see how
More informationTrigonometrically tted predictor corrector methods for IVPs with oscillating solutions
Journal of Computational and Applied Mathematics 58 (23) 35 44 www.elsevier.com/locate/cam Trigonometrically tted predictor corrector methods for IVPs with oscillating solutions G.Psihoyios a, T.E. Simos
More informationPhysics 6303 Lecture 15 October 10, Reminder of general solution in 3-dimensional cylindrical coordinates. sinh. sin
Physics 6303 Lecture 15 October 10, 2018 LAST TIME: Spherical harmonics and Bessel functions Reminder of general solution in 3-dimensional cylindrical coordinates,, sin sinh cos cosh, sin sin cos cos,
More informationResearch Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
Applied Mathematics Volume 013, Article ID 591636, 5 pages http://dx.doi.org/10.1155/013/591636 Research Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
More informationSolving Delay Differential Equations (DDEs) using Nakashima s 2 Stages 4 th Order Pseudo-Runge-Kutta Method
World Applied Sciences Journal (Special Issue of Applied Math): 8-86, 3 ISSN 88-495; IDOSI Publications, 3 DOI:.589/idosi.wasj.3..am.43 Solving Delay Differential Equations (DDEs) using Naashima s Stages
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationCurriculum Vitae Hassan Eltayeb Gadain Salaim Institute for Mathematics Research University Putra Malaysia, UPM, Serdang, Selangor
Curriculum Vitae Hassan Eltayeb Gadain Salaim Institute for Mathematics Research University Putra Malaysia, 43400 UPM, Serdang, Selangor Personal Information Name: Hassan Eltayeb Gadaub Salaim Date of
More informationDeferred correction based on exponentially-fitted mono-implicit Runge-Kutta methods
Deferred correction based on exponentially-fitted mono-implicit Runge-Kutta methods M. Van Daele Department of Applied Mathematics, Computer Science and Statistics Ghent University 3th International Conference
More informationStrauss PDEs 2e: Section Exercise 4 Page 1 of 6
Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where
More informationWhat we ll do: Lecture 21. Ordinary Differential Equations (ODEs) Differential Equations. Ordinary Differential Equations
What we ll do: Lecture Ordinary Differential Equations J. Chaudhry Department of Mathematics and Statistics University of New Mexico Review ODEs Single Step Methods Euler s method (st order accurate) Runge-Kutta
More informationTest one Review Cal 2
Name: Class: Date: ID: A Test one Review Cal 2 Short Answer. Write the following expression as a logarithm of a single quantity. lnx 2ln x 2 ˆ 6 2. Write the following expression as a logarithm of a single
More informationOn Nonlinear Methods for Stiff and Singular First Order Initial Value Problems
Nonlinear Analysis and Differential Equations, Vol. 6, 08, no., 5-64 HIKARI Ltd, www.m-hikari.com https://doi.org/0.988/nade.08.8 On Nonlinear Methods for Stiff and Singular First Order Initial Value Problems
More informationRevision Checklist. Unit FP3: Further Pure Mathematics 3. Assessment information
Revision Checklist Unit FP3: Further Pure Mathematics 3 Unit description Further matrix algebra; vectors, hyperbolic functions; differentiation; integration, further coordinate systems Assessment information
More informationFurther Pure Mathematics 3 GCE Further Mathematics GCE Pure Mathematics and Further Mathematics (Additional) A2 optional unit
Unit FP3 Further Pure Mathematics 3 GCE Further Mathematics GCE Pure Mathematics and Further Mathematics (Additional) A optional unit FP3.1 Unit description Further matrix algebra; vectors, hyperbolic
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Multistep Methods Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 1. Runge Kutta methods 2. Embedded RK methods
More informationCOSC 3361 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods
COSC 336 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods Fall 2005 Repetition from the last lecture (I) Initial value problems: dy = f ( t, y) dt y ( a) = y 0 a t b Goal:
More information6.2. The Hyperbolic Functions. Introduction. Prerequisites. Learning Outcomes
The Hyperbolic Functions 6. Introduction The hyperbolic functions cosh x, sinh x, tanh x etc are certain combinations of the exponential functions e x and e x. The notation implies a close relationship
More information= (, ) V λ (1) λ λ ( + + ) P = [ ( ), (1)] ( ) ( ) = ( ) ( ) ( 0 ) ( 0 ) = ( 0 ) ( 0 ) 0 ( 0 ) ( ( 0 )) ( ( 0 )) = ( ( 0 )) ( ( 0 )) ( + ( 0 )) ( + ( 0 )) = ( + ( 0 )) ( ( 0 )) P V V V V V P V P V V V
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationPursuit Differential Game Described by Infinite First Order 2-Systems of Differential Equations
Malaysian Journal of Mathematical Sciences 11(2): 181 19 (217) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal Pursuit Differential Game Described by Infinite
More informationPONDI CHERRY UNI VERSI TY
B.Sc. ALLIED MATHEMATICS SYLLABUS 2009-2010 onwards PONDI CHERRY UNI VERSI TY PUDUCHERRY 605 014 B.Sc. ALLIED MATHEMATICS Syllabus for Allied Mathematics for B.Sc. Physics Main/Chemistry Main/Electronics
More informationMonday, 6 th October 2008
MA211 Lecture 9: 2nd order differential eqns Monday, 6 th October 2008 MA211 Lecture 9: 2nd order differential eqns 1/19 Class test next week... MA211 Lecture 9: 2nd order differential eqns 2/19 This morning
More informationExercises, module A (ODEs, numerical integration etc)
FYTN HT18 Dept. of Astronomy and Theoretical Physics Lund University, Sweden Exercises, module A (ODEs, numerical integration etc) 1. Use Euler s method to solve y (x) = y(x), y() = 1. (a) Determine y
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationFour Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential Equations
International journal of scientific and technical research in engineering (IJSTRE) www.ijstre.com Volume Issue ǁ July 206. Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential
More informationLecture 3. Lecturer: Prof. Sergei Fedotov Calculus and Vectors. Exponential and logarithmic functions
Lecture 3 Lecturer: Prof. Sergei Fedotov 10131 - Calculus and Vectors Exponential and logarithmic functions Sergei Fedotov (University of Manchester) MATH10131 2011 1 / 7 Lecture 3 1 Inverse functions
More informationTHE CONVERGENCE AND ORDER OF THE 3-POINT BLOCK EXTENDED BACKWARD DIFFERENTIATION FORMULA
VOL 7 NO DEEMBER ISSN 89-668 6- Asian Research Publishing Networ (ARPN) All rights reserved wwwarpnournalscom THE ONVERGENE AND ORDER OF THE -POINT BLOK EXTENDED BAKWARD DIFFERENTIATION FORMULA H Musa
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationHyperbolic functions
Roberto s Notes on Differential Calculus Chapter 5: Derivatives of transcendental functions Section Derivatives of Hyperbolic functions What you need to know already: Basic rules of differentiation, including
More informationChapter 5 Exercises. (a) Determine the best possible Lipschitz constant for this function over 2 u <. u (t) = log(u(t)), u(0) = 2.
Chapter 5 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rjl/fdmbook Exercise 5. (Uniqueness for
More informationLinear Differential Equations. Problems
Chapter 1 Linear Differential Equations. Problems 1.1 Introduction 1.1.1 Show that the function ϕ : R R, given by the expression ϕ(t) = 2e 3t for all t R, is a solution of the Initial Value Problem x =
More informationFUNCTIONS OF ONE VARIABLE FUNCTION DEFINITIONS
Page of 6 FUNCTIONS OF ONE VARIABLE FUNCTION DEFINITIONS 6. HYPERBOLIC FUNCTIONS These functions which are defined in terms of e will be seen later to be related to the trigonometic functions via comple
More informationPractice Differentiation Math 120 Calculus I Fall 2015
. x. Hint.. (4x 9) 4x + 9. Hint. Practice Differentiation Math 0 Calculus I Fall 0 The rules of differentiation are straightforward, but knowing when to use them and in what order takes practice. Although
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationConsidering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.
Lecture 3 Usual complex functions MATH-GA 245.00 Complex Variables Polynomials. Construction f : z z is analytic on all of C since its real and imaginary parts satisfy the Cauchy-Riemann relations and
More informationChapter 13: General Solutions to Homogeneous Linear Differential Equations
Worked Solutions 1 Chapter 13: General Solutions to Homogeneous Linear Differential Equations 13.2 a. Verifying that {y 1, y 2 } is a fundamental solution set: We have y 1 (x) = cos(2x) y 1 (x) = 2 sin(2x)
More informationThe Method of Undetermined Coefficients.
The Method of Undetermined Coefficients. James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University May 24, 2017 Outline 1 Annihilators 2 Finding The
More information) # ( 281). Give units with your answer.
Math 120 Winter 2009 Handout 17: In-Class Review for Exam 2 The topics covered by Exam 2 in the course include the following: Implicit differentiation. Finding formulas for tangent lines using implicit
More informationHigher Order Linear ODEs
im03.qxd 9/21/05 11:04 AM Page 59 CHAPTER 3 Higher Order Linear ODEs This chapter is new. Its material is a rearranged and somewhat extended version of material previously contained in some of the sections
More informationAn Adaptive Hierarchical Matrix on Point Iterative Poisson Solver
Malaysian Journal of Mathematical Sciences 10(3): 369 382 (2016) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal An Adaptive Hierarchical Matrix on Point
More informationMath 322. Spring 2015 Review Problems for Midterm 2
Linear Algebra: Topic: Linear Independence of vectors. Question. Math 3. Spring Review Problems for Midterm Explain why if A is not square, then either the row vectors or the column vectors of A are linearly
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationKink, singular soliton and periodic solutions to class of nonlinear equations
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 10 Issue 1 (June 015 pp. 1 - Applications and Applied Mathematics: An International Journal (AAM Kink singular soliton and periodic
More informationA Special High Order Runge-Kutta Type Method for the Solution of the Schrödinger Equation
Appl. Math. Inf. Sci. 9, No. 5, 559-577 015 559 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.1785/amis/090541 A Special High Order Runge-Kutta Type Method for
More informationAccelerated Shift-and-Add algorithms
c,, 1 12 () Kluwer Academic Publishers, Boston. Manufactured in The Netherlands. Accelerated Shift-and-Add algorithms N. REVOL AND J.-C. YAKOUBSOHN nathalie.revol@univ-lille1.fr, yak@cict.fr Lab. ANO,
More information[ ] is a vector of size p.
Lecture 11 Copyright by Hongyun Wang, UCSC Recap: General form of explicit Runger-Kutta methods for solving y = F( y, t) k i = hfy n + i1 j =1 c ij k j, t n + d i h, i = 1,, p + p j =1 b j k j A Runge-Kutta
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More informationChapter 5 Logarithmic, Exponential, and Other Transcendental Functions
Chapter 5 Logarithmic, Exponential, an Other Transcenental Functions 5.1 The Natural Logarithmic Function: Differentiation 5.2 The Natural Logarithmic Function: Integration 5.3 Inverse Functions 5.4 Exponential
More informationENGI 3425 Review of Calculus Page then
ENGI 345 Review of Calculus Page 1.01 1. Review of Calculus We begin this course with a refresher on ifferentiation an integration from MATH 1000 an MATH 1001. 1.1 Reminer of some Derivatives (review from
More informationODE Runge-Kutta methods
ODE Runge-Kutta methods The theory (very short excerpts from lectures) First-order initial value problem We want to approximate the solution Y(x) of a system of first-order ordinary differential equations
More informationLaplace Theory Examples
Laplace Theory Examples Harmonic oscillator s-differentiation Rule First shifting rule Trigonometric formulas Exponentials Hyperbolic functions s-differentiation Rule First Shifting Rule I and II Damped
More informationDr. Back. Nov. 3, 2009
Dr. Back Nov. 3, 2009 Please Don t Rely on this File! In you re expected to work on these topics mostly without computer aid. But seeing a few better pictures can help understanding the concepts. A copy
More informationSecond-Order Linear ODEs
Chap. 2 Second-Order Linear ODEs Sec. 2.1 Homogeneous Linear ODEs of Second Order On pp. 45-46 we extend concepts defined in Chap. 1, notably solution and homogeneous and nonhomogeneous, to second-order
More informationMITIGATION OF POWER SYSTEM SMALL SIGNAL OSCILLATION USING POSICAST CONTROLLER AND EVOLUTIONARY PROGRAMMING
Journal of Engineering Science and Technology Special Issue on Applied Engineering and Sciences, October (2014) 39-50 School of Engineering, Taylor s University MITIGATION OF POWER SYSTEM SMALL SIGNAL
More informationMA 201 Complex Analysis Lecture 6: Elementary functions
MA 201 Complex Analysis : The Exponential Function Recall: Euler s Formula: For y R, e iy = cos y + i sin y and for any x, y R, e x+y = e x e y. Definition: If z = x + iy, then e z or exp(z) is defined
More informationResearch Article The Spectral Method for Solving Sine-Gordon Equation Using a New Orthogonal Polynomial
International Scholarly Research Network ISRN Applied Mathematics Volume, Article ID 4673, pages doi:.54//4673 Research Article The Spectral Method for Solving Sine-Gordon Equation Using a New Orthogonal
More informationThe Solution of Weakly Nonlinear Oscillatory Problems with No Damping Using MAPLE
World Applied Sciences Journal (): 64-69, 0 ISSN 88-495 IDOSI Publications, 0 DOI: 0.589/idosi.wasj.0..0.09 The Solution of Weakly Nonlinear Oscillatory Problems with No Damping Using MAPLE N. Hashim,
More information10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections )
c Dr. Igor Zelenko, Fall 2017 1 10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections 7.2-7.4) 1. When each of the functions F 1, F 2,..., F n in right-hand side
More informationBASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA
1 BASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA This part of the Basic Exam covers topics at the undergraduate level, most of which might be encountered in courses here such as Math 233, 235, 425, 523, 545.
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationOrdinary Differential Equations (ODEs)
Chapter 13 Ordinary Differential Equations (ODEs) We briefly review how to solve some of the most standard ODEs. 13.1 First Order Equations 13.1.1 Separable Equations A first-order ordinary differential
More informationCalculus Differentiation Norhafizah Md Sarif Faculty of Industrial Science & Technology
Calculus Differentiation By Norhafizah Md Sarif Faculty of Industrial Science & Technology norhafizah@ump.edu.my Description Aims This chapter is aimed to : 1. introduce the concept of integration. evaluate
More informationRegression and Nonlinear Axes
Introduction to Chemical Engineering Calculations Lecture 2. What is regression analysis? A technique for modeling and analyzing the relationship between 2 or more variables. Usually, 1 variable is designated
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt =
More information1 Holomorphic functions
Robert Oeckl CA NOTES 1 15/09/2009 1 1 Holomorphic functions 11 The complex derivative The basic objects of complex analysis are the holomorphic functions These are functions that posses a complex derivative
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More information2.1 Limits, Rates of Change and Slopes of Tangent Lines
2.1 Limits, Rates of Change and Slopes of Tangent Lines (1) Average rate of change of y f x over an interval x 0,x 1 : f x 1 f x 0 x 1 x 0 Instantaneous rate of change of f x at x x 0 : f x lim 1 f x 0
More informationNumerical solution of stiff ODEs using second derivative general linear methods
Numerical solution of stiff ODEs using second derivative general linear methods A. Abdi and G. Hojjati University of Tabriz, Tabriz, Iran SciCADE 2011, Toronto 0 5 10 15 20 25 30 35 40 1 Contents 1 Introduction
More informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationModule 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 1 : Introduction
Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 1 : Introduction 1 Introduction In this module, we develop solution techniques for numerically solving ordinary
More information