Abstract
Since obtaining an analytic solution to some mathematical and physical problems is often very difficult, academics in recent years have focused their efforts on treating these problems using numerical methods. In science and engineering, systems of integral differential equations and their solutions are extremely important. The Taylor collocation method is described as a matrix approach for solving numerically Linear Differential Equations (LDE) by using truncated Taylor series. Integral equations are used to solve problems such as radiative transmission and the oscillation of a string, membrane, or axle. Differential equations can be used to tackle oscillating difficulties. To discover approximate solutions for linear systems of integral differential equations with variable coefficients in terms of Taylor polynomials, the collocation approach, which is offered for differential and integral equation solutions, will be developed. A system of LDE will be translated into matrix equations, and a new matrix equation will be generated in terms of the Taylor coefficients matrix by employing Taylor collocation points. The needed system will be converted to a linear algebraic equation system. Finding the Taylor coefficients will lead to the Taylor series technique.
1 Introduction
Indeed, many occurrences in practically all engineering analysis and applied science concerns, such as physical implementations, potential theory, and electrodynamics, may be simplified to solving integral equations. Because these equations are rarely solved openly, approximate solutions must be obtained. There are several numerical approaches that have been developed to solve integral equations. According to their bounds, integral equations are divided into two types: Fredholm and Volterra. The distinction between these two equations is that the integral bounds in the Fredholm equation are constant, whereas the upper bound in the Volterra equation is changeable. The Taylor sorting method will be discussed in this part to acquire the Fredholm integral equation system’s finite Taylor series solutions.
Highorder linear systems in science and engineering. The F–V integrodifferential equations and their solution are extremely important. Many system components, such as biomedical application in genetics and population dynamics, wherein impulses occur naturally, systems of integrodifferential equations can be used to model systems of integrodifferential equations that are formed via control. Integral systems are also linked to some controllability difficulties for impulsive systems and integrodifferential theory conclusions, but they are yet to be thoroughly investigated. As a result, solutions to these and many other equations are critical in a variety of applications. Infinite series techniques can be used to solve a variety of physical problems. These are used to approximate various functions, solve the system of integrodifferential equations and evaluate challenging integrals. Taylor’s expansion is the most often used infinite series expansion method for this purpose. It is a strategy for modeling a function near to a start or initial point that is used locally. It does not, however, work effectively when a significant number of control factors must always be taken into account. The Taylor collocation method is the emphasis of this study, and it is applied to a variety of fascinating scenarios, including systems of integrodifferential equations. Integral equations were also influenced by mathematical physics models such as diffraction issues, quantum mechanics scattering, conformal mapping, and water waves. Integral equations are used to solve problems like radiative transfer and the oscillation of a string, membrane, or axle. Differential equations may mathematically characterize engineering difficulties and as a result, differential equations play a critical part in the solution of practical problems. Newton’s law, which states that the rate of change of a particle’s momentum is equal to the magnitude acting on it, can be expressed mathematically as a differential equation. Issues occurring in electric circuits, chemical kinetics, and heat transfer in a medium can all be theoretically expressed as differential equations. These differential equations can be turned into Volterra and Fredholmtype integral equations.
Kanwal and Liu studied Taylor expansion approaches for resolving integral equations that have already been given for Volterra integral equation [1]. Sezer probed Taylor’s polynomial solution of Volterra integral equations [2]. Sezer probed it for some differential equations [3], Nas et al. for integrodifferential equations [4], Karamete and Sezer for linear integrodifferential equations [5], and Yalçınbas and Sezer for F–V linear integrodifferential equations [6]. Akyüz and Sezer studied Chebyshev methods that have been used to solve highorder linear differential equation systems [7]. AkyüzDaşcıoǧlu explained the linear integral equations [8], which is defined as
Yusufoğlu probed the numerical solution of Duffing equation by the Laplace decomposition algorithm [9]. Gülsu and Sezer explained Taylor polynomials, a Taylor collocation approach is provided for numerically solving the system of highorder linear Fredholm–Volterra integrodifferential equations [10]. Yuldashev mixed the value problem for the nonlinear integrodifferential equation with a parabolic operator of higher power [11]. Toutounian et al. solved highorder linear differential and difference equations with constant coefficients, and Fourier operational matrices of differentiation and transmission are introduced. They also used Legendre–Gauss collocation nodes to enhance our approaches for extended Pantograph equations with variable coefficients [12].
Tohidi and Kılıçman discussed how to solve a variety of onedimensional parabolic partial differential equations (PDEs) with given initial and nonlocal boundary conditions. The central concept is based on direct collocation and translating the PDEs under consideration into corresponding algebraic equations [13]. Mirzaee et al. presented a unique method that is under mixed beginning conditions, a method for solving linear complex differential equations is developed by using complex operational matrices of Euler functions and their fascinating properties [14]. Tohidi presented an effective numerical approach with piecewise intervals for solving linear twodimensional Fredholm integral equations that produce an approximate polynomial solution [15]. Zogheib and Tohidi explained numerical solution of nonlinear one and twodimensional heat transfer problems with provided beginning conditions and linear Robin boundary conditions is considered [16]. Zogheib and Tohidi presented a computational solution for multidimensional hyperbolic telegraph formulas with suitable initial time and boundary space conditions. In both the spatial and temporal variables, the truncated Hermite series with unknown coefficients are employed to approximate the answer [17]. Zézé et al. found multipoint Taylor series to solve differential equations [18].
Hadadian and Yousefi derived a new and efficient numerical approach for solving advectiondiffusion equations (ADE). The method is based on the combination of the ADE under consideration. The ADE is transformed into a corresponding integral equation via integration, and the integral equation comprises initial and boundary conditions [19]. Yuldashev examined a boundaryvalue problem for Boussinesq type nonlinear integrodifferential equation with reflecting argument [20]. Singh et al. proposed an efficient matrix approach for numerically solving the twodimensional diffusion and telegraph equations subject to Dirichlet boundary conditions using Euler approximation [21]. Yang et al. contributed to the development of a highly accurate computational approach for solving twodimensional mass transport equations during apple slice convective air drying. Numerical technique is based on interpolating the solution of the above equations in a nodal form across the roots of orthogonal Jacobi polynomials [22].
The collocation approach, which is provided for the solution of differential and integral equations, is expanded to find approximate solutions for linear systems of integral differential equations with variable coefficients in terms of Taylor polynomials. A set of linear differential equations is converted into matrix equations, and by using Taylor collocation points, a new matrix equation is created in terms of the Taylor coefficient matrix. The required system is transformed into a system of linear algebraic equations. The Taylor series method is reached by locating the Taylor coefficients. Integrals are utilized in a variety of industries, including engineering, where engineers use integrals to determine the geometry of a building. In physics, it is utilized in the center of gravity, for example. Threedimensional designs are exhibited in the realm of schematic diagram.
2 Fundamental relations
Linear Fredholm integral equation consisting of k equations with k unknown systems is defined as
In short, this system is
where kernel functions K(x, t), P(x), and f(x) are known functions in Eq. (1) problem.
The finite Taylor approach will be examined in the form. Taylor series expansions will be presumed for the time being of K(x,t) functions that exist. Writing the system of Eq. (1) in matrix form:
The matrices are
Now to an integral part of this system
If it is said that system (3) takes shape.
Now, we add Taylor sorting points defined to this system
Matrices form
Including
can be written in this format
Let us try to write the matrix I here in terms of A. The rows of the matrix
Now, Taylor sort order
Let us place the points:
where K(x, t) function now depends only on t. Therefore,
can open a univariate Taylor expression. Thus, in the matrix form of Eq. (7) spelling
included
is in the form. If K is (x _{ s }) elements,
If we use the previously determined matrix notation for y(t), we get from Eq. (6)
If the integral is included in the above sum, then
The integral here
The elements of this matrix are
Thus, the Eq. (6) is calculated as
Let us write this clearly:
Thus, the matrix representation of Eq. (9) is possible.
The matrices are
The matrix I in Eq. (4) can specified as the terms of Taylor coefficient matrix.
To find the expression, Eq. (10) is substituted in the matrix I;
and if the derived expression is placed in Eq. (5), the result is obtained as
When this expression is edited
The matrix values in the brackets in Eq. (11) are known, and as a result of the operations, a square matrix of size k(N + 1) is obtained as
This is a linear system consisting of k(N + 1) with k(N + 1) unknowns. The Taylor coefficients are calculated by solving this system. When these coefficients are placed in Eq. (2), the Fredholm integral system’s approximate solution in terms of Taylor polynomials is found.
3 Illustrations
Example 1:
Solution:
Fundamental equation
For p(x) = 0, then
For N = 2, then the points are given by
Now
The exact solution is given by
Example 2:
Solution:
Taylor polynomial approach
Comparison tables
Similarly, the results are obtained for different values of N as presented in Figure 1 and Tables 1–3.
N = 2  N = 2  N = 3  N = 3  Chebyshev  Method  

