SHIFTED JACOBI COLLOCATION METHOD BASED ON OPERATIONAL MATRIX FOR SOLVING THE SYSTEMS OF FREDHOLM AND VOLTERRA INTEGRAL EQUATIONS

Similar documents
Jacobi Operational Matrix Approach for Solving Systems of Linear and Nonlinear Integro Differential Equations

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions

A new Approach for Solving Linear Ordinary Differential Equations

Haar Wavelet Collocation Method for the Numerical Solution of Nonlinear Volterra-Fredholm-Hammerstein Integral Equations

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

Septic B-Spline Collocation Method for the Numerical Solution of the Modified Equal Width Wave Equation

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

International Conference on Advanced Computer Science and Electronics Information (ICACSEI 2013) equation. E. M. E. Zayed and S. A.

Fixed point method and its improvement for the system of Volterra-Fredholm integral equations of the second kind

New Exact Traveling Wave Solutions for Two Nonlinear Evolution Equations

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

arxiv: v1 [math.co] 12 Sep 2014

EXACT TRAVELLING WAVE SOLUTIONS FOR THREE NONLINEAR EVOLUTION EQUATIONS BY A BERNOULLI SUB-ODE METHOD

A Hybrid Variational Iteration Method for Blasius Equation

Research Article Collocation Method Based on Genocchi Operational Matrix for Solving Generalized Fractional Pantograph Equations

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

SPECTRAL-COLLOCATION METHOD FOR FRACTIONAL FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS

Solution for singularly perturbed problems via cubic spline in tension

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

THE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Inexact Newton Methods for Inverse Eigenvalue Problems

NEW EXACT ANALYTICAL SOLUTIONS FOR THE GENERAL KDV EQUATION WITH VARIABLE COEFFICIENTS. Jiangsu, PR China

2.29 Numerical Fluid Mechanics

Numerical Heat and Mass Transfer

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method

Convexity preserving interpolation by splines of arbitrary degree

The binomial transforms of the generalized (s, t )-Jacobsthal matrix sequence

Research Article A Jacobi Dual-Petrov-Galerkin Method for Solving Some Odd-Order Ordinary Differential Equations

Numerical Solution of Nonlinear Singular Ordinary Differential Equations Arising in Biology Via Operational Matrix of Shifted Legendre Polynomials

Module 3: Element Properties Lecture 1: Natural Coordinates

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Haar wavelet collocation method to solve problems arising in induction motor

FUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM

Grid Generation around a Cylinder by Complex Potential Functions

Numerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline

High resolution entropy stable scheme for shallow water equations

Complex Solutions for the Fisher Equation and the Benjamin-Bona-Mahony Equation

1 Introduction We consider a class of singularly perturbed two point singular boundary value problems of the form: k x with boundary conditions

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

Solving a Class of Nonlinear Delay Integro differential Equations by Using Differential Transformation Method

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Numerical Solutions of Fractional Partial Differential Equations by Using Legendre Wavelets

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Deriving the X-Z Identity from Auxiliary Space Method

Copyright 2014 Tech Science Press CMC, vol.43, no.2, pp.87-95, 2014

ME 501A Seminar in Engineering Analysis Page 1

Significance of Dirichlet Series Solution for a Boundary Value Problem

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

Integral Transforms and Dual Integral Equations to Solve Heat Equation with Mixed Conditions

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

The Analytical Solution of a System of Nonlinear Differential Equations

Introduction to Density Functional Theory. Jeremie Zaffran 2 nd year-msc. (Nanochemistry)

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

A Compact Finite Difference Method for Solving the General Rosenau RLW Equation

Solution of Singularly Perturbed Differential Difference Equations Using Higher Order Finite Differences

Combined Wronskian solutions to the 2D Toda molecule equation

POLYNOMIAL BASED DIFFERENTIAL QUADRATURE FOR NUMERICAL SOLUTIONS OF KURAMOTO-SIVASHINSKY EQUATION

An efficient method to set up a Lanczos based preconditioner for discrete ill-posed problems

Haar Wavelet Based Numerical Solution of Elasto-hydrodynamic Lubrication with Line Contact Problems

On the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros

Solving 2D-BKDV Equation by a Sub-ODE Method

12. The Hamilton-Jacobi Equation Michael Fowler

A property of the elementary symmetric functions

Approximate Technique for Solving Class of Fractional Variational Problems

4DVAR, according to the name, is a four-dimensional variational method.

Computers and Mathematics with Applications

A FORMULA FOR COMPUTING INTEGER POWERS FOR ONE TYPE OF TRIDIAGONAL MATRIX

Chapter 3 Differentiation and Integration

Restricted divisor sums

Global Sensitivity. Tuesday 20 th February, 2018

