Development of a New Scheme for the Solution of Initial Value Problems in Ordinary Differential Equations

Similar documents
Chapter Direct Method of Interpolation

Inventory Management Models with Variable Holding Cost and Salvage value

Solutions to Problems from Chapter 2

ON THE STABILITY OF DELAY POPULATION DYNAMICS RELATED WITH ALLEE EFFECTS. O. A. Gumus and H. Kose

On the Pseudo-Spectral Method of Solving Linear Ordinary Differential Equations

5.1-The Initial-Value Problems For Ordinary Differential Equations

Contraction Mapping Principle Approach to Differential Equations

e t dt e t dt = lim e t dt T (1 e T ) = 1

September 20 Homework Solutions

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function

15/03/1439. Lecture 4: Linear Time Invariant (LTI) systems

Spectral Galerkin Method for Optimal Control Problems Governed by Integral and Integro- Differential Equations

CALDERON S REPRODUCING FORMULA FOR DUNKL CONVOLUTION

Solutions for Nonlinear Partial Differential Equations By Tan-Cot Method

0 for t < 0 1 for t > 0

Average & instantaneous velocity and acceleration Motion with constant acceleration

A new model for solving fuzzy linear fractional programming problem with ranking function

PART V. Wavelets & Multiresolution Analysis

Additional Exercises for Chapter What is the slope-intercept form of the equation of the line given by 3x + 5y + 2 = 0?

INTEGRALS. Exercise 1. Let f : [a, b] R be bounded, and let P and Q be partitions of [a, b]. Prove that if P Q then U(P ) U(Q) and L(P ) L(Q).

f t f a f x dx By Lin McMullin f x dx= f b f a. 2

4.8 Improper Integrals

CS537. Numerical Analysis

Convergence of Singular Integral Operators in Weighted Lebesgue Spaces

Topic 6b Finite Difference Approximations

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 9: The High Beta Tokamak

Distribution of Mass and Energy in Five General Cosmic Models

1.0 Electrical Systems

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

1. Introduction. 1 b b

IX.2 THE FOURIER TRANSFORM

Fuzzy Laplace Transforms for Derivatives of Higher Orders

A Simple Method to Solve Quartic Equations. Key words: Polynomials, Quartics, Equations of the Fourth Degree INTRODUCTION

Systems Variables and Structural Controllability: An Inverted Pendulum Case

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6.

Lecture #6: Continuous-Time Signals

Advances in Intelligent Systems Research, volume 136 4th International Conference on Sensors, Mechatronics and Automation (ICSMA 2016)

Procedia Computer Science

Abstract. W.W. Memudu 1 and O.A. Taiwo, 2

3. Renewal Limit Theorems

Asymptotic relationship between trajectories of nominal and uncertain nonlinear systems on time scales

RESPONSE UNDER A GENERAL PERIODIC FORCE. When the external force F(t) is periodic with periodτ = 2π

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE BOUNDARY-VALUE PROBLEM

GENERALIZATION OF SOME INEQUALITIES VIA RIEMANN-LIOUVILLE FRACTIONAL CALCULUS

MAT 266 Calculus for Engineers II Notes on Chapter 6 Professor: John Quigg Semester: spring 2017

Method For Solving Fuzzy Integro-Differential Equation By Using Fuzzy Laplace Transformation

CONSIDERATIONS REGARDING THE OPTIMUM DESIGN OF PRESTRESSED ELEMENTS

Mathematics Paper- II

Hermite-Hadamard-Fejér type inequalities for convex functions via fractional integrals

Location is relative. Coordinate Systems. Which of the following can be described with vectors??

An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples.

û s L u t 0 s a ; i.e., û s 0

TMA4329 Intro til vitensk. beregn. V2017

A Kalman filtering simulation

ON THE OSTROWSKI-GRÜSS TYPE INEQUALITY FOR TWICE DIFFERENTIABLE FUNCTIONS

Neural assembly binding in linguistic representation

SOME USEFUL MATHEMATICS

Math Week 5 concepts and homework, due Friday February 10

Chapter 2. Numerical Integration also called quadrature. 2.2 Trapezoidal Rule. 2.1 A basic principle Extending the Trapezoidal Rule DRAWINGS

MTH 146 Class 11 Notes

Temperature Rise of the Earth

Some Inequalities variations on a common theme Lecture I, UL 2007

FM Applications of Integration 1.Centroid of Area

Physics 2A HW #3 Solutions

Green s Functions and Comparison Theorems for Differential Equations on Measure Chains

2k 1. . And when n is odd number, ) The conclusion is when n is even number, an. ( 1) ( 2 1) ( k 0,1,2 L )

NEW FRACTIONAL DERIVATIVES WITH NON-LOCAL AND NON-SINGULAR KERNEL Theory and Application to Heat Transfer Model

Section P.1 Notes Page 1 Section P.1 Precalculus and Trigonometry Review

