The Derivation of Implicit Second Derivative Method for Solving Second-Order Stiff Ordinary Differential Equations Odes.

Similar documents
CURVE FITTING LEAST SQUARES METHOD

UNIVERSITI KEBANGSAAN MALAYSIA PEPERIKSAAN AKHIR SEMESTER I SESI AKADEMIK 2007/2008 IJAZAH SARJANAMUDA DENGAN KEPUJIAN NOVEMBER 2007 MASA : 3 JAM

Differentiation and Numerical Integral of the Cubic Spline Interpolation

FORMULAE FOR FINITE DIFFERENCE APPROXIMATIONS, QUADRATURE AND LINEAR MULTISTEP METHODS

CS321. Introduction to Numerical Methods

Outline. Finite Difference Grids. Numerical Analysis. Finite Difference Grids II. Finite Difference Grids III

Mathematically, integration is just finding the area under a curve from one point to another. It is b

Chapter Simpson s 1/3 Rule of Integration. ( x)

Basic Concepts in Numerical Analysis November 6, 2017

Review of Linear Algebra

Outline. Remaining Course Schedule. Review Systems of ODEs. Example. Example Continued. Other Algorithms for Ordinary Differential Equations

Some Unbiased Classes of Estimators of Finite Population Mean

Differential Entropy 吳家麟教授

12 Iterative Methods. Linear Systems: Gauss-Seidel Nonlinear Systems Case Study: Chemical Reactions

4. Runge-Kutta Formula For Differential Equations. A. Euler Formula B. Runge-Kutta Formula C. An Example for Fourth-Order Runge-Kutta Formula

Chapter Runge-Kutta 2nd Order Method for Ordinary Differential Equations

APPLICATION OF THE CHEBYSHEV POLYNOMIALS TO APPROXIMATION AND CONSTRUCTION OF MAP PROJECTIONS

CS321. Numerical Analysis

More Regression Lecture Notes CE 311K - McKinney Introduction to Computer Methods Department of Civil Engineering The University of Texas at Austin

Integration by Parts for D K

4. Runge-Kutta Formula For Differential Equations

An Alternative Method to Find the Solution of Zero One Integer Linear Fractional Programming Problem with the Help of -Matrix

INTERPOLATION(2) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek

Spectrum Estimation by Several Interpolation Methods

CS321. Numerical Analysis

Chapter 12-b Integral Calculus - Extra

DATA FITTING. Intensive Computation 2013/2014. Annalisa Massini

6.6 Moments and Centers of Mass

CHAPTER 6 CURVE FITTINGS

A New Method for Solving Fuzzy Linear. Programming by Solving Linear Programming

Fibonacci and Lucas Numbers as Tridiagonal Matrix Determinants

An Extended Mixture Inverse Gaussian Distribution

Numerical Differentiation and Integration

ICS141: Discrete Mathematics for Computer Science I

A Comparison of Chebyshev polynomials and Legendre polynomials in order to solving Fredholm integral equations

In Calculus I you learned an approximation method using a Riemann sum. Recall that the Riemann sum is

The Distribution of Minimizing Maximum Entropy: Alternative to Weibull distribution for wind speed

A Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions

MTH 146 Class 7 Notes

. The set of these sums. be a partition of [ ab, ]. Consider the sum f( x) f( x 1)

Numerical Solution of Higher Order Linear Fredholm Integro Differential Equations.

On Solution of Min-Max Composition Fuzzy Relational Equation

Numerical Analysis Topic 4: Least Squares Curve Fitting

Numerical Analysis Third Class Chemical Engineering Department University of Technology Assoc. Prof. Dr. Zaidoon Mohsin shakor

Fourth Order Four-Stage Diagonally Implicit Runge-Kutta Method for Linear Ordinary Differential Equations ABSTRACT INTRODUCTION

Analytical Approach for the Solution of Thermodynamic Identities with Relativistic General Equation of State in a Mixture of Gases

Algorithms behind the Correlation Setting Window

Observations on the transcendental Equation

POWERS OF COMPLEX PERSYMMETRIC ANTI-TRIDIAGONAL MATRICES WITH CONSTANT ANTI-DIAGONALS

NUMERICAL SOLUTIONS OF SOME KINDS OF FRACTIONAL DIFFERENTIAL EQUATIONS

The z-transform. LTI System description. Prof. Siripong Potisuk

Chapter Gauss-Seidel Method

Numerical Differentiation

