Wednesday 11 May 2022

NUMERICAL SOLUTION OF NONLINEAR VOLTERRA-FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS VIA DIRECT METHOD USING TRIANGULAR FUNCTIONS

 Dr.K.R.Salini[1] and Dr.S.Bamini[2]

[1] Guest Lecturer, PG & Department of Mathematics, Government Arts & Science College, Kariyampatti, Tirupattur, Tamilnadu, India.

Email: shalinimanoj1985@gmail.com

[2] Research Coordinator, Assistant Professor, PG & Research Department of Mathematics, Marudhar Kesari Jain College For Women, Vaniyambadi, Tamilnadu, India.

Email: saransham07@gmail.com

Abstract

We determine the numerical solution of the specific nonlinear Voltera-Fredholm Integro-Differential Equation is proposed. The method is based on new vector forms for representation of triangular functions and its operational matrix.

Introduction

A set of basic functions and an appropriate projection method or a direct method. These methods often transform an integro-differential equation to a linear or nonlinear system of algebraic equations which can be solved by direct or iterative methods.

Consider a Volterra-Fredhlom Integro-differential equations of the form π‘₯′(𝑠)+π‘ž(𝑠)π‘₯(𝑠)+πœ†1∫π‘˜1(𝑠,𝑑)𝐹(π‘₯(𝑑))𝑑𝑑10+πœ†2∫π‘˜2(𝑠,𝑑)𝐺(π‘₯(𝑑))𝑑𝑑10=𝑦(𝑠)

π‘₯(0)=π‘₯0

Where the functions 𝐹(π‘₯(𝑑)) and 𝐺(π‘₯(𝑑)) are polynomials of π‘₯(𝑑) with constant coefficients.

For convenience, we put 𝐹(π‘₯(𝑑))=[π‘₯(𝑑)]𝑛1 and 𝐺(π‘₯(𝑑))=[π‘₯(𝑑)]𝑛2, where 𝑛1 & 𝑛2 are positive integers.

For 𝑛1,𝑛2=1, equation (1.1) is a linear integro- differential equation.

We present new vector forms of triangular functions (TFs), operational matrix of integration, expansion of functions of one and two variables with respect to TFs and other TFs properties.

By using new representations, a nonlinear integro-differential equation can be easily reduced to a nonlinear system of algebraic equations.

Journal of the Maharaja Sayajirao University of Baroda ISSN: 0025-0422

68

Volume-56, No.1 (IV) 2022

Triangular Functions

Definition