Numerical Analysis: Trapezoidal and Simpson s Rule

Approximation and numerical methods for Volterra and Fredholm integral equations for functions with values in L-spaces

Comparison between the Discrete and Continuous Time Models

Some algorthim for solving system of linear volterra integral equation of second kind by using MATLAB 7 ALAN JALAL ABD ALKADER

4. Runge-Kutta Formula For Differential Equations

12 Basic Integration in R

Minimum Squared Error

A Time Truncated Improved Group Sampling Plans for Rayleigh and Log - Logistic Distributions

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

Minimum Squared Error

How to Prove the Riemann Hypothesis Author: Fayez Fok Al Adeh.

A NUMERICAL SOLUTION OF THE URYSOHN-TYPE FREDHOLM INTEGRAL EQUATIONS

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

S Radio transmission and network access Exercise 1-2

Chapter 2. Motion along a straight line. 9/9/2015 Physics 218

FRACTIONAL-order differential equations (FDEs) are

Dipartimento di Elettronica Informazione e Bioingegneria Robotics

f(x) dx with An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples dx x x 2

Part I: Basic Concepts of Thermodynamics

How to prove the Riemann Hypothesis

ENGI 3424 Engineering Mathematics Five Tutorial Examples of Partial Fractions

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Application on Inner Product Space with. Fixed Point Theorem in Probabilistic

Review of Calculus, cont d

SOME PROPERTIES OF GENERALIZED STRONGLY HARMONIC CONVEX FUNCTIONS MUHAMMAD ASLAM NOOR, KHALIDA INAYAT NOOR, SABAH IFTIKHAR AND FARHAT SAFDAR

Scientific Research of the Institute of Mathematics and Computer Science DIFFERENT VARIANTS OF THE BOUNDARY ELEMENT METHOD FOR PARABOLIC EQUATIONS

Estimating the population parameter, r, q and K based on surplus production model. Wang, Chien-Hsiung

( ) ( ) ( ) ( ) ( ) ( y )

COSC 3361 Numerical Analysis I Numerical Integration and Differentiation (III) - Gauss Quadrature and Adaptive Quadrature

Tax Audit and Vertical Externalities

FRACTIONAL EULER-LAGRANGE EQUATION OF CALDIROLA-KANAI OSCILLATOR

Principles of Geometric Optics

Transcription:

IOS Journl o Memics IOSJM ISSN: 78-78 Volume Issue July-Aug PP -9 www.iosrjournls.org Developmen o New Sceme or e Soluion o Iniil Vlue Problems in Ordinry Dierenil Equions Ogunrinde. B. dugb S. E. Deprmen o Memicl Sciences Ekii Se Universiy Ado Ekii Nigeri. Absrc: Tis pper is on e developmen o new sceme or e soluion o iniil vlue problems in ordinry dierenil equions. We presen some bsic conceps nd undmenl eories wic re very vil o e developmen o e sceme. Te new sceme is o order seven nd is derivion is bsed on e represenion o e eoreicl soluion by perurbed polynomil uncion o degree our. Tis sceme is suible or problems ssocied wi e sysems o ordinry dierenil equions ving oscillory or eponenil soluion. Key words: Convergence Eponenil Soluion Ordinry Dierenil Equion I. Inroducion In e yers ps lrge number o meods suible or solving ordinry dierenil equions ve been proposed. A mjor impeus o developing numericl procedures ws e invenion o e clculus by Newon nd Leibniz s is led o ccure memicl models or pysicl reliy suc s Sciences Engineering Medicine nd Buess. Tese memicl models cnno be usully solved eplicily nd numericl meod o obin pproime soluions is needed. Anoer imporn spec o e developmen o numericl meods ws e creion o logrims by Npier nd oers giving muc simpler mnner o crrying ou e rimeic operions o muliplicion division nd eponeniion. Up o e le 8 s i ppers mos memicins were quie brod in eir ineres. Generlly e eiciency o ny o ese meods depends on e sbiliy nd ccurcy properies. Te ccurcy properies o dieren meods re usully compred by considering e order o convergence s well s e runcion error coeicien o e vrious meods. Te sources o moivion or is work re ose o wo proposed numericl inegrion sceme o order our nd considered o order ive suied or e soluion o iniil vlue problem in ordinry dierenil equion o e orm b rele ving eponenil soluion were nd b re rel undeermined coeiciens nd re comple prmeers nd inroduced new sceme o order si wic is n improvemen over. Tere re mny numericl inegring scemes o genere e numericl soluions developed by severl uors suc s nd 7 jus o menion ew. In is pper we sll consider e iniil vlue problem o e orm y y y nd develop n lgorim o order seven wic cn eecively solve iniil vlue problem in ordinry dierenil equion. II. Some Bsic Conceps We sll consider e ollowing bsic conceps wic re very vil o e developmen o e new sceme.. Sbiliy A numericl meod is sid o be sble i e dierence beween e numericl soluion nd e ec soluion cn be mde s smll s possible is i ere eiss wo posiive numbers e nd K suc e ollowing olds. y y K e n n. Consisency A numericl sceme wi n incremen uncion n y is sid o be consisen wi e iniil vlue problem under considerion i y y n. Convergence A numericl meod is sid o be convergen i or ll iniil vlue problem sisying e ypoesis o Lipsciz condiion given by www.iosrjournls.org Pge

