JAMPJournal of Applied Mathematics and Physics2327-4352Scientific Research Publishing10.4236/jamp.2020.83037JAMP-98812ArticlesPhysics&Mathematics The Numerical Solutions of Systems of Nonlinear Integral Equations with the Spline Functions XiaoyanLiu1*YurongPan2University of La Verne, La Verne, USABengbu University, Bengbu, China2602202008034704808, February 20209, March 2020 12, March 2020© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

The main goal of this work is to develop an effective technique for solving nonlinear systems of Volterra integral equations. The main tools are the cardinal spline functions on small compact supports. We solve a system of algebra equations to approximate the solution of the system of integral equations. Since the matrix for the algebraic system is nearly triangular, It is relatively painless to solve for the unknowns and an approximation of the original solution with high precision is accomplished. In order to enhance the accuracy, several cardinal splines are employed in the paper. Our schemes were compared with other techniques proposed in recent papers and the advantage of our method was exhibited with several numerical examples.

System of Integral Equations Nonlinear Integral Equations Numerical Solutions Spline Functions
1. Introduction

Integral equations appear in many fields, including dynamic systems, mathematical applications in economics, communication theory, optimization and optimal control systems, biology and population growth, continuum and quantum mechanics, kinetic theory of gases, electricity and magnetism, potential theory, geophysics, etc. Many differential equations with boundary-value can be reformulated as integral equations. There are also some problems that can be expressed only in terms of integral equations. Abundant papers have appeared on solving integral equations, for example, Polyanin summarized different solutions of integral equations in  and   published in 2013 and 2016. In   and , we discussed numerical methods using cardinal splines in solving systems of linear integral equations. In this paper we are going to explore the applications of cardinal splines in solving nonlinear systems of integral equations.

We are interested in the systems of Volterra integral equations of the second kind

y ( x ) = F ( y T ( x ) ) + g ( x ) + ∫ a x   K ( x , t , y T ( t ) ) d t , (1.1)

where the kernel K ( x , t , y T ) = [ K 1 ( x , t , y T ) , K 2 ( x , t , y T ) , ⋯ , K m ( x , t , y T ) ] T and g ( x ) = [ g 1 ( x ) , g 2 ( x ) , ⋯ , g m ( x ) ] T are known functions, and y T ( x ) = [ y i ( x ) ] i = 1 m is to be determined; F ( y T ( x ) ) = [ F 1 ( y T ( x ) ) , F 2 ( y T ( x ) ) , ⋯ , F m ( y T ( x ) ) ] are known functions in y T ( x ) .

This paper is divided into six sections. In Section 2 and 3, two univariate cardinal continuous splines on small compact supports are presented. In Section 4, the applications of cardinal splines on solving integral equations are explored. The unknown functions are expressed as linear combinations of horizontal translations of a cardinal spline function. Then a system of equations on the coefficients is deducted. We can solve the system and a good approximation of the original solution is obtained. The sufficient condition for the existence of the inverse matrix is discussed and the convergence is investigated. In Section 5, the numerical examples are given. The non-linear system on unknowns is solved and an accurate approximation of the original solution is obtained in each case. Section 6 contains the concluding remarks.

2. Cardinal Splines with Small Compact Supports

Since the paper  by Schoenberg published in 1946, spline functions have been studied by many scholars. Spline functions have excellent properties and applications are endless (for example, cf. ). The spline functions on uniform partitions are simple to construct and easy to apply, and are sufficient for a variety of applications.

The starting point is frequently the zero degree polynomial B-spline, with the integral iteration formula

