The Differential Transform Method for Solving Volterra s Population Model

Similar documents
Application of Homotopy Analysis Method for Solving various types of Problems of Ordinary Differential Equations

Solution of Differential Equation from the Transform Technique

Numerical Solution of Nonlinear Integro-Differential Equations with Initial Conditions by Bernstein Operational Matrix of Derivative

Modified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations

Lebesgue Constant Minimizing Bivariate Barycentric Rational Interpolation

Analysis of Analytical and Numerical Methods of Epidemic Models

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

wavelet collocation method for solving integro-differential equation.

The Hypergeometric Coupon Collection Problem and its Dual

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

Binomial transform of products

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

FUZZY RELIABILITY ANALYSIS OF COMPOUND SYSTEM BASED ON WEIBULL DISTRIBUTION

Research Article Extension of Wolfe Method for Solving Quadratic Programming with Interval Coefficients

DISTANCE BETWEEN UNCERTAIN RANDOM VARIABLES

International Journal of Mathematical Archive-4(9), 2013, 1-5 Available online through ISSN

Solving Nonlinear Stochastic Diffusion Models with Nonlinear Losses Using the Homotopy Analysis Method

A New Type of q-szász-mirakjan Operators

Acoustic Field inside a Rigid Cylinder with a Point Source

Exact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method

Chapter 4 : Laplace Transform

Lecture 11. Solution of Nonlinear Equations - III

PAijpam.eu ON TENSOR PRODUCT DECOMPOSITION

μ are complex parameters. Other

A New Family of Multipoint Iterative Methods for Finding Multiple Roots of Nonlinear Equations

Research Article Two Expanding Integrable Models of the Geng-Cao Hierarchy

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

On Bivariate Haar Functions and Interpolation Polynomial

RAYLEIGH'S METHOD Revision D

Lecture 10: Bounded Linear Operators and Orthogonality in Hilbert Spaces

AN EFFICIENT ESTIMATION METHOD FOR THE PARETO DISTRIBUTION

Stability Analysis of the Euler Discretization for SIR Epidemic Model

The Sumudu transform and its application to fractional differential equations

