Jornal of Kerbala University, Vol. 5 No.4 Scientific.Decembar 2007

Similar documents
ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

Comparison between Fourier and Corrected Fourier Series Methods

Available online at J. Math. Comput. Sci. 4 (2014), No. 4, ISSN:

Application of the Adomian Decomposition Method (ADM) and the SOME BLAISE ABBO (SBA) method to solving the diffusion-reaction equations

Research Article A Generalized Nonlinear Sum-Difference Inequality of Product Form

The Solution of the One Species Lotka-Volterra Equation Using Variational Iteration Method ABSTRACT INTRODUCTION

EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D. S. Palimkar

BIBECHANA A Multidisciplinary Journal of Science, Technology and Mathematics

1 Notes on Little s Law (l = λw)

Application of Fixed Point Theorem of Convex-power Operators to Nonlinear Volterra Type Integral Equations

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

Supplement for SADAGRAD: Strongly Adaptive Stochastic Gradient Methods"

Basic Results in Functional Analysis

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE

Application of Homotopy Perturbation Method to Biological Population Model

Two Implicit Runge-Kutta Methods for Stochastic Differential Equation

TAKA KUSANO. laculty of Science Hrosh tlnlersty 1982) (n-l) + + Pn(t)x 0, (n-l) + + Pn(t)Y f(t,y), XR R are continuous functions.

Homotopy Analysis Method for Solving Fractional Sturm-Liouville Problems

S n. = n. Sum of first n terms of an A. P is

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract

K3 p K2 p Kp 0 p 2 p 3 p

Existence Of Solutions For Nonlinear Fractional Differential Equation With Integral Boundary Conditions

Common Fixed Point Theorem in Intuitionistic Fuzzy Metric Space via Compatible Mappings of Type (K)

Approximating Solutions for Ginzburg Landau Equation by HPM and ADM

Review Exercises for Chapter 9

Extended Laguerre Polynomials

Department of Mathematical and Statistical Sciences University of Alberta

Prakash Chandra Rautaray 1, Ellipse 2

Numerical Method for Ordinary Differential Equation

On The Generalized Type and Generalized Lower Type of Entire Function in Several Complex Variables With Index Pair (p, q)

APPROXIMATE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS WITH UNCERTAINTY

Numerical Solution of Parabolic Volterra Integro-Differential Equations via Backward-Euler Scheme

INVESTMENT PROJECT EFFICIENCY EVALUATION

A Note on Random k-sat for Moderately Growing k

DIFFERENTIAL EQUATIONS

Extremal graph theory II: K t and K t,t

Some Properties of Semi-E-Convex Function and Semi-E-Convex Programming*

A Generalization of Hermite Polynomials

CLOSED FORM EVALUATION OF RESTRICTED SUMS CONTAINING SQUARES OF FIBONOMIAL COEFFICIENTS

A Study On (H, 1)(E, q) Product Summability Of Fourier Series And Its Conjugate Series

STK4080/9080 Survival and event history analysis

2 f(x) dx = 1, 0. 2f(x 1) dx d) 1 4t t6 t. t 2 dt i)

Fermat Numbers in Multinomial Coefficients

MODIFIED ADOMIAN DECOMPOSITION METHOD FOR SOLVING RICCATI DIFFERENTIAL EQUATIONS

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Boundary Value Problem for the Higher Order Equation with Fractional Derivative

Approximately Quasi Inner Generalized Dynamics on Modules. { } t t R

SOLVING OF THE FRACTIONAL NON-LINEAR AND LINEAR SCHRÖDINGER EQUATIONS BY HOMOTOPY PERTURBATION METHOD

Some inequalities for q-polygamma function and ζ q -Riemann zeta functions

The Moment Approximation of the First Passage Time for the Birth Death Diffusion Process with Immigraton to a Moving Linear Barrier

On Stability of Quintic Functional Equations in Random Normed Spaces

A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY

