Skip to main content

Retracted Article: Jacobi spectral solution for integral algebraic equations of index 1

This article was retracted on 02 June 2015

Abstract

This paper is concerned with obtaining the approximate solution of a class of semi-explicit integral algebraic equations (IAEs) of index 1. A Jacobi collocation method is proposed for the IAEs of index 1. A rigorous analysis of the error bound in weighted L 2 -norm is also provided, which theoretically justifies the spectral rate of convergence while the kernels and the source functions are sufficiently smooth. The results of several numerical experiments are presented which support the theoretical results.

MSC:65R20, 45F15.

1 Introduction

Consider the following system of integral equations:

{ y ( t ) = f 1 ( t ) + ∫ 0 t K 11 ( t , s ) y ( s ) d s + ∫ 0 t K 12 ( t , s ) z ( s ) d s , 0 = f 2 ( t ) + ∫ 0 t K 21 ( t , s ) y ( s ) d s + ∫ 0 t K 22 ( t , s ) z ( s ) d s , t ∈ I = [ 0 , T ] ,
(1.1)

where K i j (⋅,⋅), i,j=1,2 are d i × d j matrices, f 1 , f 2 are given d 1 , d 2 dimensional vector functions, and (y(t),z(t)) is a solution to be determined. Here, we assume that the data functions f i , K i j (i,j=1,2) are sufficiently smooth such that f 2 (0)=0 and |det( K 22 (t,t))|≥ k 0 >0 for all t∈I. The existence and uniqueness results for the solution of the system (1.1) have been discussed in [1].

The system (1.1) is a particular case of the general form of the integral algebraic equations (IAEs)

A(t)X(t)=G(t)+ ∫ 0 t K ( t , s , X ( s ) ) ds

which has been introduced in [1], where det(A(t))=0, Rank(A)≥1. An initial investigation of these equations indicates that they have properties very similar to differential algebraic equations (DAEs). In analogy with the theory of DAEs (see, e.g., [2]), Kauthen [3] in 2000 has called the system (1.1) the semi-explicit IAEs of index 1.

Coupled systems of integral algebraic equations (IAEs) consisting of the first and second kind Volterra integral equations (VIEs) arise in many mathematical modeling processes, e.g., the controlled heat equation which represents a boundary reaction in diffusion of chemicals [4], the two dimensional biharmonic equation in a semi-infinite strip [5, 6], dynamic processes in chemical reactors [7], and the deformed Pohlmeyer equation [8]. Also a good source of information (including numerous additional references) on applications of IAEs system is the monograph by Brunner [1].

As far as we know, there have been a few works available in the literature which has considered the theory of IAEs system. The existence and uniqueness results of continuous solution to linear IAEs system have been investigated by Chistyakov [9]. Gear in [10] defined the index notion of IAEs system by considering the effect of perturbation of the equations on the solutions. He has also introduced the index reduction procedure for IAEs system in [10]. Bulatov and Chistyakov [11] gave the existence and uniqueness conditions of the solution for IAEs systems with convolution kernels and defined the index notion in analogy to Gear’s approach. Further details of their investigation may be found in [12, 13]. Kauthen in [3] analyzed the spline collocation method and its convergence properties for the semi-explicit IAEs system (1.1). Brunner [1] defined the index-1 tractability for the IAEs system (1.1) and he also investigated the existence of a unique solution for this type of systems. Recently, the authors in [4] have defined the index-2 tractability for a class of IAEs and presented the Jacobi collocation method including the matrix-vector multiplication representation of the equation and its convergence analysis. The authors in [14] proposed the Legendre collocation method for IAEs of index 1 and obtained the posteriori error estimation.

It is well known that the classical Jacobi polynomials have been used extensively in mathematical analysis and practical applications, and play an important role in the analysis and implementation of the spectral methods. The main purpose of this work is to use the Jacobi collocation method to numerically solve the IAEs (1.1). We will provide a posteriori error estimate in the weighted L 2 -norm which theoretically justifies the spectral rate of convergence. To do so, we use some well-known results of the approximation theory from [15–17] including the Jacobi polynomials, the Gronwall inequality and the Lebesgue constant for the Lagrange interpolation polynomials.

This paper is organized as follows. In Section 2, we carry out the Jacobi collocation approach for the IAEs system (1.1). A posteriori error estimation of the method in the weighted L 2 -norm as a main result of the paper is given in Section 3. Some numerical experiments are reported in Section 4 to verify the theoretical results obtained in Section 3. The final section contains conclusions and remarks.

2 The Jacobi collocation method

This section is devoted to applying the Jacobi collocation method to numerically solve the IAEs system (1.1). We first use some variable transformations to change the equation into a new system of integral equations defined on the standard interval [−1,1], so that the Jacobi orthogonal polynomial theory can be applied conveniently. For simplicity, we will consider the IAEs system (1.1) with d 1 = d 2 =1.

Let ω α , β (x)= ( 1 − x ) α ( 1 + x ) β (α,β>−1) be a weight function in the usual sense. As defined in [15, 18, 19], the set of Jacobi polynomials { J n α , β ( x ) } n = 0 ∞ forms a complete L ω α , β 2 (−1,1) orthogonal system, where L ω α , β 2 (−1,1) is the space of functions f:[−1,1]→R with ∥ f ∥ L ω α , β 2 ( − 1 , 1 ) 2 <∞, and

∥ f ∥ L ω α , β 2 ( − 1 , 1 ) 2 = 〈 f , f 〉 L ω α , β 2 ( − 1 , 1 ) = ∫ − 1 1 |f(x) | 2 ω α , β (x)dx.

For the sake of applying the theory of orthogonal polynomials, we use the change of variables:

η= 2 T s−1,−1≤η≤τ,τ= 2 T t−1,−1≤τ≤1,
(2.1)

to rewrite the IAEs system (1.1) as

{ y ˆ ( τ ) = f ˆ 1 ( τ ) + ∫ − 1 τ K ˆ 11 ( τ , η ) y ˆ ( η ) d η + ∫ − 1 τ K ˆ 12 ( τ , η ) z ˆ ( η ) d η , 0 = f ˆ 2 ( τ ) + ∫ − 1 τ K ˆ 21 ( τ , η ) y ˆ ( η ) d η + ∫ − 1 τ K ˆ 22 ( τ , η ) z ˆ ( η ) d η ,
(2.2)

where

f ˆ i ( τ ) = f i ( T 2 ( τ + 1 ) ) , K ˆ i j ( τ , η ) = T 2 K i j ( T 2 ( τ + 1 ) , T 2 ( η + 1 ) ) , i , j = 1 , 2 , y ˆ ( τ ) = y ( T 2 ( τ + 1 ) ) , z ˆ ( τ ) = z ( T 2 ( τ + 1 ) ) .

We consider the discrete expansion of K ˆ i j (τ,η) as follows:

I N ( K ˆ i j ( τ n , η ) ) = ∑ k = 0 N ( K ˆ i j ) k J k α , β (η)(i,j=1,2),
(2.3)

where I N is a projection to the finite dimensional space B N =span { J n α , β ( x ) } n = 0 N and J n α , β (x) is the Jacobi polynomial such that

( K ˆ i j ) k = 1 γ k ∑ l = 0 N ω l K ˆ i j ( τ n , τ l ) J k α , β ( τ l ),
(2.4)

where the quadrature points τ l are the Jacobi Gauss quadrature points, i.e., the zeros of J N + 1 α , β (x), the normalization constant γ k and the weights ω l are given in [[18], p.231].

In the Jacobi collocation method, we seek a solution of the form

I N ( y ˆ ( η ) ) = y ˆ N (η)= ∑ k = 0 N y ˆ k J k α , β (η), I N ( z ˆ ( η ) ) = z ˆ N (η)= ∑ k = 0 N z ˆ k J k α , β (η).
(2.5)

Inserting the discrete expansions (2.3) and (2.5) into (2.2), we obtain

{ y ˆ N ( τ ) = f ˆ 1 ( τ ) + ∑ k = 0 N ∑ l = k N ( c k l + c k l ′ ) V k l ( τ ) , 0 = f ˆ 2 ( τ ) + ∑ k = 0 N ∑ l = k N ( q k l + q k l ′ ) V k l ( τ ) ,
(2.6)

where V k l (τ)= ∫ − 1 τ J k α , β (η) J l α , β (η)dη, and

c k l = { y ˆ k ( K ˆ 11 ) k , k = l , y ˆ k ( K ˆ 11 ) l + y ˆ l ( K ˆ 11 ) k , k ≠ l , c k l ′ = { z ˆ k ( K ˆ 12 ) k , k = l , z ˆ k ( K ˆ 12 ) l + z ˆ l ( K ˆ 12 ) k , k ≠ l , q k l = { y ˆ k ( K ˆ 21 ) k , k = l , y ˆ k ( K ˆ 21 ) l + y ˆ l ( K ˆ 21 ) k , k ≠ l , q k l ′ = { z ˆ k ( K ˆ 22 ) k , k = l , z ˆ k ( K ˆ 22 ) l + z ˆ l ( K ˆ 22 ) k , k ≠ l .

By substituting the collocation points Ï„ n in the system (2.6) and employing the discrete representation (2.5), we obtain

{ y ˆ N ( τ n ) = f ˆ 1 ( τ n ) + ∑ k = 0 N ∑ l = k N ( c k l + c k l ′ ) V k l ( τ n ) , 0 = f ˆ 2 ( τ n ) + ∑ k = 0 N ∑ l = k N ( q k l + q k l ′ ) V k l ( τ n ) (n=0,1,…,N)
(2.7)

and we let

J = ( J i α , β ( τ j ) ) ( i , j = 0 , 1 , … , N ) , Y = ( y ˆ 0 , y ˆ 1 , … , y ˆ N ) T , Z = ( z ˆ 0 , z ˆ 1 , … , z ˆ N ) T , F 1 = ( f ˆ 1 ( τ 0 ) , f ˆ 1 ( τ 1 ) , … , f ˆ 1 ( τ N ) ) T , F 2 = ( f ˆ 2 ( τ 0 ) , f ˆ 2 ( τ 1 ) , … , f ˆ 2 ( τ N ) ) T , A = ( a l n ) , a l n = ∑ k = 0 N ( K ˆ 11 ) k V k l ( τ n ) , B = ( b l n ) , b l n = ∑ k = 0 N ( K ˆ 12 ) k V k l ( τ n ) ( l , n = 0 , 1 , … , N ) , C = ( c l n ) , c l n = ∑ k = 0 N ( K ˆ 21 ) k V k l ( τ n ) , D = ( d l n ) , d l n = ∑ k = 0 N ( K ˆ 22 ) k V k l ( τ n ) ( l , n = 0 , 1 , … , N ) ,

then system (2.7) can be rewritten as

( A − J B C D ) ( Y Z ) = ( − F 1 − F 2 ) .
(2.8)

Now, the unknown coefficients y ˆ k and z ˆ k , k=0,1,…,N are obtained by solving the linear system (2.8) and finally the approximate solutions y ˆ N (η) and z ˆ N (η) will be computed by substituting these coefficients into (2.5). Because the matrix of the coefficients of the linear system (2.8) is non-symmetric, when N is large, we can use some basic iterative methods to solve it, such as the generalized minimal residual (GMRES) method (see [[19], pp.48-60]), which is popular for solving non-symmetric linear systems.

3 Error estimation

In this section, we present a posteriori error estimate for the proposed scheme in the weighted L 2 -norm. Firstly, we recall some preliminaries and useful lemmas from [15].

Let P N be the space of all polynomials of degree not exceeding N, and Λ=(−1,1). From [15], the inverse inequality concerning differentiability of the algebraic polynomials on the interval can be expressed in terms of L p -norm.

Lemma 3.1 (see [15])

Let φ∈ P N , then for any integer r≥1, 2≤p≤∞, there exists a positive constant C independent of N such that

∥ φ ( r ) ∥ L ω α , β p ( Λ ) ≤C N 2 r ∥ φ ∥ L ω α , β p ( Λ ) .
(3.1)

From now on, to simplify the notation, we denote ∥ ⋅ ∥ L ω α , β 2 ( Λ ) by ∥⋅∥ and J k α , β (x) by J k (x). We also give some error bounds for the Jacobi system in terms of the Sobolev norms. The Sobolev norm and semi-norm of order m≥0, considered in this section, are given by

∥ u ( x ) ∥ H m ( Λ ) = ( ∑ k = 0 m ∥ u ( m ) ( x ) ∥ ) 1 2 ,
(3.2)
|u(x) | H m , N ( Λ ) = ( ∑ j = min { m , N + 1 } m ∥ u ( j ) ( x ) ∥ ) 1 2 .
(3.3)

Lemma 3.2 (see [15])

Let u∈ H m (Λ), I N u= ∑ k = 0 N u ˆ k J k is the truncated orthogonal Jacobi series of u. The truncation error u− I N u can be estimated as follows:

∥u− I N u∥≤C N − m |u | H m , N ( Λ ) ,∀u∈ H m (Λ),
(3.4)
∥ u − I N u ∥ H l ( Λ ) ≤C N 2 l − 1 2 − m |u | H m , N ( Λ ) ,∀u∈ H m (Λ),1≤l≤m.
(3.5)

The following main theorem reveals the convergence results of the presented scheme in the weighted L 2 -norm.

Theorem 3.1 Consider the system of integral algebraic equation (2.2) where y ˆ , z ˆ ∈ H m (Λ) and the data functions f ˆ i , K ˆ i j , i,j=1,2, are sufficiently smooth and | K ˆ 22 (τ,τ)|≥ k 0 >0 for all τ∈Λ. Let ( y ˆ N , z ˆ N ) be the Jacobi collocation approximation of ( y ˆ , z ˆ ) which is defined by (2.5). Then the following estimates hold:

∥ y ˆ N − y ˆ ∥ ≤ C N − m ( | f 1 | H m , N ( Λ ) + N − 1 ( | y ˆ | H m , N ( Λ ) + | z ˆ | H m , N ( Λ ) ) ) + C N − 1 ( ∥ y ˆ ∥ + ∥ z ˆ ∥ ) ∥ y ˆ N − y ˆ ∥ ≤ + { C N − 2 m log N ( Ω 11 | y ˆ | + Ω 12 | z ˆ | ) + C N − m log N ( Ω 11 ∥ y ˆ ∥ + Ω 12 ∥ z ˆ ∥ ) , − 1 < α , β ≤ − 1 2 , C N 1 2 + γ − 2 m ( Ω 11 | y ˆ | + Ω 12 | z ˆ | ) + C N 1 2 + γ − m ( Ω 11 ∥ y ˆ ∥ + Ω 12 ∥ z ˆ ∥ ) , otherwise ,
(3.6)
∥ z ˆ N − z ˆ ∥ ≤ C N − m ( | f 1 | H m , N ( Λ ) + ( N − 1 + N 1 2 ) ( | y ˆ | H m , N ( Λ ) + | z ˆ | H m , N ( Λ ) ) ) ∥ z ˆ N − z ˆ ∥ ≤ + C N − 1 ( ∥ y ˆ ∥ + ∥ z ˆ ∥ ) ∥ z ˆ N − z ˆ ∥ ≤ + { C N 2 − 2 m log N ( ( Ω 11 + Ω 21 ) | y ˆ | + ( Ω 12 + Ω 22 ) | z ˆ | ) + C N 2 − m log N ( ( Ω 11 + Ω 21 ) ∥ y ˆ ∥ + ( Ω 12 + Ω 22 ) ∥ z ˆ ∥ ) , − 1 < α , β ≤ − 1 2 , C N 5 2 + γ − 2 m ( ( Ω 11 + Ω 21 ) | y ˆ | + ( Ω 12 + Ω 22 ) | z ˆ | ) + C N 5 2 + γ − m ( ( Ω 11 + Ω 21 ) ∥ y ˆ ∥ + ( Ω 12 + Ω 22 ) ∥ z ˆ ∥ ) , otherwise
(3.7)

if N is sufficiently large, where Ω i j = max 0 ≤ n ≤ N | K i j ( τ n ,η) | H m , N ( Λ ) (i,j=1,2), γ=max(α,β).

Proof Similar to the idea in [14], we let

e ( τ ) = y ˆ N ( τ ) − y ˆ ( τ ) , ε ( τ ) = z ˆ N ( τ ) − z ˆ ( τ ) , S 1 ( τ n ) = ∑ k = 0 N ∑ l = k N c k l V k l ( τ n ) − ∫ − 1 τ n K ˆ 11 ( τ n , η ) y ˆ N ( η ) d η , S 2 ( τ n ) = ∑ k = 0 N ∑ l = k N c k l ′ V k l ( τ n ) − ∫ − 1 τ n K ˆ 12 ( τ n , η ) z ˆ N ( η ) d η , S 3 ( τ n ) = ∑ k = 0 N ∑ l = k N q k l V k l ( τ n ) − ∫ − 1 τ n K ˆ 21 ( τ n , η ) y ˆ N ( η ) d η , S 4 ( τ n ) = ∑ k = 0 N ∑ l = k N q k l ′ V k l ( τ n ) − ∫ − 1 τ n K ˆ 22 ( τ n , η ) z ˆ N ( η ) d η , F 1 = I N ( ∫ − 1 τ K ˆ 11 ( τ , η ) e ( η ) d η ) − ∫ − 1 τ K ˆ 11 ( τ , η ) e ( η ) d η , F 2 = I N ( ∫ − 1 τ K ˆ 12 ( τ , η ) ε ( η ) d η ) − ∫ − 1 τ K ˆ 12 ( τ , η ) ε ( η ) d η , F 3 = I N ( ∫ − 1 τ K ˆ 11 ( τ , η ) y ˆ ( η ) d η ) − ∫ − 1 τ K ˆ 11 ( τ , η ) y ˆ ( η ) d η , F 4 = I N ( ∫ − 1 τ K ˆ 12 ( τ , η ) z ˆ ( η ) d η ) − ∫ − 1 τ K ˆ 12 ( τ , η ) z ˆ ( η ) d η , F 5 = I N ( ∫ − 1 τ K ˆ 21 ( τ , η ) e ( η ) d η ) − ∫ − 1 τ K ˆ 21 ( τ , η ) e ( η ) d η , F 6 = I N ( ∫ − 1 τ K ˆ 22 ( τ , η ) ε ( η ) d η ) − ∫ − 1 τ K ˆ 22 ( τ , η ) ε ( η ) d η .

From systems (2.2) and (2.7), we obtain

e ( τ ) = ∫ − 1 τ K ˆ 11 ( τ , η ) e ( η ) d η + ∫ − 1 τ K ˆ 12 ( τ , η ) ε ( η ) d η + ( I N ( f ˆ 1 ( τ ) ) − f ˆ 1 ( τ ) ) e ( τ ) = + I N ( S 1 ( τ ) ) + I N ( S 2 ( τ ) ) + F 1 + F 2 + F 3 + F 4 ,
(3.8)
− K ˆ 21 ( τ , τ ) e ( τ ) − K ˆ 22 ( τ , τ ) ε ( τ ) = ∫ − 1 τ ∂ K ˆ 21 ( τ , η ) ∂ τ e ( η ) d η + ∫ − 1 τ ∂ K ˆ 22 ( τ , η ) ∂ τ ε ( η ) d η − K ˆ 21 ( τ , τ ) e ( τ ) − K ˆ 22 ( τ , τ ) ε ( τ ) = + I N ′ ( S 3 ( τ ) ) + I N ′ ( S 4 ( τ ) ) + F 5 ′ + F 6 ′ .
(3.9)

Equations (3.8) and (3.9) can be written in the compact matrix representation:

A(τ)E(τ)= ∫ − 1 τ K(τ,η)E(η)dη+B
(3.10)

with

A ( τ ) = [ 1 0 − K ˆ 21 ( τ , τ ) − K ˆ 22 ( τ , τ ) ] , K ( τ , η ) = [ K ˆ 11 ( τ , η ) K ˆ 12 ( τ , η ) ∂ K ˆ 21 ( τ , η ) ∂ τ ∂ K ˆ 21 ( τ , η ) ∂ τ ] , E ( τ ) = ( e ( τ ) ε ( τ ) )

and

B= ( ( I N ( f ˆ 1 ( τ ) ) − f ˆ 1 ( τ ) ) + I N ( S 1 ( τ ) ) + I N ( S 2 ( τ ) ) + F 1 + F 2 + F 3 + F 4 I N ′ ( S 3 ( τ ) ) + I N ′ ( S 4 ( τ ) ) + F 5 ′ + F 6 ′ ) .

Since | k 22 (τ,τ)|≥ k 0 >0, the inverse of the matrix A(τ) exists and

A − 1 (τ)= [ 1 0 − K ˆ 21 ( τ , τ ) K ˆ 22 ( τ , τ ) − 1 K ˆ 22 ( τ , τ ) ] .

Using the Gronwall inequality (see, e.g., [[20], Lemma 3.4]) on (3.10), we have

∥E∥≤C∥F∥.
(3.11)

Here F= A − 1 B.

It follows from (3.4) that

∥ ( I N ( f ˆ 1 ( τ ) ) − f ˆ 1 ( τ ) ) ∥ ≤ C N − m | f ˆ 1 | H m , N ( Λ ) , ∥ e ( τ ) ∥ ≤ C N − m | y ˆ | H m , N ( Λ ) , ∥ ε ( τ ) ∥ ≤ C N − m | z ˆ | H m , N ( Λ ) .
(3.12)

Using (3.4), (3.3) for m=1 and the Hardy inequality [[21], Lemma 3.8], we obtain

∥ F 1 ∥ ≤ C N − 1 ∥ ( ∫ − 1 τ K ˆ 11 ( τ , η ) e ( η ) d η ) τ ′ ∥ ≤ C N − 1 ∥ − K ˆ 11 ( τ , τ ) e ( τ ) + ∫ − 1 τ ∂ K ˆ 11 ( τ , η ) ∂ τ e ( η ) d η ∥ ≤ C N − 1 ∥ K ˆ 11 ( τ , τ ) + C ∥ ∥ e ( τ ) ∥ ≤ C N − 1 − m | y ˆ | H m , N ( Λ ) ,
(3.13)

and consequently

∥ F 2 ∥≤C N − 1 − m | z ˆ | H m , N ( Λ ) ,
(3.14)
∥ F 3 ∥≤C N − 1 ∥ y ˆ ∥,
(3.15)
∥ F 4 ∥≤C N − 1 ∥ z ˆ ∥.
(3.16)

Also,

I N ( S 1 ( τ ) ) = ∑ n = 0 N S 1 ( τ n ) L n (τ),
(3.17)

where L n (Ï„) is the Lagrange interpolation polynomial based on the Gauss quadrature nodes (see [19]). Therefore, we have

∥ I N ( S 1 ( τ ) ) ∥ ≤ max 0 ≤ n ≤ N | S 1 ( τ n )| max − 1 ≤ τ ≤ 1 ∑ n = 0 N L n (τ).

Moreover, using the Cauchy-Schwarz inequality [15], we have

| S 1 ( τ n ) | ≤ | ∫ − 1 τ n ( I N ( K ˆ 11 ( τ n , η ) ) − K ˆ 11 ( τ n , η ) ) y ˆ N ( η ) d η | ≤ ∥ I N ( K ˆ 11 ( τ n , η ) ) − K ˆ 11 ( τ n , η ) ∥ ∥ y ˆ N ∥ .

By (3.4), we have

| S 1 ( τ n )|≤C N − m | K ˆ 11 ( τ n ,η) | H m , N ( Λ ) ( ∥ e ∥ + ∥ y ˆ ∥ ) .
(3.18)

Now, we will make use of the result of Chen and Tang [[21], Lemma 3.4], which gives the Lebesgue constant for the Lagrange interpolating polynomials with the nodes of the Jacobi polynomials. Actually, the following relation holds:

max − 1 ≤ τ ≤ 1 ∑ n = 0 N L n (τ)= { O ( log N ) , − 1 < α , β ≤ − 1 2 , O ( N γ + 1 2 ) , γ = max ( α , β ) , otherwise .
(3.19)

So we obtain

∥ I N ( S 1 ( τ ) ) ∥ ≤ { C N − m log N Ω 11 ( ∥ e ∥ + ∥ y ˆ ∥ ) , − 1 < α , β ≤ − 1 2 , C N 1 2 + γ − m Ω 11 ( ∥ e ∥ + ∥ y ˆ ∥ ) , otherwise .
(3.20)

Similarly,

∥ I N ( S 2 ( τ ) ) ∥ ≤ { C N − m log N Ω 12 ( ∥ ε ∥ + ∥ z ˆ ∥ ) , − 1 < α , β ≤ − 1 2 , C N 1 2 + γ − m Ω 12 ( ∥ ε ∥ + ∥ z ˆ ∥ ) , otherwise .
(3.21)

It then follows from (3.1), (3.20), and (3.21) that

∥ I N ′ ( S 3 ( τ ) ) ∥ ≤ C N 2 ∥ I N ( S 3 ( τ ) ) ∥ ≤ { C N − m log N Ω 21 ( ∥ e ∥ + ∥ y ˆ ∥ ) , − 1 < α , β ≤ − 1 2 , C N 1 2 + γ − m Ω 21 ( ∥ e ∥ + ∥ y ˆ ∥ ) , otherwise .
(3.22)

Similarly,

∥ I N ′ ( S 4 ( τ ) ) ∥ ≤ C N 2 ∥ I N ( S 4 ( τ ) ) ∥ ≤ { C N − m log N Ω 22 ( ∥ ε ∥ + ∥ z ˆ ∥ ) , − 1 < α , β ≤ − 1 2 , C N 1 2 + γ − m Ω 22 ( ∥ ε ∥ + ∥ z ˆ ∥ ) , otherwise .
(3.23)

Using (3.2) and letting l=1 in (3.5), we have

∥ F 5 ′ ∥ ≤ ∥ F 5 ∥ H l ( Λ ) ≤C N 3 2 − m | ∫ − 1 τ K ˆ 21 (τ,η)e(η)dη | H m , N ( Λ ) .

Applying (3.3) for m=1, we have

∥ F 5 ′ ∥ ≤C N 1 2 ∥ K ˆ 21 ( τ , τ ) e ( τ ) + ∫ − 1 τ ∂ K ˆ 21 ( τ , η ) ∂ τ e ( η ) d η ∥ ≤C N 1 2 − m | y ˆ | H m , N ( Λ ) .
(3.24)

Similarly,

∥ F 6 ′ ∥ ≤C N 1 2 − m | z ˆ | H m , N ( Λ ) .
(3.25)

Finally, combining the above estimates and (3.11), the desired error estimates (3.6) and (3.7) are obtained. □

4 Numerical experiments

In this section, we consider some numerical examples in order to illustrate the validity of the proposed Jacobi collocation method. These problems are solved using the Jacobi collocation method for α=1/2, β=1/3. All the computations were performed using software Matlab®. To examine the accuracy of the results, L ω α , β 2 (−1,1) errors are employed to verify the efficiency of the method. To find the numerical convergence order in the numerical examples, we assume that E N =C ( 1 N ) p , where E N denotes the L ω α , β 2 (−1,1) errors between the numerical and exact solution of the system and C denotes some positive constant. Thus, the convergence order p can be computed by p= log 2 ( E N / 2 ) log 2 ( E N ) .

Example 4.1 Consider the following linear system of IAEs of index 1:

A(t)X(t)=g(t)+ ∫ 0 t K(t,s)X(s)ds,t∈[0,1],

where

A(t)= ( 1 0 0 0 ) ,X(t)= ( y ( t ) z ( t ) ) ,K(t,s)= ( t 3 + s + 1 cos ( 3 s ) + 1 t + s + 2 sin ( 3 s ) + 2 ) ,

and g(t)= ( f 1 ( t ) , f 2 ( t ) ) T are chosen so that the exact solutions of this system are

y(t)=cost,z(t)=sin3t.

Let ( y ˆ N , z ˆ N ) and ( y ˆ , z ˆ ) be the approximate and the exact solution of the system, respectively. The L ω α , β 2 (−1,1)-norm of the errors and orders of convergence are given in Tables 1 and 2.

Table 1 L ω α , β 2 (−1,1) errors for Example 4.1
Table 2 Convergence orders for Example 4.1

Example 4.2 Consider the following linear system of IAEs of index 1:

A(t)X(t)=g(t)+ ∫ 0 t K(t,s)X(s)ds,t∈[0,1],

where

A(t)= ( 1 0 0 0 ) ,X(t)= ( y ( t ) z ( t ) ) ,K(t,s)= ( t 2 + s 2 + 2 s + t + 1 e 2 t + s s + t 2 + 2 ) ,

and g(t)= ( f 1 ( t ) , f 2 ( t ) ) T are chosen so that the exact solutions of this system are

y(t)= e 2 t ,z(t)= 1 1 + t 2 .

The computational results have been reported in Tables 3 and 4.

Table 3 L ω α , β 2 (−1,1) errors for Example 4.2
Table 4 Convergence orders for Example 4.2

It can be seen from Tables 1 and 3 that the errors decay exponentially and the rates of convergence for y ˆ are larger than those for z ˆ . Tables 2 and 4 also shows that the orders of convergence for y ˆ are higher than those for z ˆ . In fact, it is noted that the differences of the convergence orders for y ˆ and z ˆ are about 2 from Tables 2 and 4. Although the exact solutions of Examples 4.1 and 4.2 are infinitely differentiable functions, it can be implied, from Theorem 3.1, that a conservative estimate of the numerical convergence order for system (1.1) is only 1. However, from numerical experiments, the numerical convergence order is much higher than 1.

Example 4.3 Consider the following linear system of IAEs:

A(t)X(t)=g(t)+ ∫ 0 t K(t,s)X(s)ds,t∈[0,1],

where

A(t)= ( 1 0 0 0 1 0 0 0 0 ) ,X(t)= ( x 1 ( t ) x 2 ( t ) x 3 ( t ) ) ,K(t,s)= ( 3 − 2 s 2 − s 3 − s 2 − s 2 ( 2 − s ) − 1 s − 2 − 1 1 s + 2 s 2 − 4 2 − s ) ,

and g(t)= ( f 1 ( t ) , f 2 ( t ) , f 3 ( t ) ) T are chosen so that the exact solutions of this system are

x 1 (t)= x 2 (t)= e t , x 3 (t)= e t t − 2 .

Let u, v, w be the approximations of the exact solutions x 1 , x 2 , x 3 , respectively. The errors and orders of convergence for the proposed method with several values of N are reported in Tables 5 and 6. The results show that the methods are convergent with a good accuracy.

Table 5 L ω α , β 2 (−1,1) errors for Example 4.3
Table 6 Convergence orders for Example 4.3

5 Conclusions

This paper studies the Jacobi collocation method for the semi-explicit IAEs system of index 1. The scheme consists of finding an explicit expression for the integral terms of the equations associated with the Jacobi collocation method. A posteriori error estimation of the method in the weighted L 2 -norm was obtained. With the availability of this methodology, it will now be possible to investigate the approximate solution of other classes of IAEs systems. Although our convergence theory does not cover the nonlinear case, it contains some complications and restrictions for establishing a convergent result similar to Theorem 3.1 which will be the subject of our future work.

References

  1. Brunner H: Collocation Methods for Volterra Integral and Related Functional Equations. Cambridge University Press, Cambridge; 2004.

    Book  MATH  Google Scholar 

  2. Hairer E, Wanner G: Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic Problems. Springer, New York; 1996.

    Book  MATH  Google Scholar 

  3. Kauthen JP: The numerical solution of integral-algebraic equations of index 1 by polynomial spline collocation methods. Math. Comput. 2000, 236: 1503–1514.

    Article  MathSciNet  MATH  Google Scholar 

  4. Hadizadeh M, Ghoreishi F, Pishbin S: Jacobi spectral solution for integral-algebraic equations of index-2. Appl. Numer. Math. 2011, 61: 131–148. 10.1016/j.apnum.2010.08.009

    Article  MathSciNet  MATH  Google Scholar 

  5. Cannon JR: The One Dimensional Heat Equation. Cambridge University Press, Cambridge; 1984.

    Book  MATH  Google Scholar 

  6. Gomilko AM: A Dirichlet problem for the biharmonic equation in a semi-infinite strip. J. Eng. Math. 2003, 46: 253–268. 10.1023/A:1025065714786

    Article  MathSciNet  MATH  Google Scholar 

  7. Kafarov VV, Mayorga B, Dallos C: Mathematical method for analysis of dynamic processes in chemical reactors. Chem. Eng. Sci. 1999, 54: 4669–4678. 10.1016/S0009-2509(99)00277-8

    Article  Google Scholar 

  8. Zenchuk AI: Combination of inverse spectral transform method and method of characteristics: deformed Pohlmeyer equation. J. Nonlinear Math. Phys. 2008, 15: 437–448. 10.2991/jnmp.2008.15.s3.42

    Article  MathSciNet  Google Scholar 

  9. Chistyakov VF: On singular systems of ordinary differential equations and their analogues. In Lyapunov Functions and Their Applications. Nauka, Novosibirsk; 1986:231–240. (Russian)

    Google Scholar 

  10. Gear CW: Differential-algebraic equations, indices, and integral-algebraic equations. SIAM J. Numer. Anal. 1990, 27: 1527–1534. 10.1137/0727089

    Article  MathSciNet  MATH  Google Scholar 

  11. Bulatov, MV, Chistyakov, VF: The properties of differential-algebraic systems and their integral analogs. Department of Mathematics and Statistics, Memorial University of Newfoundlande (1997)

    Google Scholar 

  12. Bulatov MV, Chistyakov VF: Numerical methods of investigation and solution for integral algebraic equations. Sci-CADE09 2009, 25–29.

    Google Scholar 

  13. Bulatov MV: Regularization of singular system of Volterra integral equations. Comput. Math. Math. Phys. 2002, 42: 315–320.

    MathSciNet  Google Scholar 

  14. Pishbin S, Ghoreishi F, Hadizadeh M: A posteriori error estimation for the Legendre collocation method applied to integral-algebraic Volterra equations. Electron. Trans. Numer. Anal. 2011, 38: 327–346.

    MathSciNet  MATH  Google Scholar 

  15. Canuto C, Hussaini MY, Quarteroni A, Zang TA: Spectral Methods: Fundamentals in Single Domains. Springer, New York; 2006.

    MATH  Google Scholar 

  16. Hesthaven JS, Gottlieb S, Gottlieb D: Spectral Methods for Time-Dependent Problems. Cambridge University Press, Cambridge; 2007.

    Book  MATH  Google Scholar 

  17. Qu CK, Wong R: Szego’s conjecture on Lebesgue constants for Legendre series. Pac. J. Math. 1988, 135: 157–188. 10.2140/pjm.1988.135.157

    Article  MathSciNet  MATH  Google Scholar 

  18. Guo B, Wang L: Jacobi interpolation approximations and their applications to singular differential equations. Adv. Comput. Math. 2001, 14: 227–276. 10.1023/A:1016681018268

    Article  MathSciNet  MATH  Google Scholar 

  19. Shen J, Tang T: Spectral and High-Order Methods with Applications. Science Press, Beijing; 2006.

    MATH  Google Scholar 

  20. Tang T, Xu X, Cheng J: On spectral methods for Volterra integral equations and the convergence analysis. J. Comput. Math. 2008, 26: 825–837.

    MathSciNet  MATH  Google Scholar 

  21. Chen Y, Tang T: Convergence analysis of the Jacobi spectral collocation methods for Volterra integral equations with a weakly singular kernel. Math. Comput. 2010, 79: 147–167. 10.1090/S0025-5718-09-02269-8

    Article  MathSciNet  MATH  Google Scholar 

Download references

Acknowledgements

This work is supported by the National Natural Science Foundation of China (11101109, 11271102), the Natural Science Foundation of Hei-long-jiang Province of China (A201107) and SRF for ROCS, SEM. The authors would like to acknowledge two anonymous referees for their careful reading of the manuscript and constructive comments.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jingjun Zhao.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

An erratum to this article can be found at http://dx.doi.org/10.1186/s13660-015-0700-x.

This article has been retracted by Professor Ravi P Agarwal, Editor-in-Chief of Journal of Inequalities and Applications. Following publication of this article, it was brought to the attention of the editorial and publishing staff that this article has a substantial overlap with two articles by Hadizadeh, Ghoreishi and Pishbin published in 2011 in Applied Numerical Mathematics and in Electronic Transactions on Numerical Analysis.

This is a violation of publication ethics which, in accordance with the Committee on Publication Ethics guidelines, warrants a retraction of the article and a notice to this effect to be published in the journal.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhao, J., Wang, S. Retracted Article: Jacobi spectral solution for integral algebraic equations of index 1. J Inequal Appl 2014, 156 (2014). https://doi.org/10.1186/1029-242X-2014-156

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1029-242X-2014-156

Keywords