Exact Solutions for Nonlinear D-S Equation by Two Known Sub-ODE Methods

Lecture 12: Discrete Laplacian

New Exact Solutions of the Kawahara Equation using Generalized F-expansion Method

New Method for Solving Poisson Equation. on Irregular Domains

Numerical Solution of Ordinary Differential Equations

Preconditioning techniques in Chebyshev collocation method for elliptic equations

Research Article Cubic B-Spline Collocation Method for One-Dimensional Heat and Advection-Diffusion Equations

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Finite Element Method: A Short Introduction

Implicit Integration Henyey Method

Spectral method for fractional quadratic Riccati differential equation

PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two

Numerical solution for systems of two dimensional integral equations by using Jacobi operational collocation method

Chebyshev Gauss Lobatto Pseudo spectral Method for One Dimensional Advection Diffusion Equation with Variable coefficients

The lower and upper bounds on Perron root of nonnegative irreducible matrices

1 GSW Iterative Techniques for y = Ax

A Note on the Numerical Solution for Fredholm Integral Equation of the Second Kind with Cauchy kernel

A NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT

Binomial transforms of the modified k-fibonacci-like sequence

Finite Element Modelling of truss/cable structures

arxiv: v1 [math.ho] 18 May 2008

Perfect Fluid Cosmological Model in the Frame Work Lyra s Manifold

Least-Squares Solutions of Generalized Sylvester Equation with Xi Satisfies Different Linear Constraint

Transcription:

Mathematcal and Computatonal Applcatons, Vol., o., pp. 76-93, 5 http://d.do.org/.9/mca-5-7 SHIFED JACOBI COLLOCAIO MEHOD BASED O OPERAIOAL MARIX FOR SOLVIG HE SYSEMS OF FREDHOLM AD VOLERRA IEGRAL EQUAIOS Abdollah Borhanfar and Khadeh Sadr Department of Appled Mathematcs, Faculty of Mathematcal Scences Mohaghegh Ardabl Unversty, Ardabl, Iran borhan@uma.ac.r, kh.sadr@uma.ac.r Abstract- hs paper ams to construct a general formulaton for the shfted Jacob operatonal matrces of ntegraton and product. he man am s to generalze the Jacob ntegral and product operatonal matrces to the solvng system of Fredholm and Volterra equatons. hese matrces together wth the collocaton method are appled to reduce the soluton of these problems to the soluton of a system of algebrac equatons. he method s appled to solve system of lnear and nonlnear Fredholm and Volterra equatons. Illustratve eamples are ncluded to demonstrate the valdty and effcency of the presented method. Also, several theorems, whch are related to the convergence of the proposed method, wll be presented. Key ords- Collocaton Method, Shfted Jacob Polynomals, System of Fredholm and Volterra Integral Equatons; Convergence; Operatonal Integral and Product Matrces. IRODUCIO Fndng the analytcal solutons of functonal equatons has been devoted of attenton of mathematcans nterest n recent years. Several methods are proposed to acheve ths purpose, such as [-]. Mathematcal modelng for many problems n dfferent felds, such as engneerng, chemstry, physcs and bology, leads to ntegral equatons or system of ntegral equatons. Several methods have been proposed to solve these problems. For eample, Varatonal teraton method [], dfferental transform method [], ystrom method [], Haar functons method [], Homotopy perturbaton method [3], Chebyshev wavelet method [4] and many others. Between of present methods, spectral methods have been used to solve dfferent functonal equatons, because of ther hgh accuracy and easy applyng. Specfc types of spectral methods that more applcable and wdely used, are the Galerkn, collocaton, and tau methods [5-9]. Saadatmand and Dehghan ntroduced shfted Legendre operatonal matr for fractonal dfferental equatons [3], Doha derved a new eplct formula for shfted Chebyshev polynomals for fractonal dfferental equatons [3], Bhrawy used a quadrature shfted Legendre tau method for fractonal dfferental equatons [3]. Recently, Doha ntroduced shfted Chebyshev operatonal matr and appled t wth spectral methods for solvng problems to ntal and boundary condtons [33]. he mportance of Sturm-Louvlle problems for spectral methods les n the fact that the spectral appromaton of the soluton of a functonal equaton s usually regarded as a fnte epanson of egenfunctons of a sutable Sturm-Louvlle problem.