A New Efficient Approach to Solve Multi-Objective Transportation Problem in the Fuzzy Environment (Product approach)

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ. Численные методы. Учебно-методическое пособие

Numerical Differentiation

ON NILPOTENCY IN NONASSOCIATIVE ALGEBRAS

i+1 by A and imposes Ax

this is the indefinite integral Since integration is the reverse of differentiation we can check the previous by [ ]

3/20/2013. Splines There are cases where polynomial interpolation is bad overshoot oscillations. Examplef x. Interpolation at -4,-3,-2,-1,0,1,2,3,4

SUM PROPERTIES FOR THE K-LUCAS NUMBERS WITH ARITHMETIC INDEXES

Chapter 1. Infinite Sequences and Series. 1.1 Sequences. A sequence is a set of numbers written in a definite order

Accelerated Solution of High Order Non-linear ODEs using Chebyshev Spectral Method Comparing with Adomian Decomposition Method

Numerical differentiation (Finite difference)

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

Derivation of the Modified Bi-quintic B-spline Base Functions: An Application to Poisson Equation

Construction of Canonical Polynomial Basis Functions for Solving Special N th -Order Linear Integro-Differential Equations

Calculating the Values of Multiple Integrals by Replacing the Integrands Interpolation by Interpolation Polynomial

2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission

Chapter Unary Matrix Operations

ME 501A Seminar in Engineering Analysis Page 1

Chapter Linear Regression

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

Initial-Value Problems for ODEs. numerical errors (round-off and truncation errors) Consider a perturbed system: dz dt

Design of Bayesian MDS Sampling Plan Based on the Process Capability Index

International Journal of Scientific and Research Publications, Volume 3, Issue 5, May ISSN

Methods for solving the radiative transfer equation. Part 3: Discreteordinate. 1. Discrete-ordinate method for the case of isotropic scattering.

International Journal of Scientific & Engineering Research, Volume 6, Issue 12, December ISSN

On Stability of a Class of Fractional Differential Equations

Systems of second order ordinary differential equations

ADAPTIVE CLUSTER SAMPLING USING AUXILIARY VARIABLE

Introduction to Numerical Differentiation and Interpolation March 10, !=1 1!=1 2!=2 3!=6 4!=24 5!= 120

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

MATLAB PROGRAM FOR THE NUMERICAL SOLUTION OF DUHAMEL CONVOLUTION INTEGRAL

European Journal of Mathematics and Computer Science Vol. 3 No. 1, 2016 ISSN ISSN

Chapter 2 Intro to Math Techniques for Quantum Mechanics

4 Round-Off and Truncation Errors

Zjavka: Recognition of Generalized Patterns by a Differential Polynomial eural etwork

Soo King Lim Figure 1: Figure 2: Figure 3: Figure 4: Figure 5: Figure 6: Figure 7: Figure 8: Figure 9: Figure 10: Figure 11:

A METHOD FOR THE RAPID NUMERICAL CALCULATION OF PARTIAL SUMS OF GENERALIZED HARMONICAL SERIES WITH PRESCRIBED ACCURACY

Sequences and summations

ME/BioE C117 April 13, 2006 Professor Lisa Pruitt Exam #2. This is a closed notes/closed book exam. Show all work. Write your name on every page.

Numerical Methods for Eng [ENGR 391] [Lyes KADEM 2007] Direct Method; Newton s Divided Difference; Lagrangian Interpolation; Spline Interpolation.

Schrödinger Equation Via Laplace-Beltrami Operator

Neural Networks for Nonlinear Fractional Programming

PubH 7405: REGRESSION ANALYSIS REGRESSION IN MATRIX TERMS

Assignment 7/MATH 247/Winter, 2010 Due: Friday, March 19. Powers of a square matrix

7.0 Equality Contraints: Lagrange Multipliers

Numerical Analysis Formulae Booklet

Numerical Simulations of the Complex Modied Korteweg-de Vries Equation. Thiab R. Taha. The University of Georgia. Abstract

Transcription:

