Solution of nonlinear weakly singular Volterra integral equations using the fractional‐order Legendre functions and pseudospectral method

In this article, our main goal is to render an idea to convert a nonlinear weakly singular Volterra integral equation to a non‐singular one by new fractional‐order Legendre functions. The fractional‐order Legendre functions are generated by change of variable on well‐known shifted Legendre polynomials. We consider a general form of singular Volterra integral equation of the second kind. Then the fractional Legendre–Gauss–Lobatto quadratures formula eliminates the singularity of the kernel of the integral equation. Finally, the Legendre pseudospectral method reduces the solution of this problem to the solution of a system of algebraic equations. This method also can be utilized on fractional differential equations as well. The comparison of results of the presented method and other numerical solutions shows the efficiency and accuracy of this method. Also, the obtained maximum error between the results and exact solutions shows that using the present method leads to accurate results and fast convergence for solving nonlinear weakly singular Volterra integral equations. Copyright © 2015 John Wiley & Sons, Ltd.     

The transmutation operator method for efficient solution of the inverse Sturm‐Liouville problem on a half‐line

The inverse Sturm‐Liouville problem on a half‐line is considered. With the aid of a Fourier‐Legendre series representation of the transmutation integral kernel and the Gel'fand‐Levitan equation, the numerical solution of the problem is reduced to a system of linear algebraic equations. The potential q is recovered from the first coefficient of the Fourier‐Legendre series. The resulting numerical method is direct and simple. The results of the numerical experiments are presented.     

Wavelets collocation methods for the numerical solution of elliptic BV problems

Based on collocation with Haar and Legendre wavelets, two efficient and new numerical methods are being proposed for the numerical solution of elliptic partial differential equations having oscillatory and non-oscillatory behavior. The present methods are developed in two stages. In the initial stage, they are developed for Haar wavelets. In order to obtain higher accuracy, Haar wavelets are replaced by Legendre wavelets at the second stage. A comparative analysis of the performance of Haar wavelets collocation method and Legendre wavelets collocation method is carried out. In addition to this, comparative studies of performance of Legendre wavelets collocation method and quadratic spline collocation method, and meshless methods and Sinc–Galerkin method are also done. The analysis indicates that there is a higher accuracy obtained by Legendre wavelets decomposition, which is in the form of a multi-resolution analysis of the function. The solution is first found on the coarse grid points, and then it is refined by obtaining higher accuracy with help of increasing the level of wavelets. The accurate implementation of the classical numerical methods on Neumann’s boundary conditions has been found to involve some difficulty. It has been shown here that the present methods can be easily implemented on Neumann’s boundary conditions and the results obtained are accurate; the present methods, thus, have a clear advantage over the classical numerical methods. A distinct feature of the proposed methods is their simple applicability for a variety of boundary conditions. Numerical order of convergence of the proposed methods is calculated. The results of numerical tests show better accuracy of the proposed method based on Legendre wavelets for a variety of benchmark problems.     

The comparison of two reliable methods for accurate solution of time‐fractional Kaup–Kupershmidt equation arising in capillary gravity waves

In this paper, we compared two different methods, one numerical technique, viz Legendre multiwavelet method, and the other analytical technique, viz optimal homotopy asymptotic method (OHAM), for solving fractional‐order Kaup–Kupershmidt (KK) equation. Two‐dimensional Legendre multiwavelet expansion together with operational matrices of fractional integration and derivative of wavelet functions is used to compute the numerical solution of nonlinear time‐fractional KK equation. The approximate solutions of time fractional Kaup–Kupershmidt equation thus obtained by Legendre multiwavelet method are compared with the exact solutions as well as with OHAM. The present numerical scheme is quite simple, effective, and expedient for obtaining numerical solution of fractional KK equation in comparison to analytical approach of OHAM. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical solution of two-dimensional time fractional cable equation with Mittag-Leffler kernel

The main motive of this article is to study the recently developed Atangana-Baleanu Caputo (ABC) fractional operator that is obtained by replacing the classical singular kernel by Mittag-Leffler kernel in the definition of the fractional differential operator. We investigate a novel numerical method for the nonlinear two-dimensional cable equation in which time-fractional derivative is of Mittag-Leffler kernel type. First, we derive an approximation formula of the fractional-order ABC derivative of a function t^k using a numerical integration scheme. Using this approximation formula and some properties of shifted Legendre polynomials, we derived the operational matrix of ABC derivative. In the author of knowledge, this operational matrix of ABC derivative is derived the first time. We have shown the efficiency of this newly derived operational matrix by taking one example. Then we solved a new class of fractional partial differential equations (FPDEs) by the implementation of this ABC operational matrix. The two-dimensional model of the time-fractional model of the cable equation is solved and investigated by this method. We have shown the effectiveness and validity of our proposed method by giving the solution of some numerical examples of the two-dimensional fractional cable equation. We compare our obtained numerical results with the analytical results, and we conclude that our proposed numerical method is feasible and the accuracy can be seen by error tables. We see that the accuracy is so good. This method will be very useful to investigate a different type of model that have Mittag-Leffler fractional derivative.     