B 0 ( x ) = { 1, − 1 2 < x < 1 2 0, elsewhere , (2.1)

B n ( x ) = ∫ − 1 2 1 2   B n − 1 ( x + t ) d t ,   n = 1,2,3, ⋯ (2.2)

we could construct higher order polynomial spline functions with higher degree of smoothness. More specifically, B 1 ( x ) has the expression

B 1 ( x ) = 1 2 | x − 1 | + 1 2 | x + 1 | − | x | . (2.3)

B n ( x ) are called one dimensional B-splines, which are polynomial splines and have small supports ( − n + 1 2 , n + 1 2 ) , i.e. B n ( x ) = 0 for x > n + 1 2 or x < − n + 1 2 , and excellent traits (cf. ). In my previous papers  and , low degree orthonormal

spline and cardinal splines functions with small compact supports were applied in solving the second kind of Volterra integral equations. In this paper we use the notation ‖ u ‖ A = m a x x ∈ A | u ( x ) | .

Let B 1, h ( x ) = B 1 ( x h ) . It is proved that

∑ k = − ∞ ∞ ( a + b k h ) B 1, h ( x − k h ) = ( a + b x ) .

Notice that this particular B-spline is also a cardinal spline, therefore it is straightforward to apply it in interpolations. As far as the convergence rate of interpolation is concerned, we have the following proposition (cf.   and ).

Proposition 1. Given that f ( x ) ∈ C [ a , b ] , f ′ ( x ) exists and is bounded in [ a , b ] . Let n be an integer, h = b − a n , let x i = a + i h , f i = f ( x i ) , i = 0 , 1 , 2 , ⋅ ⋅ ⋅ , n ,

Ω f ( x ) = ∑ k = 0 n     f k B 1 , h ( x − k h ) ;

then

‖ Ω f ( x ) − f ( x ) ‖ [ a , b ] ≤ 3 ‖ f ′ ( x ) ‖ [ a , b ] h 2 .

If f ( x ) ∈ C ( − ∞ , ∞ ) and f ′ ( x ) exists and is bounded, let h be a real number, let f i = f ( i h ) , i = 0 , 1 , 2 , ⋅ ⋅ ⋅ ,

Ω f ( x ) = ∑ k = − ∞ ∞ f k B 1 , h ( x − k h ) ,

then

‖ Ω f ( x ) − f ( x ) ‖ ( − ∞ , ∞ ) ≤ 3 ‖ f ′ ( x ) ‖ ( − ∞ , ∞ ) h 2 . (2.4)

3. A Univariate C<sup>2</sup> Cardinal Spline

By cardinal conditions (cf. ), we mean, let L ( x ) be a function, { x i } , i = 0 , ± 1 , ± 2 , ± 3 , ⋯ be interpolation points, then

L ( x i ) = { 1 , i = 0 0 , i ≠ 0 ,     i = 0 , ± 1 , 2 , ± 3 , ± 4 , ⋅ ⋅ ⋅ . (3.1)

The cardinal spline that was originally given in  is based on B 3 ( x ) from (1.1) using the similar process as in Section 2. Let

L 3 ( x ) = 6 B 3 ( x ) − 60 ∫ − 1 2 1 2     t 2 B 2 ( x + t ) d t = ∫ − 1 2 1 2 ( 6 − 60 t 2 ) B 2 ( x + t ) d t . (3.2)

Then L 3 ( x h ) satisfies the above cardinal condition when x i = i h . Notice that by the construction, L 3 ( x ) ∈ C 2 ( − ∞ , ∞ ) . L 3 ( x ) = 0 for x ∈ ( − ∞ , − 2 ) ∪ ( 2, ∞ ) . L 3 ( x ) is a polynomial of degree ≤ 5 in each subinterval ( − 2, − 1 ) , ( − 1,0 ) , ( 0,1 ) , ( 1,2 ) of its support. Furthermore, from direct calculation we deduct the following two propositions (cf. ).

Proposition 2. Let L 3 ( x ) be the cardinal spline constructed above, then

∑ k = − ∞ ∞ L 3 ( x h − k ) ( α ( i h ) 2 + β ( i h ) + γ ) = α x 2 + β x + γ ,

‖ ∑ k = − ∞ ∞ L 3 ( x h − k ) ( i h ) 3 − x 3 ‖ ≤ 1 6 h 3 ,

where α , β , γ are any complex numbers.

Proposition 3. If f ( x ) ∈ C 3 [ a , b ] , let h = b − a n , n be an integer, let f i = f ( a + i h ) , i = 0 , 1 , 2 , ⋅ ⋅ ⋅ ,

Ω f ( x ) = ( 3 f 0 − 3 f 1 + f 2 ) L 3 ( x h + 1 ) + ∑ k = 0 n     f i L 3 ( x h − k )     + ( 3 f n − 3 f n − 1 + f n − 2 ) L 3 ( x h − n − 1 ) ;

then

‖ Ω f ( x ) − f ( x ) ‖ [ a , b ] ≤ 7 ‖ f ‴ ( x ) ‖ [ a , b ] h 3 .

If f ( x ) ∈ C 3 ( − ∞ , ∞ ) and be bounded, let h be a real number, let f i = f ( i h ) , i = 0 , 1 , 2 , ⋅ ⋅ ⋅ ,

Ω f ( x ) = ∑ k = − ∞ ∞ f i L 3 ( x h − k ) ,

then

‖ Ω f ( x ) − f ( x ) ‖ ( − ∞ , ∞ ) ≤ 7 ‖ f ‴ ( x ) ‖ ( − ∞ , ∞ ) h 3 .

4. Numerical Methods Solving Systems of Integral Equations

Method 1-V for solving the system of nonlinear Volterra integral equations

As for the Volterra integral Equations (1.1) we solve it in an interval [ a , b ] . Again we let h = b − a n , x i = a + i h , i = 0,1, ⋯ , n . Furthermore, by pluging in y s ( x ) = ∑ k = 0 n     c s , k B 1, h ( x − x k ) , K s ( x , t , y T ) = ∑ i = 0 n ∑ j = 0 n   K s ( x i , x j , c j ) B 1, h ( x − x i ) B 1, h ( t − x j ) , g s ( x ) = ∑ k = 0 n     g s ( x k ) B 1 , h ( x − x k ) , s = 1 , 2 , ⋯ , m , we get

∑ k = 0 n     F s ( c k ) B 1, h ( x − x k ) − ∑ i = 0 n ∑ j = 0 n   B 1, h ( x − x i ) ∫ a x K ( x i , x j , c j ) B 1, h ( t − x j ) d t = ∑ k = 0 n     g s ( x k ) B 1 , h ( x − x k ) ,     s = 1 , 2 , ⋯ , m

Let x = x l , we arrive at

F s ( c l ) − ∑ j = 0 n   K s ( x l , x j , c j ) ∫ a x l B 1, h ( t − x j ) d t = g s ( x l ) , (4.1)

l = 0 , 1 , 2 , 3 , 4 , ⋯ , n ; s = 1 , 2 , ⋯ , m (4.2)

which is a simple system of ( n + 1 ) m nonlinear equations of unknowns { c s , 0 , c s , 1 , ⋯ , c s , n } s = 1 , 2 , ⋯ , m . Notice that this is a nearly triangular system and it is solvable (the solution may not be unique because it is not linear):

F s ( c 0 ) = g s ( x 0 ) ,

F s ( c 1 ) − h 2 ∑ j = 0 1   K s ( x 1 , x j , c j ) = g s ( x 1 ) ,

F s ( c l ) − h 2 ( K s ( x l , x 0 , c 0 ) + 2 ∑ j = 1 l − 1   K s ( x l , x j , c j ) + K s ( x l , x l , c l ) ) = g s (xl)

l = 2,3,4, ⋯ , n ,

where c j = [ c 1, j , c 2, j , ⋯ , c m , j ] T , j = 0,1,2, ⋯ , n .

Proposition 4. Given the Equation (1) and that f s ( x ) , g s ( x ) ∈ C [ a , b ] , f ′ s ( x ) and g ′ s ( x ) exists and is bounded in [ a , b ] , K ( x , y , u ) ∈ C 1 ( [ a , b ] × [ a , b ] × [ α 1 , β 1 ] × ⋅ ⋅ ⋅ × [ α s , β s ] ) , F ( u ) ∈ C 1 ( [ α 1 , β 1 ] × ⋅ ⋅ ⋅ × [ α s , β s ] ) . Furthermore, K ( x , y , u ) and F ( u ) satisfies the condition:

| ∫ a b ( K ( x , t , f ( t ) ) − K ( x , t , u ( t ) ) ) d t | < L M m a x x ∈ [ a , b ] | f ( x ) − u ( x ) | .

‖ F ( f ( t ) ) − F ( u ) ‖ < χ m a x x ∈ [ a , b ] | f ( x ) − u ( x ) |

where L M + χ < 1 . Let n be an integer, h = b − a n , let x i = a + i h , f s i = f s ( x i ) , i = 0 , 1 , 2 , ⋅ ⋅ ⋅ , n , { c s ,0 , c s ,1 , ⋯ , c s , n } s = 1 m satisfy the nonlinear system (2.2)

f s ∗ ( x ) = ∑ k = 0 n     c s , i B 1 , h ( x − k h ) ;

then

‖ f s ∗ ( x ) − f s ( x ) ‖ [ a , b ] = O ( h 2 ) .

where [ f 1 ( x ) , f 2 ( x ) , ⋯ , f m ( x ) ] is the exact solution of Equation (1.1).

Method 2-V for solving the Volterra integral equation

To improve the approximation rate, we apply the spline function L 3, h ( x ) . Again we let h = b − a n , x i = a + i h , i = 0 , 1 , ⋯ , n . Furthermore, let L 3, h ( x ) = L 3 ( x h ) be the cardinal spline given in Section 3, and

f s ( x ) = ∑ k = − 1 n + 1   c s , k L 3 , h ( x − x k ) ,

K s ( x , t , y ) = ∑ i = 0 n ∑ j = 0 n   K s ( x i , x j , c j ) L 3, h ( x − x i ) L 3, h ( t − x j ) ,

g s ( x ) = ∑ k = 0 n     g s ( x k ) L 3 , h ( x − x k ) ,

s = 1,2, ⋯ , m , where

K ( x − 1 , x j , c j ) = 3 K ( x 0 , x j , c j ) − 3 K ( x 1 , x j , c j ) + K ( x 2 , x j , c j ) ,

K ( x i , x − 1 , c j ) = 3 K ( x i , x 0 , c j ) − 3 K ( x i , x 1 , c j ) + K ( x i , x 2 , c j ) ,

K ( x − 1 , x − 1 , c j ) = 3 ( 3 K ( x 0 , x 0 , c j ) − 3 K ( x 1 , x 0 , c j ) + K ( x 2 , x 0 , c j ) )     − 3 ( 3 K ( x 0 , x 1 , c j ) − 3 K ( x 1 , x 1 , c j ) + K ( x 2 , x 1 , c j ) )     + ( 3 K ( x 0 , x 2 , c j ) − 3 K ( x 1 , x 2 , c j ) + K ( x 2 , x 2 , c j ) ) ,

g s ( x − 1 ) = 3 g s ( x 0 ) − 3 g s ( x 1 ) + g s (x2)

K ( x n + 1 , x j , c j ) = 3 K ( x n , x j , c j ) − 3 K ( x n − 1 , x j , c j ) + K ( x n − 2 , x j , c j ) ,

K ( x i , x n + 1 , c j ) = 3 K ( x i , x n , c j ) − 3 K ( x i , x n − 1 , c j ) + K ( x i , x 2 , c j ) ,

K ( x n + 1 , x n + 1 , c j ) = 3 ( 3 K ( x n , x n , c j ) − 3 K ( x n − 1 , x n , c j ) + K ( x n − 2 , x n , c j ) )     − 3 ( 3 K ( x n , x n − 1 , c j ) − 3 K ( x n − 1 , x n − 1 , c j ) + K ( x n − 2 , x n − 1 , c j ) )     + ( 3 K ( x n , x n − 2 , c j ) − 3 K ( x n − 1 , x n − 2 , c j ) + K ( x n − 2 , x n − 2 , c j ) ) ,

g s ( x n + 1 ) = 3 g s ( x n ) − 3 g s ( x n − 1 ) + g s ( x n − 2 )

j = 0,1, ⋯ , n ,

where c j = [ c 1, j , c 2, j , ⋯ , c m , j ] T , j = 0,1,2, ⋯ , n . Let { c s ,0 , c s ,1 , ⋯ , c s , n } s = 1,2, ⋯ , m be unknown coefficients to be determined and c s , − 1 = 3 c s ,0 − 3 c s ,1 + c s ,2 , c s , n + 1 = 3 c s , n − 3 c s , n − 1 + c s , n − 2 . Plug into the integral Equation (1.1), then we have

∑ k = − 1 n + 1     c s , k L 3, h ( x − x k ) − ∫ a x ∑ i = − 1 n + 1 ∑ j = − 1 n + 1   K s ( x i , x j , c j ) L 3, h ( x − x i ) L 3, h ( t − x j ) d t = ∑ k = − 1 n + 1     g s ( x k ) L 3, h ( x − x k ) ,       x ∈ [ a , b ] .

Let x = x l , we arrive at ( l = 0,1,2,3,4, ⋯ , n )

c s , l − ∑ j = − 1 n + 1   K s ( x l , x j , c j ) ∫ a x l L 3 , h ( x − x i ) L 3 , h ( t − x j ) d t = g s ( x l ) , (4.3)

which is still a relatively simple system of equations. For the convergency rate of solution of the Volterra integral Equations (1.1), we have a similar result.

Proposition 5. Given that f s ( x ) , g s ( x ) ∈ C 3 [ a , b ] , f s ( 4 ) ( x ) and g s ( 4 ) ( x ) exists and is bounded in [ a , b ] , K ( x , y , u ) ∈ C 3 ( [ a , b ] × [ a , b ] × [ α 1 , β 1 ] × ⋅ ⋅ ⋅ × [ α s , β s ] ) , F ( u ) ∈ C 3 ( [ α 1 , β 1 ] × ⋅ ⋅ ⋅ × [ α s , β s ] ) . Furthermore, K ( x , y , u ) satisfies the condition:

‖ ∫ a b ( K ( x , t , f ( t ) ) − K ( x , t , u ( t ) ) ) d t ‖ < L M m a x x ∈ [ a , b ] ‖ f ( x ) − u ( x ) ‖ ,

‖ F ( f ( t ) ) − F ( u ) ‖ < χ m a x x ∈ [ a , b ] | f ( x ) − u ( x ) |

where L M + χ < 1 . Let n be an integer, h = b − a n , let x i = a + i h , f s i = f s ( x i ) , i = 0 , 1 , 2 , ⋅ ⋅ ⋅ , n ; { c s ,0 , c s ,1 , ⋯ , c s , n } s = 1 m satisfy the system (4.2)

f s ∗ ( x ) = ∑ k = − 1 n + 1   c s , k L 3 , h ( x − x k ) ,

then

‖ f s ∗ ( x ) − f s ( x ) ‖ [ a , b ] = O ( h 3 ) .

where [ f 1 ( x ) , f 2 ( x ) , ⋯ , f m ( x ) ] is the exact solution of Equation (1.1).

5. Numerical Examples

Example 1. Given the system of integral equations

u ( x ) + v 3 ( x ) − 2 x ∫ 0 x ( u ( t ) + v 3 ( t ) ) d t = a (x)

v ( x ) − u 3 ( x ) − 1 5 ∫ 0 x ( u 3 ( t ) − v ( t ) ) d t = b (x)

where a ( x ) = − 1 / 2 x 5 − 2 / 3 x 4 + x 3 + x 2 , b ( x ) = x + 1 / 10 x 2 − x 6 − 1 / 35 x 7 , and [ u ( x ) , v ( x ) ] are unknown functions.

Let h = 0.1 , x k = k h , k = 0 , 1 , 2 , ⋯ , 10 , u ( x ) = ∑ k = 0 10     c k B 1 , h ( x − x k ) , v ( x ) = ∑ k = 0 10     d k B 1 , h ( x − x k ) , a ( x ) = ∑ k = 0 10     a ( x k ) B 1 , h ( x − x k ) , b ( x ) = ∑ k = 0 10     b ( x k ) B 1 , h ( x − x k ) ,

K s ( x , t , u ( t ) , v ( t ) ) = ∑ i = 0 10 ∑ j = 0 10     K s ( x i , x j , c j , d j ) B 1 , h ( x − x i ) B 1 , h ( t − x j ) ,

we get

∑ k = 0 n     F s ( c k ) B 1, h ( x − x k ) − ∑ i = 0 n ∑ j = 0 n     B 1, h ( x − x i ) ∫ a x K ( x i , x j , c j ) B 1, h ( t − x j ) d t = ∑ k = 0 n     g s ( x k ) B 1 , h ( x − x k ) ,     s = 1 , 2 , ⋯ , m

where c j = [ c 1, j , c 2, j ] = [ c j , d j ] , [ g 1 ( x ) , g 2 ( x ) ] = [ a ( x ) , b ( x ) ] . Let x = x l , we arrive at

[ c 0 = a ( 0 ) , d 0 = b ( 0 ) ]

[ c 1 + d 1 3 − h 2 ∑ j = 0 1     2 h ( c j + d j 3 ) − a ( h ) = 0 ; d 1 − c 1 3 − h 2 ∑ j = 0 1 1 5 ( − d j + c j 3 ) − b ( h ) = 0 ]

[ c l + d l 3 − h 2 ( 2 ( l h ) ( c 0 + d 0 3 ) + 2 ∑ j = 1 l − 1   2 ( l h ) ( c j + d j 3 ) + 2 ( l h ) ( c l + d l 3 ) ) − a ( l h ) = 0 ; d l − c l 3 − h 2 ( 1 5 ( − d 0 + c 0 3 ) + 2 ∑ j = 1 l − 1 1 5 ( − d j + c j 3 ) + 1 5 ( − d l + c l 3 ) ) − b ( l h ) = 0 ]

l = 2,3,4, ⋯ ,10.

Solve the above nonlinear system and we obtain

[ c 0 , c 1 , c 2 , c 3 , c 4 , c 5 , c 6 , c 7 , c 8 , c 9 , c 10 ] = [ 0 , 0.01004 , 0.04018 , 0.09046 , 0.16091 , 0.2515 ,       0.36195 , 0.4921 , 0.6419 , 0.81156 , 1.00125 ]

[ d 0 , d 1 , d 2 , d 3 , d 4 , d 5 , d 6 , d 7 , d 8 , d 9 , d 10 ] = [ 0 , 0.10000 , 0.20000 , 0.30001 , 0.40008 , 0.50031 ,         0.60084 , 0.70167 , 0.80264 , 0.90365 , 1.0047 ]

Error < 5 × 10 − 3 .

(In the paper , the error is 0.02579 for c 10 )

Method 2-V, let f ( x ) = u ( x ) + v 3 ( x ) , g ( x ) = v ( x ) − u 3 ( x ) . The original system becomes

f ( x ) − 2 x ∫ 0 x f ( t ) d t = a (x)

g ( x ) + 1 5 ∫ 0 x   g ( t ) d t = b (x)

Applying L3, we still Let h = 0.1 , n = 10 , x k = k h , k = − 1,0,1,2, ⋯ ,10,11 ,

K s ( x , t , u ( t ) ) = ∑ k = − 1 11 ∑ j = − 1 11 K s ( x k , x j , c j , d j ) L 3 , h ( x − x k ) ⋅ L 3 , h ( t − x j ) ,

a ( x ) = ∑ k = − 1 11     a ( x k ) L 3 , h ( x − x k ) , b ( x ) = ∑ k = − 1 11   b ( x k ) L 3 , h ( x − x k ) , f ( x ) = ∑ k = − 1 11     c k L 3 , h ( x − x k ) , g ( x ) = ∑ k = − 1 11     d k L 3 , h ( x − x k ) , where c − 1 = f ( 0 ) − h f ′ ( 0 ) + h 2 f ″ ( 0 ) / 2 − h 3 f ‴ ( 0 ) / 6 = 0.01 − 0.001 = 0.009 , (we use the given integral equation to find that f ( 0 ) = a ( 0 ) = 0 , f ′ ( 0 ) = a ′ ( 0 ) = 0 , f ″ ( 0 ) = a ″ ( 0 ) = 2 , f ‴ ( 0 ) = a ‴ ( 0 ) = 6 ). Similarly, d − 1 = g ( 0 ) − h g ′ ( 0 ) + h 2 g ″ ( 0 ) / 2 − h 3 g ‴ ( 0 ) / 6 = − 0.1 . In addition, we let c n + 1 = 3 c n − 3 c n − 1 + c n − 2 and d n + 1 = 3 d n − 3 d n − 1 + d n − 2 . Plug into the integral equations, we obtain the nonlinear system:

∑ k = − 1 11   c k L 3 , h ( x − x k ) − ∫ 0 x ∑ k = − 1 11 ∑ j = − 1 11   2 x k c j L 3 , h ( x − x k ) L 3 , h ( t − x j ) d t = ∑ k = − 1 11   a ( x k ) L 3 , h ( x − x k ) ,

∑ k = − 1 11   d k L 3 , h ( x − x k ) + 1 5 ∫ 0 x ∑ k = − 1 11 ∑ j = − 1 11   d j L 3 , h ( x − x k ) L 3 , h ( t − x j ) d t = ∑ k = − 1 11   b ( x k ) L 3 , h ( x − x k )

the solution of the system is

[ c 0 , c 1 , c 2 , c 3 , c 4 , c 5 , c 6 , c 7 , c 8 , c 9 , c 10 ] = [ 0 , 0.011 , 0.048 , 0.117 , 0.224 , 0.375 , 0.576 ,         0.833 , 1.152000 , 1.5389996 , 2.000054 ] ;

[ d 0 , d 1 , d 2 , d 3 , d 4 , d 5 , d 6 , d 7 , d 8 , d 9 , d 10 ] = [ 0 , 0.099999 , 0.199936 , 0.299270 , 0.395902 , 0.484370 ,         0.553336 , 0.582339 , 0.537838 , 0.368533 , 0.0000508 ] .

Furthermore, we solve u ( i h ) + v 3 ( i h ) = c i , v ( i h ) − u 3 ( i h ) = d i , i = 0 , 1 , 2 , ⋯ , 10 . The results are

[ u ( 0 ) , u ( h ) , u ( 2 h ) , u ( 3 h ) , u ( 4 h ) , u ( 5 h ) , u ( 6 h ) , u ( 7 h ) , u ( 8 h ) , u ( 9 h ) , u ( 10 h ) ] = [ 0,0.010000,0.040000,0.090000,0.160001,0.250003,         0.360006,0.490009,0.640010,0.810010,0.999990 ] ,

[ v ( 0 ) , v ( h ) , v ( 2 h ) , v ( 3 h ) , v ( 4 h ) , v ( 5 h ) , v ( 6 h ) , v ( 7 h ) , v ( 8 h ) , v ( 9 h ) , v ( 10 h ) ] = [ 0.09999996,0.1999997,0.299999,0.3999978,0.499996,       0.5999946,0.699994,0.799995,0.899995,1.000021 ] .

The error < 2 × 10−5. Our method is much better.

Example 2 Given the system of integral equations

f ( x ) = a ( x ) − 1 2 g 4 ( x ) − 1 2 ∫ 0 x ( ( − 2 x t + 3 x 2 ) f ( t ) − t g 4 ( t ) ) d t

g ( x ) = b ( x ) − f 6 ( x ) − ∫ 0 x ( ( t − 2 x ) f 6 ( t ) + g ( t ) ) d t

where a ( x ) = x + 11 12 x 4 − 1 12 x 6 , b ( x ) = x + 1 2 x 2 + x 6 − 9 56 x 8 . Let h = 0.05 , x k = k h , k = 0,1,2, ⋯ ,10 . f ( x ) = ∑ k = 0 20     c k B 1 , h ( x − x k ) , g ( x ) = ∑ k = 0 20     d k B 1 , h ( x − x k ) , K s ( x , t , u ( t ) , v ( t ) ) = ∑ i = 0 20 ∑ j = 0 20     K s ( x i , x j , c j , d j ) B 1 , h ( x − x i ) B 1 , h ( t − x j ) , a ( x ) = ∑ k = 0 20     a ( x k ) B 1 , h ( x − x k ) , b ( x ) = ∑ k = 0 20     b ( x k ) B 1 , h ( x − x k ) , where [ f ( x ) , g ( x ) ] are unknown functions. we get

∑ k = 0 n   F s ( c k ) B 1, h ( x − x k ) − ∑ i = 0 n ∑ j = 0 n   B 1, h ( x − x i ) ∫ a x K ( x i , x j , c j ) B 1, h ( t − x j ) d t = ∑ k = 0 n     g s ( x k ) B 1, h ( x − x k ) ,     s = 1,2, ⋯ , m

where c j = [ c 1, j , c 2, j ] = [ c j , d j ] , [ g 1 ( x ) , g 2 ( x ) ] = [ a ( x ) , b ( x ) ] . Let x = x s , we arrive at

c 0 = a ( 0 ) , d 0 = b (0)

c 1 + ( 1 / 2 ) d 1 4 − a ( h ) + ( h / 2 ) ∑ j = 0 1 ( 1 / 2 ) ( ( − 2 ( 1 h ) ( j h ) + 3 ( 1 h ) 2 ) c j − ( j h ) d j 4 ) = 0

d 1 + c 1 6 − b ( h ) + ( h / 2 ) ∑ j = 0 1 ( d j + ( ( j h ) − 2 ( 1 h ) ) c j 6 )

c s + ( 1 / 2 ) d s 4 + ( h / 2 ) ( ( 1 / 2 ) ( ( − 2 ( s h ) ( 0 h ) + 3 ( s h ) 2 ) c 0 − ( 0 h ) d 0 4 )   + 2 ∑ j = 1 s − 1 ( 1 / 2 ) ( ( − 2 ( s h ) ( j h ) + 3 ( s h ) 2 ) c j − ( j h ) d j 4 )   + ( 1 / 2 ) ( ( − 2 ( s h ) ( s h ) + 3 ( s h ) 2 ) c s − ( s h ) d s 4 ) ) − a ( s h ) = 0 ,

d s + c s 6 + ( h / 2 ) ( ( d 0 + ( ( 0 h ) − 2 ( s h ) ) c 0 6 ) + 2 ∑ j = 1 s − 1 ( d j + ( ( j h ) − 2 ( s h ) ) c j 6 )     + ( d s + ( ( s h ) − 2 ( s h ) ) c s 6 ) ) − b ( s h ) = 0

s = 2,3,4, ⋯ ,20.

Solve the nonlinear system and we obtain:

[ c 0 , c 1 , c 2 , c 3 , c 4 , c 5 , c 6 , c 7 , c 8 , c 9 , c 10 , c 11 , c 12 , c 13 , c 14 , c 15 , c 16 , c 17 , c 18 , c 19 , c 20 ] = [ 0 , 0.5000 , 0.1000 , 0.1500 , 0.2000 , 0.2500 , 0.3000 , 0.3501 ,         0.4001 , 0.4501 , 0.5001 , 0.5502 , 0.6002 , 0.6503 , 0.7006 ,         0.7488 , 0.7996 , 0.8498 , 0.8999 , 0.94998 , 1.0000 ] ;

[ d 0 , d 1 , d 2 , d 3 , d 4 , d 5 , d 6 , d 7 , d 8 , d 9 , d 10 , d 11 , d 12 , d 13 , d 14 , d 15 , d 16 , d 17 , d 18 , d 19 , d 20 ] = [ 0 , 0.0500 , 0.1000 , 0.1500 , 0.2000 , 0.2500 , 0.3000 , 0.3500 ,         0.4000 , 0.4500 , 0.49999 , 0.54998 , 0.59995 , 0.6499 , 0.6996 ,         0.7518 , 0.8009 , 0.8507 , 0.9006 , 0.9506 , 1.0006 ] ,

error < 2 × 10 − 3

6. Conclusion

The proposed method is a simple and effective procedure for solving nonlinear Volterra integral equations of the second kind. The methods can be adapted easily to the Volterra integral equations of the first kind, which have the form

g s ( x ) = ∫ A   K ( x , t ) y ( t ) d t , where the upper limit of the integration is a variable.

The methods can also be extended to the Fredholm and Volterra integral equations of the first kind or the second kind, where the integral is on an infinite set. The higher degree cardinal splines could also be applied to non-linear integral equations; the resulting system of coefficients will be a little more complicated non-linear systems, which takes more time and effort to solve. Compared with the recent paper , our method is more effective.

Acknowledgements

The research is partially supported by the research funds of the University of La Verne, and by key natural science foundation of Anhui education commission (NO.KJ2017A568), natural science foundation of Bengbu University (NO.2018CXY045; 00009134).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Liu, X.Y. and Pan, Y.R. (2020) The Numerical Solutions of Systems of Nonlinear Integral Equations with the Spline Functions. Journal of Applied Mathematics and Physics, 8, 470-480. https://doi.org/10.4236/jamp.2020.83037

ReferencesPolyanin, A.D. (1998) Handbook of Integral Equations. CRC Press LLC, New York.https://doi.org/10.1201/9781420050066Singh, I. and Kumar, S. (2016) Haar Wavelet Method for Some Nonlinear Volterra Integral Equations of the First Kind. Journal of Computational and Applied Mathematics, 292, 541-552. https://doi.org/10.1016/j.cam.2015.07.022Aziz, I. and Siraj-ul-Islam (2013) New Algorithms for the Numerical Solution of Nonlinear Fredholm and Volterra Integral Equations Using Haar Wavelets. Journal of Computational and Applied Mathematics, 239, 333-345. https://doi.org/10.1016/j.cam.2012.08.031Liu, X. and Xie, J. (2014) Numerical Methods for Solving Systems of Fredholm Integral Equations with Cardinal Splines. AIP Conference Proceedings, 1637, 590.https://doi.org/10.1063/1.4904628Liu, X., Liu, Z. and Xie, J. (2015) Solving Systems of Volterra Integral Equations with Cardinal Splines. Journal of Applied Mathematics and Physics, 3, 1422-1430.https://doi.org/10.4236/jamp.2015.311170Liu, X., Xie, J., Liu, Z. and Huang, J. (2020) The Cardinal Spline Methods for the Numerical Solution of Nonlinear Integral Equations. Journal of Chemistry, 2020, Article ID: 3236813. https://doi.org/10.1155/2020/3236813Schoenberg, I.J. (1946) Contributions to the Problem of Approximation of Equidistant Data by Analytic Functions. Quarterly of Applied Mathematics, 4.Chui, C. K. (1988) Multivariate Splines. SIAM, Philadelphia, PA.https://doi.org/10.1137/1.9781611970173Liu, X. (2006) Univariate and Bivariate Orthornormal Splines and Cardinal Splines on the Compact Supports. Journal of Computational and Applied Mathematics, 195, 93-105. https://doi.org/10.1016/j.cam.2005.04.070Liu X. ,et al. (2007)Interpolation by Cardinal Exponential Splines The Journal of Information and Computational Science 4, 179-194.Liu X. ,et al. (2007)Interpolation by Cardinal Trigonometric Splines International Journal of Pure and Applied Mathematics 40, 115-122.Biazar, J., Babolian, E. and Islam, R. (2003) Solution of a System of Volterra Integral Equations of the First Kind by Adomian Method. Applied Mathematics and Computation, 139, 249-258. https://doi.org/10.1016/S0096-3003(02)00173-X