Two π‘š-sets of triangular functions are defined over the intercal [0,𝑑) as 𝑇1𝑖(𝑑)={1−𝑑−π‘–β„Žβ„Žπ‘–β„Ž≤𝑑<(𝑖+1)β„Ž0π‘œπ‘‘β„Žπ‘’π‘Ÿπ‘€π‘–π‘ π‘’ 𝑇2𝑖(𝑑)={𝑑−π‘–β„Žβ„Žπ‘–β„Ž≤𝑑<(𝑖+1)β„Ž0π‘œπ‘‘β„Žπ‘’π‘Ÿπ‘€π‘–π‘ π‘’

Where 𝑖=0,1,2…π‘š−1 with a positive integer value for m.

Also consider β„Ž=𝑇/π‘š, and 𝑇1𝑖 as the π‘–π‘‘β„Ž left-handed triangular function and 𝑇2𝑖 as

the π‘–π‘‘β„Ž right-handed triangular function.

It is assumed that 𝑇=1, so TFs are defined over [0,1) and β„Ž=1/π‘š.

From the definition of TFs, it is clear that triangular functions are disjoint, orthogonal

and complete.

We can write ∫𝑇1𝑖10(𝑑)𝑇1𝑗(𝑑)𝑑𝑑=∫𝑇2𝑖10(𝑑)𝑇2𝑗(𝑑)𝑑𝑑={β„Ž3𝑖=𝑗0𝑖≠𝑗

Also, ∅𝑖(𝑑)=𝑇1𝑖(𝑑)+𝑇2𝑖,𝑖=0,1,…π‘š−1

Where ∅𝑖(𝑑) is the 𝑖th block-pulse function defined as ∅𝑖(𝑑)={1π‘–β„Ž≤𝑑<(𝑖+1)β„Ž0π‘œπ‘‘β„Žπ‘’π‘Ÿπ‘€π‘–π‘ π‘’

Where 𝑖=0,1,…π‘š−1.

Vector Forms

Consider the first π‘š terms of the left-handed triangular functions and the first π‘š terms of the right- handed triangular functions and write them concisely as π‘š-vectors. 𝑇1(𝑑)=[𝑇10(𝑑),𝑇11(𝑑),…𝑇1π‘š−1(𝑑)]𝑇

𝑇2(𝑑)=[𝑇20(𝑑),𝑇21(𝑑),…𝑇2π‘š−1(𝑑)]𝑇

Where 𝑇1(𝑑) and 𝑇2(𝑑) are called left-handed triangular functions (LHTF) and right-handed triangular functions (RHTF) vector, respectively.

The product of two TFs vectors are

Journal of the Maharaja Sayajirao University of Baroda ISSN: 0025-0422

69

Volume-56, No.1 (IV) 2022

𝑇1(𝑑)𝑇1𝑇(𝑑)≅( 𝑇10(𝑑)0⋯⋯⋯00𝑇11(𝑑)⋯⋯⋯0⋮⋮0⋮⋮0⋮⋮⋯⋮⋮⋯⋮⋮⋯⋮⋮𝑇1π‘š−1(𝑑))

𝑇2(𝑑)𝑇2𝑇(𝑑)≅( 𝑇20(𝑑)0⋯⋯⋯00𝑇21(𝑑)⋯⋯⋯0⋮⋮0⋮⋮0⋮⋮⋯⋮⋮⋯⋮⋮⋯⋮⋮𝑇2π‘š−1(𝑑)) 𝑇1(𝑑)𝑇2𝑇(𝑑)=0

𝑇2(𝑑)𝑇1𝑇(𝑑)=0

Where 0 is the π‘§π‘’π‘Ÿπ‘œ π‘š×π‘š matrix.

Also ∫𝑇110(𝑑)𝑇1𝑇(𝑑)𝑑𝑑=∫𝑇2(𝑑)𝑇1𝑇(𝑑)𝑑𝑑10≅β„Ž3𝐼

∫𝑇110(𝑑)𝑇2𝑇(𝑑)𝑑𝑑=∫𝑇2(𝑑)𝑇1𝑇(𝑑)𝑑𝑑10≅β„Ž6𝐼

In which 𝐼 is an π‘š×π‘š identity matrix.

TFs Expansion

The expansion of a function 𝑓(𝑑) over [0,1) with respect to TFs written as 𝑓(𝑑)≅Ξ£π‘π‘–π‘š−1𝑖=0𝑇1𝑖(𝑑)+Ξ£π‘‘π‘–π‘š−1𝑖=0𝑇2𝑖(𝑑) 𝑓(𝑑)=𝑐𝑇𝑇1(𝑑)+𝑑𝑇𝑇2(𝑑)

Where we may put 𝑐𝑖=𝑓(π‘–β„Ž) and 𝑑𝑖=𝑓((𝑖+𝑖)β„Ž) for 𝑖=0,1,2,..π‘š−1.

Operational Matrix of Integration

Expressing ∫𝑇1(𝜏)π‘‘πœ≅𝑃1𝑇1(𝑠)+𝑃2𝑇2(𝑠)𝑠0

∫𝑇2(𝜏)π‘‘πœ≅𝑃1𝑇1(𝑠)+𝑃2𝑇2(𝑠)𝑠0

Where 𝑃1π‘š×π‘š and 𝑃2π‘š×π‘š are called operational matrices of integration in TFs domain and represented as follows 𝑃1=β„Ž2( 011⋯1001⋯10⋮00⋮00⋮0⋯⋱⋯1⋮0)

Journal of the Maharaja Sayajirao University of Baroda ISSN: 0025-0422

70

Volume-56, No.1 (IV) 2022

𝑃2=β„Ž2( 111⋯1011⋯10⋮00⋮00⋮0⋯⋱⋯1⋮1)

So the integral of any function 𝑓(𝑑) can be approximated as ∫𝑓(𝜏)π‘‘πœ≅𝑐𝑇𝑇1(𝑠)+𝑑𝑇𝑇2(𝑠)π‘‘πœπ‘ 0 ∫𝑓(𝜏)π‘‘πœ≅(𝑐+𝑑)𝑇𝑃1𝑇1(𝑠)+(𝑐+𝑑)𝑇𝑃2𝑇2(𝑠)𝑠0

New Representation of TFs Vector Forms and Other Properties

We define a new representation of TFs vector forms. Then, some characteristics of TFs are presented using the new definition.

Definition and Expansion

Let 𝑇(𝑑) be a 2π‘š-vector defined as 𝑇(𝑑)=(𝑇1(𝑑)𝑇2(𝑑)),0<𝑑<1

Where 𝑇1(𝑑) and 𝑇2(𝑑) have been defined in (1.2).

Now, the expansion of 𝑓(𝑑) with respect to TFs can be written as 𝑓(𝑑)≅𝐹1𝑇𝑇1(𝑑)+𝐹2𝑇𝑇2(𝑑)

≅𝐹𝑇𝑇(𝑑)

Where 𝐹1 &𝐹2 are TFs co-efficients with

𝐹1𝑖=𝑓(π‘–β„Ž)& 𝐹2𝑖=𝑓((𝑖+1)β„Ž),

For 𝑖=0,1,…π‘š−1.

Also, 2π‘š-vector 𝐹 is defined as 𝐹=(𝐹1𝐹2)

Now, assume that π‘˜(𝑠,𝑑) is a function of two variables. It can be expanded with respect to TFs as follows

π‘˜(𝑠,𝑑)≅𝑇𝑇(𝑠)𝐾𝑇(𝑑)

Where 𝑇(𝑠)& 𝑇(𝑑) are 2π‘š1&2π‘š2 dimensional triangular functions and 𝐾 is a

2π‘š1×2π‘š2 TFs coefficient matrix.

For convenience, we put π‘š1=π‘š2=π‘š. So matrix 𝐾 can be written as

Journal of the Maharaja Sayajirao University of Baroda ISSN: 0025-0422

71

Volume-56, No.1 (IV) 2022

𝐾=((𝐾11)π‘š×π‘š(𝐾12)π‘š×π‘š(𝐾21)π‘š×π‘š(𝐾22)π‘š×π‘š)

Where 𝐾11,𝐾12,𝐾21,𝐾22 can be compted by sampling the function π‘˜(𝑠,𝑑) at points 𝑠𝑖 & 𝑑𝑗 such that 𝑠𝑖=π‘–β„Ž & 𝑑𝑗=π‘—β„Ž, for 𝑖,𝑗=0,1,…π‘š. (𝐾11)𝑖𝑗=π‘˜(𝑠𝑖,𝑑𝑗),𝑖=0,1,…π‘š−1,𝑗=0,1,…π‘š−1 (𝐾12)𝑖𝑗=π‘˜(𝑠𝑖,𝑑𝑗),𝑖=0,1,…π‘š−1,𝑗=0,1,…π‘š (𝐾21)𝑖𝑗=π‘˜(𝑠𝑖,𝑑𝑗),𝑖=0,1,…π‘š,𝑗=0,1,…π‘š−1 (𝐾22)𝑖𝑗=π‘˜(𝑠𝑖,𝑑𝑗),𝑖=0,1,…π‘š,𝑗=0,1,…π‘š

Product Properties

Let 𝑋 be an 2π‘š-vertor which can be written as 𝑋𝑇=(𝑋1𝑇 𝑋2𝑇) such that 𝑋1 & 𝑋2 are π‘š-vectors.

It can be concluded from (1.3) & (1.4) that 𝑇(𝑑)𝑇𝑇(𝑑)𝑋=(𝑇1(𝑑)𝑇2(𝑑))(𝑇𝑇(𝑑)𝑇2𝑇(𝑑))(𝑋1𝑋2)

≅(𝑑𝑖𝑔(𝑇1(𝑑))π‘œπ‘š×π‘šπ‘œπ‘š×π‘šπ‘‘π‘–π‘”(𝑇2(𝑑)))(𝑋1𝑋2)

=𝑑𝑖𝑔 (𝑇(𝑑))𝑋

=𝑑𝑖𝑔(𝑋)𝑇(𝑑)

Therefore,

𝑇(𝑑)𝑇𝑇(𝑑)𝑋≅𝑋̅ 𝑇(𝑑)

Where 𝑋̅=𝑑𝑖𝑔(𝑋) is an 2π‘š×2π‘š diagonal matrix.

Let 𝐡 be a 2π‘š×2π‘š matrix as 𝐡=((𝐡11)π‘š×π‘š(𝐡12)π‘š×π‘š(𝐡21)π‘š×π‘š(𝐡22)π‘š×π‘š)

So, it can be similarly concluded from equations (1.3) and (1.4) that 𝑇(𝑑)𝐡𝑇(𝑑)=𝑇1𝑇(𝑑)𝑇2𝑇(𝑑)(𝐡11𝐡12𝐡21𝐡22)(𝑇1(𝑑)𝑇2(𝑑)) ≅𝑇1𝑇(𝑑)𝐡11𝑇1(𝑑)+𝑇2𝑇(𝑑)𝐡22𝑇2(𝑑) ≅𝐡̂11𝑇𝑇1(𝑑)+𝐡̂22𝑇𝑇2(𝑑)

Where 𝐡̂11 and 𝐡̂22 are π‘š−vectors with elements equal to the diagonal entries of matrices 𝐡11 and 𝐡22 respectively.

Therefore,

𝑇𝑇(𝑑)𝐡𝑇(𝑑)≅𝐡̂ 𝑇𝑇(𝑑)

Journal of the Maharaja Sayajirao University of Baroda ISSN: 0025-0422

72

Volume-56, No.1 (IV) 2022

In which 𝐡̂ is a 2π‘š- vector with elements equal to the diagonal entries of matrix 𝐡.

It is immediately concluded from equation (1.5) that ∫𝑇10(𝑑)𝑇𝑇(𝑑)𝑑𝑑=∫(𝑇1(𝑑)𝑇2(𝑑))𝑇1𝑇(𝑑)𝑇2𝑇(𝑑)𝑑𝑑10

=∫(𝑇1(𝑑)𝑇1𝑇(𝑑)𝑇1(𝑑)𝑇2𝑇(𝑑)𝑇2(𝑑)𝑇1𝑇(𝑑)𝑇2(𝑑)𝑇2𝑇(𝑑))𝑑𝑑10

≅(β„Ž3πΌπ‘š×π‘šβ„Ž6πΌπ‘š×π‘šβ„Ž6πΌπ‘š×π‘šβ„Ž3πΌπ‘š×π‘š)

Therefore,

∫𝑇10(𝑑)𝑇𝑇(𝑑)𝑑𝑑≅𝐷

Where 𝐷 is the following 2π‘š×2π‘š matrix. 𝐷=β„Ž3( 10⋯012⁄0⋯001⋯0012⁄⋯0⋮012⁄0⋮0⋮0012⁄⋮0⋱⋯⋯⋯⋱⋯⋮100⋮12⁄⋮010⋮0⋮001⋮0⋱…⋯⋯⋱⋯⋮12⁄00⋮1)

Operational Matrix

Expressing ∫𝑇(𝜏)π‘‘πœπ‘ 0 in terms of 𝑇(𝑠), we can write ∫𝑇(𝜏)π‘‘πœπ‘ 0=∫(𝑇1(𝜏)𝑇2(𝜏))π‘‘πœπ‘ 0

≅(𝑃1𝑇1(𝑠)+𝑃2𝑇2(𝑠)𝑃1𝑇1(𝑠)+𝑃2𝑇2(𝑠))

=(𝑃1𝑃2𝑃1𝑃2)(𝑇1(𝑠)𝑇2(𝑠))

So, ∫𝑇(𝜏)π‘‘πœπ‘ 0≅𝑃𝑇(𝑠)

Where 𝑃2π‘š×2π‘š, operational matrix of 𝑇(𝑠) is

𝑃=(𝑃1𝑃2𝑃1𝑃2) (

Journal of the Maharaja Sayajirao University of Baroda ISSN: 0025-0422

73

Volume-56, No.1 (IV) 2022

Where 𝑃1 & 𝑃2 are given by

The integral of any function 𝑓(𝑑) can be approximated as ∫𝑓(𝜏)π‘‘πœπ‘ 0≅∫𝐹𝑇𝑇(𝜏)π‘‘πœπ‘ 0

≅𝐹𝑇𝑃𝑇(𝑠)

Solving Nonlinear Integro-Differential Equation

Consider the following nonlinear Volterra-Fredhlom integro-differential equation

{π‘₯′(𝑠)+π‘ž(𝑠)π‘₯(𝑠)+πœ†1∫π‘˜1(𝑠,𝑑)[π‘₯(𝑑)]𝑛1𝑑𝑑+𝑠0πœ†2∫π‘˜2(𝑠,𝑑)[π‘₯(𝑑)]𝑛2𝑑𝑑=𝑦(𝑠)𝑠0π‘₯(0)=π‘₯0, 0≤𝑠<1,𝑛1,𝑛2≥1

Where the parameters πœ†1 and πœ†2 and β„’2 functions π‘ž(𝑠),𝑦(𝑠),π‘˜1(𝑠,𝑑),π‘˜2(𝑠,𝑑) are known but π‘₯(𝑠) is not.

The appearance of initial condition equation. This is necessary to ensure the existence of a solution.

Approximating functions

π‘₯(𝑠),π‘₯′(𝑠),π‘ž(𝑠),𝑦(𝑠),[π‘₯(𝑑)]𝑛1,[π‘₯(𝑑)]𝑛2,π‘˜1(𝑠,𝑑),π‘˜2(𝑠,𝑑)

with respect to TFs, π‘₯(𝑠)≅π‘₯𝑇(𝑠)𝑋 π‘₯′(𝑠)≅𝑋′𝑇𝑇(𝑠)=𝑇𝑇𝑋′

π‘ž(𝑠)≅ 𝑄𝑇𝑇(𝑠)=𝑇𝑇(𝑠)𝑄

𝑦(𝑠)≅ π‘Œπ‘‡π‘‡(𝑠)=𝑇𝑇(𝑠)π‘Œ [π‘₯(𝑑)]𝑛1≅𝑋𝑛1𝑇𝑇(𝑠)=𝑇𝑇(𝑠)𝑋𝑛1 [π‘₯(𝑑)]𝑛2≅𝑋𝑛2𝑇𝑇(𝑠)=𝑇𝑇(𝑠)𝑋𝑛2 π‘˜1(𝑠,𝑑)≅𝑇𝑇𝐾1𝑇(𝑑)

π‘˜2(𝑠,𝑑)≅𝑇𝑇𝐾2𝑇(𝑑)

Where 2π‘š-vectors 𝑋,𝑋′,𝑄,π‘Œ,𝑋𝑛1,𝑋𝑛2 and 2π‘š×2π‘š matrices 𝐾1 & 𝐾2 are TFs coefficients of

π‘₯(𝑠),π‘₯′(𝑠),π‘ž(𝑠),𝑦(𝑠),[π‘₯(𝑑)]𝑛1,[π‘₯(𝑑)]𝑛2,π‘˜1(𝑠,𝑑),π‘˜2(𝑠,𝑑) respectively.

Lemma

Let 2π‘š-vectors 𝑋 and 𝑋𝑛 be TFs coefficients of π‘₯(𝑠) and [π‘₯(𝑠)]𝑛 respectively. If

𝑋=(𝑋1𝑇𝑋2𝑇)𝑇=(𝑋10,𝑋11,…,𝑋1π‘š−1,𝑋20,𝑋21…𝑋2π‘š−1)𝑇, then

𝑋𝑛=(𝑋10,𝑋11,…,𝑋1π‘š−1,𝑋20,𝑋21…𝑋2π‘š−1)𝑇

Where 𝑛≥1 is a positive integer.

Proof:

Journal of the Maharaja Sayajirao University of Baroda ISSN: 0025-0422

74

Volume-56, No.1 (IV) 2022

When 𝑛=1, , follows at once from [π‘₯(𝑠)]𝑛=π‘₯(𝑠).

Suppose that (4.16), holds for 𝑛, we shall deduce it for 𝑛+1.

Since [π‘₯(𝑠)]𝑛+1≅(𝑋𝑇𝑇(𝑠)).(𝑋1𝑇𝑇(𝑠))

=𝑋𝑇𝑇(𝑠)𝑇𝑇(𝑠)𝑋𝑛 =𝑋𝑇𝑋𝑛̃𝑇(𝑠)

we obtain 𝑋𝑇𝑋𝑛̃=(𝑋10𝑛+1,𝑋11𝑛+1,……,𝑋1π‘š−1𝑛+1,𝑋20𝑛+1,𝑋21𝑛+1,……,𝑋2π‘š−1𝑛+1)𝑇

Equation holds for 𝑛+1.

Components of 𝑿𝒏 can be Computed in terms of Components of Unknown Vector 𝑿

We substitute (1.15) into (1.14), we have π‘Œπ‘‡π‘‡(𝑠)≅𝑋′𝑇𝑇(𝑠)+𝑄𝑇𝑇(𝑠)𝑇𝑇(𝑠)𝑋+πœ†1𝑇𝑇(𝑠)𝐾1∫𝑇(𝑑)𝑇𝑇(𝑑)𝑋𝑛1𝑑𝑑+𝑠0πœ†2𝑇𝑇(𝑠)𝐾2∫𝑇(𝑑)𝑇𝑇(𝑑)𝑋𝑛2𝑑𝑑𝑠0

it follows that π‘Œπ‘‡π‘‡(𝑠)≅𝑋′𝑇𝑇(𝑠)+(𝑄̃ 𝑇(𝑠))𝑇𝑋+πœ†1𝑇𝑇(𝑠)𝐾1𝑋̃𝑛1∫𝑇(𝑑)𝑑𝑑+πœ†2𝑇𝑇(𝑠)𝐾2𝐷𝑋𝑛2𝑠0

Using Operational matrix 𝑃,

π‘Œπ‘‡π‘‡(𝑠)≅𝑋′𝑇𝑇(𝑠)+𝑋𝑇𝑄̃ 𝑇(𝑠)+πœ†1𝑇𝑇(𝑠)𝐾1𝑋̃𝑛1𝑃𝑇(𝑠)+πœ†2(𝐾2𝐷𝑋𝑛2)𝑇𝑇(𝑠)

In which πœ†1𝐾1𝑋̃𝑛1𝑃 is a 2π‘š×2π‘š matrix

𝑇𝑇(𝑠)πœ†1𝐾1𝑋̃𝑛1𝑃𝑇(𝑠)≅𝑋𝑛1𝑇̂𝑇(𝑠)

Where 𝑋̃𝑛1 is an 2π‘š-vector with components equal to the diagonal entries of the matrix πœ†1𝐾1𝑋̃𝑛1𝑃.

π‘Œπ‘‡π‘‡(𝑠)≅𝑋′𝑇𝑇(𝑠)+𝑋𝑇𝑄̃ 𝑇(𝑠)+𝑋𝑛1𝑇̂𝑇(𝑠)+πœ†2(𝐾2𝐷𝑋𝑛2)𝑇𝑇(+𝑄̃ 𝑋+𝑋𝑛1̂+ πœ†2𝐾2𝐷𝑋𝑛2≅π‘Œ

Where 𝑄̃ is a diagonal matrix, so 𝑄̃ 𝑇=𝑄̃.

𝑋′ must be computed in terms of 𝑋. π‘₯(𝑠)−π‘₯(0)=∫π‘₯′(𝜏)π‘‘πœπ‘ 0

≅∫𝑋′𝑇𝑇(𝜏)π‘‘πœπ‘ 0 ≅𝑋′𝑇𝑃𝑇(𝑠) π‘₯(𝑠)≅𝑋′𝑇𝑃𝑇(𝑠)+𝑋0𝑇𝑇(𝑠)

Where 𝑋0 is the 2π‘š- vector of the form

Journal of the Maharaja Sayajirao University of Baroda ISSN: 0025-0422

75

Volume-56, No.1 (IV) 2022

𝑋0=[π‘₯0,π‘₯0,…π‘₯0]𝑇

Consequently,

𝑋≅𝑃𝑇𝑋′+𝑋0

Now, combining (1.19) and (1.20) and replacing ≅ with =, it follows that

(𝐼+𝑃𝑇𝑄̃)𝑋+𝑃𝑇𝑋𝑛1̂+πœ†2𝑃𝑇𝐾2𝐷𝑋𝑛2≅π‘ƒπ‘‡π‘Œ+𝑋0

Equation is a nonlinear system of 2π‘š algebraic equations for the 2π‘š unknowns 𝑋10,𝑋11,…𝑋1π‘š−1,𝑋20,𝑋21,…𝑋2π‘š−1.

Components of 𝑋𝑇=(𝑋1𝑇 𝑋2𝑇) can be obtainbed by an iterative method.

Hence, an approximate solution π‘₯(𝑠)≅𝑋𝑇𝑇(𝑠)

π‘₯(𝑠)≅𝑋1𝑇𝑇1(𝑠)+𝑋2𝑇𝑇2(𝑠)

Can be computed for equation (1.14) without using any projection method.

Conclusion

We investigated the numerical solution of nonlinear and singularly perturbed for Vulture integro-differential equations. Also, numerical solution of high-order and fractional nonlinear Volterra-Fredholm integro-differential equations and its error analysis are discussed. Finally, some numerical results of integro-differential equations are presented to illustrated the efficiency and accuracy of the proposed methods.

Reference

1. G.Yu.Mekhtieva, V.R.Ibrahimov “On one numerical method for solution of integro-differential equations”, Vestnik BSU, ser.phys.-math. Scien., 2003, No3, pp.21-27.(Russian).

2. H.Brunner, “Collection method for Volterra integral and Relation Functional Equations”, Cambridge University Press, Cambridge,2004.

3. P.Darania, E.Abadian, “A method for the numerical solution of the integro-differential equations”, Appl.Math. and Comput. 188(2007) 657-668.

4. B. Ahmad, S. Siva Sundaram, Existence Results For Nonlinear Impulsive Hybrid Boundary Value Problems involving Fractional Differential Equations, Non linear Analysis, hybrid System, 3(2009).

5. Z.B. Bai, H.S. Lii, Positive Solutions Of Boundary Valuve problems Of Nonlinear Fractional Differential Equations, J. Math Anal.Appl.311(2005)\

6. v. Lakshmikantham, D.D. Bainov And P.S. Simeonov, Theory Of Impulsive Differential Equations, modern Applied Mathematics (1989).

7. K.S. Miller, B.Ross, An Introduction To The Fractional Caculus And Fractional Differential Equations, Wiley, New york, 1993