Operational matrix approach for solution of integro-differential equations arising in theory of anomalous relaxation processes in vicinity of singular point

In this paper we study the numerical solution of singular Abel–Volterra integro-differential equations, which are typical for the theory of anomalous diffusion and viscoelastic delayed stresses. The proposed method is based on application of the operational and almost operational matrices to derivatives and integrals in a vicinity of the kernel’s singular point. As examples, two orthonormal systems are considered: Bernstein polynomials and Legendre wavelets. The methods convert the singular integro-differential equation in to a system of algebraic equations that implies two advantages: (i) one does not need to introduce artificial smoothing factors into the singular integrand and (ii) the direct estimation of computational error around singular point is possible via the obtained explicit expression. The examples of numerical solution and their discussion are presented.     

Legendre spectral collocation method for volterra-hammerstein integral equation of the second kind

This paper is concerned with obtaining the approximate solution for VolterraHammerstein integral equation with a regular kernel. We choose the Gauss points associated with the Legendre weight function ω(x) = 1 as the collocation points. The Legendre collocation discretization is proposed for Volterra-Hammerstein integral equation. We provide an error analysis which justifies that the errors of approximate solution decay exponentially in L~2 norm and L~∞ norm. We give two numerical examples in order to illustrate the validity of the proposed Legendre spectral collocation method.     

再生核移位勒让德基函数法求解分数阶微分方程

该文以再生核理论为基础,用移位Legendre多项式作为基函数构造了一个新的再生核空间,并给出了该空间下的再生核函数.与经典的再生核函数有所不同的是该空间下的再生核函数不再是分段函数,因此可以减小分数阶算子作用在核函数上时的计算量,使近似解更为精确.数值算例表明该方法的有效性.     

Cardinal functions for Legendre pseudo-spectral method for solving the integro-differential equations

In this article, we present a new numerical method to solve the integro-differential equations (IDEs). The proposed method uses the Legendre cardinal functions to express the approximate solution as a finite series. In our method the operational matrix of derivatives is used to reduce IDEs to a system of algebraic equations. To demonstrate the validity and applicability of the proposed method, we present some numerical examples. We compare the obtained numerical results from the proposed method with some other methods. The results show that the proposed algorithm is of high accuracy, more simple and effective.     

Numerical solution of one-dimensional Burgers’ equation using reproducing kernel function

This paper presents a numerical method for one-dimensional Burgers’ equation by the Hopf–Cole transformation and a reproducing kernel function, abbreviated as RKF. The numerical solution is given as explicit integral expressions with the RKF at each time step, so that the computation is fully parallel. The stability and error estimates are derived. Numerical results for some test problems are presented and compared with the exact solutions. Some numerical results are also compared with the results obtained by other methods. The present method is easily implemented and effective.     

Computing Stability of Differential Equations with Bounded Distributed Delays

This paper deals with the stability analysis of scalar delay integro-differential equations (DIDEs). We propose a numerical scheme for computing the stability determining characteristic roots of DIDEs which involves a linear multistep method as time integration scheme and a quadrature method based on Lagrange interpolation and a Gauss–Legendre quadrature rule. We investigate to which extent the proposed scheme preserves the stability properties of the original equation. We derive and prove a sufficient condition for (asymptotic) stability of a DIDE (with a constant kernel) which we call RHP-stability. Conditions are obtained under which the proposed scheme preserves RHP-stability. We compare the obtained results with corresponding ones using Newton–Cotes formulas. Results of numerical experiments on computing the stability of DIDEs with constant and nonconstant kernel functions are presented.     

Numerical solution for the weakly singular Fredholm integro-differential equations using Legendre multiwavelets

An effective method based upon Legendre multiwavelets is proposed for the solution of Fredholm weakly singular integro-differential equations. The properties of Legendre multiwavelets are first given and their operational matrices of integral are constructed. These wavelets are utilized to reduce the solution of the given integro-differential equation to the solution of a sparse linear system of algebraic equations. In order to save memory requirement and computational time, a threshold procedure is applied to obtain the solution to this system of algebraic equations. Through numerical examples, performance of the present method is investigated concerning the convergence and the sparseness of the resulted matrix equation.     

