On The Stability and Accuracy of Some Runge-Kutta Methods of Solving Second Order Ordinary Differential Equations

Size: px
Start display at page:

Download "On The Stability and Accuracy of Some Runge-Kutta Methods of Solving Second Order Ordinary Differential Equations"

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

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 information

NUMERICAL DIFFERENTIAL 1

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

More information

1. Introduction. 2. Numerical Methods

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

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

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

More information

An Improved Self-Starting Implicit Hybrid Method

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

Solving third order boundary value problem with fifth order block method

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

A NUMERICAL METHOD OF SOLVING CAUCHY PROBLEM FOR DIFFERENTIAL EQUATIONS BASED ON A LINEAR APPROXIMATION

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

Computational Methods CMSC/AMSC/MAPL 460. Quadrature: Integration

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

On Exact Finite-Difference Scheme for Numerical Solution of Initial Value Problems in Ordinary Differential Equations.

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

Chapter 2: Numerical Methods

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

After 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

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

Solving Third Order Boundary Value Problem Using. Fourth Order Block Method

Solving 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

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

On the convergence, consistence and stability of a standard finite difference scheme

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

On the Derivation and Implementation of a Four Stage Harmonic Explicit Runge-Kutta Method *

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

x x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula

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

A Class of Blended Block Second Derivative Multistep Methods for Stiff Systems

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

Solutions 3.2-Page 215

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

A Self-Starting Hybrid Linear Multistep Method for a Direct Solution of the General Second-Order Initial Value Problem

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

PRELIMINARY MATHEMATICS LECTURE 5

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

The Jordan Normal Form: A General Approach to Solving Homogeneous Linear Systems. Mike Raugh. March 20, 2005

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

A New Hybrid in the Nonlinear Part of Adomian Decomposition Method for Initial Value Problem of Ordinary Differential Equation

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

Bernoulli Numbers. n(n+1) = n(n+1)(2n+1) = n(n 1) 2

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

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t

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

Summation Method for Some Special Series Exactly

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

Famurewa O. K. E*, Ademiluyi R. A. and Awoyemi D. O.

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

Numerical Derivatives by Symbolic Tools in MATLAB

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

ME 501A Seminar in Engineering Analysis Page 1

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

ALLOCATING SAMPLE TO STRATA PROPORTIONAL TO AGGREGATE MEASURE OF SIZE WITH BOTH UPPER AND LOWER BOUNDS ON THE NUMBER OF UNITS IN EACH STRATUM

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

CS321. Numerical Analysis and Computing

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

CS537. Numerical Analysis and Computing

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

A Pseudo Spline Methods for Solving an Initial Value Problem of Ordinary Differential Equation

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

Computation Sessional. Numerical Differentiation and Integration

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

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 4 Solutions [Numerical Methods]

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

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

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

Newton s Method. Video

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

ENGI 9420 Lecture Notes 3 - Numerical Methods Page 3.01

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

More Elementary Aspects of Numerical Solutions of PDEs!

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

2.3 Warmup. Graph the derivative of the following functions. Where necessary, approximate the derivative.

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

Computation of Hahn Moments for Large Size Images

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

Third-order Composite Runge Kutta Method for Solving Fuzzy Differential Equations

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

Stability analysis of numerical methods for stochastic systems with additive noise

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

Modification of Weerakoon-Fernando s Method with Fourth-Order of Convergence for Solving Nonlinear Equation

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

Explicit and closed formed solution of a differential equation. Closed form: since finite algebraic combination of. converges for x x0

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

The Advection-Diffusion equation!

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

Topic 9 - Taylor and MacLaurin Series

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

Numerical Integration Formulas

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

Taylor Polynomials and Approximations - Classwork

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

Some Variants of Newton's Method with Fifth-Order and Fourth-Order Convergence for Solving Nonlinear Equations

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

NUMERICAL METHODS FOR SOLVING EQUATIONS

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

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

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

CHAPTER 6d. NUMERICAL INTERPOLATION

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

On error estimation in almost Runge-Kutta (ARK) methods

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

g () 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!

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

wavelet collocation method for solving integro-differential equation.

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

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

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

SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road QUESTION BANK (DESCRIPTIVE)

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

Implicit function theorem

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

ANOTHER PROOF FOR FERMAT S LAST THEOREM 1. INTRODUCTION

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

Where do eigenvalues/eigenvectors/eigenfunctions come from, and why are they important anyway?

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

Chapter 4. Fourier Series

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

Mathematical Expressions for Estimation of Errors in the Formulas which are used to obtaining intermediate values of Biological Activity in QSAR

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

Partial Differential Equations

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

Fluids Lecture 2 Notes

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

REPRESENTING MARKOV CHAINS WITH TRANSITION DIAGRAMS

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

DEGENERACY AND ALL THAT

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

Convergence of the FEM. by Hui Zhang Jorida Kushova Ruwen Jung

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

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

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

Continuous Random Variables: Conditioning, Expectation and Independence

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

5. DIFFERENTIAL EQUATIONS

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

POWER SERIES METHODS CHAPTER 8 SECTION 8.1 INTRODUCTION AND REVIEW OF POWER SERIES

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

Thomas J. Osler Mathematics Department Rowan University Glassboro NJ Introduction

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

Numerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets

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

Numerical Methods Control Test

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

Numerical Solution of Non-Linear Ordinary Differential Equations via Collocation Method (Finite Elements) and Genetic Algorithms

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

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

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

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

Using Stirling s Interpolation to Find Gauss and Mean Curvature for the Surface

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

Chapter 9: Numerical Differentiation

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

Polynomial Generalizations and Combinatorial Interpretations for Sequences Including the Fibonacci and Pell Numbers

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

Estimation of Population Ratio in Post-Stratified Sampling Using Variable Transformation

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

Phys. 201 Mathematical Physics 1 Dr. Nidal M. Ershaidat Doc. 12

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

IJITE Vol.2 Issue-11, (November 2014) ISSN: Impact Factor

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

PAPER : IIT-JAM 2010

PAPER : 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 information

MATH CALCULUS II Objectives and Notes for Test 4

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

Numerical Methods for Ordinary Differential Equations

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

Blind Adaptive Filters & Equalization

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

Construction of Control Chart for Random Queue Length for (M / M / c): ( / FCFS) Queueing Model Using Skewness

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

When dealing with series, n is always a positive integer. Remember at every, sine has a value of zero, which means

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

C/CS/Phys C191 Deutsch and Deutsch-Josza algorithms 10/20/07 Fall 2007 Lecture 17

C/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 information

Taylor Series and the Mean Value Theorem of Derivatives

Taylor Series and the Mean Value Theorem of Derivatives 1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential

More information

Subject: Differential Equations & Mathematical Modeling -III. Lesson: Power series solutions of Differential Equations. about ordinary points

Subject: 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 information

Further Methods for Advanced Mathematics (FP2) WEDNESDAY 9 JANUARY 2008

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

Determination of Energy Involved In a Stepwise Size Reduction of Maize, Using a Numerical Approach.

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

Castiel, Supernatural, Season 6, Episode 18

Castiel, 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 information

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

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

Accuracy. Computational Fluid Dynamics. Computational Fluid Dynamics. Computational Fluid Dynamics

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

Maclaurin and Taylor series

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

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

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Similarity 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