A fast numerical algorithm based on the Taylor wavelets for solving the fractional integro‐differential equations with weakly singular kernels

In this paper, a fast numerical algorithm based on the Taylor wavelets is proposed for finding the numerical solutions of the fractional integro‐differential equations with weakly singular kernels. The properties of Taylor wavelets are given, and the operational matrix of fractional integration is constructed. These wavelets are utilized to reduce the solution of the given fractional integro‐differential equation to the solution of a linear system of algebraic equations. Also, convergence of the proposed method is studied. Illustrative examples are included to demonstrate the validity and applicability of the technique.     

Numerical solution of delay integro‐differential equations by using Taylor collocation method

This paper is concerned with the numerical solution of delay integro‐differential equations. The main purpose of this work is to provide a new numerical approach based on the use of continuous collocation Taylor polynomials for the numerical solution of delay integro‐differential equations. It is shown that this method is convergent. Numerical illustrations confirm our theoretical analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

Solution of high‐order linear Fredholm integro‐differential equations with piecewise intervals

In this study, a practical matrix method is presented to find an approximate solution for high‐order linear Fredholm integro‐differential equations with piecewise intervals under the initial boundary conditions in terms of Taylor polynomials. The method converts the integro differential equation to a matrix equation, which corresponds to a system of linear algebraic equations. Error analysis and illustrative examples are included to demonstrate the validity and applicability of the technique. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010 27: 1327–1339, 2011     

Approximations of the nonlinear Volterra's population model by an efficient numerical method

In this paper, we apply the new homotopy perturbation method to solve the Volterra's model for population growth of a species in a closed system. This technique is extended to give solution for nonlinear integro‐differential equation in which the integral term represents the total metabolism accumulated fromtime zero. The approximate analytical procedure only depends on two components. The newhomotopy perturbationmethodwas applied to nonlinear integro‐differential equations directly and by converting the problem into nonlinear ordinary differential equation. We also compare this method with some other numerical results and show that the present approach is less computational and is applicable for solving nonlinear integro‐differential equations and ordinary differential equations as well. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical solution of nonlinear partial quadratic integro‐differential equations of fractional order via hybrid of block‐pulse and parabolic functions

In this paper, an effective numerical approach based on a new two‐dimensional hybrid of parabolic and block‐pulse functions (2D‐PBPFs) is presented for solving nonlinear partial quadratic integro‐differential equations of fractional order. Our approach is based on 2D‐PBPFs operational matrix method together with the fractional integral operator, described in the Riemann–Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations, which greatly simplifies the problem. By using Newton's iterative method, this system is solved, and the solution of fractional nonlinear partial quadratic integro‐differential equations is achieved. Convergence analysis and an error estimate associated with the proposed method is obtained, and it is proved that the numerical convergence order of the suggested numerical method is O(h³) . The validity and applicability of the method are demonstrated by solving three numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the exact solutions much easier.     

Regularity theory for fully nonlinear integro‐differential equations

We consider nonlinear integro‐differential equations like the ones that arise from stochastic control problems with purely jump Lévy processes. We obtain a nonlocal version of the ABP estimate, Harnack inequality, and interior C^{1, α} regularity for general fully nonlinear integro‐differential equations. Our estimates remain uniform as the degree of the equation approaches 2, so they can be seen as a natural extension of the regularity theory for elliptic partial differential equations. © 2008 Wiley Periodicals, Inc.     

Two‐dimensional shifted Legendre polynomial collocation method for electromagnetic waves in dielectric media via almost operational matrices

In this paper, a numerical solution of fractional partial differential equations (FPDEs) for electromagnetic waves in dielectric media will be discussed. For the solution of FPDEs, we developed a numerical collocation method using an algorithm based on two‐dimensional shifted Legendre polynomials approximation, which is proposed for electromagnetic waves in dielectric media. By implementing the partial Riemann–Liouville fractional derivative operators, two‐dimensional shifted Legendre polynomials approximation and its operational matrix along with collocation method are used to convert FPDEs first into weakly singular fractional partial integro‐differential equations and then converted weakly singular fractional partial integro‐differential equations into system of algebraic equation. Some results concerning the convergence analysis and error analysis are obtained. Illustrative examples are included to demonstrate the validity and applicability of the technique. Copyright © 2017 John Wiley & Sons, Ltd.     