IOSR Jourl o Mtetcs (IOSR-JM) e-issn: - p-issn: 9-6X. Volue Issue Ver. I (Mr. - Apr. ) PP - www.osrourls.org Te Dervto o Iplct Secod Dervtve Metod or Solvg Secod-Order St Ordr Deretl Equtos Odes. Y. Skwe Deprtet o Mtetc Adw Stte Uverst Mub-Nger. Correspodg Autor: Y. Skwe Abstrct: A sgle-step brd block etod or tl vlue probles o geerl secod order Ordr Deretl Equtos s bee studed ts pper. I te dervto o te etod power seres s dopted s bss ucto to obt te cotuous scee troug collocto d terpoltos pproc. Tlor etod s lso used togeter wt ew etod to geerte te o-overlppg uercl results. Te ewl costructed etod s te ppled to solve te sste o secod-order st ordr deretl equtos d te ccurc s better we copred wt te estg etods ters o error. Kewords: Power Seres Collocto d Iterpolto Metod Hbrd Block Metod St ODEs Sste o Secod Order ODEs. ----------------------------------------------------------------------------------------------------------------------------- ---------- Dte o Subsso: -- Dte o cceptce: -- ----------------------------------------------------------------------------------------------------------------------------- ---------- I. Itroducto Most rel le probles tt rse vrous elds o stud be t egeerg or scece re odeled s tetcl odels beore te re solved. Tese odels ote led to deretl equtos. Nuerous probles suc s cecl ketcs orbtl dcs crcut d cotrol teor d Newto s secod lw pplctos volve secod-order ODEs [-]. Ordr deretl equtos (ODEs) re cool used or tetcl odelg dverse elds suc s egeerg operto reserc dustrl tetcs bevorl sceces rtcl tellgece geet d socolog. Ts tetcl odelg s te rt o trsltgproble ro pplcto re to trctble tetcl orultos wose teoretcl d uercl lss provde sgt swers d gudce useul or te orgtg pplcto []. Ts tpe o proble c be orulted eter ters o rst-order or ger order ODEs. I ts rtcle te sste o secod-order ODEs o te ollowg or s cosdered. " " " b b b Te etod o solvg ger-order ODEs b reducg te to sste o rst-order pproc volves ore uctos to evlute te d te leds to coputtol burde s etoed []-[]. Te ultstep etods or solvg ger-order ODEs drectl ve bee developed b scolrs suc s [6]- []. Te o ts pper s to develop ew uercl etod or solvg sstes o secod-order st ODEs. II. Dervto O Te Metod I ts secto oe-step brd block etod wt two o-step pots or solvg Equto () s derved. Let te power seres o te or d () v.... DOI:.99/- www.osrourls.org Pge ()

Te Dervto o Iplct Secod dervtve etod or solvg Secod-order St Ordr Deretl DOI:.99/- www.osrourls.org Pge be te pprote soluto to Equto () or were s N re te rel coecets to be detered v s te uber o collocto pots s te uber o terpolto pots d s costt step sze o te prtto o tervl b wc s gve b b N. Derettg Equto () oce d twce elds:. v (). " v () Iterpoltg Equto () t te selected tervls.e. d colloctg Equto () d () t ll pots te selected tervl.e. d gves te ollowg equtos wc c be wrtte tr or: 6 6 6 6 6 9 6 6 66 6 6 6 6 6 6 96 6 6 6 6 96 6. () Applg te Guss elto etod o Equto () gves te coecet or s. Tese vlues re te substtuted to Equto () to gve te plct cotuous brd etod o te or:. (6) Derettg Equto (6) oce elds:

Te Dervto o Iplct Secod dervtve etod or solvg Secod-order St Ordr Deretl DOI:.99/- www.osrourls.org Pge d d d d. () Were te cotuous scees re

I(z) Te Dervto o Iplct Secod dervtve etod or solvg Secod-order St Ordr Deretl III. Covergecelss. Order d error Costts o te Metods Accordg to [] te order o te ew etod Equto () s obted b usg te Tlor seres d t s oud tt te developed etod s uorl order Te wt error costts vector o: C..9. 9. T. Cosstec Deto.:Te brd block etod () s sd to be cosstet t s order ore t or equl to oe.e. P. Tereore te etod s cosstet.. Zero Stblt Deto.: Te brd block etod ()sd to be zero stble te rst crcterstc polol r vg roots suc tt r z d rz te te ultplct o rz ust ot greter t two. I order to d te zero-stblt o brd block etod () we ol cosder te rst crcterstc polol o te etod ccordg to deto (.) s ollows r r B r r z Wc ples r.hece te etod s zero-stble sce r z d rz.. Covergece Teore (.): Cosstec d zero stblt re sucet codto or ler ultstep etod to be coverget.sce te etod () s cosstet d zero stble t ples te etod s coverget or ll pot.. Regos o Absolute Stblt (RAS) Usg te MATLAB pckge we were ble to plot te stblt regos o te block etod (see g. below). Ts s doe b reorultg te block etod s geerl ler etod to obt te vlues o te trces ccordg to [] []. Te trces re substtuted to te stblt tr d usg MATLAB sotwre te bsolute stblt regos o te ew etods re plotted s sow g. below...6.. -. -. -.6 -. -.6 -. -. - -. -.6 -. -.. Re(z) Fgure: Rego o Absolute Stblt. DOI:.99/- www.osrourls.org Pge