t _{ i }  y _{1}(t _{ i })  y _{2}(t _{ i })  y _{1}(t _{ i })  y _{2}(t _{ i })  y _{1}(t _{ i })  y _{2}(t _{ i }) 
0  −2.858 × 10^{−4}  1  −8.549 × 10^{−4}  1  0  1 
0.1  0.1  1  0.0999  1  0.1  1 
0.2  0.2  1  0.1999  1  0.2  1 
0.3  0.3  1  0.2999  1  0.3  1 
0.4  0.4  1  0.3999  1  0.4  1 
0.5  0.5  1  0.5  1.001  0.5  1 
0.6  0.6  1  0.6  1.001  0.6  1 
0.7  0.7  1  0.7  1.001  0.7  1 
0.8  0.8  1  0.8  1.001  0.8  1 
0.9  0.9  1  0.9  1  0.9  1 
1  1  1  1  1  1  1 

Present  Chebyshev  Exact 

N = 6  N = 6  N = 6  
−1  0.476  0.498  0.4987 
−1/2  0.5987  0.5876  0.5321 
0  0.8765  0.8543  0.86554 
1/2  0.465  0.4987  0.4987 
1  0.5876  0.567  0.5434 

Present  Chebyshev  Exact 

N = 6  N = 6  N = 6  
–1  0.6557  0.6574  0.6572 
–1/2  0.70435  0.7314  0.7288 
0  0.56543  0.64375  0.5743 
1/2  0.5389  0.5763  0.5373 
1  0.5876  0.567  0.5434 
Example 3:
Solution:
For N = 2
After simplification, the Taylor coefficient of the matrix is obtained
4 Conclusion
Many scientific and engineering issues use integral equations. Volterra or Fredholm integral equations can be used to solve a wide range of starting and boundary value issues. Potential theory assisted more than any other discipline to the development of integral equations. Integral equations were also influenced by mathematical physics models such as dispersion difficulties, scattered in quantum physics, conformal mapping, and water waves. Solving integraldifferential equalities analytically is frequently difficult. Obtaining approximate solutions is usually required. For this purpose, the procedure mentioned can be recommended. The coefficients in Taylor expansion of the solution of a linear integrodifferential equation are calculated using the approach outlined in this study, and it works whenever the function is a linear integrodifferential equation. P _{ k }(x) is the kernel function of K(x,t). In this range, there is a Taylor series expansion. Furthermore, it appears that the approach works best when the functions are combined. It is possible to extend the Taylor series of K(x,t), f(x), and P _{ k }(x), which converges quickly. We select several terms from the Taylor expansion to get the best approximation solution of the equation; that is, the truncation limit N must have been large enough. An estimate for N, the degree of the approximation polynomial, for computing efficiency (the truncation limit of the Taylor series) to y(x) should also be available. The Taylor collocation method can be used to solve both differential and integral problems. This is demonstrated in example 3. Both differential and integral formulations are the foundation of algorithms for the numerical solution of continuous electromagnetic field issues. The former strategy, utilizing the finite element method, has been the most popular, at least for lowfrequency issues. The distinctive benefits of integral equations over differential equations are explored, along with some of the associated challenges. It is suggested that the integral equation technique may be particularly successful when applied to more complex situations, such as moving systems and optimization.
This method can also be used to find an exact solution N or less than N if the equation has an exact solution which is a polynomial of degree. Partial integrodifferential equations and ordinary differential equations with variable coefficients can also be solved using this method.
Acknowledgments
This research was funded by the Princess Nourahbint Abdulrahman University Researchers Supporting Project number (PNURSP2022R152), Princess Nourahbint Abdulrahman University, Riyadh, Saudi Arabia. The authors express their gratitude.