x+ 2 + c p () x c p () x is an arbitrary function. ( sin x ) dx p f() d d f() dx = x dx p cosx = cos x+ 2 d p () x + x-a r (1.

Exact Solutions of the Generalized Benjamin Equation and (3 + 1)- Dimensional Gkp Equation by the Extended Tanh Method

A comparative study of a system of Lotka-Voltera type of PDEs through perturbation methods

Crank-Nicolson Implicit Method For The Nonlinear Schrodinger Equation With Variable Coefficient

The Adomian Polynomials and the New Modified Decomposition Method for BVPs of nonlinear ODEs

A new sequence convergent to Euler Mascheroni constant

Numerical Solutions of Second Order Boundary Value Problems by Galerkin Residual Method on Using Legendre Polynomials

Orthogonal Function Solution of Differential Equations

HOMOTOPY PERTURBATION METHOD FOR VISCOUS HEATING IN PLANE COUETTE FLOW

Reliability Equivalence Analysis of a Parallel-Series System Subject to Degradation Facility

Evaluation of Bessel Functions Using a Computer Program

Application of Homotopy Perturbation Method for the Large Angle period of Nonlinear Oscillator

A Taylor Series Based Method for Solving a Two-dimensional Second-order Equation

NTMSCI 5, No. 1, (2017) 26

Automated Proofs for Some Stirling Number Identities

A collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation

AVERAGE MARKS SCALING

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

FREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES

Wavelet Transform Theory. Prof. Mark Fowler Department of Electrical Engineering State University of New York at Binghamton

Energy Identities of ADI-FDTD Method with Periodic Structure

Solving a Nonlinear Equation Using a New Two-Step Derivative Free Iterative Methods

A Generalization of Ince s Equation

Generating Functions for Laguerre Type Polynomials. Group Theoretic method

Observations on Derived K-Fibonacci and Derived K- Lucas Sequences

Int. Journal of Math. Analysis, Vol. 6, 2012, no. 31, S. Panayappan

APPLICATION OF UNCERTAIN NONLINEAR SYSTEMS PARTIAL STATE VARIABLES CONTROL TO A CLASS OF PENDULUM SYSTEMS

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

The Phi Power Series

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan

Stopping oscillations of a simple harmonic oscillator using an impulse force

Direct Solution of Initial Value Problems of Fourth Order Ordinary Differential Equations Using Modified Implicit Hybrid Block Method

Stability of fractional positive nonlinear systems

Estimation of Backward Perturbation Bounds For Linear Least Squares Problem

Adv. Studies Theor. Phys., Vol. 5, 2011, no. 4, A Numerical Computation of Non-Dimensional. Form of a Mathematical Model of Soil Salinity

SphericalHarmonicY. Notations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation

Uniform Strict Practical Stability Criteria for Impulsive Functional Differential Equations

GENERALIZED KERNEL AND MIXED INTEGRAL EQUATION OF FREDHOLM - VOLTERRA TYPE R. T. Matoog

Name Period ALGEBRA II Chapter 1B and 2A Notes Solving Inequalities and Absolute Value / Numbers and Functions

FUZZY TRANSPORTATION PROBLEM WITH ADDITIONAL RESTRICTIONS

Application 10.5D Spherical Harmonic Waves

Average Number of Real Zeros of Random Fractional Polynomial-II

Bernoulli Polynomials Talks given at LSBU, October and November 2015 Tony Forbes

Orthogonal Functions

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS

A Tabu Search Method for Finding Minimal Multi-Homogeneous Bézout Number

ASYMPTOTIC STABILITY OF POSITIVE FRACTIONAL 2D LINEAR SYSTEMS WITH DELAYS

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

International Journal of Multidisciplinary Research and Modern Education (IJMRME) ISSN (Online): (

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Fourth Order Positively Smoothed Padé Schemes. for Parabolic Partial Differential Equations with. Nonlocal Boundary Conditions

You may work in pairs or purely individually for this assignment.

Some results on the Apostol-Bernoulli and Apostol-Euler polynomials

Castiel, Supernatural, Season 6, Episode 18

A Steady State Heat Conduction Problem in. a Thick Annular Disc Due to Arbitrary. Axisymmetric Heat Flux

MULTIPLE TIME SCALES SOLUTION OF AN EQUATION WITH QUADRATIC AND CUBIC NONLINEARITIES HAVING FRAC- TIONAL-ORDER DERIVATIVE

Research Article Numerical Solution of Fractional Integro-Differential Equations by Least Squares Method and Shifted Chebyshev Polynomial

New Aspects for Non-linear Semigroups in L1 Applied to Heat Equation s Analysis

The Binomial Multi-Section Transformer

Multinomial. Notations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

Higher-order iterative methods by using Householder's method for solving certain nonlinear equations

1 6 = 1 6 = + Factorials and Euler s Gamma function

MATHEMATICAL MODELS - Vol. I - Controllability, Observability, and Stability of Mathematical Models - Abderrahman Iggidr

Some New Iterative Methods for Solving Nonlinear Equations

BINOMIAL COEFFICIENT HARMONIC SUM IDENTITIES ASSOCIATED TO SUPERCONGRUENCES

FIR Filter Design: Part II

Transcription:

AASCIT Couicatios Volue, Issue 6 Septeber, 15 olie ISSN: 375-383 The Differetial Trasfor Method for Solvig Volterra s Populatio Model Khatereh Tabatabaei Departet of Matheatics, Faculty of Sciece, Kafas Uiversity, Kars, Turey Era Guerha Departet of Coputer, Faculty of Egieerig, Kafas Uiversity, Kars, Turey I Keywords Volterra s Populatio Model, Itegro-differetial equatio, Differetial Trasfor Method this article, Differetial trasfor ethod is preseted for solvig Volterra s populatio odel for populatio growth of a species i a closed syste. This odel is a oliear itegro-differetial where the itegral ter represets the effect of toi. This powerful ethod catches the eact solutio. First Volterra s populatio odel has bee coverted to power series by oe-diesioal differetial trasforatio. Thus we obtaied uerical solutio Volterra s populatio odel. Itroductio This ote deals with a atheatical odel of the acculated effect of tois o a populatio livig i a closed syste. Scudo idicates i his review that volterra proposed this odel for a populatio q( ) of idetical idividuals wich ehibits crowdig ad sesitivity to total etabolis : dq t = aq bq cq q( ) d, dt q() = q, (1.1) Where a > is the birth rate coefficiet, b > is the crowdig coefficiet ad c > is the toicity coefficiet. The coefficiet c idicates the essetial behavior of the populatio evolutio before its level falls to zero i the log ter. q is the iitial populatio ad q = q() t deotes the populatio at tie t. This odel is a first-order itegro-ordiary t differetial equatio where the ter cq q( ) d, represets the effect of toi accuulatio o the species. We apply scale tie ad populatio by itroducig the odiesioal variables tc qb λ =, u =, b a to obtai the odiesioal proble: du λ u u u u( ) d dλ = u() = u, (1.) Where u( λ ) is the scaled populatio of idetical idividuals at tie λ ad is a prescribed o- diesioal paraeter. The

63 15; (6): 6-67 odel is characterized by the oliear Volterra itegro-differetial equatio du λ u u u u( ) d dλ = u() =.1, (1.3) c ad is a prescribed paraeter. The o diesioal paraeter is =. ba I recet years, uerous wors have bee focusig o the developet of ore advaced ad efficiet ethods for Volterra's populatio odel such as sigular perturbatio ethod [5]. For eaple hootopy perturbatio ethod [6], quasilierizatio approach ethod [4] ad Adoia decopositio ethod ad Sic-Galeri ethod copared for the solutio of soe atheatical populatio growth odels [3]. I [1], the series solutio ethod ad The decopositio ethod for Volterra's populatio odel is cosidered. I 1986, Zhou [] first itroduced the differetial trasfor ethod (DTM) i solvig liear ad oliear iitial value probles i the electrical circuit aalysis. The differetial trasfor ethod obtais a aalytical solutio i the for of a polyoial. It is differet fro the traditioal high order Taylor s series ethod, which requires sybolic copetitio of the ecessary derivatives of the data fuctios. the results of applyig differetial trasforatio ethod to the Volterra s populatio odel will be preseted. Oe-Diesioal Differetial Trasfor Differetial trasfor of fuctio y( ) is defied as follows: 1 d y( ) Y( ) =,! d (.1) I equatio (.1), y( ) is the origial fuctio ad Y( ) is the trasfored fuctio, which is called the T-fuctio. Differetial iverse trasfor of Y( ) is defied as Fro equatio (.1) ad (.), we obtai y( ) = Y( ), (.) = d y( ) y( ) =,! = d (.3) Equatio (.3) iplies that the cocept of differetial trasfor is derived fro Taylor series epasio, but the ethod does ot evaluate the derivatives sybolically. However, relative derivatives are calculated by a iterative way which are described by the trasfored equatios of the origial fuctios. I this study we use the lower case letter to represet the origial fuctio ad upper case letter represet the trasfored fuctio. Fro the defiitios of equatios (1.) ad (.), it is easily prove that the trasfored fuctios coply with the basic atheatics operatios show i Table 1. I actual applicatios, the fuctio y( ) is epressed by a fiite series ad equatio (.) ca be writte as: = + y( ) = Y( ), (.4) = Equatio (.3) iplies that 1 Y( ) is egligibly sall. I fact, is decided by the covergece of atural frequecy i this study.

ISSN: 375-383 64 Table 1. The fudaetal operatios of oe-diesioal differetial trasfor ethod. Origial fuctio Trasfored fuctio y( ) = u( ) ± v( ) Y( ) = U( ) ± V( ) y( ) = cw( ) Y( ) = cw( ) y( ) = dw/ d Y( ) = ( + 1) W( + 1) j j y( ) = d w/ d Y( ) = ( + 1)( + ) ( + j) W( + j) y( ) u( ) v( ) = r= Theore 1. if y ( ) = ep ( ) the W ( ) 1 =.! Y( ) = U( r) V( r) Proof: By usig equatio (.1), we have ( ) W ( e ) 1 1 1 = = =. e!!! t= Theore. if u( ) ( + )! y( ) = the Y( ) = U( + ).! Proof: By usig equatio (.1), we have + ( ) ( ) ( ) 1 u 1 u +! Y ( ) = ( ) = = U ( + ).!!! + t= t= Theore 3. if y( ) = the Y( ) = δ( ) = 1 = Proof: By usig equatio (.1), we have 1 ( ) W( ) = =! t= 1 ( )! = = 1 =!! 1 ( ) =! Theore 4. if w( ) = si ( w + α ) the W ( ) Proof: By usig equatio (.1), we have ( ) W( ) = δ = w π. = si( + α )! 1 =. 1 si( w + α) 1 = 1 W 1 = = wcos( w + ) 1! 1! ( ) α 1 π 1 π = wsi(( + α) + w) = wsi( + α), 1! 1!

65 15; (6): 6-67 1 si( w + α) 1 π = W ( ) = = w cos( + + w)! t! I the geeral for we have α 1 π π 1 π 1 π = w si( + α + + w) = w si( + α + w) = w si( + α),!!! 1 si( w + α) w π = W ( ) = = si( + α + w)!! w π = si( + α).! w π! Theore 5. if w( ) = cos( w + α ) the W ( ) = cos( + α ). Proof: By usig equatio (.1), we have I the geeral for we have 1 cos( w + α) 1 = 1 W 1 = = wsi( w + ) 1! t 1! ( ) α 1 π 1 π = wcos(( + α) + w) = wcos( + α), 1! 1! 1 cos( w + α) 1 π = W ( ) = = w si( + α + w)!! t 1 π π 1 π 1 π = w cos + α + + w == w cos( + α + w) = w cos( + α),!!! 1 si( w + α) w π = W ( ) = = cos( + α + w)!! w π w π = cos( + α) = cos( + α).!! Nuerical Eaple I this sectio, we use oe diesioal differetial trasfor ethod for solvig the populatio growth odel characterized by the oliear Volterra itegro-differetial equatio With the iitial coditio du 1 u() t u () t 1 u() t t u( ) d, dt = (3.1) Whe taig the oe diesioal differetial trasfor of (3.1), we ca obtia: u () =.1, (3.) U 1 1 ( h h h 1) 1 ( ) 1 ( ) ( ) 1 ( ) ( 1), ( h 1) U h U s U h s U h s U s + = + s= s= 1 s (3.3)

ISSN: 375-383 66 Fro the iitial coditio (3.), we ca write U () =.1, (3.4) For each h substitutig ito Eq. (3.3) ad by recursive ethod, the values U( h ) ca be evaluated as follows: U() =.1, U(1) =.9, U() = 3.55, U(3) = 6.316666667, U(4) = 5.5375499998, U(5) = 63.74166667, U(6) = 156.1375,..., (3.5) by usig the iverse trasforatio rule for oe diesioal i Eq. (.1), the followig solutio ca be obtaied: h= h 1 3 4 5 6 u() t = U( h) = U() + U(1) t + U() t + U(3) t + U(4) t + U(5) t + U(6) t +..., 3 4 5 6 7 =.1+.9t + 3.55t + 6.316666667t 5.5375499998 t 63.74166667t 156.1375 t + O( t ). (3.6) I a coplete agreeet with the results previously obtaied i the previous sectios. Coclusio I this study, the differetial trasfor ethod for the solutio of the Volterra s populatio odel is successfully epaded. It is observed that the ethod is robust. The ethod gives rapidly covergig series solutios.the uerical results show that differetial trasfor ethod is a accurate ad reliable uerical techique for the solutio of the Volterra s populatio odel. Khatereh Tabatabaei Kh. Tabatabaei received the PhD degree i Matheatics of Atatur Uiversity. Her research iterests are i the areas of applied atheatics icludig the uerical ethods for special odels of differetial equatios ad differetial-algebraic equatios ad itegral equatios. Eail address: Kh_tabatabaey@yahoo.co Era Güerha Lecturer, Kafas Uiversity Faculty of Egieerig, Departet of Coputer, Kars- Turey His research iterests are i the areas of atheatics Egieerig ad coputer sciece ad software prograig ad web desiger. Eail address: Era.guerha@gail.co Refereces [1] A. M. Wazwaz, Aalytical approiatios ad Padé approiats for Volterra s populatio odel, Applied Matheatics ad Coputatio, 1 (1999) 13-5. [] J. K. Zhou, Huazhog Uiversity Press, Wuha, Chia, 1986. [3] K. Al-Khaled, Nuerical Approiatios for Populatio Growth Models. Appl. Math. Coput, 16 (5) 865-873. [4] K. Parad, M. Ghasei ASEMI, S. Rezazadeh, A. Peiravi, A. Ghorbapour, A. Tavaoli golpaygai, quasilierizatio approach for solvig Volterra s populatio odel, Appl. Coput. Math., (1)95-13. [5] R. D. Sall, Populatio growth i a closed syste ad Matheatical Modelig, i: Classroo Notes i Applied Matheatics, SIAM, Philadelphia, PA, (1989) 317-3.

67 15; (6): 6-67 [6] S. T. Mohyud-Di, A. Yıldırı, Y. Gūlaat, Aalytical solutio of Volterra s populatio odel, Joural of Kig Saud Uiversity (1) 47 5. [7] Arioglu A, Ozol I. Solutio of differetial?differece equatios by usig differetial trasfor ethod. Appl. Math. Coput, 181 (6) 153 16. [8] Arioglu A, Ozol I, Solutio of fractioal differetial equatios by usig differetial trasfor ethod. Chaos Solito. Fract, 34 (7) 1473 1481. [9] E. Celi, Kh. Tabatabaei, Solvig a class of volterra itegral equatio systes by the Differetial Trasfor Method, 9 (1) 1-5. [1] V. S. Ertür, Differetial trasforatio ethod for solvig differetial equatios of Lae-Ede type, Math. Coput. Appl, 1 (7) 135-139.