Shfted Jacob Collocaton Method Based on Operatonal Matr 77 he Jacob polynomalsp ( )(,, ) play mportant roles n mathematcal analyss and ts applcatons [34]. It s proven that Jacob polynomals are precsely the only polynomals arsng as egenfunctons of a sngular Sturm-Louvlle problem [35-36]. hs class of polynomals comprses all the polynomal soluton to sngular Sturm-Louvlle problems on [,]. Chebyshev, Legendre, and ultrasphercal polynomals are partcular cases of the Jacob polynomals. In ths paper, the shfted Jacob operatonal matrces of ntegraton and product s ntroduced, whch s based on Jacob collocaton method for solvng numercally the systems of the lnear and nonlnear Fredholm and Volterra ntegral equatons on the nterval [,], to fnd the appromate soluton u ( ). he each of equaton of the systems resulted are collocated at ( ) nodes of the shfted Jacob- Gauss nterpolaton on(,). hese equatons generate n ( ) lnear or nonlnear algebrac equatons. he nonlnear systems resulted can be solved usng ewton teratve method. he remander of ths paper s organzed as follows: he Jacob polynomals and ther ntegral and product operatonal matrces to ntegral equatons are obtaned n Secton. Secton 3 s devoted to applyng the Jacob operatonal matrces for solvng system of ntegral equatons. In Secton 4, the proposed method s appled to several eamples. A concluson s presented n Secton 5.. JACOBI POLYOMIALS AD HEIR OPERAIOAL MARICES.. Propertes of shfted Jacob polynomals he Jacob polynomals, assocated wth the real parameters (, ) are a sequence of polynomals P ( t),(,,,...), each of degree, are orthogonal wth Jacob weghted functon, w( ) ( ) ( ) over I [,], and P ( t) P ( t) w( t) dt h, n m n nm where s Kroneker functon and h nm n ( n ) ( n ) (n ) n! ( n ) hese polynomals can be generated wth the followng recurrence formula; ( ) t( )( ) P ( t) P ( t) ( )( ) ( )( )( ) P ( ),, 3,..., t ( )( ) where, P ( t ), and P ( t) t..

7 A. Borhanfar and K. Sadr In order to use these polynomals on the nterval [,], shfted Jacob polynomals are defned by ntroducng the change of varable t. Let the shfted Jacob polynomalsp ( ) be denoted by from followng formula; R ( ), then R ( ) can be generated ( ) ( )( )( ) R ( ) R ( t) ( )( ) () ( )( )( ) R ( ),, 3,..., t ( )( ) where, R ( ), and R ( ) ( ). Remark. Of ths polynomals, the most commonly used are the shfted Gegenbauer polynomals, C ( ), the shfted Chebyshev polynomals of the frst knd, ( ), the S, S, shfted Legendre polynomals, P ( ), the shfted Chebyshev polynomals of the second S, knd, U ( ). hese orthogonal polynomals are related to the shfted Jacob S, polynomals by the followng relatons.! ( )! ( ) C ( ) R ( ), ( ) R ( ), S, S, ( ) ( ) ( )! ( ) (, ) P ( ) R ( ), U ( ) R ( ). S, S, 3 ( ) he analytc form of the shfted Jacob polynomals, R R (), s gven by k k ( ) ( ) ( k ) ( ), ( k ) ( )( k)! k! k Some propertes of the shfted Jacob polynomals are as follows, () R () ( ), () R () ( ), d ( n ) (, ) (3) ( ) ( ). R n R n d ( n ) he orthogonalty condton of shfted Jacob polynomals s R ( ) R ( ) ( ) d, (3) k k k ()

Shfted Jacob Collocaton Method Based on Operatonal Matr 79 where ( ), shfted weghted functon, s as follows, hk ( ) ( ), and. k Lemma. he shfted Jacob polynomal where Proof. p are ( n ) p ( n) p can be obtaned as, ( n) n Rn n ( n) n R ( ) p, n n ( ). n ( ) can be obtaned n the form of: ow, usng propertes () and (3) n above, the lemma can be proved. Lemma. For m, one has ( ) l l m R ( ) ( ) d p B( m l, ), where B( s, t ) s the Beta functon and s defned as Proof. Usng Lemma and ( n) d p R ( ). n! d s t ( s) ( t) ( ). B s t v v dv ( ) ( ) one has ( s t) m ( ) m l l l ( ) m l l l ( ) l l R ( ) ( ) d p ( ) d p ( ) d p B( m l, )... he appromaton of functons Let (,), and for r ( s the set of all non-negatve ntegers), the r weghted Sobolov space H ( ) s defned n the usual way and s denoted nner product, sem-norm and norm by ( uv, ), v and v, respectvely. In r, r, r, partcular, L ( ) H ( ), ( u, v ) ( u, v ) and v v,,,