Funding information: This research was funded by the Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R152), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Author contributions: M.G. and L.W. framed the issue. M.G., K.S., W.J., and R.S. resolved the problem. M.G., L.W., K.S., R.S., S.E.A., M.M.A.L., W.J., M.R.E., and S.A. computed and analyzed the results. All the authors equally contributed to the writing and proofreading of the article. All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

Conflict of interest: The authors state no conflict of interest.
References
[1] Kanwal RP, Liu KC. A Taylor expansion approach for solving integral equation. Int J Math Edu Sci Tech. 1989;20:411–4.10.1080/0020739890200310Search in Google Scholar
[2] Sezer M. Taylor polynomial solution of Volterra integral equations. Int J Math Edu Sci Tech. 1994;25:625–33.10.1080/0020739940250501Search in Google Scholar
[3] Sezer M. A method for approximate solution of the second order linear differential equations in terms of Taylor polynomials. Int J Math Edu Sci Tech. 1996;27:821–34.10.1080/0020739960270606Search in Google Scholar
[4] Nas S, Yalçınba S, Sezer M. A Taylor polynomial approach for solving highorder linear Fredholm integrodifferential equations. Int J Math Edu Sci Tech. 2000;31:213–25.10.1080/0020739X.2022.12131593Search in Google Scholar
[5] Karamete A, Sezer M. A Taylor collocation method for the solution of linear integrodifferential equations. Int J Comput Math. 2002;79:987–1000.10.1080/00207160212116Search in Google Scholar
[6] Yalçınbas S, Sezer M. The approximate solution of highorder linear Volterra–Fredholm integrodifferential equations in terms of Taylor polynomials. Appl Math Comput. 2002;112:291–308.10.1016/S00963003(99)000594Search in Google Scholar
[7] Akyüz A, Sezer M. Chebyshev polynomial solutions of systems of highorder linear differential equations with variable coefficients. Appl Math Comput. 2003;144:237–47.10.1016/S00963003(02)004034Search in Google Scholar
[8] AkyüzDaşcıoǧlu A. Chebyshev polynomial solutions of systems of linear integral equations. Appl Math Comput. 2004;151:221–32.10.1016/S00963003(03)003345Search in Google Scholar
[9] Yusufoğlu E. Numerical solution of Duffing equation by the Laplace decomposition algorithm. Appl Math Comput. 2006;177:572–80.10.1016/j.amc.2005.07.072Search in Google Scholar
[10] Gülsu M, Sezer M. A fast and efficient numerical approach for solving advectiondiffusion equations by using hybrid functions. Comput Appl Math. 2006;38(4):1–19.Search in Google Scholar
[11] Yuldashev TK. Mixed value problem for nonlinear integrodifferential equation with parabolic operator of higher power. Comput Math Math Phys. 2012;52:105.10.1134/S0965542512010150Search in Google Scholar
[12] Toutounian F, Tohidi E, Kilicman A. Fourier operational matrices of differentiation and transmission: introduction and applications. Abs Appl Anal. 2013;2013:198926.10.1155/2013/198926Search in Google Scholar
[13] Tohidi E, Kılıçman A. An efficient spectral approximation for solving several types of parabolic PDEs with nonlocal boundary conditions. Math Prob Eng. 2014;2014:369029.10.1155/2014/369029Search in Google Scholar
[14] Mirzaee F, Bimesl S, Tohidi E. A new complexvalued method and its applications in solving differential equations. Sci Iran. 2015;22(6):2424–31.Search in Google Scholar
[15] Tohidi E. Taylor matrix method for solving linear twodimensional Fredholm integral equations with piecewise intervals. Comput Appl Math. 2015;34(3):1117–30.10.1007/s4031401401663Search in Google Scholar
[16] Zogheib B, Tohidi E. An accurate spacetime pseudospectral method for solving nonlinear multidimensional heat transfer problems. Mediterr J Math. 2017;14(1):1–19.10.1007/s0000901608270Search in Google Scholar
[17] Zogheib B, Tohidi E. Modal Hermite spectral collocation method for solving multidimensional hyperbolic telegraph equations. Comput Math Appl. 2018;75(10):3571–88.10.1016/j.camwa.2018.02.018Search in Google Scholar
[18] Zézé DS, PotierFerry M, Tampango Y. Multipoint Taylor series to solve differential equations. Discrete Cont DynS. 2019;12:1791–8.10.3934/dcdss.2019118Search in Google Scholar
[19] Hadadian M, Yousefi N. A new spectral integral equation method for solving twodimensional unsteady advectiondiffusion equations via Chebyshev polynomials. Eng Comput. 2019;36(7):2327–68.10.1108/EC0220180063Search in Google Scholar
[20] Yuldashev TK. On a boundaryvalue problem for Boussinesq type nonlinear integradifferential equation with reflecting argument. Lobachevskii J Math. 2020;53:111–23.10.1134/S1995080220010151Search in Google Scholar
[21] Singh S, Devi V, Tohidi E, Singh VK. An efficient matrix approach for twodimensional diffusion and telegraph equations with Dirichlet boundary conditions. Phys A Stat Mech Appl. 2020;545:123784.10.1016/j.physa.2019.123784Search in Google Scholar
[22] Yang Y, Rządkowski G, Pasban A, Tohidi E. A high accurate scheme for numerical simulation of twodimensional mass transfer processes in food engineering. Alex Eng J. 2021;60(2):2629–39.10.1016/j.aej.2020.12.055Search in Google Scholar
© 2022 Madiha Ghamkhar et al., published by De Gruyter
This work is licensed under the Creative Commons Attribution 4.0 International License.