Numerical Approximate Methods for Solving Linear and Nonlinear Integral Equations

Integral equation has been one of the essential tools for various area of applied mathematics. In this work, we employed different numerical methods for solving both linear and nonlinear Fredholm integral equations. A goal is to categorize the selected methods and assess their accuracy and efficiency. We discuss challenges faced by researchers in this field, and we emphasize the importance of interdisciplinary effort for advancing the study on numerical methods for solving integral equations. Integral equations can be viewed as equations which are results of transformation of points in a given vector spaces of integrable functions by the use of certain specific integral operators to points in the same space. If, in particular, one is concerned with function spaces spanned by polynomials for which the kernel of the corresponding transforming integral operator is separable being comprised of polynomial functions only, then several approximate methods of solution of integral equations can be developed. This work, specially, deals with the development of different wavelet methods for solving integral and intgro-differential equations. Wavelets theory is a relatively new and emerging area in mathematical research. It has been applied in a wide range of engineering disciplines; particularly, wavelets are very successfully used in signal analysis for waveform representations and segmentations, time frequency analysis, and fast algorithms for easy implementation. Wavelets permit the accurate representation of a variety of functions and operators. Moreover, wavelets establish a connection with fast numerical algorithms. Wavelets can be separated into two distinct types, orthogonal and semi-orthogonal. The preliminary concept of integral equations and wavelets are first presented in Chapter 1. Classification of integral equations, construction of wavelets and multi-resolution analysis (MRA) have been briefly discussed and provided in this chapter. In Chapter 2, different wavelet methods are constructed and function approximation by these methods with convergence analysis have been presented. In Chapter 3, linear semi-orthogonal compactly supported B-spline wavelets together with their dual wavelets have been applied to approximate the solutions of Fredholm integral equations (both linear and nonlinear) of the second kind and their systems. Properties of these wavelets are first presented; these properties are then utilized to reduce the computation of integral equations to some algebraic equations. Convergence analysis of B-spline method has been discussed in this chapter. Again, in Chapter 4, system of nonlinear Fredholm integral equations have been solved by using hybrid Legendre Block-Pulse functions and xiii Bernstein collocation method. In Chapter 5, two practical problems arising from chemical phenomenon, have been modeled as Fredholm- Hammerstein integral equations and solved numerically by different numerical techniques. First, COSMO-RS model has been solved by Bernstein collocation method, Haar wavelet method and Sinc collocation method. Second, Hammerstein integral equation arising from chemical reactor theory has been solved by B-spline wavelet method. Comparison of results have been demonstrated through illustrative examples. In Chapter 6, Legendre wavelet method and Bernoulli wavelet method have been developed to solve system of integro-differential equations. Legendre wavelets along with their operational matrices are developed to approximate the solutions of system of nonlinear Volterra integro-differential equations. Also, nonlinear Volterra weakly singular integro-differential equations system has been solved by Bernoulli wavelet method. The properties of these wavelets are used to reduce the system of integral equations to a system of algebraic equations which can be solved numerically by Newton's method. Rigorous convergence analysis has been done for these wavelet methods. Illustrative examples have been included to demonstrate the validity and applicability of the proposed techniques. In Chapter 7, we have solved the second order Lane-Emden type singular differential equation. First, the second order differential equation is transformed into integro-differential equation and then solved by Legendre multi-wavelet method and Chebyshev wavelet method. Convergence of these wavelet methods have been discussed in this chapter. In Chapter 8, we have developed a efficient collocation technique called Legendre spectral collocation method to solve the Fredholm integro-differential-difference equations with variable coefficients and system of two nonlinear integro-differential equations which arise in biological model. The proposed method is based on the Gauss-Legendre points with the basis functions of Lagrange polynomials. The present method reduces this model to a system of nonlinear algebraic equations and again this algebraic system has been solved numerically by Newton's method. The study of fuzzy integral equations and fuzzy differential equations is an emerging area of research for many authors. In Chapter 9, we have proposed some numerical techniques for solving fuzzy integral equations and fuzzy integro-differential equations. Fundamentals of fuzzy calculus have been discussed in this chapter. Nonlinear fuzzy Hammerstein integral equation has been solved by Bernstein polynomials and Legendre wavelets, and then compared with homotopy analysis method. We have solved nonlinear fuzzy Hammerstein Volterra integral equations with constant delay by Bernoulli wavelet method and then compared with B-spline wavelet method. Finally, fuzzy integro-differential equation has been solved by Legendre wavelet method and compared with homotopy analysis method. In fuzzy case, we have applied two-dimensional numerical methods which are discussed in chapter 2. Convergence analysis and error estimate have been also provided for Bernoulli wavelet method. xiv The study of fractional calculus, fractional differential equations and fractional integral equations has a great importance in the field of science and engineering. Most of the physical phenomenon can be best modeled by using fractional calculus. Applications of fractional differential equations and fractional integral equations create a wide area of research for many researchers. This motivates to work on fractional integral equations, which results in the form of Chapter 10. First, the preliminary definitions and theorems of fractional calculus have been presented in this chapter. The nonlinear fractional mixed Volterra-Fredholm integro-differential equations along with mixed boundary conditions have been solved by Legendre wavelet method. A numerical scheme has been developed by using Petrov-Galerkin method where the trial and test functions are Legendre wavelets basis functions. Also, this method has been applied to solve fractional Volterra integro-differential equations. Uniqueness and existence of the problem have been discussed and the error estimate of the proposed method has been presented in this work. Sinc Galerkin method is developed to approximate the solution of fractional Volterra-Fredholm integro-differential equations with weakly singular kernels. The proposed method is based on the Sinc function approximation. Uniqueness and existence of the problem have been discussed and the error analysis of the proposed method have been presented in this chapte