A. Borhanfar and K. Sadr H ( ) f f s ntegrable and v, r r k ( k, k) r, k v v v r r, (, ) A functon u( ) H r ( ) can be epanded n (,, ). r, P span R ( ), R ( ),..., R ( ) as the below formula,, where the coeffcents c are gven by u( ) c R ( ), (4) c R ( ) u( ) ( ) d,,,,.... (5) By notng n practce, only the frst ( ) terms shfted polynomals are consdered, then one has where Snce from u( ) u ( ) c R ( ) ( ) C, (6) C [ c, c,..., c ], ( ) [ R ( ), R ( ),..., R ( )]. (7) (,, ) P s a fnte dmensonal vector space, u ( ) has a unque best appromaton (,, ) P, say (,, ) u () P, that s (,, ) y P, u( ) u ( ) u( ) y. In [37] s shown that for any u( ) H r ( ), r and r, a postve constant C ndependent of any functon,, and est that r, r, u ( ) u ( ) c ( ( )) u. () Let u ( ) s tmes contnuously dfferentable. he followng heorem can present an upper bound for estmatng the error. heorem. Let u( ) : [,] s tmes contnuously dfferentable for, and (,, ) P span R ( ), R ( ),..., R ( ). If u ( ) ( ) A s the best appromaton to u ( ) from follows: (,, ) P then the error bound s presented as

Shfted Jacob Collocaton Method Based on Operatonal Matr MS u( ) u ( ) B(, ), ( )! ( ) where M ma u ( ) and S ma,. [,] Proof. Let consder the aylor epanson ( ) ( ) y( ) u( ) u ( )( )... u ( )!. herefore u u u ( )! ( ) ( ) ( ) ( ), (,), Snce ( ) A s the best appromaton to u ( ) from P (,, ), and has (,, ) y () P one M ( ) ( ) ( ) ( ) ( ) ( ) ( ). (( )!) u u u y d () Snce s always postve n (,), by choosng S ma, one has ( ) ( ) M S M S ( ) ( ) ( ) (, ). (( )!) (( )!) u u d B hs error bound shows appromaton of polynomals converges to u ( ) as..3. he Jacob ntegral operatonal matr In ths subsecton, Jacob operatonal matr of the ntegraton s derved. Let ( t) dt P ( ), (9) where matr P s called the Jacob operatonal matr of the ntegraton. heorem 3. Let P s ( ) ( ) operatonal matr of ntegral. hen the elements of ths matr are obtaned as P p p B( m n, ),,,,.... ( ) ( ) m n m n m Proof. Usng Eq. (9) and orthogonalty property of Jacob polynomals one has, where ( ( t) dt, ( )) and follows, ( ( ), ( )), P t dt are two( ) ( ) matrces defned as

A. Borhanfar and K. Sadr Set t dt R t dt R, ( ( ), ( )) ( ( ), ( )), dag. ( R ( t) dt, R ( )) R ( t) dt R ( ) ( ) d. R () t dt and R ( ) by usng Lemma can be obtaned as m R ( t) dt p, m ( ) n n () m m n R ( ) p,,,,...,. herefore, by usng Lemma can be obtaned as follows, ( ) ( ) m n m n m n m ( ) ( ) p p B m n m n m n m So, the elements of matr P s obtaned as ( ) ( ) m n m n m p p ( ) d (, ). P p p B( m n, ),,,,.... ow the followng theorem can present an upper bound for estmatng the error of ntegral operator. he error vector E s defned as, where E ( t) dt P ( ) [ E, E,..., E ], k k k E R ( t) dt P R ( ), k,,...,. r heorem 4. If E R ( t) dt P R ( ) H ( ), then an error bound of ntegral operator of vector k k k can be epressed by k k r ( k) ( k) k, E c ( ( )) B( r 3, r ).

Shfted Jacob Collocaton Method Based on Operatonal Matr 3 Proof. By usng nequalty (), Lemma, and settng u( ) R ( t) dt one has where r, r () k ( r, r) k r ( k ) D p ( r, r) k! ( k ) r p ( r ) ( r, r) k k ( k ) r ( k ) r ( r, r) k k ( k ) ( k ) r r ( ) d k k ( k) ( k) B r r u D R t dt! p ( r ) ( k) ( k) ( )( ) ( ) d ( 3, ), and the theorem can be proved..4 he product operatonal matr he followng property of the product of two Jacob functon vector wll also be appled to solve the Volterra ntegral equatons. ( ) ( ) Y Y ( ) () where Y s a ( ) ( ) product operatonal matr and t`s elements are determned n terms of the vector Y `s elements. Usng Eq.() and by the orthogonalty property of Jacob polynomals the elements Y can be calculated as follows, Y Y ( ( )) ( ( )) ( ( )) ( ) d k k k Y h Y R ( ) R ( ) R ( ) ( ) d k k k k where h R ( ) R ( ) R ( ) ( ) d. k k k, k