Solution of Volterra's population model via block‐pulse functions and Lagrange‐interpolating polynomials

A numerical method for solving Volterra's population model for population growth of a species in a closed system is proposed. Volterra's model is a nonlinear integro‐differential equation where the integral term represents the effects of toxin. The approach is based on hybrid function approximations. The properties of hybrid functions that consist of block‐pulse and Lagrange‐interpolating polynomials are presented. The associated operational matrices of integration and product are then utilized to reduce the solution of Volterra's model to the solution of a system of algebraic equations. The method is easy to implement and computationally very attractive. Applications are demonstrated through an illustrative example. Copyright © 2008 John Wiley & Sons, Ltd.     

On integro‐differential inclusions with operator‐valued kernels

We study integro‐differential inclusions in Hilbert spaces with operator‐valued kernels and give sufficient conditions for the well‐posedness. We show that several types of integro‐differential equations and inclusions are covered by the class of evolutionary inclusions, and we therefore give criteria for the well‐posedness within this framework. As an example, we apply our results to the equations of visco‐elasticity and to a class of nonlinear integro‐differential inclusions describing phase transition phenomena in materials with memory. Copyright © 2014 John Wiley & Sons, Ltd.     

Boundary‐value problems for the third‐order loaded equation with noncharacteristic type‐change boundaries

The purpose of this paper is to establish unique solvability for a certain generalized boundary‐value problem for a loaded third‐order integro‐differential equation with variable coefficients. Moreover, the method of integral equations is applied to obtain an equation related to the Riemann‐Liouville operators.     

A Taylor polynomial approach for solving the most general linear Fredholm integro‐differential‐difference equations

In this study, a matrix method is developed to solve approximately the most general higher order linear Fredholm integro‐differential‐difference equations with variable coefficients under the mixed conditions in terms of Taylor polynomials. This technique reduces the problem into the linear algebraic system. The method is valid for any combination of differential, difference and integral equations. An initial value problem and a boundary value problem are also presented to illustrate the accuracy and efficiency of the method. Copyright © 2012 John Wiley & Sons, Ltd.     

Spectral collocation methods for nonlinear weakly singular Volterra integro‐differential equations

In this paper, a spectral collocation approximation is proposed for neutral and nonlinear weakly singular Volterra integro‐differential equations (VIDEs) with non‐smooth solutions. We use some suitable variable transformations to change the original equation into a new equation, so that the solution of the resulting equation possesses better regularity, and the the Jacobi orthogonal polynomial theory can be applied conveniently. Under reasonable assumptions on the nonlinearity, we carry out a rigorous error analysis in L^∞ norm and weighted L² norm. To perform the numerical simulations, some test examples (linear and nonlinear) are considered with nonsmooth solutions, and numerical results are presented. Further more, the comparative study of the proposed methods with some existing numerical methods is provided.     

Numerical approximations for population growth model by rational Chebyshev and Hermite functions collocation approach: A comparison

This paper aims to compare rational Chebyshev (RC) and Hermite functions (HF) collocation approach to solve Volterra's model for population growth of a species within a closed system. This model is a nonlinear integro‐differential equation where the integral term represents the effect of toxin. This approach is based on orthogonal functions, which will be defined. The collocation method reduces the solution of this problem to the solution of a system of algebraic equations. We also compare these methods with some other numerical results and show that the present approach is applicable for solving nonlinear integro‐differential equations. Copyright © 2010 John Wiley & Sons, Ltd.     

H1‐Galerkin expanded mixed finite element methods for nonlinear pseudo‐parabolic integro‐differential equations

H¹‐Galerkin mixed finite element method combined with expanded mixed element method is discussed for nonlinear pseudo‐parabolic integro‐differential equations. We conduct theoretical analysis to study the existence and uniqueness of numerical solutions to the discrete scheme. A priori error estimates are derived for the unknown function, gradient function, and flux. Numerical example is presented to illustrate the effectiveness of the proposed scheme. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Computational methods of solving the boundary value problems for the loaded differential and Fredholm integro‐differential equations