Status of the differential transformation method

Further to a recent controversy on whether the differential transformation method (DTM) for solving a differential equation is purely and solely the traditional Taylor series method, it is emphasized that the DTM is currently used, often only, as a technique for (analytically) calculating the power series of the solution (in terms of the initial value parameters). Sometimes, a piecewise analytic continuation process is implemented either in a numerical routine (e.g., within a shooting method) or in a semi-analytical procedure (e.g., to solve a boundary value problem). Emphasized also is the fact that, at the time of its invention, the currently-used basic ingredients of the DTM (that transform a differential equation into a difference equation of same order that is iteratively solvable) were already known for a long time by the "traditional"-Taylor-method users (notably in the elaboration of software packages --numerical routines-- for automatically solving ordinary differential equations). At now, the defenders of the DTM still ignore the, though much better developed, studies of the "traditional"-Taylor-method users who, in turn, seem to ignore similarly the existence of the DTM. The DTM has been given an apparent strong formalization (set on the same footing as the Fourier, Laplace or Mellin transformations). Though often used trivially, it is easily attainable and easily adaptable to different kinds of differentiation procedures. That has made it very attractive. Hence applications to various problems of the Taylor method, and more generally of the power series method (including noninteger powers) has been sketched. It seems that its potential has not been exploited as it could be. After a discussion on the reasons of the "misunderstandings" which have caused the controversy, the preceding topics are concretely illustrated.Comment: To appear in Applied Mathematics and Computation, 29 pages, references and further considerations adde

Lp-solution to the random linear delay differential equation with stochastic forcing term