4 A. Borhanfar and K. Sadr 3. APPLICAIOS OF HE OPERAIOAL MARICES OF IEGRAIO AD PRODUC In ths secton, the presented operatonal matrces are appled to solve the system of lnear and nonlnear Fredholm and Volterra ntegral equatons. 3.. he system of Fredholm ntegral equatons A system of Fredholm ntegral equatons can be presented as follows; m u ( ) f ( ) k (, t) G ( u ( t), u ( t),..., u ( t)) dt,,,,..., n, () n where k (, t) L [,] [,], f are known functons, and G are lnear or nonlnear functons n terms of unknown functons u ( ), u ( ),..., and u ( ). o solve system n (), the functons u ( ), G () t and k (, t) can be appromated as follows, u ( ) c R ( ) ( ) C, G ( t) ( t) Y, k (, t) ( ) K ( t), where K and vectors respectvely and Y are ( ) ( ),,..., n,,,..., m, () known matrces and ( ) unknown C [ c, c,..., c ], ( ) [ R ( ), R ( ),..., R ( )]. th substtutng appromatons (3) n system () one has m ( ) C f ( ) ( ) K ( t) ( t) Y dt m f ( ) ( ) K ( t) ( t) dt Y m f ( ) ( ) K DY,,,..., n, where D s the followng ( ) ( ) known matr, D ( t) ( t) dt. he system (3) have n ( ) unknown coeffcents (3) c. For collocatng, ( ) roots of Jacob polynomals R ( ) are appled and the equatons are collocated at ( ) them. Unknown coeffcents are determned wth solvng the resulted system of lnear or nonlnear algebrac equatons. 3.. System of Volterra ntegral equatons A system of Volterra ntegral equatons of the frst knd can be presented as follows,

Shfted Jacob Collocaton Method Based on Operatonal Matr 5 m f ( ) k (, t) G ( u ( t), u ( t),..., u ( t)) dt,,,,..., n. (4) n By usng the appromate relatons () one has m ( ) ( ) ( ) ( ) m ( ) ( ) ( ) m ( ) K Y ( t) dt m f K t t Y dt K t t Y dt ( ) K Y P ( ),,,..., n, where Y and P are product and ntegral operatonal matrces, respectvely. Also, a system of Volterra ntegral equatons of the second knd can be presented as m u ( ) f ( ) k (, t) G ( u ( t), u ( t),..., u ( t)) dt,,,,..., n. (6) n In the same way, the system of followng equatons s resulted. m (5) ( ) C f ( ) ( ) K Y P ( ),,,..., n, (7) By usng the frst ( ) roots of Jacob polynomals (7), unknown coeffcents c are determned. R ( ) and collocated system ( ) 4. ILLUSRAIVE EXAMPLES In ths secton, the presented method s appled to solve some eamples. Comparson between the results of present method wth the correspondng analytc solutons s gven. For ths purpose, the mamum of absolute error s computed. Eample. he followng system of lnear Volterra ntegral equatons of the second knd s consdered, 3 3 3 3 t 3 3 3 ( 3) u ( ) f ( ) ( t) u ( t) dt ( t ) u ( t) dt t u ( t) dt, ( 3 ) u ( ) f ( ) t( ) u ( t) dt t( ) u ( t) dt ( t ) u ( t) dt, () (3 6) u ( ) f ( ) ( t ) u ( t) dt t( t ) u ( t) dt (t ) u ( t) dt,

6 A. Borhanfar and K. Sadr he eact solutons are u ( ) 5, u ( ) 5 and u ( ). By 3 the applyng the technque descrbed n pervous secton wth 4, solutons and kernels are appromated as: u ( ) ( ) C, u ( ) ( ) C, u ( ) ( ) C, t ( ) K ( t ), 3 3 t K t t K t t K t 3 4 ( ) ( ), ( ) ( ), ( ) ( ) ( ), t( ) ( ) K ( t), t ( ) K ( t), t ( ) K ( t ), 3 5 6 7 t( t ) ( ) K ( t), t t ( ) K ( t ). 3 9 he system () by usng above equatons s rewrtten as, ( 3) ( ) C ( ) K C K C K C P ( ), 3 3 ( 3 ) ( ) C ( ) K C K C K C P ( ), 4 5 6 3 (3 6) ( ) C ( ) K C K C K C P ( ), 3 7 9 3 where C, C and C are the operatonal matrces of product correspondng to unknown 3 vectors C, C and C. ow usng the roots of R ( ) and collocatng the system 3 5 (), reduces the problem to solve a system of algebrac equatons. Unknown coeffcents are obtaned for some values of parameters and. Mamum absolute error for 4 and dfferent values of and has been lsted n able. () able. Mamum absolute error for 4 and dfferent values of and for Eample. Error( ) Error( ) Error( u ) 4.33 / / 7.7 / /.7 / / 6. / / 4.63 / 4 / 4 4.435. /. / 3/ 3. u u 3 5.4.335.769 3.7.937.95 3.6.9 4. 3.9353 6.3.3 5..75.333.7659.996 6. able shows that a good appromaton can be acheved for the eact solutons by usng a few terms of shfted Jacob polynomals for varous values of parameters and.

