On The Stability and Accuracy of Some Runge-Kutta Methods of Solving Second Order Ordinary Differential Equations
|
|
- Stella Atkinson
- 5 years ago
- Views:
Transcription
1 Iteratioal Joural o Computatioal Egieerig Resear Vol Issue 7 O Te Stabilit ad Aura o Some Ruge-Kutta Metods o Solvig Seod Order Ordiar Dieretial Euatios S.O. Salawu R.A. Kareem ad O.T. Arowolo Departmet o Matematis Uiversit o Ilori Ilori Nigeria Departmet o Matematis Lagos State Poltei Iorodu Nigeria ABSTRACT. Tis paper sees umerial solutios to seod order dieretial euatios o te orm wit iitial value usig dieret Ruge-Kutta metods o order two. Two ases o Epliit Ruge-Kutta metod were derived ad teir stabilit was determied tis is te implemeted. Te results were ompared wit te Euler s metod or aura. KEYWORDS. Euler s Metod Ordiar Dieretial Euatio Ruge-Kutta metod Stabilit ad Talor series I. INTRODUCTION. Solutios to Dieretial euatio ave over te ears bee a ous to Applied Matematis. Te uestio te beomes ow to id te solutios to tose euatios. Te rage o Dieretial Euatios tat a be solved b straigt orward aaltial metod is relativel restrited []. Eve solutio i series ma ot alwas be satisator eiter beause o te slow overgee o te resultig series or beause o te ivolved maipulatio i repeated stages o dieretiatio [9].Ruge-Kutta metods are amog te most popular ODE solvers. Te were irst studied b Carle Ruge ad Marti Kutta aroud 9. Moder developmets are mostl due to Jo Buter i te 9s []. Te Ruge-Kutta metod is ot restrited to solvig ol irstorder dieretial euatios but Ruge-Kutta metods are also used to solve iger order ordiar dieretial euatios or oupled (simultaeous) dieretial euatios []. Te iger order euatios a be solved b osiderig a euivalet sstem o irst order euatios. However it is also possible to develop diret sigle steps metods to solve iger order euatio. Deiitio: We deie a epliit Ruge-Kutta metod wit slopes b te ollowig euatios []. b r i i i j! j i ij Were a i j r r... ij r are arbitrar. Iss 5-5 Jul Page
2 O Te Stabilit Ad Aura O Some II. DERIVATION OF SECOND ORDER RUNGE-KUTTA METHOD. To derive te Ruge-Kutta metods or seod order ordiar dieretial euatio o te orm Let us deie wit te iitial oditio..!!. r r. r r. Were b r r r ad r are arbitrar ostats to be determied. Usig Talor series epasio orm euatio. ad. gives. iv....5!!! iv....!!.7.8 iv [ ].9 We re-write euatios.7.8 ad.9 as D iv D D Were D So euatios.5 ad. beomes Iss 5-5 Jul Page
3 O Te Stabilit Ad Aura O Some D D D....!!! D D D.... Simpliig! we ave!! O D D O. Were we ave used Te substitutio o ad. i euatio. ad. ield 5 r r r D r D r O. r r r D r D r O.5 We ompare te oeiiet o euatios r r r r. ad.5 wit euatios. ad. r r Te oeiiet o i ad i orrespodig oeiiets i euatio. truatio error is O i ad.to obtai.. o euatios. ad.5 ad. will ot mat wit te or a oie o r ad r. Tus te loal O i. A solutio o euatio. ad. osidered or dieret values o r. Here we osider two ases o r. ma be Case I: r Te rom euatio. we get r b r r I te utio is idepedet o te we a ostrut a Ruge-Kutta metod i wi te loal truatio error i ad is O. Here we obtai Iss 5-5 Jul Page
4 O Te Stabilit Ad Aura O Some r r r r wit te solutio r r r r r r r Tus te Ruge-utta metod or te seod order iitial value problem Tus rom euatio wit iitial value... ad. we ave! ( ) Case II: r From euatio. we ave! r r r b 9 Tereore!!.9 9. Te two metods derived are epliit seod order Ruge-Kutta. III. STABILITY ANALYSIS. Wile umeriall solvig a iitial value problem or ordiar dieretial euatios a error is itrodued at ea itegratio step due to te iaura o te ormula. Eve we te loal error at ea step is small te total error ma beome large due to aumulatio ad ampliiatio o tese loal errors. Tis growt peomeo is alled umerial istabilit []. We sall disuss te stabilit o te Ruge-Kutta metod i.8 ad. Let us osider te dieretial euatio. Subjet to iitial oditio b Were is a real umber Iss 5-5 Jul Page
5 We sall osider te ase order ordiar dieretial euatio is o te orm O Te Stabilit Ad Aura O Some ad wile give a trival solutio. Reall tat te seod Te rom euatio. Substitute tis. i euatio.7 we ave 9 Substitutig te euatio. ito euatio.8 we ave. 8. Let us ow osider te ase we. Te solutio i tis ase is osillatig. We tereore osider te eigevalues o te matri wi is give b I te. To determie te value o ad we osider.5 We a re-write euatio.5 as. we we ave Iss 5-5 Jul Page
6 O Te Stabilit Ad Aura O Some Iss 5-5 Jul Page Substitutig euatio.7 i. we ave Were ad 578. Computig ad as utios o we id tat te roots ave uit modulus.. Tus te stabilit iterval o te Ruge-Kutta metod.8 is.. wi also sow tat is o order. For te solutios o. are epoetial i ature ad a be writte i te matri orm as e b Ad b or poit tis solutio beomes Ad te 7 lim e b Te relative error or te metod uder disussio i ase o a large umber o itegratio iterval (large small ) are to be osidered. Te maimum eigevalue o te matri.
7 O Te Stabilit Ad Aura O Some Iss 5-5 Jul Page 8 Is obtaied as... Ad tereore te relative error is give b 7 log log log log e e F We sall ow osider te stabilit o te Ruge-Kutta metod. usig dieretial euatio.. Euatio.9 beomes.8 Substitutig te euatio.9 ito euatio. we ave.9 We ow osider te ase we. Te eigevalues is give b Were ad. Computig ad as utios o te roots ave uit modulus 9 5. Tus te stabilit iterval o te Ruge-Kutta metod. is 9 5. ad is o order. For te solutios o. are epoetial i ature as above i ase I.
8 O Te Stabilit Ad Aura O Some Te maimum eigevalue o te matri. Is obtaied as 7... Te relative error is give b F log log e log e log 8 IV. IMPLEMENTATION OF THE METHOD. Cosider te iitial value problem [] Subjeted to iitial oditio.. Te omplemetar solutio is iitial oditio we ave A ad B Te geeral solutio beome A Be ad te partiular solutio gives p Usig te e Numerial solutios are preerred to te derived ases I ad ase II o epliit seod order Ruge-Kutta metod o te IVP i. obtaiig umerial solutios or values o up to ad iludig wit a step size o. as oud i Table ad Table. Table: Solutios o ase I seod ordered Ruge-Kutta metods wit. X Numerial solutio () Aaltial solutio Y() Absolute Error Iss 5-5 Jul Page 7
9 O Te Stabilit Ad Aura O Some Figure :.. Table: Solutios o ase II seod ordered Ruge-Kutta metods wit. X Numerial solutio () Aaltial solutio Y() Absolute Error Figure : Euler s metod Euler s metod is oe o ma metods or geeratig umerial solutios to dieretial euatios. Besides most o te oter metods tat migt be disussed are reiemets o Euler s metod. Tis metod is implemeted ad ompared its aura ad te error wit te metod i setio. Te Euler s metod geeralized i te orm [] Cosider te iitial value problem i euatio. Subjeted to iitial oditio. Table: Solutios o Euler s metods wit. X Numerial solutio () Aaltial solutio Y() Absolute Error Iss 5-5 Jul Page 8
10 O Te Stabilit Ad Aura O Some Figure : V. CONCLUSION. Ivestigatio arried out o some epliit seod ordered Ruge-Kutta metod i tis paper as sow tat te stabilit iterval o te Ruge-Kutta metod i ase I is. ad te stabilit iterval o te Ruge-Kutta metod i ase II is It is lear tat ase I is more stable ta ase II o te derived Ruge-Kutta metods. Te two metods are sow to be aurate eiiet ad geeral i appliatio or suiietl solutio o. Te result obtaied i te preset wor demostrate te eetiveess ad superiorit or te solutio o seod order ordiar dieretial euatio wi gave a ver ig aura we ompared wit eat solutio. REFERENCE. [] Autar Kaw Ruge-Kutta Seod Order Metod or Ordiar Dieretial Euatios. Uiversit o Sout Florida; Holisti Numerial metods Istitute. [] Bu R.C. ad Bu E.F. Itrodutio to Dieretial Euatio; Hogto Mili Compa 97 [] Buter J.C. Numerial Metods or Ordiar Dieretial Euatios; Jo Wile ad sos New Yor. [] [Buter Jo Ruge-Kutta Metods or Ordiar Dieretial Euatios; COE Worsop o Numerial Aalsis Kusu Uiversit 5. [5] Cristoper P. Grat Teor o Ordiar Dieretial Euatios; Brigam Youg Uiversit. [] Corliss G. ad Kirliger O Impliit Talor Series metods or Sti Ordiar Dieretial Euatio; Cetre or Resear o Parallel Computatio Rie Uiversit Housto 99. [7] Curtis F.G. ad Patri O.W. Applied Numerial Aalsis; Addiso Wesle Publisig Compa 989. [8] Greespa D. Teor ad Solutios o Ordiar Dieretial Euatios; Te MaMilla Compa New Yor 9. [9] Iserles A. Numerial Aalsis or Ordiar Dieretial Euatios; Cambridge Uiversit Press Cambridge 99. [] Jai M.K. Numerial Solutio o Dieretial Euatios Wile Easter Limited. [] Joe D. Homa Numerial Metods or Egieers ad Sietists; marel Deer I. [] Retewald G. Numerial Itegratio o ODEs or Iitial Value Problems. Departmet o Meaial egieerig Portlad State Uiversit.. [] Ouuga S.A. ad Aabi M.A. Computatioal Matematis; A First Course. Wim Publiatio Lagos Nigeria. [] Patil P.B. ad Verma U.P. Numerial Computatioal Metods Narosa Publisig House New Deli. Iss 5-5 Jul Page 9
Study on Solution of Non-homogeneous Linear Equation based on Ordinary Differential Equation Driving Jing Zhang
Iteratioal Coeree o Automatio Meaial Cotrol ad Computatioal Egieerig AMCCE 05 Stud o Solutio o No-omogeeous Liear Equatio based o Ordiar Dieretial Equatio Drivig Jig Zag Meaial ad Eletrial Egieerig College
More informationNUMERICAL 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 information1. Introduction. 2. Numerical Methods
America Joural o Computatioal ad Applied Matematics, (5: 9- DOI:.59/j.ajcam.5. A Stud o Numerical Solutios o Secod Order Iitial Value Problems (IVP or Ordiar Dieretial Equatios wit Fourt Order ad Butcer
More informationd 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 informationAn Improved Self-Starting Implicit Hybrid Method
IOSR Joural o Matematics (IOSR-JM e-issn: 78-78, p-issn:9-76x. Volume 0, Issue Ver. II (Mar-Apr. 04, PP 8-6 www.iosrourals.org A Improved Sel-Startig Implicit Hbrid Metod E. O. Adeea Departmet o Matematics/Statistics,
More informationSolving third order boundary value problem with fifth order block method
Matematical Metods i Egieerig ad Ecoomics Solvig tird order boudary value problem wit it order bloc metod A. S. Abdulla, Z. A. Majid, ad N. Seu Abstract We develop a it order two poit bloc metod or te
More informationA NUMERICAL METHOD OF SOLVING CAUCHY PROBLEM FOR DIFFERENTIAL EQUATIONS BASED ON A LINEAR APPROXIMATION
U.P.B. Sci. Bull., Series A, Vol. 79, Iss. 4, 7 ISSN -77 A NUMERICAL METHOD OF SOLVING CAUCHY PROBLEM FOR DIFFERENTIAL EQUATIONS BASED ON A LINEAR APPROXIMATION Cristia ŞERBĂNESCU, Marius BREBENEL A alterate
More informationComputational Methods CMSC/AMSC/MAPL 460. Quadrature: Integration
Computatioal Metods CMSC/AMSC/MAPL 6 Quadrature: Itegratio Ramai Duraiswami, Dept. o Computer Siee Some material adapted rom te olie slides o Eri Sadt ad Diae O Leary Numerial Itegratio Idea is to do itegral
More informationOn Exact Finite-Difference Scheme for Numerical Solution of Initial Value Problems in Ordinary Differential Equations.
O Exact Fiite-Differece Sceme for Numerical Solutio of Iitial Value Problems i Ordiar Differetial Equatios. Josua Suda, M.Sc. Departmet of Matematical Scieces, Adamawa State Uiversit, Mubi, Nigeria. E-mail:
More informationChapter 2: Numerical Methods
Chapter : Numerical Methods. Some Numerical Methods for st Order ODEs I this sectio, a summar of essetial features of umerical methods related to solutios of ordiar differetial equatios is give. I geeral,
More informationAfter the completion of this section the student. V.4.2. Power Series Solution. V.4.3. The Method of Frobenius. V.4.4. Taylor Series Solution
Chapter V ODE V.4 Power Series Solutio Otober, 8 385 V.4 Power Series Solutio Objetives: After the ompletio of this setio the studet - should reall the power series solutio of a liear ODE with variable
More informationSolving Third Order Boundary Value Problem Using. Fourth Order Block Method
Applied Matematical Scieces, Vol. 7,, o. 5, 69-65 HIKARI Ltd, www.m-ikari.com Solvig Tird Order Boudar Value Problem Usig Fourt Order Block Metod Amad Sa Abdulla, *Zaaria Abdul Maid, ad Norazak Seu, Istitute
More information(8) 1f = f. can be viewed as a real vector space where addition is defined by ( a1+ bi
Geeral Liear Spaes (Vetor Spaes) ad Solutios o ODEs Deiitio: A vetor spae V is a set, with additio ad salig o elemet deied or all elemets o the set, that is losed uder additio ad salig, otais a zero elemet
More informationOn the convergence, consistence and stability of a standard finite difference scheme
AMERICAN JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH 2, Sciece Huβ, ttp://www.sciub.org/ajsir ISSN: 253-649X, doi:.525/ajsir.2.2.2.74.78 O te covergece, cosistece ad stabilit of a stadard fiite differece
More informationOn the Derivation and Implementation of a Four Stage Harmonic Explicit Runge-Kutta Method *
Applied Mathematics, 05, 6, 694-699 Published Olie April 05 i SciRes. http://www.scirp.org/joural/am http://dx.doi.org/0.46/am.05.64064 O the Derivatio ad Implemetatio of a Four Stage Harmoic Explicit
More informationx x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula
NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS. If g( is cotiuous i [a,b], te uder wat coditio te iterative (or iteratio metod = g( as a uique solutio i [a,b]? '( i [a,b].. Wat
More informationA Class of Blended Block Second Derivative Multistep Methods for Stiff Systems
Iteratioal Joural of Iovative Mathematics, Statistics & Eerg Policies ():-6, Ja.-Mar. 7 SEAHI PUBLICATIONS, 7 www.seahipa.org ISSN: 67-8X A Class of Bleded Bloc Secod Derivative Multistep Methods for Stiff
More informationSolutions 3.2-Page 215
Solutios.-Page Problem Fid the geeral solutios i powers of of the differetial equatios. State the reurree relatios ad the guarateed radius of overgee i eah ase. ) Substitutig,, ad ito the differetial equatio
More informationA Self-Starting Hybrid Linear Multistep Method for a Direct Solution of the General Second-Order Initial Value Problem
IOS Joural o Matematics (IOS-JM) ISSN: 78-578. Volume 4 Issue 6 (Ja. - eb. ) PP 7- www.iosrourals.org Sel-Startig Hbrid Liear Multistep Metod or a Direct Solutio o te Geeral Secod-Order Iitial Value Problem
More informationPRELIMINARY MATHEMATICS LECTURE 5
SCHOOL OF ORIENTAL AND AFRICAN STUDIES UNIVERSITY OF LONDON DEPARTMENT OF ECONOMICS 5 / - 6 5 MSc Ecoomics PRELIMINARY MATHEMATICS LECTURE 5 Course website: http://mercur.soas.ac.uk/users/sm97/teachig_msc_premath.htm
More informationThe Jordan Normal Form: A General Approach to Solving Homogeneous Linear Systems. Mike Raugh. March 20, 2005
The Jorda Normal Form: A Geeral Approach to Solvig Homogeeous Liear Sstems Mike Raugh March 2, 25 What are we doig here? I this ote, we describe the Jorda ormal form of a matrix ad show how it ma be used
More informationA New Hybrid in the Nonlinear Part of Adomian Decomposition Method for Initial Value Problem of Ordinary Differential Equation
Joural of Matematics Researc; Vol No ; ISSN - E-ISSN - Publised b Caadia Ceter of Sciece ad Educatio A New Hbrid i te Noliear Part of Adomia Decompositio Metod for Iitial Value Problem of Ordiar Differetial
More informationBernoulli Numbers. n(n+1) = n(n+1)(2n+1) = n(n 1) 2
Beroulli Numbers Beroulli umbers are amed after the great Swiss mathematiia Jaob Beroulli5-705 who used these umbers i the power-sum problem. The power-sum problem is to fid a formula for the sum of the
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
More informationSummation Method for Some Special Series Exactly
The Iteratioal Joural of Mathematis, Siee, Tehology ad Maagemet (ISSN : 39-85) Vol. Issue Summatio Method for Some Speial Series Eatly D.A.Gismalla Deptt. Of Mathematis & omputer Studies Faulty of Siee
More informationFamurewa O. K. E*, Ademiluyi R. A. and Awoyemi D. O.
Africa Joural of Matematics ad omputer Sciece Researc Vol. (), pp. -, Marc Available olie at ttp://www.academicourals.org/ajmsr ISSN 6-97 Academic Jourals Full Legt Researc Paper A comparative stud of
More informationNumerical Derivatives by Symbolic Tools in MATLAB
Numerical Derivatives by Symbolic ools i MALAB Mig-Gog Lee * ad Rei-Wei Sog ad Hsuag-Ci Cag mglee@cu.edu.tw * Departmet o Applied Matematics Cug Hua Uiversity Hsicu, aiwa Abstract: Numerical approaces
More informationME 501A Seminar in Engineering Analysis Page 1
Accurac, Stabilit ad Sstems of Equatios November 0, 07 Numerical Solutios of Ordiar Differetial Equatios Accurac, Stabilit ad Sstems of Equatios Larr Caretto Mecaical Egieerig 0AB Semiar i Egieerig Aalsis
More informationALLOCATING SAMPLE TO STRATA PROPORTIONAL TO AGGREGATE MEASURE OF SIZE WITH BOTH UPPER AND LOWER BOUNDS ON THE NUMBER OF UNITS IN EACH STRATUM
ALLOCATING SAPLE TO STRATA PROPORTIONAL TO AGGREGATE EASURE OF SIZE WIT BOT UPPER AND LOWER BOUNDS ON TE NUBER OF UNITS IN EAC STRATU Lawrece R. Erst ad Cristoper J. Guciardo Erst_L@bls.gov, Guciardo_C@bls.gov
More informationCS321. Numerical Analysis and Computing
CS Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 September 8 5 What is the Root May physical system ca
More informationCS537. Numerical Analysis and Computing
CS57 Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 Jauary 9 9 What is the Root May physical system ca be
More informationA Pseudo Spline Methods for Solving an Initial Value Problem of Ordinary Differential Equation
Joural of Matematics ad Statistics 4 (: 7-, 008 ISSN 549-3644 008 Sciece Publicatios A Pseudo Splie Metods for Solvig a Iitial Value Problem of Ordiary Differetial Equatio B.S. Ogudare ad G.E. Okeca Departmet
More informationComputation Sessional. Numerical Differentiation and Integration
CE 6: Egieerig Computatio Sessioal Numerical Dieretiatio ad Itegratio ti di ad gradiet commad di() Returs te dierece betwee adjacet elemets i. Typically used or uequally spaced itervals = gradiet(, ) Determies
More informationENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 4 Solutions [Numerical Methods]
ENGI 3 Advaced Calculus or Egieerig Facult o Egieerig ad Applied Sciece Problem Set Solutios [Numerical Methods]. Use Simpso s rule with our itervals to estimate I si d a, b, h a si si.889 si 3 si.889
More informationNewton s Method. y f x 1 x x 1 f x 1. Letting y 0 and solving for x produces. x x 1 f x 1. x 1. x 2 x 1 f x 1. f x 1. x 3 x 2 f x 2 f x 2.
460_008.qd //04 :7 PM Page 9 SECTION.8 Newto s Method 9 (a) a a Sectio.8 (, ( )) (, ( )) Taget lie c Taget lie c b (b) The -itercept o the taget lie approimates the zero o. Figure.60 b Newto s Method Approimate
More informationNewton s Method. Video
SECTION 8 Newto s Method 9 (a) a a Sectio 8 (, ( )) (, ( )) Taget lie c Taget lie c b (b) The -itercept o the taget lie approimates the zero o Figure 60 b Newto s Method Approimate a zero o a uctio usig
More informationENGI 9420 Lecture Notes 3 - Numerical Methods Page 3.01
ENGI 940 Lecture Notes 3 - Numerical Methods Page 3.0 3. Numerical Methods The majority of equatios of iterest i actual practice do ot admit ay aalytic solutio. Eve equatios as simple as = e ad I = e d
More informationMore Elementary Aspects of Numerical Solutions of PDEs!
ttp://www.d.edu/~gtryggva/cfd-course/ Outlie More Elemetary Aspects o Numerical Solutios o PDEs I tis lecture we cotiue to examie te elemetary aspects o umerical solutios o partial dieretial equatios.
More information2.3 Warmup. Graph the derivative of the following functions. Where necessary, approximate the derivative.
. Warmup Grap te erivative of te followig fuctios. Were ecessar, approimate te erivative. Differetiabilit Must a fuctio ave a erivative at eac poit were te fuctio is efie? Or If f a is efie, must f ( a)
More informationComputation of Hahn Moments for Large Size Images
Joural of Computer Sciece 6 (9): 37-4, ISSN 549-3636 Sciece Publicatios Computatio of Ha Momets for Large Size Images A. Vekataramaa ad P. Aat Raj Departmet of Electroics ad Commuicatio Egieerig, Quli
More informationThird-order Composite Runge Kutta Method for Solving Fuzzy Differential Equations
Global Joural of Pure ad Applied Mathematics. ISSN 097-768 Volume Number (06) pp. 7-76 Research Idia Publicatios http://www.ripublicatio.com/gjpam.htm Third-order Composite Ruge Kutta Method for Solvig
More informationStability analysis of numerical methods for stochastic systems with additive noise
Stability aalysis of umerical metods for stoctic systems wit additive oise Yosiiro SAITO Abstract Stoctic differetial equatios (SDEs) represet pysical peomea domiated by stoctic processes As for determiistic
More informationModification of Weerakoon-Fernando s Method with Fourth-Order of Convergence for Solving Nonlinear Equation
ISSN: 50-08 Iteratioal Joural o AdvacedResearch i Sciece, Egieerig ad Techology Vol. 5, Issue 8, August 018 Modiicatio o Weerakoo-Ferado s Method with Fourth-Order o Covergece or Solvig Noliear Equatio
More informationExplicit and closed formed solution of a differential equation. Closed form: since finite algebraic combination of. converges for x x0
Chapter 4 Series Solutios Epliit ad losed formed solutio of a differetial equatio y' y ; y() 3 ( ) ( 5 e ) y Closed form: sie fiite algebrai ombiatio of elemetary futios Series solutio: givig y ( ) as
More informationThe Advection-Diffusion equation!
ttp://www.d.edu/~gtryggva/cf-course/! Te Advectio-iffusio equatio! Grétar Tryggvaso! Sprig 3! Navier-Stokes equatios! Summary! u t + u u x + v u y = P ρ x + µ u + u ρ y Hyperbolic part! u x + v y = Elliptic
More informationTopic 9 - Taylor and MacLaurin Series
Topic 9 - Taylor ad MacLauri Series A. Taylors Theorem. The use o power series is very commo i uctioal aalysis i act may useul ad commoly used uctios ca be writte as a power series ad this remarkable result
More informationNumerical Integration Formulas
Numerical Itegratio Formulas Berli Che Departmet o Computer Sciece & Iormatio Egieerig Natioal Taiwa Normal Uiversity Reerece: 1. Applied Numerical Methods with MATLAB or Egieers, Chapter 19 & Teachig
More informationTaylor Polynomials and Approximations - Classwork
Taylor Polyomials ad Approimatios - Classwork Suppose you were asked to id si 37 o. You have o calculator other tha oe that ca do simple additio, subtractio, multiplicatio, or divisio. Fareched\ Not really.
More informationSome Variants of Newton's Method with Fifth-Order and Fourth-Order Convergence for Solving Nonlinear Equations
Copyright, Darbose Iteratioal Joural o Applied Mathematics ad Computatio Volume (), pp -6, 9 http//: ijamc.darbose.com Some Variats o Newto's Method with Fith-Order ad Fourth-Order Covergece or Solvig
More informationNUMERICAL METHODS FOR SOLVING EQUATIONS
Mathematics Revisio Guides Numerical Methods for Solvig Equatios Page 1 of 11 M.K. HOME TUITION Mathematics Revisio Guides Level: GCSE Higher Tier NUMERICAL METHODS FOR SOLVING EQUATIONS Versio:. Date:
More informationA NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS
Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet
More informationCHAPTER 6d. NUMERICAL INTERPOLATION
CHAPER 6d. NUMERICAL INERPOLAION A. J. Clark School o Egieerig Departmet o Civil ad Evirometal Egieerig by Dr. Ibrahim A. Assakka Sprig ENCE - Computatio Methods i Civil Egieerig II Departmet o Civil ad
More informationOn error estimation in almost Runge-Kutta (ARK) methods
Joural of Applied Mathematics & Bioiformatics, vol.3, o., 3, 5-35 ISSN: 79-66 prit), 79-6939 olie) Sciepress Ltd, 3 O error estimatio i almost Ruge-Kutta ARK) methods Ochoche Abraham, Kaode R. Adeboe ad
More informationg () n = g () n () f, f n = f () n () x ( n =1,2,3, ) j 1 + j 2 + +nj n = n +2j j n = r & j 1 j 1, j 2, j 3, j 4 = ( 4, 0, 0, 0) f 4 f 3 3!
Higher Derivative o Compositio. Formulas o Higher Derivative o Compositio.. Faà di Bruo's Formula About the ormula o the higher derivative o compositio, the oe by a mathematicia Faà di Bruo i Italy o about
More informationwavelet collocation method for solving integro-differential equation.
IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 3 (arch. 5), V3 PP -7 www.iosrje.org wavelet collocatio method for solvig itegro-differetial equatio. Asmaa Abdalelah Abdalrehma
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationSIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road QUESTION BANK (DESCRIPTIVE)
QUESTION BANK 8 SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayaavaam Road 5758 QUESTION BANK (DESCRIPTIVE) Subject with Code : (6HS6) Course & Brach: B.Tech AG Year & Sem: II-B.Tech&
More informationImplicit function theorem
Jovo Jaric Implicit fuctio theorem The reader kows that the equatio of a curve i the x - plae ca be expressed F x, =., this does ot ecessaril represet a fuctio. Take, for example F x, = 2x x =. (1 either
More informationANOTHER PROOF FOR FERMAT S LAST THEOREM 1. INTRODUCTION
ANOTHER PROOF FOR FERMAT S LAST THEOREM Mugur B. RĂUŢ Correspodig author: Mugur B. RĂUŢ, E-mail: m_b_raut@yahoo.om Abstrat I this paper we propose aother proof for Fermat s Last Theorem (FLT). We foud
More informationWhere do eigenvalues/eigenvectors/eigenfunctions come from, and why are they important anyway?
Where do eigevalues/eigevectors/eigeuctios come rom, ad why are they importat ayway? I. Bacgroud (rom Ordiary Dieretial Equatios} Cosider the simplest example o a harmoic oscillator (thi o a vibratig strig)
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationMathematical Expressions for Estimation of Errors in the Formulas which are used to obtaining intermediate values of Biological Activity in QSAR
Iteratioal Joural o Iovatio ad Applied Studies ISSN 08-934 Vol. No. 3 Mar. 03, pp. 7-79 03 Iovative Space o Scietiic Researc Jourals ttp//www.issr-jourals.org/ijias/ Matematical Expressios or Estimatio
More informationPartial Differential Equations
EE 84 Matematical Metods i Egieerig Partial Differetial Eqatios Followig are some classical partial differetial eqatios were is assmed to be a fctio of two or more variables t (time) ad y (spatial coordiates).
More informationFluids Lecture 2 Notes
Fluids Leture Notes. Airfoil orte Sheet Models. Thi-Airfoil Aalysis Problem Readig: Aderso.,.7 Airfoil orte Sheet Models Surfae orte Sheet Model A aurate meas of represetig the flow about a airfoil i a
More informationREPRESENTING MARKOV CHAINS WITH TRANSITION DIAGRAMS
Joural o Mathematics ad Statistics, 9 (3): 49-54, 3 ISSN 549-36 3 Sciece Publicatios doi:.38/jmssp.3.49.54 Published Olie 9 (3) 3 (http://www.thescipub.com/jmss.toc) REPRESENTING MARKOV CHAINS WITH TRANSITION
More informationDEGENERACY AND ALL THAT
DEGENERACY AND ALL THAT Te Nature of Termodyamics, Statistical Mecaics ad Classical Mecaics Termodyamics Te study of te equilibrium bulk properties of matter witi te cotext of four laws or facts of experiece
More informationConvergence of the FEM. by Hui Zhang Jorida Kushova Ruwen Jung
Covergece o te FEM by Hi Zg Jorid Ksov Rwe Jg I order to proo FEM soltios to be coverget, mesremet or teir qlity is reqired. A simple pproc i ect soltio is ccessible is to qtiy te error betwee FEMd te
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 informationContinuous Random Variables: Conditioning, Expectation and Independence
Cotiuous Radom Variables: Coditioig, Expectatio ad Idepedece Berli Che Departmet o Computer ciece & Iormatio Egieerig Natioal Taiwa Normal Uiversit Reerece: - D.. Bertsekas, J. N. Tsitsiklis, Itroductio
More information5. DIFFERENTIAL EQUATIONS
5-5. DIFFERENTIAL EQUATIONS The most commo mathematical structure emploed i mathematical models of chemical egieerig professio ivolve differetial equatios. These equatios describe the rate of chage of
More informationPOWER SERIES METHODS CHAPTER 8 SECTION 8.1 INTRODUCTION AND REVIEW OF POWER SERIES
CHAPTER 8 POWER SERIES METHODS SECTION 8. INTRODUCTION AND REVIEW OF POWER SERIES The power series method osists of substitutig a series y = ito a give differetial equatio i order to determie what the
More informationThomas J. Osler Mathematics Department Rowan University Glassboro NJ Introduction
Ot 0 006 Euler s little summtio formul d speil vlues of te zet futio Toms J Osler temtis Deprtmet Row Uiversity Glssboro J 0608 Osler@rowedu Itrodutio I tis ote we preset elemetry metod of determiig vlues
More informationNumerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets
Australia Joural of Basic ad Applied Scieces, 5(): 98-5, ISSN 99-878 Numerical Solutio of the Two Poit Boudary Value Problems By Usig Wavelet Bases of Hermite Cubic Splie Wavelets Mehdi Yousefi, Hesam-Aldie
More informationNumerical Methods Control Test
IGEE/Promotio M Tuesday,6/5/4 Numerial Metods Cotrol Test E 7 poits: Fid te root o = usig Bisetio metod wit a error o e. E 6 poits: Use te ormula o Newto-Rapso metod irst order to solve: log E 7 poits:
More informationNumerical Solution of Non-Linear Ordinary Differential Equations via Collocation Method (Finite Elements) and Genetic Algorithms
Proceedigs of te 6t WSEAS It. Cof. o EVOLUTIONARY COPUTING Lisbo Portugal Jue 6-8 5 pp6- Numerical Solutio of No-Liear Ordiary Differetial Equatios via Collocatio etod Fiite Elemets ad Geetic Algoritms
More information1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations
. Liearizatio of a oliear system give i the form of a system of ordiary differetial equatios We ow show how to determie a liear model which approximates the behavior of a time-ivariat oliear system i a
More informationμ are complex parameters. Other
A New Numerical Itegrator for the Solutio of Iitial Value Problems i Ordiary Differetial Equatios. J. Suday * ad M.R. Odekule Departmet of Mathematical Scieces, Adamawa State Uiversity, Mubi, Nigeria.
More informationUsing Stirling s Interpolation to Find Gauss and Mean Curvature for the Surface
IOSR Joural o Mathematics (IOSR-JM ISSN: - Volume Issue (Nov - Dec PP - wwwiosrjouralsorg Usig Stirlig s Iterolatio to Fid Gauss ad Mea Curvature or the Surace R B M Amer M A Abd El- Mageed (Deartmet o
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationPolynomial Generalizations and Combinatorial Interpretations for Sequences Including the Fibonacci and Pell Numbers
Ope Joural o Discrete Mathematics,,, - http://dxdoiorg/46/odm6 Published Olie Jauary (http://wwwscirporg/oural/odm) Polyomial Geeralizatios ad Combiatorial Iterpretatios or Seueces Icludig the Fiboacci
More informationEstimation of Population Ratio in Post-Stratified Sampling Using Variable Transformation
Ope Joural o Statistics, 05, 5, -9 Published Olie Februar 05 i SciRes. http://www.scirp.org/joural/ojs http://dx.doi.org/0.436/ojs.05.500 Estimatio o Populatio Ratio i Post-Stratiied Samplig Usig Variable
More informationPhys. 201 Mathematical Physics 1 Dr. Nidal M. Ershaidat Doc. 12
Physics Departmet, Yarmouk Uiversity, Irbid Jorda Phys. Mathematical Physics Dr. Nidal M. Ershaidat Doc. Fourier Series Deiitio A Fourier series is a expasio o a periodic uctio (x) i terms o a iiite sum
More informationIJITE Vol.2 Issue-11, (November 2014) ISSN: Impact Factor
IJITE Vol Issue-, (November 4) ISSN: 3-776 ATTRACTIVITY OF A HIGHER ORDER NONLINEAR DIFFERENCE EQUATION Guagfeg Liu School of Zhagjiagag Jiagsu Uiversit of Sciece ad Techolog, Zhagjiagag, Jiagsu 56,PR
More informationPAPER : IIT-JAM 2010
MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure
More informationMATH CALCULUS II Objectives and Notes for Test 4
MATH 44 - CALCULUS II Objectives ad Notes for Test 4 To do well o this test, ou should be able to work the followig tpes of problems. Fid a power series represetatio for a fuctio ad determie the radius
More informationNumerical Methods for Ordinary Differential Equations
Numerical Methods for Ordiary Differetial Equatios Braislav K. Nikolić Departmet of Physics ad Astroomy, Uiversity of Delaware, U.S.A. PHYS 460/660: Computatioal Methods of Physics http://www.physics.udel.edu/~bikolic/teachig/phys660/phys660.html
More informationBlind Adaptive Filters & Equalization
Blid Adaptive Filters & Equalizatio Table of otets: Itrodutio about Adaptive Filters Itrodutio about Blid Adaptive Filters Covergee aalysis of LMS algoritm Ho trasforms improve te overgee rate of LMS Wy
More informationConstruction of Control Chart for Random Queue Length for (M / M / c): ( / FCFS) Queueing Model Using Skewness
Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember ISSN 5-5 Costrutio of Cotrol Chart for Radom Queue Legth for (M / M / ): ( / FCFS) Queueig Model Usig Skewess Dr.(Mrs.) A.R. Sudamai
More informationWhen dealing with series, n is always a positive integer. Remember at every, sine has a value of zero, which means
Fourier Series Some Prelimiar Ideas: Odd/Eve Fuctios: Sie is odd, which meas si ( ) si Cosie is eve, which meas cos ( ) cos Secial values of siie a cosie at Whe dealig with series, is alwas a ositive iteger.
More informationC/CS/Phys C191 Deutsch and Deutsch-Josza algorithms 10/20/07 Fall 2007 Lecture 17
C/CS/Phs C9 Deutsch ad Deutsch-Josza algorithms 0/0/07 Fall 007 Lecture 7 Readigs Beeti et al., Ch. 3.9-3.9. Stolze ad Suter, Quatum Computig, Ch. 8. - 8..5) Nielse ad Chuag, Quatum Computatio ad Quatum
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationSubject: Differential Equations & Mathematical Modeling -III. Lesson: Power series solutions of Differential Equations. about ordinary points
Power series solutio of Differetial equatios about ordiary poits Subject: Differetial Equatios & Mathematical Modelig -III Lesso: Power series solutios of Differetial Equatios about ordiary poits Lesso
More informationFurther Methods for Advanced Mathematics (FP2) WEDNESDAY 9 JANUARY 2008
ADVANCED GCE 7/ MATHEMATICS (MEI) Furter Metods for Advaced Matematics (F) WEDNESDAY 9 JANUARY 8 Additioal materials: Aswer Booklet (8 pages) Grap paper MEI Eamiatio Formulae ad Tables (MF) Afteroo Time:
More informationDetermination of Energy Involved In a Stepwise Size Reduction of Maize, Using a Numerical Approach.
IOSR Joural of Matematics (IOSR-JM) e-issn: 78-578, p-issn: 9-765X. Volume, Issue Ver. V (Ja - Feb. 5), PP 6-65 www.iosrourals.org Determiatio of ergy Ivolved I a Stepwise Size Reductio of Maize, Usig
More informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
More informationECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More informationAccuracy. Computational Fluid Dynamics. Computational Fluid Dynamics. Computational Fluid Dynamics
http://www.d.edu/~gtryggva/cfd-course/ Computatioal Fluid Dyamics Lecture Jauary 3, 7 Grétar Tryggvaso It is clear that although the umerical solutio is qualitatively similar to the aalytical solutio,
More informationMaclaurin and Taylor series
At the ed o the previous chapter we looed at power series ad oted that these were dieret rom other iiite series as they were actually uctios o a variable R: a a + + a + a a Maclauri ad Taylor series +
More informationSOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.
SOLUTIONS TO EXAM 3 Problem Fid the sum of the followig series 2 + ( ) 5 5 2 5 3 25 2 2 This series diverges Solutio: Note that this defies two coverget geometric series with respective radii r 2/5 < ad
More informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More information