[EN] This paper aims at extending a previous contribution dealing with the random autonomous-homogeneous linear differential equation with discrete delay tau > 0, by adding a random forcing term f(t) that varies with time: x'(t) = ax(t) + bx(t-tau) + f(t), t >= 0, with initial condition x(t) = g(t), -tau <= t <= 0. The coefficients a and b are assumed to be random variables, while the forcing term f(t) and the initial condition g(t) are stochastic processes on their respective time domains. The equation is regarded in the Lebesgue space L-p of random variables with finite p-th moment. The deterministic solution constructed with the method of steps and the method of variation of constants, which involves the delayed exponential function, is proved to be an L-p-solution, under certain assumptions on the random data. This proof requires the extension of the deterministic Leibniz's integral rule for differentiation to the random scenario. Finally, we also prove that, when the delay tau tends to 0, the random delay equation tends in L-p to a random equation with no delay. Numerical experiments illustrate how our methodology permits determining the main statistics of the solution process, thereby allowing for uncertainty quantification.This work has been supported by the Spanish Ministerio de Economia, Industria y Competitividad (MINECO), the Agencia Estatal de Investigacion (AEI) and Fondo Europeo de Desarrollo Regional (FEDER UE) grant MTM2017-89664-P.Cortés, J.; Jornet, M. (2020). Lp-solution to the random linear delay differential equation with stochastic forcing term. Mathematics. 8(6):1-16. https://doi.org/10.3390/math8061013S11686Xiu, D., & Karniadakis, G. E. (2004). Supersensitivity due to uncertain boundary conditions. International Journal for Numerical Methods in Engineering, 61(12), 2114-2138. doi:10.1002/nme.1152Casabán, M.-C., Cortés, J.-C., Navarro-Quiles, A., Romero, J.-V., Roselló, M.-D., & Villanueva, R.-J. (2016). A comprehensive probabilistic solution of random SIS-type epidemiological models using the random variable transformation technique. Communications in Nonlinear Science and Numerical Simulation, 32, 199-210. doi:10.1016/j.cnsns.2015.08.009Strand, J. . (1970). Random ordinary differential equations. Journal of Differential Equations, 7(3), 538-553. doi:10.1016/0022-0396(70)90100-2Villafuerte, L., Braumann, C. A., Cortés, J.-C., & Jódar, L. (2010). Random differential operational calculus: Theory and applications. Computers & Mathematics with Applications, 59(1), 115-125. doi:10.1016/j.camwa.2009.08.061Cortés, J.-C., Jódar, L., Roselló, M.-D., & Villafuerte, L. (2012). Solving initial and two-point boundary value linear random differential equations: A mean square approach. Applied Mathematics and Computation, 219(4), 2204-2211. doi:10.1016/j.amc.2012.08.066Calatayud, J., Cortés, J.-C., Jornet, M., & Villafuerte, L. (2018). Random non-autonomous second order linear differential equations: mean square analytic solutions and their statistical properties. Advances in Difference Equations, 2018(1). doi:10.1186/s13662-018-1848-8Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). Improving the Approximation of the First- and Second-Order Statistics of the Response Stochastic Process to the Random Legendre Differential Equation. Mediterranean Journal of Mathematics, 16(3). doi:10.1007/s00009-019-1338-6Licea, J. A., Villafuerte, L., & Chen-Charpentier, B. M. (2013). Analytic and numerical solutions of a Riccati differential equation with random coefficients. Journal of Computational and Applied Mathematics, 239, 208-219. doi:10.1016/j.cam.2012.09.040Burgos, C., Calatayud, J., Cortés, J.-C., & Villafuerte, L. (2018). Solving a class of random non-autonomous linear fractional differential equations by means of a generalized mean square convergent power series. Applied Mathematics Letters, 78, 95-104. doi:10.1016/j.aml.2017.11.009Nouri, K., & Ranjbar, H. (2014). Mean Square Convergence of the Numerical Solution of Random Differential Equations. Mediterranean Journal of Mathematics, 12(3), 1123-1140. doi:10.1007/s00009-014-0452-8Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). Random differential equations with discrete delay. Stochastic Analysis and Applications, 37(5), 699-707. doi:10.1080/07362994.2019.1608833Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). $$\mathrm {L}^p$$-calculus Approach to the Random Autonomous Linear Differential Equation with Discrete Delay. Mediterranean Journal of Mathematics, 16(4). doi:10.1007/s00009-019-1370-6Caraballo, T., Cortés, J.-C., & Navarro-Quiles, A. (2019). Applying the random variable transformation method to solve a class of random linear differential equation with discrete delay. Applied Mathematics and Computation, 356, 198-218. doi:10.1016/j.amc.2019.03.048Zhou, T. (2014). A Stochastic Collocation Method for Delay Differential Equations with Random Input. Advances in Applied Mathematics and Mechanics, 6(4), 403-418. doi:10.4208/aamm.2012.m38Shi, W., & Zhang, C. (2017). Generalized polynomial chaos for nonlinear random delay differential equations. Applied Numerical Mathematics, 115, 16-31. doi:10.1016/j.apnum.2016.12.004Khusainov, D. Y., Ivanov, A. F., & Kovarzh, I. V. (2009). Solution of one heat equation with delay. Nonlinear Oscillations, 12(2), 260-282. doi:10.1007/s11072-009-0075-3Shaikhet, L. (2016). Stability of equilibrium states of a nonlinear delay differential equation with stochastic perturbations. International Journal of Robust and Nonlinear Control, 27(6), 915-924. doi:10.1002/rnc.3605Benhadri, M., & Zeghdoudi, H. (2018). Mean square asymptotic stability in nonlinear stochastic neutral Volterra-Levin equations with Poisson jumps and variable delays. Functiones et Approximatio Commentarii Mathematici, 58(2). doi:10.7169/facm/1657Santonja, F.-J., & Shaikhet, L. (2012). Analysing Social Epidemics by Delayed Stochastic Models. Discrete Dynamics in Nature and Society, 2012, 1-13. doi:10.1155/2012/530472Liu, L., & Caraballo, T. (2018). Analysis of a Stochastic 2D-Navier–Stokes Model with Infinite Delay. Journal of Dynamics and Differential Equations, 31(4), 2249-2274. doi:10.1007/s10884-018-9703-xLupulescu, V., & Abbas, U. (2011). Fuzzy delay differential equations. Fuzzy Optimization and Decision Making, 11(1), 99-111. doi:10.1007/s10700-011-9112-7Krapivsky, P. L., Luck, J. M., & Mallick, K. (2011). On stochastic differential equations with random delay. Journal of Statistical Mechanics: Theory and Experiment, 2011(10), P10008. doi:10.1088/1742-5468/2011/10/p10008GARRIDO-ATIENZA, M. J., OGROWSKY, A., & SCHMALFUSS, B. (2011). RANDOM DIFFERENTIAL EQUATIONS WITH RANDOM DELAYS. Stochastics and Dynamics, 11(02n03), 369-388. doi:10.1142/s0219493711003358Cortés, J.-C., Villafuerte, L., & Burgos, C. (2017). A Mean Square Chain Rule and its Application in Solving the Random Chebyshev Differential Equation. Mediterranean Journal of Mathematics, 14(1). doi:10.1007/s00009-017-0853-6Cortés, J. C., Jódar, L., & Villafuerte, L. (2007). Numerical solution of random differential equations: A mean square approach. Mathematical and Computer Modelling, 45(7-8), 757-765. doi:10.1016/j.mcm.2006.07.017Braumann, C. A., Cortés, J.-C., Jódar, L., & Villafuerte, L. (2018). On the random gamma function: Theory and computing. Journal of Computational and Applied Mathematics, 335, 142-155. doi:10.1016/j.cam.2017.11.045Khusainov, D. Y., & Pokojovy, M. (2015). Solving the Linear 1D Thermoelasticity Equations with Pure Delay. International Journal of Mathematics and Mathematical Sciences, 2015, 1-11. doi:10.1155/2015/47926