Shfted Jacob Collocaton Method Based on Operatonal Matr 7 Eample. Consder the followng system of lnear Fredholm ntegral equaton of second knd, u ( ) 3 u ( ) f ( ) ( t) u ( t) dt t u ( t) dt, 3 u ( ) 4 u ( ) f ( ) (t ) u ( t) dt ( t) u ( t) dt, he eact solutons are u () e and u (). th 5, solutons and kernels are appromated as: u ( ) ( ) C, u ( ) ( ) C, t ( ) K ( t), t ( ) K ( t ), () 3 4 t ( ) K ( t), t ( ) K ( t), D ( t) ( t) dt. he system () by usng above equatons s rewrtten as, ( ) C 3 ( ) C ( ) K DC ( ) K DC f ( ), 3 4 3 ( ) C 4 ( ) C ( ) K DC ( ) K DC f ( ). () ow usng the roots of R ( ) and collocatng the system (), reduces the problem 6 to solve a system of lnear algebrac equatons. Solvng the system (), the unknown coeffcents wll be obtaned. able dsplays the mamum absolute errors for varous and wth 5. able. Mamum absolute error for 5 and dfferent values of and for Eample. Error( ) Error( u ) 5. / / 6.96 5 / /.95 5 / / 3.97 5 / 4 / 4.46 5 / / 5.63 6.556 4.6 3 3/ 4. u 3.9674 9.36.54.9 6.759 3.45.795 4.743 4.647

A. Borhanfar and K. Sadr Eample 3. Consder the followng system of lnear Fredholm ntegral equaton of second knd, u ( ) f ( ) t cos( ) u ( t) dt c sn( t) u ( t) dt, t ( ) ( ) ( ) ( ) ( ). u f e u t dt t u t dt he eact solutons are u ( ) and u ( ) cos( ). th 5, solutons and () kernels are appromated as: u ( ) ( ) C, u ( ) ( ) C, t cos( ) ( ) K ( t), sn( t) ( ) K ( t), ep( t ) ( ) K ( t), ( t) ( ) K ( t), D ( t) ( t) dt. 3 4 he system () by usng above equatons s rewrtten as, ( ) C ( ) K DC ( ) K DC f ( ), 3 4 ( ) C ( ) K DC ( ) K DC f ( ). (3) ow usng the roots of R ( ) and collocatng the system (3), reduces the problem 6 to solve a system of lnear algebrac equatons. Solvng the system (3), the unknown coeffcents wll be obtaned. able 3 dsplays the mamum absolute errors for varous and wth 5. able 3. Mamum absolute error for 5 and dfferent values of and for Eample 3. Error( ) Error( u ).63 / /.966 / / 7.797 / /.57 / 4 / 4.5 / / 7.9 4.795 9.45 3/ 4.9567 7 u.9 7.5. 4. 4... 5. 7.69 7

Shfted Jacob Collocaton Method Based on Operatonal Matr 9 Eample 4. In ths eample, the followng nonlnear Volterra ntegral equaton of frst knd s consdered, ( u ( t) ( t) u ( t) u ( t)) dt f ( ), ( u ( t) ( t) u ( t) u ( t)) dt f ( ). he eact solutons are () u e and () u e. th 5, solutons and kernels are appromated as: u ( ) ( ) C, u ( ) ( ) C, ( ) K ( t), t ( ) K ( t), u ( ) u ( ) C C ( ), he system (4) by usng above equatons s rewrtten as, ( ) K C K E P ( ) f ( ), ( ) K C K E P ( ) f ( ), where C, C and E are the operatonal matrces of product correspondng wth the vectors C, C and E C C, respectvely. ow usng the roots of R ( ) and collocatng the each equaton of system (5), reduces the problem to solve a system of nonlnear algebrac equatons. Solvng the system (5) by ewton teratve method, the unknown coeffcents wll be obtaned. able 4 dsplays the mamum absolute errors for varous and wth. able 4. Mamum absolute error for and dfferent values of and for Eample 4. Error( ) Error( u ) 3.374 / / 5.4 / / 7.76 / /.67 / 4 / 4 4.57 / / 3.574 6.69 9.536 3/ 4.465 u 3.374 5.4 7.76.67.67 3.574 6.69 9.536.465 (4) (5) Eample 5. Ffth eample covers the system of nonlnear Volterra ntegral equaton of second knd,