ABSOLUTE INDEXED SUMMABILITY FACTOR OF AN INFINITE SERIES USING QUASI-F-POWER INCREASING SEQUENCES

Section 8 Convolution and Deconvolution

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

Mean Square Convergent Finite Difference Scheme for Stochastic Parabolic PDEs

On stability of first order linear impulsive differential equations

A Complex Neural Network Algorithm for Computing the Largest Real Part Eigenvalue and the corresponding Eigenvector of a Real Matrix

Online Supplement to Reactive Tabu Search in a Team-Learning Problem

The Connection between the Basel Problem and a Special Integral

Averaging of Fuzzy Integral Equations

Lecture 15 First Properties of the Brownian Motion

METHOD OF THE EQUIVALENT BOUNDARY CONDITIONS IN THE UNSTEADY PROBLEM FOR ELASTIC DIFFUSION LAYER

On the Existence and Uniqueness of Solutions for Nonlinear System Modeling Three-Dimensional Viscous Stratified Flows

F D D D D F. smoothed value of the data including Y t the most recent data.

On Another Type of Transform Called Rangaig Transform

Boundary-to-Displacement Asymptotic Gains for Wave Systems With Kelvin-Voigt Damping

Generalized Jacobi spectral-galerkin method for nonlinear Volterra integral equations with weakly singular

VIM for Determining Unknown Source Parameter in Parabolic Equations

The analysis of the method on the one variable function s limit Ke Wu

Lecture 15: Three-tank Mixing and Lead Poisoning

Notes 03 largely plagiarized by %khc

Inverse Heat Conduction Problem in a Semi-Infinite Circular Plate and its Thermal Deflection by Quasi-Static Approach

Fuzzy Dynamic Equations on Time Scales under Generalized Delta Derivative via Contractive-like Mapping Principles

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3

Parametric Iteration Method for Solving Linear Optimal Control Problems

REDUCED DIFFERENTIAL TRANSFORM METHOD FOR GENERALIZED KDV EQUATIONS. Yıldıray Keskin and Galip Oturanç

Paper 3A3 The Equations of Fluid Flow and Their Numerical Solution Handout 1

MATH 507a ASSIGNMENT 4 SOLUTIONS FALL 2018 Prof. Alexander. g (x) dx = g(b) g(0) = g(b),

N! AND THE GAMMA FUNCTION

Sampling Example. ( ) δ ( f 1) (1/2)cos(12πt), T 0 = 1

Calculus BC 2015 Scoring Guidelines

Bernstein Direct Method for Solving. Variational Problems

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are:

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi

Research Article Generalized Equilibrium Problem with Mixed Relaxed Monotonicity

Inference of the Second Order Autoregressive. Model with Unit Roots

The Central Limit Theorem

FRACTIONAL VARIATIONAL ITERATION METHOD FOR TIME-FRACTIONAL NON-LINEAR FUNCTIONAL PARTIAL DIFFERENTIAL EQUATION HAVING PROPORTIONAL DELAYS

ECE 570 Session 7 IC 752-E Computer Aided Engineering for Integrated Circuits. Transient analysis. Discuss time marching methods used in SPICE

Fractional Lévy Cox-Ingersoll-Ross and Jacobi processes

Fresnel Dragging Explained

LINEAR APPROXIMATION OF THE BASELINE RBC MODEL JANUARY 29, 2013

Economics 8723 Macroeconomic Theory Problem Set 2 Professor Sanjay Chugh Spring 2017

NEWTON METHOD FOR DETERMINING THE OPTIMAL REPLENISHMENT POLICY FOR EPQ MODEL WITH PRESENT VALUE

SUMMATION OF INFINITE SERIES REVISITED

Institute of Actuaries of India

Actuarial Society of India

BE.430 Tutorial: Linear Operator Theory and Eigenfunction Expansion

Problems and Solutions for Section 3.2 (3.15 through 3.25)

Transcription:

Joral of Kerbala Uiversiy, Vol. No. Scieific.Decembar 7 Soluio of Delay Fracioal Differeial Equaios by Usig Liear Mulise Mehod حل الوعادالث التفاضل ت الكسز ت التباطؤ ت باستخذام طز قت هتعذد الخطىاث الخط ت Basim K. AL-Salai Dearme of Comuer, College of Sciece, Kerbala Uiversiy, Iraq (basim78@yahoo.com Absrac: The objecive of his aer o fid he soluio of delay fracioal differeial equaios as well as, heir umerical soluio. The modificaio was mode by iroducig he ew mehod which is liear mulise mehods, which is solvig umerically ordiary differeial equaios, delay differeial equaios ad delay fracioal differeial equaios. Furhermore his aer reses he saeme ad roof of he fudameal heorem of coverge ad sabiliy of delay fracioal differeial equaios, his had bee doe wih some modificaio. الخالصة الهذف هي هذا البحث هى إل جاد ح ل الوعادالث التفاضل ت الكسىر ت التباطؤ ت باإلضافت إلى ح ل هن العذدي. التعذ ل ك اى وطا بت قذ ن الطز قت الجذ ذة الت ه طزق هتعذدة الخطىاث الخط ت الت ت ح ل هعادالث تفاضل ت عاد ت بشكل عذدي الوعادالث التفاضل ت التباطؤ ت والوعادالث التفاضل ت الكسز ت التباطؤ ت. اضافت اى هذا البحث ب ي و بزهي ال ظز ت األساس ت الوتقارب ه وإستقزارالوعادالث التفاضل ت الكسز ت التباطؤ ت هذا هع ع و ل ب ع ض التعذ ل عل ها..Iroducio The delay fracioal differeial equaios lay a imora role i he heory of fucioal differeial equaios.i rece year, he heory of his class of equaios had become a ideede red ad lieraure o his subjec comrises over iles. delay fracioal differeial equaios have bee sudied durig he las wo ceuries. The sigificace of hese equaios lies i heir abiliy o describe rocesses wih rearded ime. The imorace of hese equaios i various braches of echology, ecoomics, biology, medical sciece, hysics ad social sciece has bee recogized recely ad has caused mahemaics o sudy hem wih icreasig ieres.oe of he mos imora braches of fracioal differeial equaios is he so called delay fracioal differeial equaios, which is a fracioal differeial equaio cosis of oe or some is argumes evaluaed a a ime which is differ by ay fixed umber or variable which is called he ime or lag. Also, heory of delay fracioal differeial equaio ca be cosidered as a geeralizaio o he heory of fracioal differeial equaios ad usig Liear Mulise Mehods (LMM's o solve delay fracioal differeial equaios of he rearded ad derivig some of he well ow umerical mehods o solve delay fracioal differeial equaios..mehod of soluio delay fracioal differeial equaios Because of he iiial codiio which a fucio give for ime se ierval wih legh equals, mus fid his soluio for divided i o ses wih legh also. Hece, we mus fid he soluio for each ime se. The bes well ow mehod for solvig delay fracioal differeial equaios is he mehod of ses (soluio of fracioal differeial equaios, [Driver, 977],[Halaay, 966], [El'sgol'c ad Nori, 96] [Basim ], which is for simliciy wih ou los of geeraliy ad illusraio urose of solvig delay fracioal differeial equaios of he form: x q ( = f(, x(, x(-, x ' (-...( for each equaio he soluio is cosruced se by se as follows: Give ha a fucio ( coiues o [ -, ], herefore oe ca form he soluio i he ex ime se ierval [, +], by solvig he followig equaio: x q ( = f(, x(, (-, ' (-, for +, q< 7

Joral of Kerbala Uiversiy, Vol. No. Scieific.Decembar 7 wih iiial codiio x( =(. If we suose ha ( is he exisig by virue of coiuiy hyoheses. Similarly if ( is defied o he whole segme [, +], hece by formig he ew equaio: x q ( = f(, x(, (-, ' (-, for + +, q< oe ca fid he soluio of his equaio ( o he ime se [ +, +], wih he iiial codiio x( += ( +. I geeral, by assumig ha - ( is defied o he ierval [ +(-, +(-], =,,,, he followig equaio is cosruced: x q ( = f(, x(, - (-, ' - (-, for +(- +, q< ad cosider a soluio of his equaio wih he iiial codiio: x( +(- = - ( +(-, which is deoed by ( for all =,. Now, cosider he examle o exlai he mehod of ses. Examle (.: Cosider he rearded delay fracioal differeial equaios: x ( x(...( wih iiial codiio x(= (=, for -, ad o fid he soluio i he firs se ierval [, ], we have o solve he followig equaio: x he, we have: x ( (, for...( (, for...( by usig he soluio of fracioal differeial equaios, [Basim, ], alyig d o boh sides of d equaio (, we ge: x( c, ( 7 ( i order o fid he soluio i he secod se ierval, suose ha: ( x( c, ( 7 ( is he iiial codiio. Sice ( is defied o he whole segme [, ], hece by formig a ew equaio: x wih iiial codiio, ( (, for...( ( c, for ( 7 ( oe ca fid he soluio i he ex se ierval [, ], ad solvig equaio (, for x( we have: x ( ( ( ( c(,...(6 ( 7 ( by usig fracioal differeial equaio alyig d d o boh sides of equaio (6, we ge: 8

Joral of Kerbala Uiversiy, Vol. No. Scieific.Decembar 7 x( ( ( ( ( ( 7 ( 9 c ( c ( ( c...(7 for all,. herefore x( is he desired secod se soluio. similarly, we roceed o he ex iervals..the Exisece ad Uiqueess Theorem As i ordiary ad fracioal differeial equaios he exisece ad uiqueess delay fracioal differeial equaios could also be esablished for his urose, cosider he delay fracioal differeial equaios: y q (=f(, y(, y(-w( wih iiial codiio, y(=(, -w( ad suose ha R deoed he recagular regio defied by: R={(, y: -a-w +a, y-y b}.. Exisece heorem cosider he iiial value roblem of delay fracioal differeial equaios: y q (=f(, y(, y(-w(...(8 wih iiial codiio, y(=(, -w( f f ad suose ha f, ad are coiuous o a regio R, ad f saisfies Lischiz y( y( w( cosa codiio. The equaio (8 has a soluio o R. oe hyohesis of he exisece i he recagle R. I is follows ha f mus be boded i R. le M> be umber such ha f(, y(, y(-w( M, for every oi i R. if we ow ae = mi{a, M b }, ad defie he recagle R o R ={(, y:-, y-y b} clearly R is a subse of R. If we ow cosider he sequece of fucios: ( f (s, (s, hus, we have he sequece of soluios will coverge. (s w(sds.. The uiqueess heorem The soluio of he delay fracioal differeial equaio: y q (=f(, y(, y(-w( wih iiial codiio y(=(, -w(, uique. Proof: Suose ha ( ad ( are ay wo soluio, he: [f (s, (s, (s w(s f (s, (s, (s w(s]ds ad suose ha K is he Lischiz cosa codiio of he fucio f i a comac subse R coaiig (, (, (-w( ad (, (, (-w(, for each R. he: [f (s, (s, (s w(s f (s, 8 (s, (s w(s]ds

Joral of Kerbala Uiversiy, Vol. No. Scieific.Decembar 7 K [f (s, (s, (s (s w(s f (s, (s (s, (s w(s (s (s w(s]ds w(s ds K Su (s (s (s w(s (s w(s asb aig < K <, he we have: [ (s (s (s w(s ad hece: (- ( = herefore ( = (, foe each R. (s w(s ].Numerical Soluio of Delay Fracioal Differeial Equaios Usig Liear Mulise Mehods (LMM's. The rese secio cosiss of a geeralizaio ad modificaio of he geeral LMM o solve DFDE. For his urose, well cosider, for simliciy ad wihou lose of geeraliy he rearded DFDE., which has he form : y (q ( = f(, y(, y(-,...(9 wih iiial codiio: y( = (, -, q <, where is a fixed umber ad y( =y =(... Geeral LMM's i DFDE's: Cosider IVP give be equaio (9, where y( = (, -, ad y( =y =(. Hece he geeral form of LMM for solvig equaio (9 is give by: j y j h jf y y j...( j j where x +j =x + jh, y +j = y +jh ad y +j- =y +jh- ad j, j are cosa. Hece: j y j h jf y j j j exaded his formulae, gives: y + y + + + y + = h[ f, y, - + f +, y +, +- + + f( x +, y +, +- ] which is geeral Delay Fracioal Liear Mulise Mehod (DFLMM. As a classificaio o he DFLMM, we say ha he DFLMM is of exlici ye if = ad i is imlici ye if. Also i some cases, imlici mehod may reduce o exlici mehods deedig o he fucio f wheher i is liear or o-liear. The local rucaio error of DFLMM is defied o be, T jy jh h jf jh, y jh, y jh... ( j j or 88

Joral of Kerbala Uiversiy, Vol. No. Scieific.Decembar 7 T (q jy jh h jy jh j j afer exadig he fucio y +jh ad i's by usig soluio fracioal differeial equaio alyig d d ad q q y y (q o boh sides, ad collocaio erms i equaio (. y y jh y jh x!! y y y jh (jh...!! jh y (q y (q! jh y (q ad so o. by subsiuig i he local rucaio error, we ge: T jy j h jf y y j j j j j y j h jf y j j! (jh jh x... ax q... = y + y + + + y + - h[ f, y, - + f +, y +, +- + + f( x +, y +, +- ] (q (... y h( (...... y (q h ( (...... y (q h (! (!... 8!!!... (q q! y... ax herefore T =C y (q- + C hy (q +C h y (q+ + +C h y (q+ + +ax q-... ( where C, C,,C,a are cosas, ad a simle calculaio yields he followig formula for he cosas C i 's, i=,,,, i erms of he coefficies j ad j : C... C (...... ( C...... C! (! C! ( (... 8!...! 88!... ( (! hese formula ca be used o derived a cerai DFLMM of give srucure.!...

Joral of Kerbala Uiversiy, Vol. No. Scieific.Decembar 7. Illusraive examle Examle (.: Cosider he Delay fracioal differeial equaio: x ( x(, wih iiial codiio: x(=(=, for -, ad suose ha we wa o fid he umerical soluio i firs ime se, wih se legh h=. give by: y + =y +hf, y, y -, by usig he comuer rogram o fid he soluio o he firs ime se, ge he resuls reseed i able (. ad i's comariso wih he exac soluio. Table (.. x i Numerical soluio Exac soluio Absolue Error. -... -.7 -.66.766. -.7 -.7.. -.6 -.7.. -. -.66.66.6 -. -.6..7 -.69 -.6.9.8.8..666.9.8.8.99.6.6666.66..66.768.6 -..86.89.897 -..68.68.89 -..69.7976.988..886.97 6.776 -.6.9.78 7.8 -.7.68.6 9.679 -.8.667.77.69.9.98887.7.9.77.8.78 Refereces:. Al-Dayee, K. K., O Soluio of Delay Boudary Value Problems, M.sc. heses, dearme of mahemaics, Al_ahrai uiversiy, Baghdad, Iraq,.. Basim K., Soluio of Fracioal Differeial Equaios Usig Variaioal Aroach, Alahrai uiversiy, Baghdad, Iraq,.. Driver R. D. E., Ordiary ad Delay Differeial Equaios, Surigier-Verlag Ic., New Yor,977.. El'sgol'z, L.E. Qualiaive Mehod i Mahemaical Aalysis, Tras. Mah. Moo., Vol., Amer. Mah Soc.,96.. Halaay A., Differeial Equaios; Sabiliy, Oscillaios ad Time Lages, Academic ress, New Yor ad Lodo Ic., 966. 888