Physical states in the canonical tensor model from the perspective of random tensor networks

Tensor models, generalization of matrix models, are studied aiming for quantum gravity in dimensions larger than two. Among them, the canonical tensor model is formulated as a totally constrained system with first-class constraints, the algebra of which resembles the Dirac algebra of general relativity. When quantized, the physical states are defined to be vanished by the quantized constraints. In explicit representations, the constraint equations are a set of partial differential equations for the physical wave-functions, which do not seem straightforward to be solved due to their non-linear character. In this paper, after providing some explicit solutions for $N=2,3$, we show that certain scale-free integration of partition functions of statistical systems on random networks (or random tensor networks more generally) provides a series of solutions for general $N$. Then, by generalizing this form, we also obtain various solutions for general $N$. Moreover, we show that the solutions for the cases with a cosmological constant can be obtained from those with no cosmological constant for increased $N$. This would imply the interesting possibility that a cosmological constant can always be absorbed into the dynamics and is not an input parameter in the canonical tensor model. We also observe the possibility of symmetry enhancement in $N=3$, and comment on an extension of Airy function related to the solutions.Comment: 41 pages, 1 figure; typos correcte

Review of modern numerical methods for a simple vanilla option pricing problem

Option pricing is a very attractive issue of financial engineering and optimization. The problem of determining the fair price of an option arises from the assumptions made under a given financial market model. The increasing complexity of these market assumptions contributes to the popularity of the numerical treatment of option valuation. Therefore, the pricing and hedging of plain vanilla options under the Black–Scholes model usually serve as a bench-mark for the development of new numerical pricing approaches and methods designed for advanced option pricing models. The objective of the paper is to present and compare the methodological concepts for the valuation of simple vanilla options using the relatively modern numerical techniques in this issue which arise from the discontinuous Galerkin method, the wavelet approach and the fuzzy transform technique. A theoretical comparison is accompanied by an empirical study based on the numerical verification of simple vanilla option prices. The resulting numerical schemes represent a particularly effective option pricing tool that enables some features of options that are depend-ent on the discretization of the computational domain as well as the order of the polynomial approximation to be captured better

Quantum mechanics on non commutative spaces and squeezed states: a functional approach

We review here the quantum mechanics of some noncommutative theories in which no state saturates simultaneously all the non trivial Heisenberg uncertainty relations. We show how the difference of structure between the Poisson brackets and the commutators in these theories generically leads to a harmonic oscillator whose positions and momenta mean values are not strictly equal to the ones predicted by classical mechanics. This raises the question of the nature of quasi classical states in these models. We propose an extension based on a variational principle. The action considered is the sum of the absolute values of the expressions associated to the non trivial Heisenberg uncertainty relations. We first verify that our proposal works in the usual theory i.e we recover the known Gaussian functions. Besides them, we find other states which can be expressed as products of Gaussians with specific hyper geometrics. We illustrate our construction in two models defined on a four dimensional phase space: a model endowed with a minimal length uncertainty and the non commutative plane. Our proposal leads to second order partial differential equations. We find analytical solutions in specific cases. We briefly discuss how our proposal may be applied to the fuzzy sphere and analyze its shortcomings.Comment: 15 pages revtex. The title has been modified,the paper shortened and misprints have been corrected. Version to appear in JHE