9 A. Borhanfar and K. Sadr f ( ) u ( t) u ( t) dt u ( ), f ( ) u ( t) dt u ( t) dt u ( ). he eact solutons are u ( ) cos( ) and u ( ) sn( ). Solutons and kernels are appromated as: u ( ) ( ) C, u ( ) ( ) C, ( ) K ( t ), (6) u ( ) u ( ) C C ( ), u ( ) C C ( ), u ( ) C C ( ), he system (6) by usng above equatons s rewrtten as, ( ) C ( ) K E P ( ) f ( ), 3 ( ) C ( ) K E E P ( ) f ( ), where E, E and E are the operatonal matrces of product correspondng wth the 3 vectors, CC (7) CCand E C C, respectvely. he mamum absolute errors for and 7,,5 are lsted n able 5. Also, usng the roots of R ( ) 6 and collocatng the each equaton of system (7), reduces the problem to solve a system of nonlnear algebrac equatons. Solvng the system (7) by ewton teratve method, the unknown coeffcents wll be obtaned. able 6 dsplays the mamum absolute errors for varous and wth 5. able 5. Mamum absolute error for 7,,5 and for Eample 5. Error( ) 7.965 3.745 5. u u Error( ).679 6.637. able 6. Mamum absolute error for 5 and dfferent values of and for Eample 5. Error( ) Error( u ) / / 6.747 / / 4.366 / /.7 4.44 3.57 4 4.569 / 4 / 3.474 9 4 u.556.754.6.435.47.7 5.65

Shfted Jacob Collocaton Method Based on Operatonal Matr 9 5. COCLUIO In ths paper, the shfted Jacob collocaton method was employed to solve a class of systems of Fredholm and Volterra ntegral equatons of frst and second knds. Frst, a general formulaton for the Jacob operatonal matr of ntegral has been derved. hs matr s used to appromate numercal soluton of system of lnear and nonlnear Volterra ntegral equatons. Proposed approach was based on the shfted Jacob collocaton method. he solutons obtaned usng the proposed method shows that ths method s a powerful mathematcal tool for solvng the ntegral equatons. Provng the convergence of the method, consstency and stablty are ensured automatcally. Moreover, only a small number of shfted Jacob polynomals are needed to obtan a satsfactory result. 6. ACKOLEDGEMES e would lke to epress our grattude to the anonymous revewer for the careful readng of the manuscrpt and for hs nvaluable comments that helped us to mprove t consderably. 7. REFERECES. A. Borhanfar and M. M. Kabr, ew perodc and solton solutons by applcaton of applcaton of ep-functon method for nonlnear evoluton equaton, Computatonal and Appled Mathematcs, 9, 5-67, 9.. A. Borhanfar, M. M. Kabr and L. Maryam Vahdat, ew perodc and solton wave solutons for the generalzed zakharov system and (+)-dmensonal nzhnk-novkovveselov system, Chaos, Soltons and Fractals, 4, 646-654, 9. 3. A. Borhanfar and M. M. Kabr, Solton and perodc solutons for (3+)-dmensonal nonlnear evoluton equatons by ep-functon method, Applcaton and Appled Mathemahcs, 9, 5-67, 9. 4. A. Borhanfar and A. zamr, Applcaton of G / G epanson method for the Zhber -Shabat equaton and other related equatons, Mathematcal and Computer Modellng, 54, 9-6,. 5. A. Borhanfar, A. zamr and M. M. Kabr, Eact travellng wave solutons for the generalzed shallow water wave (GS) equaton, Mddle-East Journal of Scentfc Research, (3), 3-35,. 6. A. Borhanfar, M. M. Kabr and R. Abazar, Applcaton of G / G epanson method to regularzed (RL) equaton, Computer and Mathematcs wth Applcatons, 6, 44-47,. 7. A. Borhanfar and R. Abazar, General soluton of generalzed (+)-dmensonal Kadomtsev-Petvashvl (KP) equaton by usng the G / G epanson method, Journal of Computer Mathematcs,, -5,.. A. Borhanfar, H. Jafar and S. A. Karm, ew soltary wave solutons for generalzed and regularzed long-wave equaton, Internatonal Journal of Computer Mathematcs, 7 (3), 59-54,.

