Approximation Technique for Solving Linear Volterra Integro-Differential Equations with Boundary Conditions

Alnair, Mohamed E. A.; Khidir, Ahmed A.

doi:https://doi.org/10.1155/2022/2217882

Abstract and Applied Analysis

On this page

Abstract Introduction Results and Discussion Conclusion Data Availability Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 2217882 | https://doi.org/10.1155/2022/2217882

Approximation Technique for Solving Linear Volterra Integro-Differential Equations with Boundary Conditions

Mohamed E. A. Alnair^1,2and Ahmed A. Khidir^1,3

Academic Editor: Yufeng Xu

Received01 Mar 2022

Revised20 Mar 2022

Accepted29 Mar 2022

Published12 Apr 2022

Abstract

This paper presents a new technique for solving linear Volterra integro-differential equations with boundary conditions. The method is based on the blending of the Chebyshev spectral methods. The application of the proposed method leads the Volterra integro-differential equation to a system of algebraic equations that are easy to solve. Some examples are introduced and the obtained results are compared with exact solution as well as the methods that reported in the literature to illustrate the effectiveness and accuracy of the method. The results demonstrate that there is congruence between the numerical and the exact results to a high order of accuracy. Tables were generated to verify the accuracy convergence of the method and error. Figures are presented to show the excellent agreement between the results of this study and the results from literature.

1. Introduction

Many Mathematical models of physical phenomena usually results in linear or nonlinear Volterra integral and integro-differential equations. These equations play a crucial role in many fields of sciences such as economics, chemistry, finance, biology, and engineering. In engineering, physics, and mathematics areas, solving a problem means solving a set of differential equations and integral or integro-differential equations. In general, finding analytic solutions of Volterra integro-differential equations is usually difficult so it is required to obtain an approximate solution. Therefore, they have been of great interest by several researchers and scientists of Volterra integro-differential equations. Several studies have been carried out and developed in the literature to find analytical or numerical solutions for Volterra integro-differential equations. Loh and Phang [1] used Genocchi polynomials to solve numerically a system of Volterra integro-differential equations which is done by approximating functions using Genocchi polynomials and derivatives using Genocchi polynomials operational matrix of integer-order derivative. Loh et. al [2] obtained the approximate solution for a new class of time-fractional partial integro-differential equation of the Caputo-Volterra type. Mirzaee et. al [3] used the moving least squares and spectral collocation method to estimate the solution of nonlinear stochastic Volterra integro-differential equations. A class of block boundary value methods has been constructed by Zhou and Stynes [4] for the solution of linear neutral Volterra integro-differential equations with weakly singular kernels. Adewumi et al. [5] presented new efficient numerical methods for solving Volterra integro-differential equations and a system of nonlinear delay integro-differential equations which arises in biology. In [6], the author discussed variational iteration method for the solutions of different linear and nonlinear Volterra integral equations. For some more recent and interesting results, we refer to [7–12] and the references cited therein.

Recently, Lotfi and Alipanah [13] presented a spectral element method to solve linear Volterra integro-differential equations; Alrabaiah et al. [14] used Laplace Adomian decomposition method to find the semianalytic solutions of nonlinear Volterra fractional integro-differential equations; a numerical method was obtained to find a solution for a singularly perturbed Volterra integro-differential equation by Nana et al. [15].