Iterative reproducing kernel method for nonlinear oscillator with discontinuity

In this paper, iterative reproducing kernel method is applied to obtain the analytical approximate solution of a nonlinear oscillator with discontinuities. The solution obtained by using the method takes the form of a convergent series with easily computable components. An illustrative example is given to demonstrate the effectiveness of the present method. The results obtained using the scheme presented here show that the numerical scheme is very effective and convenient for the nonlinear oscillator with discontinuities.     

Solving high‐order linear differential equations by a Legendre matrix method based on hybrid Legendre and Taylor polynomials

A numerical method for solving the high‐order linear differential equations with variable coefficients under the mixed conditions is presented. The method is based on the hybrid Legendre and Taylor polynomials. The solution is obtained in terms of Legendre polynomials. Comparison of the present solution is made with the existing solution and excellent agreement is noted. Illustrative examples are included to demonstrate the validity and applicability of the technique. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Optimal control of lumped parameter systems via shifted Legendre polynomial approximation

Shifted Legendre polynomial functions are employed to solve the linear-quadratic optimal control problem for lumped parameter system. Using the characteristics of the shifted Legendre polynomials, the system equations and the adjoint equations of the optimal control problem are reduced to functional ordinary differential equations. The solution of the functional differential equations are obtained in a series of the shifted Legendre functions. The operational matrix for the integration of the shifted Legendre polynomial functions is also introduced in the simulation step in order to simplify the computational procedure. An illustrative example of an optimal control problem is given, and the computational results are compared with those of the exact solution. The proposed method is effective and accurate.     

The Sinc–Legendre collocation method for a class of fractional convection–diffusion equations with variable coefficients

This paper deals with the numerical solution of classes of fractional convection–diffusion equations with variable coefficients. The fractional derivatives are described based on the Caputo sense. Our approach is based on the collocation techniques. The method consists of reducing the problem to the solution of linear algebraic equations by expanding the required approximate solution as the elements of shifted Legendre polynomials in time and the Sinc functions in space with unknown coefficients. The properties of Sinc functions and shifted Legendre polynomials are then utilized to evaluate the unknown coefficients. Several examples are given and the numerical results are shown to demonstrate the efficiency of the newly proposed method.     

第二类Volterra型积分方程的Chebyshev-Legendre谱配置方法

提出了一种新的求解第二类线性Volterra型积分方程的Chebyshev谱配置方法.该方法分别对方程中积分部分的核函数和未知函数在Chebyshev-Gauss-Lobatto点上进行插值,通过Chebyshev-Legendre变换,把插值多项式表示成Legendre级数形式,从而将积分转换为内积的形式,再利用Legendre多项式的正交性进行计算.利用Chebyshev插值算子在不带权范数意义下的逼近结果,对该方法在理论上给出了L∞范数意义下的误差估计,并通过数值算例验证了算法的有效性和理论分析的正确性.     

Legendre spectral method for the fractional Bratu problem

In this paper, the Legendre spectral collocation method (LSCM) is applied for the solution of the fractional Bratu's equation. It shows the high accuracy and low computational cost of the LSCM compared with some other numerical methods. The fractional Bratu differential equation is transformed into a nonlinear system of algebraic equations for the unknown Legendre coefficients and solved with some spectral collocation methods. Some illustrative examples are also given to show the validity and applicability of this method, and the obtained results are compared with the existing studies to highlight its high efficiency and neglectable error.     

Inverse heat problem of determining time-dependent source parameter in reproducing kernel space

A new method by the reproducing kernel Hilbert space is applied to an inverse heat problem of determining a time-dependent source parameter. The problem is reduced to a system of linear equations. The exact and approximate solutions are both obtained in a reproducing kernel space. The approximate solution and its partial derivatives are proved to converge to the exact solution and its partial derivatives, respectively. The proposed method improves the existing method. Our numerical results show that the method is of high precision.     

B-Spline collocation method for linear and nonlinear Fredholm and Volterra integro-differential equations

A numerical method based on quintic B-spline has been developed to solve the linear and nonlinear Fredholm and Volterra integro-differential equations up to order 4. The solution and its derivatives are collocated by quintic B-spline and then the integral equation is approximated by the 4-points Gauss–Turán quadrature formula with respect to the weight function Legendre. The error analysis of proposed numerical method is studied theoretically. Numerical results are given to illustrate the efficiency of the proposed method which shows that our method can be applied for large values of N. The results are compared with the results obtained by other methods which show that our method is accurate.