9 A. Borhanfar and K. Sadr 9. A. Borhanfar, H. Jafar and S. A. Karm, ew soltary wave solutons for the bad boussnesq and good boussnesq equatons, umercal Methods for partal Dfferental Equatons, 5, 3-37, 9.. A. Borhanfar, H. Jafar and S. A. Karm, ew soltons and perodc solutons for the Kadomtsev-Petvashvl equaton, onlnear Scences and Applcatons, 4, 4-9,.. A. Borhanfar and R. Abazar, Eact solutons for nonlnear Schrodnger equatons by dfferental transformaton method, Appled Mathematcs and Computng, -5, 9.. A. Borhanfar and R. Abazar, umercal study of nolnear Schrodnger and cou transformaton method, Optcs Communcatons, 3, 6-3,. 3. J. Bazar, E. Babolan and R. Islam, Soluton of the system of ordnary dfferental equatons by Adoman decomposton method, Appled Mathematcs and Computaton, 47, 73-7, 4. 4. E. Yusufoglu, An mprovement to Homotopy perturbaton method for solvng system of lnear equatons, Computers and Mathematcs wth Applcatons, 5, 3-35, 9. 5. M. M. Al-Sawalha and M. S. M. oor, Applcaton of the dfferental transformaton method for the soluton of the hyperchaotc Rossler system, Communcatons n onlnear Scence and umercal Smulaton, 4, 59-54, 9. 6. A. Golbaba and M. Javd, Applcaton of Homotopy perturbaton method for solvng eght-order boundary value problems, Appled Mathematcs and Computaton,, 334-346, 7. 7. J.H. He, Varatonal teraton method for autonomous ordnary dfferental systems, Appled Mathematcs and Computaton, 4, 5-3,.. M. atar and M. Dehghan, Improvement of He's varatonal teraton method for solvng systems of dfferental equatons, Computers and Mathematcs wth Applcatons, 5, 6-66, 9.. L. Xu, Varatonal teraton method for solvng ntegral equatons, Computers and Mathematcs wth Applcatons, 54, 7-7, 7.. A. Arkoglu and I. Ozkol, Solutons of ntegral and ntegro-dfferental equaton systems by dfferental transform method, Computers and Mathematcs wth Applcatons, 56, 4-47,.. M. C. De Bons and C. Laurta, umercal treatment of second knd Fredholm ntegral equatons on bounded ntervals, Journal of Computatonal and Appled Mathematcs, 7, 64-7,.. F. Mrzaee, umercal computatonal soluton of the lnear Volterra ntegral equatons system va ratonalzed Haar functons, J. Kng. Saud. Unversty, (Scence),, 65-6,. 3. J. Bazar and M. Eslam, Modfed HPM for solvng systems of Volterra ntegral equatons of the second knd, J. Kng. Saud. Unversty. (Scence), 3, 35-39,. 4. H. Adb and P. Assar, Chebyshev avelet Method for umercal Soluton of Fredholm Integral Equaton of the Frst Knd, Mathematcal Problems n Engneerng, do:.55//34,.

Shfted Jacob Collocaton Method Based on Operatonal Matr 93 5. J. P. Boyd, Chebyshev and Fourer Spectral Methods, second ed., Dover Publcatons, Mneola,. 6. E. A. Coutsas,. Hagstrom, and D. orres, An effcent spectral method for ordnary dfferental equatons wth ratonal functon, Mathematcs of Computaton., 65, 6-635, 96. 7. D. Funaro, Polynomal Appromatons of Dfferental Equatons, Sprnger-Verlag, Berln, 9.. D. Gottleb and S. A. Orszag, umercal Analyss of Spectral Methods: heory and Applcatons, SIAM, Phladelpha, 97. 9. E. H. Doha and A. H. Bhrawy, Effcent spectral-galerkn algorthms for drect soluton of fourth-order dfferental equatons usng Jacob polynomals, Appled umercal Mathematcs, 5, 4-44,. 3. A. Saadatmand and M. Dehghan, A new operatonal matr for solvng fractonalorder dfferental equatons, Computers and Mathematcs wth Applcatons, 59, 36-336,. 3. E. H. Doha, A. H. Bhrawy and S. S. Ezz-Elden, Effcent Chebyshev spectral methods for solvng mult-term fractonal orders dfferental equatons, Appled Mathematcal Modellng, 35, 566-567,. 3. A. H. Bhrawy, A. S. Alof and S. S. Ezz-Elden, A quadrature tau method for varable coeffcents fractonal dfferental equatons, Appled Mathematcs Letters, 4, 46-5,. 33. E. H. Doha, A. H. Bhrawy and S. S. Ezz-Elden, A Chebyshev spectral method based on operatonal matr for ntal and boundary value problems of fractonal order, Computers and Mathematcs wth Applcatons, 6, 364-373,. 34. C. Canuto, M. Y. Hussan, A. Quarteron and. A. Zang, Spectral Methods n Flud Dynamcs, Sprnger-Verlag, ew York, 9. 35. E. H. Doha, On the constructon of recurrence relatons for the epanson and connecton coeffcents n seres of Jacob polynomals, Journal of Physcs A: Mathematcal and General, 37, 657-675, 4. 36. Y. Luke, he Specal Functons and her Appromatons, vol., Academc Press, ew York, 69. 37. B. Y. Guo and L. L., ang, Jacob appromaton n non-unformly Jacobweghted Sobolov spaces, Journal of Appromaton heory,, -4, 4.