Developmen o New Sceme or e Soluion o Iniil Vlue Problems in Ordinry Dierenil y y L y y Were e Lipsciz condiion is denoed by L.Te necessry nd suicien condiions or convergence re e sbiliy nd consisency.. ound o Error Tis cn be deined s e error due o compuing device. Tey rise becuse i is possible o represen ll rel numbers ecly on inie-se mcine. I cn be represened memiclly s y p n n n Were yn is e pproime soluion nd pn is e compuer oupu. Te mgniude depends in e sorge nd e rimeic operion doped. In suc cses double precision re employed o gurnee n deque pproimion. III. Developmen o New Sceme o Order Seven Tis secion presens e derivion o e new sceme o order seven or e soluion o iniil vlue problems in ordinry dierenil equions.. Te Bsic Inerpoln Le us ssume e eoreicl soluion y o e IVP problem y y y cn be represened loclly in e inervl by e polynomil inerpoling uncion. b rele Were I we pu nd b re rel undeermined coeiciens nd re comple prmeers. i i i And in we obin e ollowing inerpoling uncion be I we move urer o deine e uncions nd s ollows e 8 rom 7 nd 8 we obin 9 We sll ssume yis e numericl soluion o e eoreicl soluion y nd y. We deine mes poin s ollows:... We en move orwrd o impose e ollowing consrins on e inerpoling uncion 9.. T e inerpoling uncion mus coincide wi e eoreicl soluion nd. In order words we required y i.e. nd y. i.e. b. We lso require e irs second ird our nd i derivives wi respec o o e inerpoling uncion respecively coincide wi e dierenil equion s well s is irs second ird our nd i derivives wi respec o. In e oer words we require : 7 www.iosrjournls.org Pge

Developmen o New Sceme or e Soluion o Iniil Vlue Problems in Ordinry Dierenil www.iosrjournls.org Pge Tis implies 7 7 8 8 Noe were i denoes e i derivive o y wi respec o by seing nd. We obin direcly rom 8 nd e 9 rom 8 we ve b rom 7 en rom were

Developmen o New Sceme or e Soluion o Iniil Vlue Problems in Ordinry Dierenil www.iosrjournls.org 7 Pge rom Ten rom Tereore I we subrc rom nd pply 8 we obin e ollowing relions y y e Were = 7 8 To obin e new sceme puing e vlues o 7 nd 8 in we ve

Developmen o New Sceme or e Soluion o Iniil Vlue Problems in Ordinry Dierenil www.iosrjournls.org 8 Pge y y e 9 Equion 9 is e new sceme o order seven.

Developmen o New Sceme or e Soluion o Iniil Vlue Problems in Ordinry Dierenil IV. Conclusion Tis pper considers e developmen o new numericl sceme or e soluion o iniil vlue problem in ordinry dierenil equions.aer successul derivion o e lgorim i is obvious e sceme will be useul in solving iniil vlue problem mos especilly problems ving oscillory or eponenil soluions. Compring e sceme wi e eig sndrd scemes especilly o Ibijol 997 wic is o order ive Ogunrinde wic is o order Si in ours is o order seven wic sows i is o iger order.owing o ime nd compuer ime llocion we ve no considered e nlysis o sbiliy propery o is sceme in gre deils. eerences unl S. O. 98: Numericl Inegrors or Si nd Higly Oscillory Dierenil Equions Memics o Compuion 7-9. Gusci W. 9: Numericl Inegrion o Ordinry Dierenil Equions bsed on Trigonomeric Polynomils Numerisce Memik 8-97. Ibijol E. A. 997: New Scemes or Numericl Inegrion o Specil Iniil Vlue Problems in Ordinry Dierenil Equions P.D. Tesis Universiy o Benin Nigeri. Ibijol E. A nd Ogunrinde. B. : On New Numericl Sceme or e Soluion o Iniil Vlue Problems Acceped Ausrlin Journl o Bsic nd Applied Sciences. Josep-Loius Lgrnge. 8: On oo inding nd Polynomil Inerpolion in Ordinry Dierenil Equions. Ogunrinde. B. 9: On New Numericl Sceme or e Soluion o Iniil Vlue Problem in Ordinry Dierenil Equions. 7 Wllce C. S nd Gup G. K. 97: Generl Liner Mulisep Meods o Solve Ordinry Dierenil Equions. Te Ausrlin Compuer Journl Vol. -9. www.iosrjournls.org 9 Pge