Te Dervto o Iplct Secod dervtve etod or solvg Secod-order St Ordr Deretl.6 Nuercl Ipleetto To stud te ecec o te block brd etod or K we preset soe uercl eples b Skwe et-l []. I ts secto te perorce o te ew sgle-step brd block etod s eed usg te ollowg two sstes o secod-order tl vlue probles o ordr deretl equtos wt K d K. Tbles d sow te coprso o te uercl results o te ew etod wt te estg etods Skwe et-l []or solvg Eple d. Eple 99 999 99 Wt Ect Soluto e e 999 e e. (See Skwe et-l []) Tble : Coprso o bsolute errors or eple. Eple 9 9 99 99 Wt Ect Soluto e e (See Skwe et-l []). DOI:.99/- www.osrourls.org Pge

Te Dervto o Iplct Secod dervtve etod or solvg Secod-order St Ordr Deretl Tble : Coprso o bsolute errors or eple. IV. Coclusos I ts rtcle uorl order egt plct sgle-step block etod wt two o-step pots s derved v te terpolto d collocto pproc. Te bsolute errors rsg ro eples d usg te ew etod were copred wt te estg etod []. Skwe et-l solves eples d t s evdet ro te tbles preseted bove tt te ewl proposed etod perors better t Skwe et-l []. Te etod s lso desrble b vrtue o possessg g order o ccurc. Te developed etod s cosstet A stble coverget wt rego o bsolute stblt d s uorl order egt. Reereces []. Mud R.; Y.Y.A A st order plct brd bckwrd deretto orule (HBDF) or block soluto o ordr deretl equtos. Aer. J. Mt. Sttstcs. p. 9-9. []. M. Alksssbe; Zur O. Iplct oe-step block brd trd-dervtve etod or te drect soluto o tl vlue probles o secod order ordr deretl equtos. J. ppl.t. p. []. Or Z.; Sul M. Prllel r-pot plct block etod or solvg ger order ordr deretl equtos drectl. J. ICT 66. []. Y. Skwe; J. Sbo; P. Tub; T. Y. kg.order Te Iplct Oe-Step Hbrd Block Metod or Te Soluto o St Secodorder Ordr Deretl Equtos.IJEAS. 9-66. []. Jes A.; Ades A.; Josu S. Cotuous block etod or te soluto o secod order tl vlue probles o ordr deretl equto. It. J. Pure Appl. Mt. 6. [6]. Or Z.; Sule M.B. Prllel two-pot eplct block etod or solvg g-order ordr deretl equtos. It. J. Sul. Process Model.6.. J. Coput. Appl. []. Vgo-Agur J.; Ros H. Vrble stepsze pleetto o ultstep etods or Mt. 6 9. []. Ades A.O.; Ake T.A.; Udo O. Iproved cotuous etod or drect soluto o geerl secod order ordr deretl equtos. J. Nger. Assoc. Mt. Ps. 9 6. [9]. Awoe D.O. A P-stble ler ultstep etod or solvg geerl trd order o ordr deretl equtos. It. J. Coput. Mt. 9 99. []. Herc P. Soe Applctos o te Quotet-Derece Algort. Proc. Sp. Appl. Mt. 96 9. []. Lbert J.D. Coputtol Metods ODEs; Jo Wle d Sos: New York NY USA 9. []. Y. Skwe.; J. Sbo.; T. Y. Kg. Te costructo o plct oe-step block brd etods wt ultple o-grd pots or te soluto o st ODEs. JSRR -. []. Y. Skwe; G. M. Kuleg; I. A. Bkr.Secod dervtve Hbrd Block Bckwrd Deretto Forule or Nuercl Soluto o St Sstes. JAMCS p. (): -. Y. Skwe. " Te Dervto o Iplct Secod Dervtve Metod or Solvg Secod-Order St Ordr Deretl Equtos Odes." IOSR Jourl o Mtetcs (IOSR-JM). (): -. DOI:.99/- www.osrourls.org Pge