The main contribution of this work is to introduce an efficient numerical technique for solving the linear Volterra integro-differential equations together with boundary conditions. The technique is based on the Chebyshev spectral collocation method. The Chebyshev spectral collocation methods have been applied successfully in different fields of sciences and engineering because of their ability to give high accuracy of boundary value problems (see [16–20]. The application of the current method leads the Volterra integro-differential equation to a system of algebraic equations that are easy to solve when compared to a system of Volterra integro-differential equation. The main benefit of this method is that we used our proposed technique without any assumptions on large or small parameters. The investigation is mainly targeted to device this consistent method to integro-differential equations with boundary conditions. The main advantages of this approach over the standard same approaches are (i) the technique suggests a standard way of choosing the linear operator of the integro-differential equation whereas the other methods are choosing a linear operator to be simple in order to ensure that the integro-differential equation can be easily solved, and (ii) this algorithm transforms the integro-differential equation into a system of linear algebraic equations that is easier and faster to solve when compared to a system of integral equations.

This paper is organized as follows: in Section 2, we introduce a description of the proposed method. In Section 3, the method is implemented using some examples with known analytical solution. The results are discussed and investigated in Section 4. Finally, the conclusions are given in Section 5.

2. Description of the New Technique

Consider the following linear second order Volterra integro-differential equation of the second kind given by under the boundary conditions where is an unknown function to be determined; , and are known real valued functions defined on .

To illustrate the idea of the new algorithm, we assume that the kernel of Equation (1) can be separated into a product of two functions, namely, and ; consequently, Equation (1) can be obtained as follows:

Now, the integral term of Volterra integro-differential equation above is expressed as

One approach is to note that the integral equation (4) is an initial value problem expressed as

This differential equation is easy to solve using any method; here, we used the Chebyshev spectral collocation method; the functions and are approximated as a truncated series of Chebyshev polynomials given by the form [21–23]. where is the th Chebyshev polynomial and are the Chebyshev coefficients, and are the Gauss-Lobatto collocation points (see [21]) defined by where is the number of collocation points. The derivatives at the collocation points are represented as where is the order of differentiation and is the Chebyshev spectral differentiation matrix whose entries (see [22]) are given by

Here, and with .

Substituting the above assumptions in Equation (5) yields a system of algebraic equations expressed as the following matrix equation where , and . The solution of this system for is obtained by where we observe that the first row of Chebyshev differential matrix in equation above is replaced by the row , and we add the vector ; this is caused by imposing the condition into Equation (10). According to the integral Equation (4) and Equation (11), we introduce the following integral operator defined as where and

The operator is square matrix of size

2.1. The Linearity of the Operator

It is clear that the above operator is a linear operator because

Therefore, the integral operator is a linear operator.

2.2. Existence and Uniqueness of the Operator

In this subsection, we show that the operator exists and is unique.

Theorem 1. (i)Any square matrix is invertible (nonsingular) if and only if det(ii)If is an invertible matrix, then its inverse is exist and unique

Proof. See Larson and David [24].

According to the theorem above, the existence and uniqueness of the linear operator depends on the existence and uniqueness of inverse of matrix

The determinants of matrix above have been computed on the domain for . In Figure 1, we plotted the determinants for different values of varied . We observe that from the figure, all the computed determinants are not equal to zero, and it is noticed that all the determinants are greater than or equal to , i.e., the inverse of the operator exists and is unique.

Now, replace the derivative operator of the integro-differential Equation (1) by the Chebyshev spectral differentiation matrix defined in (9) and also replace the integral term of the integro-differential Equation (1) by the operator , which gives

Here, the superscript denotes matrix transpose, is the diagonal matrix of size that puts the vector of size on the main diagonal, and , , , , , and are vectors of sizes . The functions and are obtained from the kernel function which can be written as as we mentioned before.

Thus, Equation (16) can be written as where is a square matrix and is column vector of sizes defined by

The boundary conditions and are implemented into the system (17) by modifying the first and last rows of matrices and in such a way that the the system of linear Equation (17) takes the form

After modifying the matrix system (17) to incorporate boundary conditions, the final approximate solution of the linear second order Volterra integro-differential Equation (1) is obtained as where and are the the modified matrices of and , respectively.

The general idea underpinning the use of the proposed method is to convert the linear Volterra integro-differential equation by a system of linear algebraic equations that replace the derivative and integral parts in the integral equation by the differential matrix and the introduced integral operator, respectively. The obtained linear algebraic equations can easily be solved with the help of symbolic computation software such as Maple, Mathematica, MATLAB, and other symbolic computer packages.

3. Numerical Examples

In this section, some numerical examples are carried out to illustrate the effectiveness and accuracy of the proposed method. In addition, the numerical results are compared with exact solution, and all of them were performed on the computer using a program written in MATLAB v7.5.0 (R2007b). We also showed in figures and tables the approximate solutions, exact solutions, and the absolute error function at the selected points of the given interval.

Example 1. Consider the following linear second-order Volterra integro-differential Equation [25] together with the boundary conditions and . The exact solution is . The technique of the solution starts by rewriting the integral term in the above equation as

Applying the proposed algorithm on Equation (21) and according to the assumptions (8) and (12), one can write the integro-differential Equation (21) in a system of linear equations: where , , and is an identity matrix.

To implement the boundary conditions to the systems (23), we replace all the entries of the first and last rows of the matrix by zeros and then set , and . After modifying the matrix system (23) to incorporate the boundary conditions, the solution is obtained as where and are the modified matrices of and , respectively.

Example 2. Consider the following Volterra integro-differential equation [26]: subject to the boundary conditions . The exact solution is .

The integral term of integro-differential Equation (25) can be expressed as

According to the assumptions (8) and (12), we can apply the proposed algorithm on the integro-differential Equation (25); the equation is transformed into the following system of linear equations where and . After modifying the matrix system (27) to incorporate the boundary conditions of the integro-differential Equation (25), the solution is obtained as where and are the modified matrices of and , respectively.

Example 3. Consider the linear Volterra integro-differential equation given by [26, 27] with the initial conditions and the exact solution .

According to the assumptions (8) and (12), one can express the integral term of integro-differential Equation (29) as

Applying the proposed algorithm on the integro-differential Equation (29), the equation is transformed into the following system of linear equations where and .

Now, incorporate the initial condition into the system (31) to obtain the final solution of integro-differential Equation (25) as where and are the the modified matrices of and , respectively, after imposing the initial condition.

Example 4. Consider the following system of linear Volterra integral equations subject to the boundary conditions where . The exact solutions are and .

Applying the proposed algorithm, one can express the integral parts of integro-differential Equations (33) and (34) as

Applying the proposed algorithm on the system of integro-differential Equations (33) and (34), the equation is transformed into the following system of linear equations where is a square matrix while and are column vectors defined by where

In the above definitions, is a column vector whose interies are s. The boundary conditions (35) are imposed on Equation (38) by modifying the first and last rows of and in such a way that the modified matrices and take the form

After modifying the matrix system (38) to incorporate boundary conditions, the solution is obtained as where and are the the modified matrices of and , respectively.

4. Results and Discussion

In this section, we give the results obtained by applying the proposed algorithm on various linear Volterra integro-differential equations. Implementation of the numerical schemes was performed using personal computer of 2.5 GHz CPU speed including Matlab2017 package to perform the simulation results. The accuracy of the method is demonstrated by presenting infinity error between exact and approximate solutions computed as where and are the exact and numerical solutions, respectively. To check the validity and accuracy of the current algorithm, we made comparisons between the obtained results with the exact solutions and also with those other methods reported in the literature. The results of our solutions for all introduced examples in this report are showed and discussed in tables and graphs.

Table 1 shows absolute errors of Example 1 for , and at selected values of between and . Also, the time taken for the computation has been presented at the bottom of the table. The absolute errors for various values of number of grid points are shown in Figure 2. From Table 1 and Figure 2, clearly that the maximum absolute errors are generally very small and results obtained by the current algorithm are almost the same as the results of the exact solution, it is also noted that when increasing the number of grid points does not result in a significant improvement in the accuracy of the current approximation. It is clear that from Table 1, the present algorithm is computationally fast as accurate results are generated in a fraction of a second as it is shown in bottom of the table.

In Table 2, the exact errors of Example 2 have been calculated for various in to measure the extent of agreement between the exact and approximate solution. We also presented the time taken for the computation at the end of the table. From the table, it can be seen that the current algorithm provides us with the accurate approximate solution of Example 2. We remark that increasing the number of nodes (i.e., increasing ) does not result in a significant improvement in the accuracy of the approximation.

We also compare the relative error estimates with those obtained using the homotopy analysis method (HAM) reported by El-Ajou et al. [28] for in in Table 3. Full convergence upto decimal places accurate results reported by HAM is achieved for , and an increase in the values of results in an increase in the maximum error form, while high convergence is achieved to decimal places accurate results for the present method for all values of ; also, increasing in nodes () leads to decrease in the accuracy of HAM for Example 2 while in increasing does not result in a significant improvement in the accuracy using the current approximation. One can note that the present method gives much better numerical results compared to HAM. Figure 3 shows the errors between the exact and numerical solutions for Example 2 varied . Very small errors are achieved () for the values , i.e., that accuracy does not improved when for large values of .

In Table 4, the absolute error is obtained for Example 3 using current method, Tau method investigated by Hosseini and Shahmorad [26], and practical matrix method obtained by Yüzbaşı et al. [27]. It is seen from the table that the results obtained by the present method is better than those obtained by Tau and practical matrix methods; also, it is noticed from Table 4 that when increased, the error is increased for Tau and practical matrix methods while in the present method the error is approximately fixed for all values of in the domain. Figure 4 showed the error of the numerical results against the nodes of Example 3. It can be seen from this graph that when the number of nodes () increase, the error becomes smaller, and then, no effect when . In Table 5, we give the maximum errors between the exact and present results of Example 3; the results were computed in the domain . To give a sense of the computational efficiency of the method, the computational time to generate the results is also given. The accuracy is seen to improve with an increase in the number of collocation points . It is observed that accurate results with errors of order up to are obtained using very few collocation points. We remark, also, that this algorithm is computationally fast as accurate results are generated in a fraction of a second.

Table 6 and Figure 5 are showed the maximum absolute errors of the approximate solutions for Example 4 at selected values of and . The maximum absolute errors are very small for and and further decrease with an increase in the values of upto , and then, increasing does not result in a significant improvement in the accuracy of the approximate solution. We also observed that the current algorithm is computationally fast as accurate results are generated in a fraction of a second even for large values on . Figure 5 shows a comparison between the exact and approximate solutions for Example 4 of and . We observe that there is good agreement between the two sets of results, proving the efficiency of the proposed approach.

(a)

(b)

In Figure 6, the maximum absolute errors for Example 4 is presented. The figures display a variation of the errors at a various values of . It can be seen that, in the problem considered, the algorithm needs small values of to achieve good currency and to converge fully between the exact and approximate results. Also, we notice that there is no significant improvement in the accuracy when the number of nodes () is increasing.

(a)

(b)

A comparison between numerical results with the exact solution for Example 4 is shown in Figure 5. We can see that the two graphs overlap indicating the very good accuracy of the method. In Figure 6, convergence analysis graphs for Example 4 are presented. The figures present a variation of the error at various values of collocation points . It can be seen that, in almost all, the iteration scheme starts at to converge fully. It is remarkable to note accurate results with errors of order up to . It is worth remarking that the accuracy of the method is not dependent on the number of collocation points.

All the above observations are clear indication that this algorithm is a powerful method that is appropriate in solving linear Volterra integro-differential equations.

5. Conclusion

In this study, we have proposed a simple and easy algorithm for solving linear Volterra integro-differential equations with boundary conditions. This technique has been constructed by introducing a new integral operator together with Chebyshev pseudospectral method. The method allows us to choose the integral operators in the integro-differential equations in terms of the Chebyshev spectral differentiation matrix. The numerical solutions have been shown in tables and graphs and compared with the exact solutions and other methods in the literature. The accuracy of the obtained result is a point of interest; moreover, if we increase the number of nodes approximation, the numerical solution get closer to the exact solutions. The fast convergence and accuracy of the proposed algorithm are valid reasons for researcher to use this method for different problems of Volterra integro-differential equations arising in various areas of science.

The main conclusions emerging from this study may be summarized as follows: (1)The proposed method is highly accurate and efficient, converges rapidly, very easy to understand, and is sufficient to give good agreement with the exact solution(2)The suggested technique does not depend on the number of collocation points, and it requires a few numbers of collocation points to achieve high accuracy of results(3)No any iteration is required to achieve the currency while other methods required more iterations(4)The method is more flexible than other numerical methods in the literature such as HAM, Tau method, and practical matrix method

Finally, for the future study, we suggest that how is the nonseparable kernel dealt with and also how to apply this technique on nonlinear Volterra and Fredholm integral equations.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

J. R. Loh and C. Phang, “A new numerical scheme for solving system of Volterra integro-differential equation,” Alexandria Engineering Journal, vol. 57, no. 2, pp. 1117–1124, 2018.
View at: Publisher Site | Google Scholar
J. R. Loh, C. Phang, and K. G. Tay, “New method for solving fractional partial integro-differential equations by combination of Laplace transform and resolvent kernel method,” Chinese Journal of Physics, vol. 67, pp. 666–680, 2020.
View at: Publisher Site | Google Scholar
F. Mirzaee, E. Solhi, and S. Naserifar, “Approximate solution of stochastic Volterra integro-differential equations by using moving least squares scheme and spectral collocation method,” Applied Mathematics and Computation, vol. 410, article 126447, 2021.
View at: Publisher Site | Google Scholar
Y. Zhou and M. Stynes, “Block boundary value methods for solving linear neutral Volterra integro-differential equations with weakly singular kernels,” Journal of Computational and Applied Mathematics, vol. 401, article 113747, 2022.
View at: Publisher Site | Google Scholar
A. O. Adewumi, S. O. Akindeinde, and R. S. Lebelo, “Sumudu Lagrange-spectral methods for solving system of linear and nonlinear Volterra integro-differential equations,” Applied Numerical Mathematics, vol. 169, pp. 146–163, 2021.
View at: Publisher Site | Google Scholar
A. M. Wazwaz, “The variational iteration method for solving linear and nonlinear Volterra integral and integro-differential equations,” International Journal of Computer Mathematics, vol. 87, no. 5, pp. 1131–1141, 2010.
View at: Publisher Site | Google Scholar
H. Rezazadeh, M. S. Osman, M. Eslami et al., “Hyperbolic rational solutions to a variety of conformable fractional Boussinesq-like equations,” Nonlinear Engineering, vol. 8, no. 1, pp. 224–230, 2019.
View at: Publisher Site | Google Scholar
M. S. Osman, H. Rezazadeh, and M. Eslami, “Traveling wave solutions for (3+1) dimensional conformable fractional Zakharov-Kuznetsov equation with power law nonlinearity,” Nonlinear Engineering, vol. 8, no. 1, pp. 559–567, 2019.
View at: Publisher Site | Google Scholar
M. S. Osman, “New analytical study of water waves described by coupled fractional variant Boussinesq equation in fluid dynamics,” Pramana, vol. 93, no. 2, p. 26, 2019.
View at: Publisher Site | Google Scholar
M. S. Osman, D. Lu, and M. M. A. Khater, “A study of optical wave propagation in the nonautonomous Schrodinger-Hirota equation with power-law nonlinearity,” Results in Physics, vol. 13, article 102157, 2019.
View at: Google Scholar
M. S. Osman and A. M. Wazwaz, “A general bilinear form to generate different wave structures of solitons for a (3+1)-dimensional Boiti-Leon-Manna-Pempinelli equation,” Mathematical Methods in the Applied Sciences, vol. 42, no. 18, pp. 6277–6283, 2019.
View at: Publisher Site | Google Scholar
M. Higazy, S. Aggarwal, and T. A. Nofal, “Sawi decomposition method for Volterra integral equation with application,” Journal of Mathematics, vol. 2020, Article ID 6687134, 13 pages, 2020.
View at: Publisher Site | Google Scholar
M. Lotfi and A. Alipanah, “Legendre spectral element method for solving Volterra-integro differential equations,” Results in Applied Mathematics, vol. 7, pp. 100116–100116, 2020.
View at: Publisher Site | Google Scholar
H. Alrabaiah, M. Jamil, K. Shah, and R. A. Khan, “Existence theory and semi-analytical study of non-linear Volterra fractional integro-differential equations,” Alexandria Engineering Journal, vol. 59, no. 6, pp. 4677–4686, 2020.
View at: Publisher Site | Google Scholar
N. Nana, S. C. O. Noutchie, and R. Y. Massoukou, “A second order finite difference scheme for singularly perturbed Volterra integro-differential equation,” Alexandria Engineering Journal, vol. 59, no. 4, pp. 2441–2447, 2020.
View at: Google Scholar
A. A. Khidir and S. L. Alsharari, “Application of successive linearisation method on the boundary layer flow problem of heat and mass transfer with radiation effect,” International Journal of Analysis and Applications, vol. 19, no. 5, pp. 725–742, 2021.
View at: Google Scholar
A. A. Khidir, “A new numerical technique for solving Volterra integral equations using Chebyshev spectral method,” Mathematical Problems in Engineering, vol. 2021, Article ID 9230714, 11 pages, 2021.
View at: Publisher Site | Google Scholar
Y. Daoud, M. Abdalbagi, and A. A. Khidir, “On the solution of magneto-hydrodynamics three-dimensional flow due to a stretching sheet in a porous medium using the successive linearization method,” Chinese Journal of Physics, vol. 73, pp. 232–238, 2021.
View at: Publisher Site | Google Scholar
A. A. Khidir and A. F. Aljohani, “On successive linearization method for differential equations with nonlinear conditions,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 23, 2021.
View at: Publisher Site | Google Scholar
A. A. Khidir and W. I. Alfaifi, “Application of successive linearisation method on mixed convection from an exponentially stretching surface with magnetic field effect,” International Journal of Scientific Research in Mathematical and Statistical Sciences, vol. 8, no. 6, pp. 13–19, 2021.
View at: Google Scholar
C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamics, Springer-Verlag, Berlin, 1988.
View at: Publisher Site
W. S. Don and A. Solomonoff, “Accuracy and speed in computing the Chebyshev collocation derivative,” SIAM Journal on Scientific Computing, vol. 16, no. 6, pp. 1253–1268, 1995.
View at: Publisher Site | Google Scholar
L. N. Trefethen, Spectral Methods in MATLAB, SIAM, 2000.
View at: Publisher Site
R. Larson and D. C. Falve, Elementary Linear Algebra, Houghton mifflin harcourt publishing company Boston, New York, NY, USA, 6th edition, 2008.
S. Cravero, G. Pittaluga, and L. Sacripante, “An algorithm for solving linear Volterra integro-differential equations,” Numerical Algorithms, vol. 60, no. 1, pp. 101–114, 2012.
View at: Publisher Site | Google Scholar
S. M. Hosseini and S. Shahmorad, “Numerical solution of a class of integro-differential equations by the Tau method with an error estimation,” Applied Mathematics and Computation, vol. 136, no. 2-3, pp. 559–570, 2003.
View at: Publisher Site | Google Scholar
Ş. Yüzbaşı, N. Şahı̇n, and M. Sezer, “Bessel polynomial solutions of high-order linear Volterra integro-differential equations,” Computers and Mathematics with Applications, vol. 62, no. 4, pp. 1940–1956, 2011.
View at: Publisher Site | Google Scholar
A. El-Ajou, O. Abu Arqub, and S. Momani, “Homotopy analysis method for second-order boundary value problems of integrodifferential equations,” Discrete Dynamics in Nature and Society, vol. 2012, Article ID 365792, 18 pages, 2012.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Mohamed E. A. Alnair and Ahmed A. Khidir. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

643

Downloads

884

Citations