The article presents a new general solution to a loaded differential equation and describes its properties. Solving a linear boundary value problem for loaded differential equation is reduced to the solving a system of linear algebraic equations with respect to the arbitrary vectors of general solution introduced. The system's coefficients and right sides are computed by solving the Cauchy problems for ordinary differential equations. Algorithms of constructing a new general solution and solving a linear boundary value problem for loaded differential equation are offered. Linear boundary value problem for the Fredholm integro‐differential equation is approximated by the linear boundary value problem for loaded differential equation. A mutual relationship between the qualitative properties of original and approximate problems is obtained, and the estimates for differences between their solutions are given. The paper proposes numerical and approximate methods of solving a linear boundary value problem for the Fredholm integro‐differential equation and examines their convergence, stability, and accuracy.     

Schwarz problem for higher‐order complex partial differential equations in the upper half plane

Linear and nonlinear elliptic complex partial differential equations of higher‐order are considered under Schwarz conditions in the upper‐half plane. Firstly, using the integral representations for the solutions of the inhomogeneous polyanalytic equation with Schwarz conditions, a class of integral operators is introduced together with some of their properties. Then, these operators are used to transform the problem for linear equations into singular integral equations. In the case of nonlinear equations such a transformation yields a system of integro‐differential equations. Existence of the solutions of the relevant boundary value problems for linear and nonlinear equations are discussed via Fredholm theory and fixed point theorems, respectively.     

A numerical approach for solving weakly singular partial integro‐differential equations via two‐dimensional‐orthonormal Bernstein polynomials with the convergence analysis

In this paper, we develop an efficient matrix method based on two‐dimensional orthonormal Bernstein polynomials (2D‐OBPs) to provide approximate solution of linear and nonlinear weakly singular partial integro‐differential equations (PIDEs). First, we approximate all functions involved in the considerable problem via 2D‐OBPs. Then, by using the operational matrices of integration, differentiation, and product, the solution of Volterra singular PIDEs is transformed to the solution of a linear or nonlinear system of algebraic equations which can be solved via some suitable numerical methods. With a small number of bases, we can find a reasonable approximate solution. Moreover, we establish some useful theorems for discussing convergence analysis and obtaining an error estimate associated with the proposed method. Finally, we solve some illustrative examples by employing the presented method to show the validity, efficiency, high accuracy, and applicability of the proposed technique.     

The approximate solution of integro‐ differential equations

We consider a class of partial integro‐differential equations with appropriate conditions and its corresponding equation in the field of Mikusiński operators. As is usual in numerical analysis, we construct the corresponding difference equation, determine its solution, analyze its character and treat it as the approximate solution of the considered problem. We also estimate the error of approximation.     

The Taylor method in the field of Mikusiński operators

The approximate solution of a class of differential equations, in the field of Mikusiński operators is constructed by using the Taylor method. The obtained results are applied to a class of partial integro–differential equations. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

An efficient pseudo‐spectral Legendre–Galerkin method for solving a nonlinear partial integro‐differential equation arising in population dynamics

The pseudo‐spectral Legendre–Galerkin method (PS‐LGM) is applied to solve a nonlinear partial integro‐differential equation arising in population dynamics. This equation is a competition model in which similar individuals are competing for the same resources. It is a kind of reaction–diffusion equation with integral term corresponding to nonlocal consumption of resources. The proposed method is based on the Legendre–Galerkin formulation for the linear terms and interpolation operator at the Chebyshev–Gauss–Lobatto (CGL) points for the nonlinear terms. Also, the integral term, which is a kind of convolution, is directly computed by a fast and accurate method based on CGL interpolation operator, and thus, the use of any quadrature formula in its computation is avoided. The main difference of the PS‐LGM presented in the current paper with the classic LGM is in treating the nonlinear terms and imposing boundary conditions. Indeed, in the PS‐LGM, the nonlinear terms are efficiently handled using the CGL points, and also the boundary conditions are imposed strongly as collocation methods. Combination of the PS‐LGM with a semi‐implicit time integration method such as second‐order backward differentiation formula and Adams‐Bashforth method leads to reducing the complexity of computations and obtaining a linear algebraic system of equations with banded coefficient matrix. The desired equation is considered on one and two‐dimensional spatial domains. Efficiency, accuracy, and convergence of the proposed method are demonstrated numerically in both cases. Copyright © 2012 John Wiley & Sons, Ltd.